OilfieldToolsNet/IFEM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

remove code that is unreachable

Browse Source
This commit is contained in:
Arne Morten Kvarving
2021-09-30 12:06:36 +02:00
parent 88bcd66794
commit 24456cb047
1 changed files with 2 additions and 13 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
src/ASM/SplineInterpolator.C
View File

@@ -42,19 +42,13 @@ void SplineInterpolator::interpolate(const std::vector<double>& params,
int ti = 0; // index to first unused element of tangent_points
std::vector<double> tmp(2*order);
for (i = 0; i < num_points; ++i) {
bool der = false;
double par = params[i];
ki = basis.knotIntervalFuzzy(par); // knot-interval of param.
basis.computeBasisValues(params[i], &tmp[0], 1);
for (j = 0; j < order; ++j)
if ((ki-order+1+j>=0) && (ki-order+1+j<=num_coefs)) {
if ((ki-order+1+j>=0) && (ki-order+1+j<=num_coefs))
A(i+ti+1,ki-order+1+j+1) = tmp[2*j];
if (der)
A(i+ti+1,ki-order+1+j+1) = tmp[2*j+1];
}
if (der)
++ti;
}
// generating right-hand side
@@ -94,18 +88,13 @@ void SplineInterpolator::quasiinterpolate(const std::vector<double>& params,
// setting up interpolation matrix A
std::vector<double> tmp(2*order);
for (i = 0; i < num_points; ++i) {
bool der = false;
double par = params[i];
ki = basis.knotIntervalFuzzy(par);
basis.computeBasisValues(params[i], &tmp[0], 1);
for (j = 0; j < order; ++j)
if((ki-order+1+j-index>=0) && (ki-order+1+j-index<num_coefs))
{
if ((ki-order+1+j-index>=0) && (ki-order+1+j-index<num_coefs))
A(i+1,ki-order+1+j-index+1) = tmp[2*j];
if (der)
A(i+1+1,ki-order+1+j+1) = tmp[2*j+1];
}
}
// generating right-hand side