1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
|
! Module with routines for fuzzy comparisons.
! $Header$
#include "cctk.h"
module EHFinder_fuzzy
use EHFinder_Constants
interface fuzzy_eq
module procedure fuzzy_eq_0d, fuzzy_eq_1d, fuzzy_eq_2d, fuzzy_eq_3d
end interface
interface fuzzy_ne
module procedure fuzzy_ne_0d, fuzzy_ne_1d, fuzzy_ne_2d, fuzzy_ne_3d
end interface
contains
function fuzzy_eq_0d ( x, y, eps )
CCTK_REAL, intent(IN) :: x
CCTK_REAL, intent(IN) :: y, eps
logical :: fuzzy_eq_0d
fuzzy_eq_0d = ( ( x .gt. y-eps ) .and. ( x .lt. y+eps ) )
end function fuzzy_eq_0d
function fuzzy_eq_1d ( x, y, eps )
CCTK_REAL, intent(IN), dimension(:) :: x
CCTK_REAL, intent(IN) :: y, eps
logical, dimension(size(x,1)) :: fuzzy_eq_1d
fuzzy_eq_1d = ( ( x .gt. y-eps ) .and. ( x .lt. y+eps ) )
end function fuzzy_eq_1d
function fuzzy_eq_2d ( x, y, eps )
CCTK_REAL, intent(IN), dimension(:,:) :: x
CCTK_REAL, intent(IN) :: y, eps
logical, dimension(size(x,1),size(x,2)) :: fuzzy_eq_2d
fuzzy_eq_2d = ( ( x .gt. y-eps ) .and. ( x .lt. y+eps ) )
end function fuzzy_eq_2d
function fuzzy_eq_3d ( x, y, eps )
CCTK_REAL, intent(IN), dimension(:,:,:) :: x
CCTK_REAL, intent(IN) :: y, eps
logical, dimension(size(x,1),size(x,2),size(x,3)) :: fuzzy_eq_3d
fuzzy_eq_3d = ( ( x .gt. y-eps ) .and. ( x .lt. y+eps ) )
end function fuzzy_eq_3d
function fuzzy_ne_0d ( x, y, eps )
CCTK_REAL, intent(IN) :: x
CCTK_REAL, intent(IN) :: y, eps
logical :: fuzzy_ne_0d
fuzzy_ne_0d = ( ( x .lt. y-eps ) .or. ( x .gt. y+eps ) )
end function fuzzy_ne_0d
function fuzzy_ne_1d ( x, y, eps )
CCTK_REAL, intent(IN), dimension(:) :: x
CCTK_REAL, intent(IN) :: y, eps
logical, dimension(size(x,1)) :: fuzzy_ne_1d
fuzzy_ne_1d = ( ( x .lt. y-eps ) .or. ( x .gt. y+eps ) )
end function fuzzy_ne_1d
function fuzzy_ne_2d ( x, y, eps )
CCTK_REAL, intent(IN), dimension(:,:) :: x
CCTK_REAL, intent(IN) :: y, eps
logical, dimension(size(x,1),size(x,2)) :: fuzzy_ne_2d
fuzzy_ne_2d = ( ( x .lt. y-eps ) .or. ( x .gt. y+eps ) )
end function fuzzy_ne_2d
function fuzzy_ne_3d ( x, y, eps )
CCTK_REAL, intent(IN), dimension(:,:,:) :: x
CCTK_REAL, intent(IN) :: y, eps
logical, dimension(size(x,1),size(x,2),size(x,3)) :: fuzzy_ne_3d
fuzzy_ne_3d = ( ( x .lt. y-eps ) .or. ( x .gt. y+eps ) )
end function fuzzy_ne_3d
end module EHFinder_fuzzy
|
