Compute the Jacobian matrix of distance Rij = sqrt((ci(x)-cj(x))2 + (ci(y)-cj(y))2 + (ci(z)-cj(z))**2) Derivatives wrt ci(:) are saved in J_i and wrt cj(:) in J_j; the distance between the two points is also provided in output in Rij.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=rp), | intent(in) | :: | ci(3) | |||
real(kind=rp), | intent(in) | :: | cj(3) | |||
real(kind=rp), | intent(out) | :: | Rij | |||
real(kind=rp), | intent(out) | :: | J_i(3) | |||
real(kind=rp), | intent(out) | :: | J_j(3) |