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

Fixed: Account for small round-offs in start time for adaptive loop.

Browse Source 
Changed: Ignore pre-refinement tags when mesh already is refined.
This commit is contained in:
Knut Morten Okstad 2021-05-27 16:21:00 +02:00
parent 1561294100
commit ee17cff565
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Apps/Common/SIMSolverTS.h
View File

@ -80,7 +80,7 @@ public:
this->printHeading(heading);
// Pre-adaptive loop, solve for a certain time period on the initial mesh
while (this->tp.time.t < aStart && this->advanceStep())
while (this->tp.time.t < aStart-1.0e-12 && this->advanceStep())
if (!this->S1.solveStep(this->tp))
return 3;
else if (!this->saveState(geoBlk,nBlock))
@ -214,6 +214,8 @@ protected:
maxPred = atoi(value);
else if ((value = utl::getValue(child,"min_frac")))
minFrac = atof(value);
else if (this->S1.getRefined())
continue; // ignore pre-refinement tags if mesh already refined
else if ((value = utl::getValue(child,"prerefine")))
{
if (utl::getAttribute(child,"levels",preRef))