X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=check%2Ftime_vandermonde.cpp;h=675c1074e59643bca6af0fea8414559d43c0716e;hp=19415af99cb29d442c59cfd2bae8f287af5d03d4;hb=9f3e78a13c72ca73ade4d0829186edcdc1c1b2a0;hpb=695f6ae955ec530cded8f21efd5569df39447f76 diff --git a/check/time_vandermonde.cpp b/check/time_vandermonde.cpp index 19415af9..675c1074 100644 --- a/check/time_vandermonde.cpp +++ b/check/time_vandermonde.cpp @@ -7,7 +7,7 @@ */ /* - * GiNaC Copyright (C) 1999-2005 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2019 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 @@ -21,10 +21,16 @@ * * You should have received a copy of the GNU General Public License * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA */ -#include "times.h" +#include "ginac.h" +#include "timer.h" +using namespace GiNaC; + +#include +#include +using namespace std; static unsigned vandermonde_det(unsigned size) { @@ -52,7 +58,7 @@ static unsigned vandermonde_det(unsigned size) vanddet *= M(i,1) - M(j,1); if (expand(det - vanddet) != 0) { - clog << "Determaint of Vandermonde matrix " << endl + clog << "Determinant of Vandermonde matrix " << endl << "M==" << M << endl << "was miscalculated: det(M)==" << det << endl; ++result; @@ -66,17 +72,11 @@ unsigned time_vandermonde() unsigned result = 0; cout << "timing determinant of univariate symbolic Vandermonde matrices" << flush; - clog << "-------determinant of univariate symbolic Vandermonde matrices:" << endl; - vector sizes; + vector sizes = {8, 10, 12, 14}; vector times; timer swatch; - sizes.push_back(6); - sizes.push_back(8); - sizes.push_back(10); - sizes.push_back(12); - for (vector::iterator i=sizes.begin(); i!=sizes.end(); ++i) { int count = 1; swatch.start(); @@ -90,20 +90,23 @@ unsigned time_vandermonde() cout << '.' << flush; } - if (!result) { - cout << " passed "; - clog << "(no output)" << endl; - } else { - cout << " failed "; - } // print the report: cout << endl << " dim: "; for (vector::iterator i=sizes.begin(); i!=sizes.end(); ++i) cout << '\t' << *i << 'x' << *i; cout << endl << " time/s:"; for (vector::iterator i=times.begin(); i!=times.end(); ++i) - cout << '\t' << int(1000*(*i))*0.001; + cout << '\t' << *i; cout << endl; return result; } + +extern void randomify_symbol_serials(); + +int main(int argc, char** argv) +{ + randomify_symbol_serials(); + cout << setprecision(2) << showpoint; + return time_vandermonde(); +}