Changeset 1064 for library/tests/stresssuite/particle_stress.cpp
- Timestamp:
- 06/09/10 14:00:40 (14 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
library/tests/stresssuite/particle_stress.cpp
r887 r1064 11 11 12 12 TEST ( particle_stress ) { 13 14 15 16 17 13 RV x ( "1" ); 14 RV xm = x; 15 xm.t_plus ( -1 ); 16 const 17 RV y ( "2" ); 18 18 19 20 21 19 mat A = "1"; 20 vec vR = "1"; 21 ldmat R ( vR ); 22 22 23 24 25 26 27 28 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 31 32 33 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 36 35 cout << "p:" << S << endl; 36 cout << "w:" << v << endl; 37 37 38 39 38 ivec ind; 39 resample ( v, ind ); 40 40 41 41 cout << ind << endl; 42 42 }