Points x(3, n) into p(xyz and distance, high/low, n): the floats nearest, and the rest (the distances are left as they are).
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=wp), | intent(in) | :: | x(:,:) |
(3, n) |
||
| real(kind=c_float), | intent(inout) | :: | p(:,:,:) |
(4, 2, n) |
pure subroutine split_points(x, p) real(wp), intent(in) :: x(:,:) !! (3, n) real(c_float), intent(inout) :: p(:,:,:) !! (4, 2, n) p(1:3, 1, :) = real(x, c_float) p(1:3, 2, :) = real(x - real(p(1:3, 1, :), wp), c_float) end subroutine split_points