gemlab 2.0.0

Geometry and meshes laboratory for finite element analyses
Documentation 
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

use gemlab::mesh::Mesh;
use gemlab::mesh::Unstructured;
use gemlab::StrError;
use std::path::Path;
use std::path::PathBuf;
use structopt::StructOpt;

// TODO: Finish msh2tet implementation

const OUT_DIR: &str = "/tmp/gemlab";

/// Command line options
#[derive(Debug, StructOpt)]
#[structopt(name = "msh2tet", about = "Generates tetrahedral meshes from a MSH file.")]
struct Options {
    /// Input file
    #[structopt(parse(from_os_str))]
    input: PathBuf,

    /// Point defining the region
    #[structopt(
        long,
        value_names = &["x", "y", "z"],
        number_of_values = 3,
        allow_hyphen_values = true
    )]
    region: Vec<f64>,

    /// Generate Tet10 instead of Tet4 elements
    #[structopt(short, long)]
    tet10: bool,

    /// Apply node renumbering to reduce the matrix bandwidth
    #[structopt(short, long)]
    renumber: bool,

    #[structopt(long)]
    global_max_volume: Option<f64>,
}

fn main() -> Result<(), StrError> {
    // parse options
    let options = Options::from_args();

    // extract file stem
    let in_path = Path::new(&options.input);
    let fn_stem = in_path
        .file_stem()
        .ok_or("cannot get file stem")?
        .to_str()
        .ok_or("cannot convert file stem to str")?;

    // read MSH file as PLC
    let plc = Mesh::read(in_path)?;

    // data for Trigen
    let regions = vec![(1, options.region[0], options.region[1], options.region[2])];
    let holes = Vec::new();
    let o2 = options.tet10;
    let max_volumes = None;
    let global_min_angle = None;
    let renumber = options.renumber;

    // call Trigen
    let mesh = Unstructured::call_tetgen(
        &plc,
        &regions,
        &holes,
        o2,
        max_volumes,
        options.global_max_volume,
        global_min_angle,
        renumber,
    )?;

    // save MSH file
    let path_out = format!("{}/{}.msh", OUT_DIR, fn_stem);
    mesh.write(&path_out)?;
    println!("Generated MSH file: {}", path_out);
    Ok(())
}