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
|
module unit_ring_m
use semiring_m, only: semiring
implicit none
private
public :: &
unit_ring_only_minus, &
unit_ring_only_negate, &
unit_ring, &
derive_unit_ring_from_minus, &
derive_unit_ring_from_negate
requirement unit_ring_only_minus(T, plus, zero, mult, one, minus)
require :: semiring(T, plus, zero, mult, one)
elemental function minus(x, y) result(difference)
type(T), intent(in) :: x, y
type(T) :: difference
end function
end requirement
requirement unit_ring_only_negate(T, plus, zero, mult, one, negate)
require :: semiring(T, plus, zero, mult, one)
elemental function negate(x) result(negated)
type(T), intent(in) :: x
type(T) :: negated
end function
end requirement
requirement unit_ring(T, plus, zero, mult, one, minus, negate)
require :: unit_ring_only_minus(T, plus, zero, mult, one, minus)
require :: unit_ring_only_negate(T, plus, zero, mult, one, negate)
end requirement
template derive_unit_ring_from_minus(T, plus, zero, mult, one, minus)
require :: unit_ring_only_minus(T, plus, zero, mult, one, minus)
private
public :: negate
contains
elemental function negate(x) result(negated)
type(T), intent(in) :: x
type(T) :: negated
negated = minus(zero(), x)
end function
end template
template derive_unit_ring_from_negate(T, plus, zero, mult, one, negate)
require :: unit_ring_only_negate(T, plus, zero, mult, one, negate)
private
public :: minus
contains
elemental function minus(x, y) result(difference)
type(T), intent(in) :: x, y
type(T) :: difference
difference = plus(x, negate(y))
end function
end template
end module
|
