#include <shortestDist.hh>
Public Member Functions | |
Real | operator() (const SingularSet &singularities, concepts::Real3d point, const Hexahedron *elm) const |
Weight implemented by Monique Dauge in Melina. Her Fortran code:
XX= POINT(1) YY= POINT(2) ZZ= POINT(3) R1=SQRT (ZZ*ZZ+YY*YY) R2=SQRT (ZZ*ZZ+XX*XX) R3=SQRT (XX*XX+YY*YY) IF (XX+YY+ZZ.GE.0.) THEN R=SQRT (XX*XX+YY*YY+ZZ*ZZ) ELSE R=MIN(MIN(R1,R2),R3) ENDIF IF (R.GE.1.E-12) THEN RHO=S*EXP(ALPHA*LOG(R)) ELSE RHO=0. ENDIF
This computes a correct weight for the Fichera corner.
Definition at line 229 of file shortestDist.hh.