public function hmap_eval_Jh_aux(sf, q1, q2, xv) result(Jh)
evaluate Jacobian of mapping h: J_h=sqrt(det(G)) at q=(X^1,X^2,zeta)
INFO: default routine that can be overwritten by specific hmap class,
not using additional hmap-dependent auxiliary variables,
but calling the pointwise routine eval_Jh
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed
arrows point from an interface to procedures which implement that interface.
This could include the module procedures in a generic interface or the
implementation in a submodule of an interface in a parent module.
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed
arrows point from an interface to procedures which implement that interface.
This could include the module procedures in a generic interface or the
implementation in a submodule of an interface in a parent module.
Source Code
FUNCTION hmap_eval_Jh_aux(sf,q1,q2,xv)RESULT(Jh)IMPLICIT NONE!-----------------------------------------------------------------------------------------------------------------------------------! INPUT VARIABLESCLASS(c_hmap),INTENT(IN)::sfREAL(wp),INTENT(IN)::q1,q2CLASS(c_hmap_auxvar),INTENT(IN)::xv!-----------------------------------------------------------------------------------------------------------------------------------! OUTPUT VARIABLESREAL(wp)::Jh!===================================================================================================================================Jh=sf%eval_Jh((/q1,q2,xv%zeta/))END FUNCTION hmap_eval_Jh_aux