Changeset 357 for applications/pmsm/TR2245/unitsteps.cpp
- Timestamp:
- 06/08/09 02:15:30 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
applications/pmsm/TR2245/unitsteps.cpp
r317 r357 17 17 #include <stat/libFN.h> 18 18 19 #include <stat/loggers_ui.h>20 #include <stat/libEF_ui.h>21 19 22 20 #include "../pmsm.h" 23 21 #include "simulator.h" 24 22 #include "../sim_profiles.h" 23 #include "user_info.h" 24 #include "stat/loggers.h" 25 25 26 26 using namespace bdm; … … 30 30 if ( argc>1 ) {fname = argv[1]; } 31 31 else { fname = "unitsteps.cfg"; } 32 UI File F ( fname );32 UI_File F ( fname ); 33 33 34 int Ndat;35 int Npart;36 34 double h = 1e-6; 37 35 int Nsimstep = 125; 38 36 37 38 // Kalman filter 39 int Ndat; 40 int Npart; 41 F.lookupValue ( "ndat", Ndat ); 42 F.lookupValue ( "Npart",Npart ); 43 mpdf* evolQ = UI::build<mpdf>( F, "Qrw" ); 39 44 vec Qdiag; 40 45 vec Rdiag; 46 UI::get( Qdiag, F, "dQ" ); //( "1e-6 1e-6 0.001 0.0001" ); //zdenek: 0.01 0.01 0.0001 0.0001 47 UI::get( Rdiag, F, "dR" );// ( "1e-8 1e-8" ); //var(diff(xth)) = "0.034 0.034" 41 48 42 mpdf* evolQ ; 43 try { 44 // Kalman filter 45 F.lookupValue ( "ndat", Ndat ); 46 F.lookupValue ( "Npart",Npart ); 47 48 UIbuild ( F.lookup ( "Qrw" ),evolQ ); 49 Qdiag= getvec ( F.lookup ( "dQ" ) ); //( "1e-6 1e-6 0.001 0.0001" ); //zdenek: 0.01 0.01 0.0001 0.0001 50 Rdiag=getvec ( F.lookup ( "dR" ) );// ( "1e-8 1e-8" ); //var(diff(xth)) = "0.034 0.034" 51 } 52 catch UICATCH; 53 // internal model 49 // internal model 54 50 55 51 IMpmsm fxu; … … 81 77 M.set_rv ( concat ( rQ,rx ) ); 82 78 83 dirfilelog *L ; UIbuild ( F.lookup ( "logger" ), L);// ( "exp/mpf_test",100 );79 dirfilelog *L = UI::build<dirfilelog>( F, "logger" );// ( "exp/mpf_test",100 ); 84 80 int l_X = L->add ( rx, "xt" ); 85 81 int l_D = L->add ( concat ( ry,ru ), "" ); … … 99 95 vec xtm=zeros ( 4 ); 100 96 double Ww=0.0; 101 vec vecW=getvec ( F.lookup ( "profile" ) ); 97 vec vecW; 98 UI::get( vecW, F, "profile" ); 102 99 103 100 for ( int tK=1;tK<Ndat;tK++ ) {