Weight implemented by Monique Dauge in Melina. More...
#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.
ENDIFThis computes a correct weight for the Fichera corner.
Definition at line 227 of file shortestDist.hh.
| Real hp3D::DaugeWeight::operator() | ( | const SingularSet & | singularities, |
| concepts::Real3d | point, | ||
| const Hexahedron * | elm | ||
| ) | const |