public function hmap_eval_Jh_dq_aux(sf, q1, q2, q1_vec, q2_vec, q3_vec, xv) result(Jh_dq)
evaluate derivative of Jacobian of mapping h: sum_k dJ_h(q)/dq^k q_vec^k, k=1,2 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_dq
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_dq_aux(sf,q1,q2,q1_vec,q2_vec,q3_vec,xv)RESULT(Jh_dq)IMPLICIT NONE!-----------------------------------------------------------------------------------------------------------------------------------! INPUT VARIABLESCLASS(c_hmap),INTENT(IN)::sfREAL(wp),INTENT(IN)::q1,q2REAL(wp),INTENT(IN)::q1_vec,q2_vec,q3_vecCLASS(c_hmap_auxvar),INTENT(IN)::xv!-----------------------------------------------------------------------------------------------------------------------------------! OUTPUT VARIABLESREAL(wp)::Jh_dq!===================================================================================================================================Jh_dq=sf%eval_Jh_dq((/q1,q2,xv%zeta/),(/q1_vec,q2_vec,q3_vec/))END FUNCTION hmap_eval_Jh_dq_aux