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
//
// GENERATED FILE
//
use super::*;
use crate::SpiceContext;
use f2rust_std::*;
/// Approximate equality
///
/// Return .TRUE. if two double precision numbers are equal to
/// within some tolerance.
///
/// # Brief I/O
///
/// ```text
/// VARIABLE I/O DESCRIPTION
/// -------- --- --------------------------------------------------
/// X,
/// Y I Double precision numbers.
/// TOL I Tolerance.
///
/// The function is .TRUE. whenever |X - Y| <= TOL.
/// ```
///
/// # Detailed Input
///
/// ```text
/// X,
/// Y are arbitrary double precision numbers.
///
/// TOL is a tolerance. X and Y are considered to be equal
/// if they differ by no more than this amount. If TOL
/// is negative, X and Y are never considered equal.
/// ```
///
/// # Detailed Output
///
/// ```text
/// The function is .TRUE. whenever |X - Y| <= TOL, and is .FALSE.
/// otherwise.
/// ```
///
/// # Exceptions
///
/// ```text
/// Error free.
/// ```
///
/// # Examples
///
/// ```text
/// C
/// C If the eccentricity is near one, this a parabola.
/// C
/// IF ( APPROX ( ECC, 1.D0, 10.D-12 ) ) THEN
/// TYPE = 'PARABOLA'
///
/// ELSE IF ( ECC .LT. 1 ) THEN
/// TYPE = 'ELLIPSE'
///
/// ELSE
/// TYPE = 'HYPERBOLA'
/// END IF
/// ```
///
/// # Author and Institution
///
/// ```text
/// J. Diaz del Rio (ODC Space)
/// W.L. Taber (JPL)
/// I.M. Underwood (JPL)
/// ```
///
/// # Version
///
/// ```text
/// - SPICELIB Version 1.1.0, 26-OCT-2021 (JDR)
///
/// Added IMPLICIT NONE statement.
///
/// Edited the header to comply with NAIF standard.
///
/// - SPICELIB Version 1.0.1, 10-MAR-1992 (WLT)
///
/// Comment section for permuted index source lines was added
/// following the header.
///
/// - SPICELIB Version 1.0.0, 31-JAN-1990 (IMU)
/// ```
pub fn approx(x: f64, y: f64, tol: f64) -> bool {
let ret = APPROX(x, y, tol);
ret
}
//$Procedure APPROX ( Approximate equality )
pub fn APPROX(X: f64, Y: f64, TOL: f64) -> bool {
let mut APPROX: bool = false;
//
// Just shorthand, really.
//
APPROX = (f64::abs((X - Y)) <= TOL);
APPROX
}