runmat-meshing-core 0.6.0

Deterministic meshing preparation contracts for RunMat
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
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

use super::*;

#[test]
fn rejects_missing_boundary_edge_recovery_when_required() {
    let mesh = valid_tetrahedron_mesh();
    let err = validate_analysis_mesh_with_options(
        &mesh,
        AnalysisMeshValidationOptions {
            min_boundary_edge_recovery_ratio: 1.0,
            ..AnalysisMeshValidationOptions::default()
        },
    )
    .expect_err("missing boundary edge recovery should fail");
    assert_eq!(
        err,
        AnalysisMeshValidationError::BoundaryEdgeRecoveryFailed {
            recovery_ratio: "0.000000".to_string(),
            required_ratio: "1.000000".to_string(),
        }
    );
}

#[test]
fn rejects_unrecovered_boundary_faces_when_required() {
    let mut mesh = valid_tetrahedron_mesh();
    mesh.boundary_faces[0].adjacent_volume_element_ids.clear();
    let err = validate_analysis_mesh_with_options(
        &mesh,
        AnalysisMeshValidationOptions {
            min_boundary_face_recovery_ratio: 1.0,
            ..AnalysisMeshValidationOptions::default()
        },
    )
    .expect_err("missing boundary recovery should fail");
    assert_eq!(
        err,
        AnalysisMeshValidationError::BoundaryFaceRecoveryFailed {
            recovery_ratio: "0.000000".to_string(),
            required_ratio: "1.000000".to_string(),
        }
    );
}

#[test]
fn rejects_unrecovered_tetrahedron_components_recovery_when_policy_requires_strict_recovery() {
    let mut mesh = valid_tetrahedron_mesh();
    mesh.backend.tetrahedron_unrecovered_component_count = 1;

    let err = validate_analysis_mesh_with_options(
        &mesh,
        AnalysisMeshValidationOptions {
            require_no_unrecovered_tetrahedron_components: true,
            ..AnalysisMeshValidationOptions::default()
        },
    )
    .expect_err("strict recovery policy should reject unrecovered Tetrahedron components evidence");

    assert_eq!(
        err,
        AnalysisMeshValidationError::UnrecoveredTetrahedronComponentsPresent { component_count: 1 }
    );
    assert_eq!(
        analysis_mesh_validation_error_code(&err),
        "unrecovered_tetrahedron_components_present"
    );
}

#[test]
fn rejects_rolled_back_material_interface_partition_recovery() {
    let mut mesh = valid_tetrahedron_mesh();
    mesh.backend
        .tetrahedron_rolled_back_absent_material_partition_recovery_item_count = 1;
    mesh.backend
        .tetrahedron_rolled_back_absent_material_partition_element_count = 2;
    mesh.backend
        .tetrahedron_rolled_back_absent_material_partition_boundary_face_count = 3;

    let err = validate_analysis_mesh_with_options(&mesh, AnalysisMeshValidationOptions::default())
        .expect_err("rolled-back material-interface partition recovery should fail readiness");

    assert_eq!(
        err,
        AnalysisMeshValidationError::RolledBackMaterialInterfacePartitionRecoveryPresent {
            recovery_item_count: 1,
            element_count: 2,
            boundary_face_count: 3,
            post_insertion_audit_rejection_count: 0,
        }
    );
    assert_eq!(
        analysis_mesh_validation_error_code(&err),
        "rolled_back_material_interface_partition_recovery_present"
    );
}

#[test]
fn rejects_material_interface_partition_post_insertion_audit_rejection() {
    let mut mesh = valid_tetrahedron_mesh();
    mesh.backend
        .tetrahedron_rejected_absent_material_partition_post_insertion_audit_count = 1;

    let err = validate_analysis_mesh(&mesh, Default::default())
        .expect_err("post-insertion audit rejection should fail readiness");

    assert_eq!(
        err,
        AnalysisMeshValidationError::RolledBackMaterialInterfacePartitionRecoveryPresent {
            recovery_item_count: 0,
            element_count: 0,
            boundary_face_count: 0,
            post_insertion_audit_rejection_count: 1,
        }
    );
}