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
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
//
// GENERATED FILE
//
use super::*;
use crate::SpiceContext;
use f2rust_std::*;
struct SaveVars {
SVSTEP: f64,
SVINIT: bool,
}
impl SaveInit for SaveVars {
fn new() -> Self {
let mut SVSTEP: f64 = 0.0;
let mut SVINIT: bool = false;
SVINIT = false;
SVSTEP = -1.0;
Self { SVSTEP, SVINIT }
}
}
/// GF, step size
///
/// Return the time step set by the most recent call to GFSSTP.
///
/// # Required Reading
///
/// * [GF](crate::required_reading::gf)
/// * [TIME](crate::required_reading::time)
///
/// # Brief I/O
///
/// ```text
/// VARIABLE I/O DESCRIPTION
/// -------- --- --------------------------------------------------
/// TIME I Ignored ET value.
/// STEP O Time step to take.
/// ```
///
/// # Detailed Input
///
/// ```text
/// TIME is an ignored double precision number. This argument
/// is present so the argument list of this routine is
/// compatible with the GF step size routine argument list
/// specification.
///
/// When this routine is called from within the GF
/// root-finding system, either the initial ET value of the
/// current interval of the confinement window, or the
/// value resulting from the last search step, is passed in
/// via the TIME argument.
/// ```
///
/// # Detailed Output
///
/// ```text
/// STEP is the output step size. This is the value set by the
/// most recent call to GFSSTP. Units are TDB seconds.
///
/// STEP is used in the GF search root-bracketing process.
/// STEP indicates how far to advance TIME so that TIME and
/// TIME+STEP may bracket a state transition and definitely
/// do not bracket more than one state transition.
/// ```
///
/// # Exceptions
///
/// ```text
/// 1) If this routine is called before a step size has been
/// set via a call to GFSSTP, the error SPICE(NOTINITIALIZED)
/// is signaled.
/// ```
///
/// # Particulars
///
/// ```text
/// This routine returns the time step set by the most recent call to
/// GFSSTP.
/// ```
///
/// # Examples
///
/// ```text
/// 1) In normal usage of a high-level GF API routine, the caller
/// will pass in a constant step size STEP. The API routine will
/// then make the call
///
/// CALL GFSSTP ( STEP )
///
/// Subsequent calls to GFSTEP during the search process conducted
/// by the API routine will return STEP.
///
///
/// 2) User applications can pass GFSTEP to mid-level GF API routines
/// expecting a step size routine as an input argument. For
/// example, the GF API routine GFOCCE can be called as follows:
///
///
/// Set the step size.
///
/// CALL GFSSTP ( STEP )
///
///
/// Look for solutions. (GFSTEP is the 11th argument.)
///
/// CALL GFOCCE ( OCCTYP, FRONT, FSHAPE, FFRAME,
/// . BACK, BSHAPE, BFRAME, ABCORR,
/// . OBSRVR, CNVTOL, GFSTEP, GFREFN,
/// . RPT, GFREPI, GFREPU, GFREPF,
/// . BAIL, GFBAIL, CNFINE, RESULT )
/// ```
///
/// # Author and Institution
///
/// ```text
/// N.J. Bachman (JPL)
/// J. Diaz del Rio (ODC Space)
/// L.S. Elson (JPL)
/// W.L. Taber (JPL)
/// I.M. Underwood (JPL)
/// E.D. Wright (JPL)
/// ```
///
/// # Version
///
/// ```text
/// - SPICELIB Version 1.1.1, 03-JUN-2021 (JDR)
///
/// Edited the header to comply with NAIF standard.
///
/// - SPICELIB Version 1.1.0, 31-AUG-2010 (EDW)
///
/// Expanded error message on STEP for clarity.
///
/// Added TIME = TIME declaration to eliminate unused dummy
/// variable warning during compilation.
///
/// - SPICELIB Version 1.0.0, 05-MAR-2009 (NJB) (LSE) (IMU) (WLT) (EDW)
/// ```
pub fn gfstep(ctx: &mut SpiceContext, time: &mut f64, step: &mut f64) -> crate::Result<()> {
GFSTEP(time, step, ctx.raw_context())?;
ctx.handle_errors()?;
Ok(())
}
//$Procedure GFSTEP ( GF, step size )
pub fn GFSTEP(TIME: &mut f64, STEP: &mut f64, ctx: &mut Context) -> f2rust_std::Result<()> {
let save = ctx.get_vars::<SaveVars>();
let save = &mut *save.borrow_mut();
//
// SPICELIB functions
//
//
// Local variables
//
//
// Saved variables
//
//
// Initial values
//
//
// Discovery check-in.
//
if !save.SVINIT {
CHKIN(b"GFSTEP", ctx)?;
SETMSG(b"Step size was never initialized.", ctx);
SIGERR(b"SPICE(NOTINITIALIZED)", ctx)?;
CHKOUT(b"GFSTEP", ctx)?;
return Ok(());
}
//
// Set STEP to the saved value from the last call to GFSSTP.
//
*STEP = save.SVSTEP;
*TIME = *TIME;
Ok(())
}
/// Geometry finder set step size
///
/// Set the step size to be returned by GFSTEP.
///
/// # Required Reading
///
/// * [GF](crate::required_reading::gf)
/// * [TIME](crate::required_reading::time)
///
/// # Brief I/O
///
/// ```text
/// VARIABLE I/O DESCRIPTION
/// -------- --- --------------------------------------------------
/// STEP I Time step to take.
/// ```
///
/// # Detailed Input
///
/// ```text
/// STEP is the output step size to be returned by the next call
/// GFSTEP. Units are TDB seconds.
///
/// STEP is used in the GF search root-bracketing process.
/// STEP indicates how far to advance TIME so that TIME and
/// TIME+STEP may bracket a state transition and definitely
/// do not bracket more than one state transition.
/// ```
///
/// # Exceptions
///
/// ```text
/// 1) If the input step size is non-positive, the error
/// SPICE(INVALIDSTEP) is signaled. The stored step value
/// is not updated.
/// ```
///
/// # Particulars
///
/// ```text
/// See the header of GFSTEP above.
/// ```
///
/// # Examples
///
/// ```text
/// See the header of GFSTEP above.
/// ```
///
/// # Restrictions
///
/// ```text
/// 1) This routine must be called before the first time
/// GFSTEP is called during a program run.
/// ```
///
/// # Author and Institution
///
/// ```text
/// N.J. Bachman (JPL)
/// J. Diaz del Rio (ODC Space)
/// L.S. Elson (JPL)
/// E.D. Wright (JPL)
/// ```
///
/// # Version
///
/// ```text
/// - SPICELIB Version 1.1.1, 03-JUN-2021 (JDR)
///
/// Edited the header to comply with NAIF standard.
///
/// - SPICELIB Version 1.1.0, 31-AUG-2010 (EDW)
///
/// Expanded error message on STEP for clarity.
///
/// - SPICELIB Version 1.0.0, 15-APR-2009 (LSE) (NJB)
/// ```
pub fn gfsstp(ctx: &mut SpiceContext, step: f64) -> crate::Result<()> {
GFSSTP(step, ctx.raw_context())?;
ctx.handle_errors()?;
Ok(())
}
//$Procedure GFSSTP ( Geometry finder set step size )
pub fn GFSSTP(STEP: f64, ctx: &mut Context) -> f2rust_std::Result<()> {
let save = ctx.get_vars::<SaveVars>();
let save = &mut *save.borrow_mut();
//
// Check the step size.
//
if (STEP <= 0.0) {
CHKIN(b"GFSSTP", ctx)?;
SETMSG(b"Step has value #; step size must be positive.", ctx);
ERRDP(b"#", STEP, ctx);
SIGERR(b"SPICE(INVALIDSTEP)", ctx)?;
CHKOUT(b"GFSSTP", ctx)?;
return Ok(());
}
save.SVSTEP = STEP;
save.SVINIT = true;
Ok(())
}