Changeset 2557 in Sophya for trunk/SophyaProg/Tests
- Timestamp:
- Jul 21, 2004, 6:23:00 PM (21 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaProg/Tests/tsttminv.cc
r2555 r2557 21 21 22 22 // --- Choix de ce que doit faire Lapack 23 #ifdef USE_LAPACK24 23 #define ALSO_LAPACK_INV 25 24 #define ALSO_LAPACK_INV_SYM … … 27 26 #define ALSO_LAPACK_EV 28 27 #define ALSO_LAPACK_EV_SYM 29 #endif30 28 31 29 ////////////////////////////////////////////////// … … 161 159 <<"-- Only GausPiv --"<<endl 162 160 <<"-t 0/1/2 : data scaling 0=no, 1=global (def), 2=row-by-row"<<endl 163 <<"-d : also compute determinant "<<endl;161 <<"-d : also compute determinant (ne pas utiliser si N est grand)"<<endl; 164 162 return(-1); 165 163 } … … 257 255 258 256 259 #ifdef ALSO_LAPACK_INV 260 ///////// ///////////////////////////261 ///////// Test avec Lapack/////////262 //////////////////////////////////// 257 ////////////////////////////////////////////// 258 ///////// Test Inversion avec Lapack ///////// 259 ////////////////////////////////////////////// 260 #if defined(USE_LAPACK) && defined(ALSO_LAPACK_INV) 263 261 { 264 262 cout<<"\n=========================================="<<endl; … … 282 280 283 281 284 #ifdef ALSO_LAPACK_INV_SYM 285 ///////// ///////////////////////////////286 ///////// Test avec Lapack sym/////////287 //////////////////////////////////////// 282 ////////////////////////////////////////////////// 283 ///////// Test Inversion avec Lapack sym ///////// 284 ////////////////////////////////////////////////// 285 #if defined(USE_LAPACK) && defined(ALSO_LAPACK_INV_SYM) 288 286 { 289 287 cout<<"\n=========================================="<<endl; … … 307 305 308 306 309 #ifdef ALSO_LAPACK_INV_LSS310 307 //////////////////////////////////////////////// 311 308 ///////// Test avec Lapack LeastSquare ///////// 312 309 //////////////////////////////////////////////// 310 #if defined(USE_LAPACK) && defined(ALSO_LAPACK_INV_LSS) 313 311 { 314 312 cout<<"\n=========================================="<<endl; … … 332 330 333 331 334 #ifdef ALSO_LAPACK_EV335 332 //////////////////////////////////////////////// 336 333 ///////// Test avec Lapack pour EV ///////// 337 334 //////////////////////////////////////////////// 335 #if defined(USE_LAPACK) && defined(ALSO_LAPACK_EV) 338 336 { 339 337 cout<<"\n=========================================="<<endl; … … 347 345 PrtTim("--- End of LapackEigen decompositon ---"); 348 346 if(nprline>0) {cout<<Eval; cout<<endl; cout<<Evec; cout<<endl;} 349 #ifndef COMPLEX350 347 //-- Find the complex conjugate pairs 348 #if !defined(COMPLEX) 351 349 TVector< uint_2 > Evalconj(N); Evalconj = 0; 352 350 int_4 nconj=0; … … 369 367 for(int i=0;i<N;i++) { 370 368 complex<TYPER> Evec_il; 371 #if def COMPLEX369 #if defined(COMPLEX) 372 370 Evec_il = Evec(i,l); 373 371 #else … … 378 376 for(int j=0;j<N;j++) { 379 377 complex<TYPER> Evec_jl; 380 #if def COMPLEX378 #if defined(COMPLEX) 381 379 Evec_jl = Evec(j,l); 382 380 #else … … 398 396 399 397 400 #ifdef ALSO_LAPACK_EV_SYM401 398 //////////////////////////////////////////////// 402 399 ///////// Test avec Lapack sym pour EV ///////// 403 400 //////////////////////////////////////////////// 401 #if defined(USE_LAPACK) && defined(ALSO_LAPACK_EV_SYM) 404 402 { 405 403 cout<<"\n=========================================="<<endl; … … 428 426 429 427 430 #ifdef USE_GAUSPIV 431 ///////// ///////////////////////////432 ///////// Test avec GausPiv/////////433 //////////////////////////////////// 428 /////////////////////////////////////////////// 429 ///////// Test Inversion avec GausPiv ///////// 430 /////////////////////////////////////////////// 431 #if defined(USE_GAUSPIV) 434 432 { 435 433 cout<<"\n==========================================\n"
Note:
See TracChangeset
for help on using the changeset viewer.