* Lewis and Michael Wester. */
/*
- * GiNaC Copyright (C) 1999-2001 Johannes Gutenberg University Mainz, Germany
+ * GiNaC Copyright (C) 1999-2002 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
hilbert.set(r,c,numeric(1,r+c+1));
ex det = hilbert.determinant();
- // The closed form of the determinant of n x n Hilbert matrices is:
- //
- // n-1 / n-1 \
- // ----- | ----- |
- // | | | pow(factorial(r),2) | | hilbert(r,c) |
- // | | | | | |
- // r = 0 \ c = 0 /
+ /*
+ The closed form of the determinant of n x n Hilbert matrices is:
+
+ n-1 / n-1 \
+ ----- | ----- |
+ | | | pow(factorial(r),2) | | hilbert(r,c) |
+ | | | | | |
+ r = 0 \ c = 0 /
+ */
ex hilbdet = 1;
for (unsigned r=0; r<n; ++r) {
double time = .0;
cout << "timing Lewis-Wester test H (det of 80x80 Hilbert)" << flush;
- clog << "-------Lewis-Wester test H (det of 80x80 Hilbert)" << endl;
+ clog << "-------Lewis-Wester test H (det of 80x80 Hilbert):" << endl;
rolex.start();
// correct for very small times: