- fixed three little standard-conformance issues.
[ginac.git] / check / time_vandermonde.cpp
index 364ca5a5ec8c340302494a53f261fa0304e26169..f878b0b886ddc398355b836d97e694852838e61e 100644 (file)
@@ -7,7 +7,7 @@
  */
 
 /*
- *  GiNaC Copyright (C) 1999-2000 Johannes Gutenberg University Mainz, Germany
+ *  GiNaC Copyright (C) 1999-2001 Johannes Gutenberg University Mainz, Germany
  *
  *  This program is free software; you can redistribute it and/or modify
  *  it under the terms of the GNU General Public License as published by
 
 static unsigned vandermonde_det(unsigned size)
 {
-    unsigned result = 0;
-    symbol a("a");
-    
-    // construct Vandermonde matrix:
-    matrix M(size,size);
-    for (unsigned ro=0; ro<size; ++ro)
-        for (unsigned co=0; co<size; ++co)
-            if (ro%2)
-                M.set(ro,co,pow(-pow(a,1+ro/2),co));
-            else
-                M.set(ro,co,pow(pow(a,1+ro/2),co));
-    
-    // compute determinant:
-    ex vdet = M.determinant();
-    
-    // dirty consistency check of result:
-    if (!vdet.subs(a==1).is_zero()) {
-        clog << "Determaint of Vandermonde matrix " << endl
-             << "M==" << M << endl
-             << "was miscalculated: det(M)==" << vdet << endl;
-        ++result;
-    }
-    
-    return result;
+       unsigned result = 0;
+       symbol a("a");
+       
+       // construct Vandermonde matrix:
+       matrix M(size,size);
+       for (unsigned ro=0; ro<size; ++ro) {
+               for (unsigned co=0; co<size; ++co) {
+                       if (ro%2)
+                               M.set(ro,co,pow(-pow(a,1+ro/2),co));
+                       else
+                               M.set(ro,co,pow(pow(a,1+ro/2),co));
+               }
+       }
+       
+       // compute determinant:
+       ex vdet = M.determinant();
+       
+       // dirty consistency check of result:
+       if (!vdet.subs(a==1).is_zero()) {
+               clog << "Determaint of Vandermonde matrix " << endl
+                    << "M==" << M << endl
+                    << "was miscalculated: det(M)==" << vdet << endl;
+               ++result;
+       }
+       
+       return result;
 }
 
 unsigned time_vandermonde(void)
 {
-    unsigned result = 0;
-    
-    cout << "timing determinant of univariate symbolic Vandermonde matrices" << flush;
-    clog << "-------determinant of univariate symbolic Vandermonde matrices:" << endl;
-    
-    vector<unsigned> sizes;
-    vector<double> times;
-    timer swatch;
-    
-    sizes.push_back(4);
-    sizes.push_back(6);
-    sizes.push_back(8);
-    sizes.push_back(10);
-    
-    for (vector<unsigned>::iterator i=sizes.begin(); i!=sizes.end(); ++i) {
-        int count = 1;
-        swatch.start();
-        result += vandermonde_det(*i);
-        // correct for very small times:
-        while (swatch.read()<0.02) {
-            vandermonde_det(*i);
-            ++count;
-        }
-        times.push_back(swatch.read()/count);
-        cout << '.' << flush;
-    }
-    
-    if (!result) {
-        cout << " passed ";
-        clog << "(no output)" << endl;
-    } else {
-        cout << " failed ";
-    }
-    // print the report:
-    cout << endl << "    dim:   ";
-    for (vector<unsigned>::iterator i=sizes.begin(); i!=sizes.end(); ++i)
-        cout << '\t' << *i << 'x' << *i;
-    cout << endl << "    time/s:";
-    for (vector<double>::iterator i=times.begin(); i!=times.end(); ++i)
-        cout << '\t' << int(1000*(*i))*0.001;
-    cout << endl;
-    
-    return result;
+       unsigned result = 0;
+       
+       cout << "timing determinant of univariate symbolic Vandermonde matrices" << flush;
+       clog << "-------determinant of univariate symbolic Vandermonde matrices:" << endl;
+       
+       vector<unsigned> sizes;
+       vector<double> times;
+       timer swatch;
+       
+       sizes.push_back(4);
+       sizes.push_back(6);
+       sizes.push_back(8);
+       sizes.push_back(10);
+       
+       for (vector<unsigned>::iterator i=sizes.begin(); i!=sizes.end(); ++i) {
+               int count = 1;
+               swatch.start();
+               result += vandermonde_det(*i);
+               // correct for very small times:
+               while (swatch.read()<0.02) {
+                       vandermonde_det(*i);
+                       ++count;
+               }
+               times.push_back(swatch.read()/count);
+               cout << '.' << flush;
+       }
+       
+       if (!result) {
+               cout << " passed ";
+               clog << "(no output)" << endl;
+       } else {
+               cout << " failed ";
+       }
+       // print the report:
+       cout << endl << "       dim:   ";
+       for (vector<unsigned>::iterator i=sizes.begin(); i!=sizes.end(); ++i)
+               cout << '\t' << *i << 'x' << *i;
+       cout << endl << "       time/s:";
+       for (vector<double>::iterator i=times.begin(); i!=times.end(); ++i)
+               cout << '\t' << int(1000*(*i))*0.001;
+       cout << endl;
+       
+       return result;
 }