Changeset 1064 for library/tests/stresssuite/particle_stress.cpp
- Timestamp:
- 06/09/10 14:00:40 (15 years ago)
- Files:
-
- 1 modified
-
library/tests/stresssuite/particle_stress.cpp (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
library/tests/stresssuite/particle_stress.cpp
r887 r1064 11 11 12 12 TEST ( particle_stress ) { 13 RV x ( "1" );14 RV xm = x;15 xm.t_plus ( -1 );16 const17 RV y ( "2" );13 RV x ( "1" ); 14 RV xm = x; 15 xm.t_plus ( -1 ); 16 const 17 RV y ( "2" ); 18 18 19 mat A = "1";20 vec vR = "1";21 ldmat R ( vR );19 mat A = "1"; 20 vec vR = "1"; 21 ldmat R ( vR ); 22 22 23 eEmp emp;24 euni eun;25 eun.set_parameters ( "0", "1" );26 emp.set_statistics ( ones ( 10 ), eun );27 vec &v = emp._w();28 Array<vec> &S = emp._samples();23 eEmp emp; 24 euni eun; 25 eun.set_parameters ( "0", "1" ); 26 emp.set_statistics ( ones ( 10 ), eun ); 27 vec &v = emp._w(); 28 Array<vec> &S = emp._samples(); 29 29 30 for ( int i = 0; i < 10; i++ ) {31 v ( i ) = exp ( -0.5 * sum ( pow ( S ( i ) - 1, 2.0 ) ) * 10 );32 }33 v /= sum ( v );30 for ( int i = 0; i < 10; i++ ) { 31 v ( i ) = exp ( -0.5 * sum ( pow ( S ( i ) - 1, 2.0 ) ) * 10 ); 32 } 33 v /= sum ( v ); 34 34 35 cout << "p:" << S << endl;36 cout << "w:" << v << endl;35 cout << "p:" << S << endl; 36 cout << "w:" << v << endl; 37 37 38 ivec ind;39 resample ( v, ind );38 ivec ind; 39 resample ( v, ind ); 40 40 41 cout << ind << endl;41 cout << ind << endl; 42 42 }
