-// Copyright (C) 2007-2014 CEA/DEN, EDF R&D
+// Copyright (C) 2007-2021 CEA/DEN, EDF R&D
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
calculateDeterminant();
LOG(3, "determinant before inverse = " << _determinant);
+
+ double ni(1./INTERP_KERNEL::normInf(_linear_transform));
+ ni = ni*ni*ni;
// check that tetra is non-planar -> determinant is not zero
+ // AGY : the check to 0. must integrate the infinite norm of _linear_transform matrix.
// otherwise set _determinant to zero to signal caller that transformation did not work
- if(epsilonEqual(_determinant, 0.0))
+ if(epsilonEqual(ni*_determinant, 0.0))
{
_determinant = 0.0;
return;