Line data    Source code

       1              : module test_eos_blend
       2              : 
       3              :    use const_def, only: dp
       4              :    use auto_diff
       5              :    use eos_blend
       6              : 
       7              :    implicit none
       8              : 
       9              :    private
      10              :    public :: do_test_eos_blend
      11              : 
      12              : contains
      13              : 
      14            1 :    subroutine do_test_eos_blend()
      15           11 :       real(dp) :: coords(4, 2)
      16              :       type(auto_diff_real_2var_order1) :: p(2), d
      17              : 
      18            1 :       coords(1, 1) = -1d0
      19            1 :       coords(1, 2) = -1d0
      20            1 :       coords(2, 1) = -1d0
      21            1 :       coords(2, 2) = 1d0
      22            1 :       coords(3, 1) = 1d0
      23            1 :       coords(3, 2) = 1d0
      24            1 :       coords(4, 1) = 1d0
      25            1 :       coords(4, 2) = -1d0
      26              : 
      27            1 :       p(1) = 0d0
      28            1 :       p(2) = 1d0
      29            1 :       write (*, *) is_contained(4, coords, p)
      30            1 :       d = min_distance_to_polygon(4, coords, p)
      31            1 :       write (*, '(99(1pd26.16))') d
      32            1 :       p(1) = -1d0
      33            1 :       p(2) = 1d0
      34            1 :       write (*, *) is_contained(4, coords, p)
      35            1 :       d = min_distance_to_polygon(4, coords, p)
      36            1 :       write (*, '(99(1pd26.16))') d
      37            1 :       p(1) = 0d0
      38            1 :       p(2) = 0d0
      39            1 :       write (*, *) is_contained(4, coords, p)
      40            1 :       d = min_distance_to_polygon(4, coords, p)
      41            1 :       write (*, '(99(1pd26.16))') d
      42            1 :       p(1) = -2d0
      43            1 :       p(2) = 0d0
      44            1 :       write (*, *) is_contained(4, coords, p)
      45            1 :       d = min_distance_to_polygon(4, coords, p)
      46            1 :       write (*, '(99(1pd26.16))') d
      47            1 :       p(1) = -1d0
      48            1 :       p(2) = 3d0
      49            1 :       write (*, *) is_contained(4, coords, p)
      50            1 :       d = min_distance_to_polygon(4, coords, p)
      51            1 :       write (*, '(99(1pd26.16))') d
      52              : 
      53            1 :    end subroutine do_test_eos_blend
      54              : 
      55              : end module test_eos_blend