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