plot3d 0.1.13

Utilities for reading, writing, and manipulating NASA PLOT3D structured grids.
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

use plot3d::{
    connectivity_fast, create_rotation_matrix, read_plot3d_ascii, read_plot3d_binary,
    rotate_block_with_matrix, rotated_periodicity, translational_periodicity, FaceMatchPrinter,
    FaceRecordTraits, Float,
};

#[test]
fn rotational_periodicity_test() {
    let url = "https://nasa-public-data.s3.amazonaws.com/plot3d_utilities/VSPT_ASCII.xyz";
    let ascii_path = "VSPT_ASCII.xyz";
    if !std::path::Path::new(ascii_path).exists() {
        let bytes = reqwest::blocking::get(url).unwrap().bytes().unwrap();
        std::fs::write(ascii_path, &bytes).unwrap();
    }

    let blocks = read_plot3d_ascii(ascii_path).unwrap();

    let number_of_blades = 55;
    let rotation_angle = 360.0 / number_of_blades as Float;
    let copies: usize = 3;

    let mut rotated_blocks = Vec::with_capacity(blocks.len() * copies);
    rotated_blocks.extend(blocks.iter().cloned());

    for copy_idx in 1..copies {
        let angle_rad = (rotation_angle * copy_idx as Float).to_radians();
        let rotation_matrix = create_rotation_matrix(angle_rad, 'x');
        for block in &blocks {
            rotated_blocks.push(rotate_block_with_matrix(block, rotation_matrix));
        }
    }

    assert_eq!(rotated_blocks.len(), blocks.len() * copies);

    let (face_matches, outer_faces) = connectivity_fast(&blocks);

    // face_matches.print();

    let (periodic_faces, outer_faces_rotated) = rotated_periodicity(
        &blocks,
        &face_matches,
        &outer_faces,
        rotation_angle,
        'x',
        true,
    );
    println!("Printing Periodic faces");
    periodic_faces.print();
    outer_faces_rotated.print();
}

#[test]
fn translational_periodicity_test() {
    let url = "https://nasa-public-data.s3.amazonaws.com/plot3d_utilities/iso65_64blocks.xyz";
    let file_path = "iso65_64blocks.xyz";
    if !std::path::Path::new(file_path).exists() {
        let bytes = reqwest::blocking::get(url).unwrap().bytes().unwrap();
        std::fs::write(file_path, &bytes).unwrap();
    }

    let blocks = read_plot3d_binary(
        file_path,
        plot3d::BinaryFormat::Raw,
        plot3d::FloatPrecision::F32,
        plot3d::Endian::Little,
    )
    .unwrap();

    let (face_matches, outer_faces) = connectivity_fast(&blocks);

    let mut remaining_faces = outer_faces.clone();

    let (x_periodicity, rest) =
        translational_periodicity(&blocks, &remaining_faces, None, "x", None, 0.5, 1, 1, 1);
    remaining_faces = rest;

    let (y_periodicity, rest) =
        translational_periodicity(&blocks, &remaining_faces, None, "y", None, 0.5, 1, 1, 1);
    remaining_faces = rest;

    let (z_periodicity, _rest) =
        translational_periodicity(&blocks, &remaining_faces, None, "z", None, 0.5, 1, 1, 1);
    let mut face_matches_all = face_matches;
    face_matches_all.extend(x_periodicity);
    face_matches_all.extend(y_periodicity);
    face_matches_all.extend(z_periodicity);
}