Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=dp), | intent(in) | :: | u(:) | |||
real(kind=dp), | intent(in) | :: | v(:) |
function inner_product_real(u, v) result(inner_product)
real(dp), intent(in) :: u(:), v(:)
real(dp) :: inner_product, ddot
character(*), parameter :: this_routine = "inner_product_real"
ASSERT(size(u) == size(v))
inner_product = ddot(size(u), u, 1, v, 1)
end function inner_product_real