public function hmap_cyl_eval_gij_dq(sf, qL_in, q_G, qR_in, q_vec) result(g_ab_dq)
evaluate sum_ij (qL_i d/dq^k(G_ij(q_G)) qR_j) , k=1,2
where qL=(dX^1/dalpha,dX^2/dalpha [,dzeta/dalpha]) and qR=(dX^1/dbeta,dX^2/dbeta [,dzeta/dbeta]) and
where qL=(dX^1/dalpha,dX^2/dalpha ,dzeta/dalpha) and qR=(dX^1/dbeta,dX^2/dbeta ,dzeta/dbeta) and
dzeta_dalpha then known to be either 0.0 for ds and dtheta and 1.0 for dzeta
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_cyl_eval_gij_dq(sf,qL_in,q_G,qR_in,q_vec)RESULT(g_ab_dq)CLASS(t_hmap_cyl),INTENT(IN)::sfREAL(wp),INTENT(IN)::qL_in(3)REAL(wp),INTENT(IN)::q_G(3)REAL(wp),INTENT(IN)::qR_in(3)REAL(wp),INTENT(IN)::q_vec(3)!-----------------------------------------------------------------------------------------------------------------------------------! OUTPUT VARIABLESREAL(wp)::g_ab_dq!===================================================================================================================================g_ab_dq=0.0_wpEND FUNCTION hmap_cyl_eval_gij_dq