OilfieldToolsNet/opm-upscaling
4
0
Fork 0
mirror of https://github.com/OPM/opm-upscaling.git synced 2025-02-25 18:45:23 -06:00
Code Issues Packages Projects Releases Wiki Activity

added: warn if volume cell is degenerated

Browse Source 
also gives the collapse tolerance that will get rid of the cell
This commit is contained in:
Arne Morten Kvarving 2012-11-20 11:23:52 +01:00
parent 9c56c926f7
commit 29ec9824e9
1 changed files with 7 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
dune/elasticity/elasticity_upscale_impl.hpp
View File

@ -578,8 +578,13 @@ void ElasticityUpscale<GridType>::assemble(int loadcase, bool matrix)
it->geometry().jacobianInverseTransposed(r->position());
ctype detJ = it->geometry().integrationElement(r->position());
if (detJ <= 0)
continue;
if (detJ <= 1.e-4) {
std::cout << "cell " << m << " is (close to) degenerated, detJ " << detJ << std::endl;
double zdiff=0.0;
for (int i=0;i<4;++i)
zdiff = std::max(zdiff, it->geometry().corner(i+4)[2]-it->geometry().corner(i)[2]);
std::cout << " - Consider setting ctol larger than " << zdiff << std::endl;
}
Dune::FieldMatrix<ctype,comp,dim*bfunc> B;
E.getBmatrix(B,r->position(),jacInvTra);