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

Added: Support centered geometric refinement (refinement towards both sides)

Browse Source
This commit is contained in:
Arne Morten Kvarving 2016-10-06 12:54:17 +02:00 committed by Knut Morten Okstad
parent 5d3de9fd16
commit 04943191b2
1 changed files with 12 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
src/Utility/Utilities.C
View File

@ -45,6 +45,7 @@ bool utl::parseKnots (std::vector<Real>& xi)
{
// Geometric grading
bool biased = toupper(cline[0]) == 'B';
bool centred = toupper(cline[0]) == 'C';
int ru = atoi(strtok(nullptr," "));
double alpha = atof(strtok(nullptr," "));
double xi1 = (cline = strtok(nullptr," ")) ? atof(cline) : 0.0;
@ -52,11 +53,13 @@ bool utl::parseKnots (std::vector<Real>& xi)
if (xi1 < 0.0 || xi2 <= xi1 || xi2 > 1.0 || ru < 1)
return false;
if (biased && ru > 1 && alpha != 1.0)
alpha = pow(alpha,1.0/Real(ru));
alpha = pow(alpha,1.0/double(ru));
else if (centred && ru > 1)
ru /= 2;
double x = xi1;
double D = xi2 - xi1;
D *= (alpha == 1.0 ? 1.0/Real(ru+1) : (1.0-alpha)/(1.0-pow(alpha,ru+1)));
double D = (xi2-xi1) * (centred ? 0.5 : 1.0);
D *= (alpha == 1.0 ? 1.0/double(ru+1) : (1.0-alpha)/(1.0-pow(alpha,ru+1)));
if (xi1 > 0.0) xi.push_back(xi1);
for (int i = 0; i < ru; i++)
{
@ -64,6 +67,12 @@ bool utl::parseKnots (std::vector<Real>& xi)
D *= alpha;
xi.push_back(x);
}
if (centred)
{
xi.push_back(0.5*(xi1+xi2));
for (int i = 1; i <= ru; i++)
xi.push_back(xi2 - xi[ru-i] + xi1);
}
if (xi2 < 1.0) xi.push_back(xi2);
}
else