computeDisplacement Function

public function computeDisplacement(x1, x2, boundaryData) result(r12)

computeDistance computes the shortest distance between points x1 and x2 In doing so it takes possible periodic boundaries into account.

Arguments

Type IntentOptional Attributes Name
real(kind=rk) :: x1(3)

xyz coordinates of point 1

real(kind=rk) :: x2(3)

xyz coordinates of point 2

type(mus_particle_boundarydata_type), intent(in) :: boundaryData

Boundary data tells us the domain bounds and if we have periodic bounds

Return Value real(kind=rk), (3)