get_polynomial_3D Subroutine

public subroutine get_polynomial_3D(u, s, a, b, c)

Construct polynomial and interpolate at target position. To construct polynomial, velocity and strain rate of 4 source locations are needed: H(0,0,0), K(1,1,0), M(1,0,1), N(0,1,1)

Arguments

TypeIntentOptionalAttributesName
real(kind=rk), intent(in) :: u(3,4)
real(kind=rk), intent(in) :: s(6,4)
real(kind=rk), intent(out) :: a(10)
real(kind=rk), intent(out) :: b(10)
real(kind=rk), intent(out) :: c(10)

Contents


Variables

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: a0
real(kind=rk), private :: ax
real(kind=rk), private :: ay
real(kind=rk), private :: az
real(kind=rk), private :: axx
real(kind=rk), private :: ayy
real(kind=rk), private :: azz
real(kind=rk), private :: axy
real(kind=rk), private :: ayz
real(kind=rk), private :: axz
real(kind=rk), private :: b0
real(kind=rk), private :: bx
real(kind=rk), private :: by
real(kind=rk), private :: bz
real(kind=rk), private :: bxx
real(kind=rk), private :: byy
real(kind=rk), private :: bzz
real(kind=rk), private :: bxy
real(kind=rk), private :: byz
real(kind=rk), private :: bxz
real(kind=rk), private :: c0
real(kind=rk), private :: cx
real(kind=rk), private :: cy
real(kind=rk), private :: cz
real(kind=rk), private :: cxx
real(kind=rk), private :: cyy
real(kind=rk), private :: czz
real(kind=rk), private :: cxy
real(kind=rk), private :: cyz
real(kind=rk), private :: cxz
real(kind=rk), private :: P
real(kind=rk), private :: Q
real(kind=rk), private :: R
real(kind=rk), private :: Pa
real(kind=rk), private :: Qa
real(kind=rk), private :: Ra
real(kind=rk), private :: Pb
real(kind=rk), private :: Qb
real(kind=rk), private :: Rb
real(kind=rk), private :: Pc
real(kind=rk), private :: Qc
real(kind=rk), private :: Rc