Changeset 34 for include/fflas-ffpack/ffpack.h
- Timestamp:
- 08/27/07 10:16:36 (1 year ago)
- Files:
-
- 1 modified
-
include/fflas-ffpack/ffpack.h (modified) (12 diffs)
Legend:
- Unmodified
- Added
- Removed
-
include/fflas-ffpack/ffpack.h
r33 r34 269 269 size_t *rowP = new size_t[M]; 270 270 271 Timer t1;272 t1.clear();273 t1.start();271 // Timer t1; 272 // t1.clear(); 273 // t1.start(); 274 274 275 275 nullity = M - LUdivine( F, FflasNonUnit, FflasNoTrans, M, M, A, lda, P, rowP, FfpackLQUP); 276 276 277 t1.stop();277 // t1.stop(); 278 278 //cerr<<"LU --> "<<t1.usertime()<<endl; 279 279 … … 284 284 } else { 285 285 // Initializing X to 0 286 t1.clear();287 t1.start();286 // t1.clear(); 287 // t1.start(); 288 288 for (size_t i=0; i<M; ++i) 289 289 for (size_t j=0; j<M;++j) … … 299 299 for (size_t i=1; i<M; ++i) 300 300 fcopy (F, i, (X+i*ldx), 1, (A+i*lda), 1); 301 t1.stop();301 // t1.stop(); 302 302 //cerr<<"U^-1 --> "<<t1.usertime()<<endl; 303 303 … … 305 305 // X = Q^-1.X is not necessary since Q = Id 306 306 307 // X = U^-1.X308 t1.clear();309 t1.start();307 // X = U^-1.X 308 // t1.clear(); 309 // t1.start(); 310 310 ftrsm( F, FflasLeft, FflasUpper, FflasNoTrans, FflasNonUnit, 311 311 M, M, one, A, lda , X, ldx); 312 t1.stop();312 // t1.stop(); 313 313 //cerr<<"ftrsm --> "<<t1.usertime()<<endl; 314 314 … … 512 512 size_t r; 513 513 514 Timer t1;515 t1.clear();516 t1.start();514 // Timer t1; 515 // t1.clear(); 516 // t1.start(); 517 517 r = LUdivine (F, FflasUnit, FflasNoTrans, M, N, A, lda, P, Qt); 518 t1.stop();518 // t1.stop(); 519 519 //cerr<<"LU --> "<<t1.usertime()<<endl; 520 520 521 Timer t2;522 t2.clear();523 t2.start();521 // Timer t2; 522 // t2.clear(); 523 // t2.start(); 524 524 ftrtri (F, FflasUpper, FflasUnit, r, A, lda); 525 525 … … 528 528 mone, A, lda, A+r, lda); 529 529 530 t2.stop();530 // t2.stop(); 531 531 //cerr<<"U^-1 --> "<<t2.usertime()<<endl; 532 532 … … 553 553 size_t r; 554 554 555 Timer t1;556 t1.clear();557 t1.start();555 // Timer t1; 556 // t1.clear(); 557 // t1.start(); 558 558 r = LUdivine (F, FflasUnit, FflasTrans, M, N, A, lda, P, Qt); 559 t1.stop();559 // t1.stop(); 560 560 //cerr<<"LU --> "<<t1.usertime()<<endl; 561 561 562 Timer t2;563 t2.clear();564 t2.start();562 // Timer t2; 563 // t2.clear(); 564 // t2.start(); 565 565 ftrtri (F, FflasLower, FflasUnit, r, A, lda); 566 566 … … 569 569 mone, A, lda, A+r*lda, lda); 570 570 571 t2.stop();571 // t2.stop(); 572 572 //cerr<<"U^-1 --> "<<t2.usertime()<<endl; 573 573 … … 597 597 r = ColumnEchelonForm (F, M, N, A, lda, P, Qt); 598 598 599 Timer t1;600 t1.clear();601 t1.start();599 // Timer t1; 600 // t1.clear(); 601 // t1.start(); 602 602 603 603 // M = Q^T M … … 616 616 617 617 ftrtrm (F, FflasUnit, r, A, lda); 618 t1.stop();618 // t1.stop(); 619 619 //cerr<<"U^-1L^-1 --> "<<t1.usertime()<<endl; 620 620 … … 645 645 r = RowEchelonForm (F, M, N, A, lda, P, Qt); 646 646 647 Timer t1;648 t1.clear();649 t1.start();647 // Timer t1; 648 // t1.clear(); 649 // t1.start(); 650 650 // M = M Q^T 651 651 for (int i=0; i<r; ++i){ … … 664 664 ftrtrm (F, FflasNonUnit, r, A, lda); 665 665 666 t1.stop();666 // t1.stop(); 667 667 //cerr<<"U^-1L^-1 --> "<<t1.usertime()<<endl; 668 668
