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