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

changed: do not do first adaptive step with dense solver

Browse Source 
not a good idea if the initial model is of a non-trivial size
This commit is contained in:
Arne Morten Kvarving 2016-11-17 10:28:38 +01:00
parent 17bb95ace9
commit 18d3d976e2
1 changed files with 1 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
src/SIM/AdaptiveSIM.C
View File

@ -204,8 +204,7 @@ bool AdaptiveSIM::solveStep (const char* inputfile, int iStep)
// Assemble the linear FE equation system // Assemble the linear FE equation system
model.setMode(SIM::STATIC,true); model.setMode(SIM::STATIC,true);
model.initSystem(iStep == 1 ? SystemMatrix::DENSE : opt.solver, model.initSystem(opt.solver, 1, model.getNoRHS());
1, model.getNoRHS());
model.setQuadratureRule(opt.nGauss[0],true); model.setQuadratureRule(opt.nGauss[0],true);
if (!model.assembleSystem()) if (!model.assembleSystem())