Skip to content
Snippets Groups Projects
Commit dc82b7a9 authored by Peter Bannister's avatar Peter Bannister
Browse files

Cleaned up verbose reporting in kernel::hermite_interpolation1d

parent 0a432f69
No related branches found
No related tags found
No related merge requests found
...@@ -167,10 +167,10 @@ namespace MISCMATHS { ...@@ -167,10 +167,10 @@ namespace MISCMATHS {
float retval, r1 = 0.0, r4 = 0.0; float retval, r1 = 0.0, r4 = 0.0;
if (!in_bounds(data,p1) || !in_bounds(data,p4)) { if (!in_bounds(data,p1) || !in_bounds(data,p4)) {
cerr << "ERROR: One or more indicies lie outside the data range. Returning ZERO" << endl; cerr << "Hermite Interpolation - ERROR: One or more indicies lie outside the data range. Returning ZERO" << endl;
retval = 0.0; retval = 0.0;
} else if ((t < 0) || (t > 1)) { } else if ((t < 0) || (t > 1)) {
cerr << "ERROR: Interpolation index must lie between 0 and 1. Returning ZERO" << endl; cerr << "Hermite Interpolation - ERROR: Interpolation index must lie between 0 and 1. Returning ZERO" << endl;
retval = 0.0; retval = 0.0;
/* } else if (t == 0.0) { /* } else if (t == 0.0) {
retval = data(p1); retval = data(p1);
...@@ -185,7 +185,7 @@ namespace MISCMATHS { ...@@ -185,7 +185,7 @@ namespace MISCMATHS {
retval = (2*t3 - 3*t2 + 1)*data(p1) + (-2*t3 + 3*t2)*data(p4) + (t3 - 2*t2 + t)*r1 + (t3 - t2)*r4; retval = (2*t3 - 3*t2 + 1)*data(p1) + (-2*t3 + 3*t2)*data(p4) + (t3 - 2*t2 + t)*r1 + (t3 - t2)*r4;
} }
cerr << "p1, p4, t, r1, r4 = " << p1 << ", " << p4 << ", " << t << ", " << r1 << ", " << r4 << endl; // cerr << "p1, p4, t, r1, r4 = " << p1 << ", " << p4 << ", " << t << ", " << r1 << ", " << r4 << endl;
return retval; return retval;
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment