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
170
171
172
//! The `OpenCASCADE` geometry kernel
use super::*;
use crate::{GmshError, GmshResult, check_main_error, check_model_error};
use crate::interface::occ as factory;
/// All angle values are in radians, commonly given as fractions of π.
include!("common_occ.rs");
impl<'a> OccModel<'a> {
/// Add a box with a starting point and side lengths from that point.
#[must_use]
pub fn add_box(&mut self, start_point: (f64, f64, f64),
extents: (f64, f64, f64)) -> GmshResult<VolumeTag> {
self.set_current()?;
let mut ierr: c_int = 0;
let automatic_tag: c_int = -1;
unsafe {
let out_tag =
factory::add_box(start_point.0,
start_point.1,
start_point.2,
extents.0,
extents.1,
extents.2,
automatic_tag,
&mut ierr);
check_model_error!(ierr, VolumeTag(out_tag))
}
}
/// Add a sphere with a centroid and radius.
#[must_use]
pub fn add_sphere(&mut self, centroid: (f64, f64, f64), radius: f64) -> GmshResult<VolumeTag> {
let polar_min = -std::f64::consts::FRAC_PI_2;
let polar_max = std::f64::consts::FRAC_PI_2;
let azimuth = 2. * std::f64::consts::PI;
self.add_sphere_gen(centroid, radius, (polar_min, polar_max), azimuth)
}
/// Add an angular section of a sphere.
/// The polar angles are `(min, max)` relative to the x-axis on the xy-plane,
/// and the azimuth is relative to the z-axis.
/// ```
/// # use rgmsh::{Gmsh, GmshResult, add_points};
/// # use rgmsh::model::occ::*;
/// # fn main() -> GmshResult<()> {
/// # let gmsh = Gmsh::initialize()?;
/// # let mut geom = gmsh.create_occ_model("model")?;
/// use std::f64::consts;
///
/// // make a 1/8th sphere section
/// let sphere_section = geom.add_sphere_section(
/// (0., 0., 0.), 1.,
/// (0., consts::FRAC_PI_2),
/// consts::FRAC_PI_2)?;
///
/// // recover a basic sphere by setting polar = (-π/2, π/2) and azimuth = 2π.
/// let basic_sphere = geom.add_sphere_section(
/// (10., 10., 10.), 1.,
/// (-consts::FRAC_PI_2, consts::FRAC_PI_2),
/// 2. * consts::PI)?;
/// # Ok(())
/// # }
/// ```
#[must_use]
pub fn add_sphere_section(&mut self, centroid: (f64, f64, f64), radius: f64, polar: (f64, f64), azimuth: f64 // fx: f64, y: f64, z: f64,
) -> GmshResult<VolumeTag> {
self.add_sphere_gen(centroid, radius, polar, azimuth)
}
#[doc(hidden)]
#[must_use]
fn add_sphere_gen(
&mut self,
centroid: (f64, f64, f64),
radius: f64,
polar: (f64, f64),
azimuth: f64,
) -> GmshResult<VolumeTag> {
self.set_current()?;
unsafe {
let mut ierr: c_int = 0;
let automatic_tag: c_int = -1;
let out_tag = factory::add_sphere(
centroid.0,
centroid.1,
centroid.2,
radius,
automatic_tag,
polar.0,
polar.1,
azimuth,
&mut ierr,
);
check_model_error!(ierr, VolumeTag(out_tag))
}
}
/// Add a torus to the model.
/// ```
/// # use rgmsh::{Gmsh, GmshResult};
/// # use rgmsh::model::shapes::*;
/// # use rgmsh::model::occ::*;
/// # fn main() -> GmshResult<()> {
/// # let gmsh = Gmsh::initialize()?;
/// # let mut geom = gmsh.create_occ_model("model")?;
/// let centroid = Point { x: 1.0, y: 1.0, z: 1.0 };
/// let torus = geom.add_torus_experimental(
/// Torus {
/// centroid,
/// main_radius: 3.0,
/// pipe_radius: 0.2
/// })?;
/// # Ok(())
/// # }
/// ```
#[must_use]
pub fn add_torus_experimental(&mut self, torus: Torus) -> GmshResult<VolumeTag> {
self.add_torus((torus.centroid.x, torus.centroid.y, torus.centroid.z),
(torus.main_radius, torus.pipe_radius))
}
/// Add a torus with a centroid and radii values `(main_radius, pipe_radius)`.
#[must_use]
pub fn add_torus(
&mut self, centroid: (f64, f64, f64),
radii: (f64, f64)
) -> GmshResult<VolumeTag> {
let angle = 2. * std::f64::consts::PI;
self.add_torus_gen(centroid, radii, angle)
}
/// Add an angular torus section using the main torus radius. To recover
/// a basic torus, set the `angle` parameter to 2π.
#[must_use]
pub fn add_torus_section(
&mut self, centroid: (f64, f64, f64), radii: (f64, f64), angle: f64
) -> GmshResult<VolumeTag> {
self.add_torus_gen(centroid, radii, angle)
}
#[doc(hidden)]
#[must_use]
fn add_torus_gen(
&mut self,
centroid: (f64, f64, f64),
radii: (f64, f64),
angle: f64
) -> GmshResult<VolumeTag> {
self.set_current()?;
unsafe {
let mut ierr: c_int = 0;
let automatic_tag: c_int = -1;
let out_tag = factory::add_torus(
centroid.0,
centroid.1,
centroid.2,
radii.0,
radii.1,
automatic_tag,
angle,
&mut ierr,
);
check_model_error!(ierr, VolumeTag(out_tag))
}
}
}