Changeset 80ba530 in git
- Timestamp:
- Feb 1, 2010, 11:29:06 AM (13 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a657104b677b4c461d018cbf3204d72d34ad66a9')
- Children:
- 835604146a3a1dabecdf5653b203f926bfc98aa1
- Parents:
- fc323e3421a512eb4770e20b1957b2e139709eb9
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
kernel/matpol.cc
rfc323e r80ba530 50 50 static void mpSwapCol(matrix, int, int, int); 51 51 static void mpElimBar(matrix, matrix, poly, int, int); 52 53 #ifdef HAVE_MINOR54 unsigned long SM_MULT_Counter = 0;55 unsigned long SM_DIV_Counter = 0;56 int recursionLevel = 0;57 #endif // HAVE_MINOR58 52 59 53 /*2 … … 458 452 matrix nextLevel=mpNew(kr,kc); 459 453 460 #ifdef HAVE_MINOR461 if (recursionLevel == 0)462 {463 SM_MULT_Counter = 0; // for counting basic operations464 SM_DIV_Counter = 0; // for counting basic operations465 }466 #endif // HAVE_MINOR467 468 454 loop 469 455 { … … 479 465 if (ar>1) 480 466 { 481 #ifdef HAVE_MINOR482 recursionLevel++;483 #endif // HAVE_MINOR484 467 mpRecMin(ar-1,result,elems,nextLevel,kr,k,a->m[kr*a->ncols+k],R); 485 #ifdef HAVE_MINOR486 recursionLevel--;487 #endif // HAVE_MINOR488 468 mpPartClean(nextLevel,kr,k); 489 469 } … … 496 476 kr--; 497 477 } 498 #ifdef HAVE_MINOR499 if (recursionLevel == 0)500 {501 char h[100];502 sprintf(h, "%s%lu", "total number of calls to 'SM_MULT': ", SM_MULT_Counter);503 PrintLn(); PrintS(h);504 sprintf(h, "%s%lu", "total number of calls to 'SM_DIV': ", SM_DIV_Counter);505 PrintLn(); PrintS(h); PrintLn();506 }507 #endif // HAVE_MINOR508 478 mpFinalClean(nextLevel); 509 479 } … … 1857 1827 { 1858 1828 q1 = SM_MULT(a[j], piv, div); 1859 #ifdef HAVE_MINOR1860 SM_MULT_Counter++;1861 #endif // HAVE_MINOR1862 1829 if (ap[j] != NULL) 1863 1830 { 1864 q2 = SM_MULT(ap[j], elim, div); 1865 #ifdef HAVE_MINOR 1866 SM_MULT_Counter++; 1867 #endif // HAVE_MINOR 1831 q2 = SM_MULT(ap[j], elim, div); 1868 1832 q1 = pAdd(q1,q2); 1869 1833 } 1870 1834 } 1871 1835 else if (ap[j] != NULL) 1872 q1 = SM_MULT(ap[j], elim, div); 1873 #ifdef HAVE_MINOR 1874 SM_MULT_Counter++; 1875 #endif // HAVE_MINOR 1836 q1 = SM_MULT(ap[j], elim, div); 1876 1837 if (q1 != NULL) 1877 1838 { 1878 1839 if (div) 1879 1840 SM_DIV(q1, div); 1880 #ifdef HAVE_MINOR1881 SM_DIV_Counter++;1882 #endif // HAVE_MINOR1883 1841 q[j] = q1; 1884 1842 } … … 1891 1849 if (a[j] != NULL) 1892 1850 { 1893 q1 = SM_MULT(a[j], piv, div); 1894 #ifdef HAVE_MINOR 1895 SM_MULT_Counter++; 1896 #endif // HAVE_MINOR 1851 q1 = SM_MULT(a[j], piv, div); 1897 1852 if (div) 1898 1853 SM_DIV(q1, div); 1899 #ifdef HAVE_MINOR1900 SM_DIV_Counter++;1901 #endif // HAVE_MINOR1902 1854 q[j] = q1; 1903 1855 }
Note: See TracChangeset
for help on using the changeset viewer.