fordot_options.f90 Source File

module fordot_opts

   use kinds

   implicit none

   private
   public :: dot_opts

   interface dot_opts
      procedure :: dot_R0R1R1_rel_opts
   end interface

contains

   !> author: Seyed Ali Ghasemi
   pure function dot_R0R1R1_rel_opts(u,v, option) result(a)
      real(rk), dimension(:), intent(in), contiguous  :: u, v
      character(*),           intent(in)              :: option
      real(rk)                                        :: a

      select case (option)
       case ('m1')
         a = dot_product(u,v)
       case ('m2')
         a = md_2(u,v)
      case ('m3')
         a = md_3(u,v)
      case ('m4')
         a = md_4(u,v)
       case default
         a = dot_product(u,v)
      end select
   end function dot_R0R1R1_rel_opts

   !> author: Seyed Ali Ghasemi
   pure function md_2(u,v) result(a)
      use external_interfaces_dot
      real(rk), dimension(:), intent(in), contiguous :: u, v
      real(rk)                                       :: a

      a = dot(size(u),u,1,v,1)
   end function md_2

   !> author: Seyed Ali Ghasemi
   pure function md_3(u,v) result(a)
      real(rk), dimension(:), intent(in), contiguous :: u, v
      real(rk) :: a
      integer :: i

      a = 0.0_rk
      do i = 1, size(u)
         a = a + u(i)*v(i)
      end do
   end function md_3

   !> author: Seyed Ali Ghasemi
   pure function md_4(u,v) result(a)
      real(rk), dimension(:), intent(in), contiguous :: u, v
      real(rk) :: a
      a = sum(u*v)
   end function md_4

end module fordot_opts