rsspice 0.1.0

Pure Rust port of the SPICE Toolkit for space geometry
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
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

//
// GENERATED FILE
//

use super::*;
use f2rust_std::*;

const TIGHT: f64 = 0.000000000001;
const MAXDEG: i32 = 15;

//$Procedure F_CHBSHV ( Chebyshev expansion tests )
pub fn F_CHBSHV(OK: &mut bool, ctx: &mut Context) -> f2rust_std::Result<()> {
    let mut CP = StackArray::<f64, 15>::new(1..=MAXDEG);
    let mut DPDX: f64 = 0.0;
    let mut DPDXS = StackArray::<f64, 16>::new(0..=MAXDEG);
    let mut P: f64 = 0.0;
    let mut PARTDP = StackArray2D::<f64, 48>::new(1..=3, 0..=MAXDEG);
    let mut X: f64 = 0.0;
    let mut X2S = StackArray::<f64, 2>::new(1..=2);
    let mut DEGP: i32 = 0;
    let mut NDERIV: i32 = 0;

    //
    // SPICELIB functions
    //

    //
    // Local Parameters
    //

    //
    // Local Variables
    //

    //
    // Open the test family.
    //
    testutil::TOPEN(b"F_CHBSHV", ctx)?;

    //
    // *****************************************************************
    //
    // Error cases
    //
    // *****************************************************************
    //
    //
    // None of the routines participate in error handling.
    //
    //
    // *****************************************************************
    //
    // Normal cases
    //
    // *****************************************************************
    //
    //
    // --- Case: -------------------------------------------------------
    //
    testutil::TCASE(
        b"CHBDER normal case #1: find value of polynomial and its first 3 derivatives at X=1",
        ctx,
    )?;

    CP[1] = 1.0;
    CP[2] = 3.0;
    CP[3] = 0.5;
    CP[4] = 1.0;
    CP[5] = 0.5;
    CP[6] = -1.0;
    CP[7] = 1.0;

    X2S[1] = 0.5;
    X2S[2] = 3.0;

    DEGP = 6;
    NDERIV = 3;
    X = 1.0;

    spicelib::CHBDER(
        CP.as_slice(),
        DEGP,
        X2S.as_slice(),
        X,
        NDERIV,
        PARTDP.as_slice_mut(),
        DPDXS.as_slice_mut(),
    );
    testutil::CHCKXC(false, b" ", OK, ctx)?;

    testutil::CHCKSD(b"f(1)", DPDXS[0], b"~/", -0.34087791495199, TIGHT, OK, ctx)?;
    testutil::CHCKSD(
        b"df\'(1)",
        DPDXS[1],
        b"~/",
        0.38271604938272,
        TIGHT,
        OK,
        ctx,
    )?;
    testutil::CHCKSD(b"df\"(1)", DPDXS[2], b"~/", 4.2880658436214, TIGHT, OK, ctx)?;
    testutil::CHCKSD(
        b"df\'\"(1)",
        DPDXS[3],
        b"~/",
        -1.5144032921811,
        TIGHT,
        OK,
        ctx,
    )?;

    //
    // --- Case: -------------------------------------------------------
    //
    testutil::TCASE(
        b"CHBINT normal case #1: find value of polynomial and its first derivative at X=1",
        ctx,
    )?;

    CP[1] = 1.0;
    CP[2] = 3.0;
    CP[3] = 0.5;
    CP[4] = 1.0;
    CP[5] = 0.5;
    CP[6] = -1.0;
    CP[7] = 1.0;

    X2S[1] = 0.5;
    X2S[2] = 3.0;

    DEGP = 6;
    X = 1.0;

    spicelib::CHBINT(CP.as_slice(), DEGP, X2S.as_slice(), X, &mut P, &mut DPDX);
    testutil::CHCKXC(false, b" ", OK, ctx)?;

    testutil::CHCKSD(b"f(1)", P, b"~/", -0.34087791495199, TIGHT, OK, ctx)?;
    testutil::CHCKSD(b"df\'(1)", DPDX, b"~/", 0.38271604938272, TIGHT, OK, ctx)?;

    //
    // --- Case: -------------------------------------------------------
    //
    testutil::TCASE(
        b"CHBVAL normal case #1: find value of polynomial at X=1",
        ctx,
    )?;

    CP[1] = 1.0;
    CP[2] = 3.0;
    CP[3] = 0.5;
    CP[4] = 1.0;
    CP[5] = 0.5;
    CP[6] = -1.0;
    CP[7] = 1.0;

    X2S[1] = 0.5;
    X2S[2] = 3.0;

    DEGP = 6;
    X = 1.0;

    spicelib::CHBVAL(CP.as_slice(), DEGP, X2S.as_slice(), X, &mut P);
    testutil::CHCKXC(false, b" ", OK, ctx)?;

    testutil::CHCKSD(b"f(1)", P, b"~/", -0.34087791495199, TIGHT, OK, ctx)?;

    //
    // Close out the test family.
    //
    testutil::T_SUCCESS(OK, ctx);
    Ok(())
}