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
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
open Bigarray
include Mesh_common
let pslg (type l) ?(hole: l mat option) ?(region: l mat option)
?(point_marker: l int_vec option) (point: l mat)
?(segment_marker: l int_vec option) (segment: l int_mat) : l pslg =
match Array2.layout point with
| C_layout ->
MeshC.pslg ~hole ~region ~point_marker ~point ~segment_marker ~segment
| Fortran_layout ->
MeshF.pslg ~hole ~region ~point_marker ~point ~segment_marker ~segment
let create (type l) ?(hole: l mat option) ?(region: l mat option)
?(point_marker: l int_vec option) (point:l mat)
?(segment_marker: l int_vec option) ?(segment: l int_mat option)
?(neighbor: l int_mat option)
?(edge: l int_mat option) ?(edge_marker: l int_vec option)
(triangle: l int_mat) : l t =
match Array2.layout point with
| C_layout ->
MeshC.create ~hole ~region ~point_marker ~point
~segment_marker ~segment ~neighbor ~edge ~edge_marker ~triangle
| Fortran_layout ->
MeshF.create ~hole ~region ~point_marker ~point
~segment_marker ~segment ~neighbor ~edge ~edge_marker ~triangle
let copy_vec v =
let v' = Array1.create (Array1.kind v) (Array1.layout v) (Array1.dim v) in
Array1.blit v v';
v'
let copy_mat m =
let m' = Array2.create (Array2.kind m) (Array2.layout m)
(Array2.dim1 m) (Array2.dim2 m) in
Array2.blit m m';
m'
let copy (mesh: 'l t) =
make_mesh
~point: (copy_mat mesh#point)
~point_marker: (copy_vec mesh#point_marker)
~segment: (copy_mat mesh#segment)
~segment_marker: (copy_vec mesh#segment_marker)
~hole: (copy_mat mesh#hole)
~region: (copy_mat mesh#region)
~triangle: (copy_mat mesh#triangle)
~neighbor: (copy_mat mesh#neighbor)
~edge: (copy_mat mesh#edge)
~edge_marker: (copy_vec mesh#edge_marker)
let sub (type l) (mesh: l t) ?pos len : l t =
match layout mesh with
| C_layout -> MeshC.sub mesh ?pos len
| Fortran_layout -> MeshF.sub mesh ?pos len
let band_height_P1 (type l) ?filter (mesh: l #t) =
let mesh = (mesh :> l t) in
match layout mesh with
| C_layout -> MeshC.band_height_P1 filter mesh
| Fortran_layout -> MeshF.band_height_P1 filter mesh
let cuthill_mckee (type l) ?(rev=true) ?(perm: l int_vec option)
(mesh: l t) : l t =
match layout mesh with
| C_layout -> MeshC.cuthill_mckee ~rev perm mesh
| Fortran_layout -> MeshF.cuthill_mckee ~rev perm mesh
let permute_points (type l) (mesh: l t) ?(inv=false) (perm: l int_vec) : l t =
match layout mesh with
| C_layout -> MeshC.permute_points mesh ~inv perm
| Fortran_layout -> MeshF.permute_points mesh ~inv perm
let permute_triangles (type l) (mesh: l t) ?(inv=false)
(perm: l int_vec) : l t =
match layout mesh with
| C_layout -> MeshC.permute_triangles mesh ~inv perm
| Fortran_layout -> MeshF.permute_triangles mesh ~inv perm
module LaTeX =
struct
type color = int
let save (type l) ?edge (mesh: l #t) filename =
let mesh = (mesh :> l t) in
match layout mesh with
| C_layout -> MeshC.latex ?edge mesh filename
| Fortran_layout -> MeshF.latex ?edge mesh filename
let write (type l) ?edge (mesh: l #t) fh =
let mesh = (mesh :> l t) in
match layout mesh with
| C_layout -> MeshC.latex_write ?edge mesh fh
| Fortran_layout -> MeshF.latex_write ?edge mesh fh
let level_curves (type l) ?boundary (mesh: l #t) (z: l vec)
?level_eq levels filename =
let mesh = (mesh :> l t) in
match layout mesh with
| C_layout -> MeshC.level_curves ?boundary mesh z
?level_eq levels filename
| Fortran_layout -> MeshF.level_curves ?boundary mesh z
?level_eq levels filename
let super_level (type l) ?boundary (mesh: l #t) (z: l vec) level color
filename =
let mesh = (mesh :> l t) in
match layout mesh with
| C_layout ->
MeshC.super_level ?boundary mesh z level color filename
| Fortran_layout ->
MeshF.super_level ?boundary mesh z level color filename
let sub_level (type l) ?boundary (mesh: l #t) (z: l vec) level color
filename =
let mesh = (mesh :> l t) in
match layout mesh with
| C_layout -> MeshC.sub_level ?boundary mesh z level color filename
| Fortran_layout -> MeshF.sub_level ?boundary mesh z level color filename
end
let scilab (type l) (mesh: l #t) ?longitude ?azimuth ?mode ?box ?edgecolor
(z: l vec) filename =
let mesh = (mesh :> l t) in
match layout mesh with
| C_layout -> MeshC.scilab mesh ?longitude ?azimuth ?mode ?box
?edgecolor z filename
| Fortran_layout -> MeshF.scilab mesh ?longitude ?azimuth ?mode ?box
?edgecolor z filename
let matlab (type l) (mesh: l #t) ?edgecolor ?linestyle ?facealpha
(z: l vec) filename =
let mesh = (mesh :> l t) in
match layout mesh with
| C_layout -> MeshC.matlab mesh ?edgecolor ?linestyle ?facealpha
z filename
| Fortran_layout -> MeshF.matlab mesh ?edgecolor ?linestyle ?facealpha
z filename
let mathematica (type l) (mesh: l #t) (z: l vec) filename =
let mesh = (mesh :> l t) in
match layout mesh with
| C_layout -> MeshC.mathematica mesh z filename
| Fortran_layout -> MeshF.mathematica mesh z filename