portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-06-15 06:17:05 +01:00
Code Releases Activity

Remez tested

Browse Source
This commit is contained in:
Peter Boyle
2015-05-18 12:09:25 +01:00
parent 11cb3e9a01
commit 17835c6f42
7 changed files with 290 additions and 180 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
lib/algorithms/approx/Remez.cc
View File

@ -116,7 +116,7 @@ double AlgRemez::generateApprox(int num_degree, int den_degree,
int a_len, double *a_param, int *a_pow)
{
std::cout<<"Degree of the approximation is ("<<num_degree<<","<<den_degree<<")\n";
std::cout<<"Approximating the function x^("<<pnum<<"/"<<pden<<"\n";
std::cout<<"Approximating the function x^("<<pnum<<"/"<<pden<<")\n";
// Reallocate arrays, since degree has changed
if (num_degree != n || den_degree != d) allocate(num_degree,den_degree);
@ -149,6 +149,10 @@ double AlgRemez::generateApprox(int num_degree, int den_degree,
std::cout<<"Iteration " <<iter-1<<" spread "<<(double)spread<<" delta "<<(double)delta<<std::endl;
equations();
if (delta < tolerance) {
std::cout<<"Delta too small, try increasing precision\n";
assert(0);
};
assert( delta>= tolerance);
search(step);