Changeset 766 for library/tests/stresssuite
- Timestamp:
- 01/11/10 22:55:57 (15 years ago)
- Location:
- library/tests/stresssuite
- Files:
-
- 3 modified
Legend:
- Unmodified
- Added
- Removed
-
library/tests/stresssuite/arx_stress.cpp
r722 r766 24 24 //Test constructor 25 25 mat V0 = 0.00001 * eye ( ord + 1 ); 26 V0 ( 0 .0) = 1; //26 V0 ( 0 ) = 1; // 27 27 double nu0 = ord + 5.0; 28 28 -
library/tests/stresssuite/mixtures_stress.cpp
r750 r766 67 67 emix Simul; 68 68 Simul.set_rv ( x ); 69 Simul.set_parameters ( "0.5 0.6", Sim ); 69 Simul._w() = "0.5 0.6"; 70 Simul._Coms() = Sim; 70 71 Simul.validate(); 71 72 -
library/tests/stresssuite/square_mat_stress.cpp
r717 r766 60 60 61 61 tt.tic(); 62 TMatrix sq mat ( A );62 TMatrix sq_mat ( A ); 63 63 double elapsed = tt.toc(); 64 64 cout << "ctor(" << sz << " x " << sz << "): " << elapsed << " s" << endl; 65 65 66 66 tt.tic(); 67 mat res = sq mat.to_mat();67 mat res = sq_mat.to_mat(); 68 68 elapsed = tt.toc(); 69 69 … … 76 76 vec v = point->get_vector(); 77 77 double w = point->get_scalar(); 78 TMatrix sq mat2 = sqmat;79 80 tt.tic(); 81 sq mat2.opupdt ( v, w );78 TMatrix sq_mat2 = sq_mat; 79 80 tt.tic(); 81 sq_mat2.opupdt ( v, w ); 82 82 elapsed = tt.toc(); 83 83 84 84 if ( !fast ) { 85 85 mat expA = A + w * outer_product ( v, v ); 86 CHECK_CLOSE ( expA, sq mat2.to_mat(), epsilon );86 CHECK_CLOSE ( expA, sq_mat2.to_mat(), epsilon ); 87 87 } 88 88 … … 92 92 93 93 tt.tic(); 94 sq mat.inv ( invmat );94 sq_mat.inv ( invmat ); 95 95 elapsed = tt.toc(); 96 96 … … 104 104 105 105 tt.tic(); 106 double ld = sq mat.logdet();106 double ld = sq_mat.logdet(); 107 107 elapsed = tt.toc(); 108 108 … … 115 115 116 116 tt.tic(); 117 double q = sq mat.qform ( ones ( sz ) );117 double q = sq_mat.qform ( ones ( sz ) ); 118 118 elapsed = tt.toc(); 119 119 … … 125 125 126 126 tt.tic(); 127 q = sq mat.qform ( v );127 q = sq_mat.qform ( v ); 128 128 elapsed = tt.toc(); 129 129 … … 136 136 137 137 tt.tic(); 138 q = sq mat.invqform ( v );138 q = sq_mat.invqform ( v ); 139 139 elapsed = tt.toc(); 140 140 … … 146 146 cout << "invqform: " << elapsed << " s" << endl; 147 147 148 TMatrix twice = sq mat;149 150 tt.tic(); 151 twice += sq mat;148 TMatrix twice = sq_mat; 149 150 tt.tic(); 151 twice += sq_mat; 152 152 elapsed = tt.toc(); 153 153 … … 159 159 cout << "+=: " << elapsed << " s" << endl; 160 160 161 sq mat2 = sqmat;162 163 tt.tic(); 164 sq mat2.mult_sym ( A );161 sq_mat2 = sq_mat; 162 163 tt.tic(); 164 sq_mat2.mult_sym ( A ); 165 165 elapsed = tt.toc(); 166 166 167 167 if ( !fast ) { 168 168 res = ( A * A ) * A.T(); 169 CHECK_CLOSE ( res, sq mat2.to_mat(), epsilon );169 CHECK_CLOSE ( res, sq_mat2.to_mat(), epsilon ); 170 170 } 171 171