scalpr Subroutine

public subroutine scalpr(x, y, z)

Calculates the scalar product between two vectors

Note

Jonathan Chico: This can be generalized to vectors of any dimmension by using the call sum(x(:)*y(:)). This should be more efficient for ifort as it vectorizes the sum call automatically

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in), dimension(*) :: x
real(kind=dp), intent(in), dimension(*) :: y
real(kind=dp), intent(out) :: z