]> git.lizzy.rs Git - rust.git/blob - library/core/src/num/f64.rs
Add test for eval order for a+=b
[rust.git] / library / core / src / num / f64.rs
1 //! This module provides constants which are specific to the implementation
2 //! of the `f64` floating point data type.
3 //!
4 //! *[See also the `f64` primitive type](../../std/primitive.f64.html).*
5 //!
6 //! Mathematically significant numbers are provided in the `consts` sub-module.
7 //!
8 //! Although using these constants won’t cause compilation warnings,
9 //! new code should use the associated constants directly on the primitive type.
10
11 #![stable(feature = "rust1", since = "1.0.0")]
12
13 use crate::convert::FloatToInt;
14 #[cfg(not(test))]
15 use crate::intrinsics;
16 use crate::mem;
17 use crate::num::FpCategory;
18
19 /// The radix or base of the internal representation of `f64`.
20 /// Use [`f64::RADIX`](../../std/primitive.f64.html#associatedconstant.RADIX) instead.
21 ///
22 /// # Examples
23 ///
24 /// ```rust
25 /// // deprecated way
26 /// let r = std::f64::RADIX;
27 ///
28 /// // intended way
29 /// let r = f64::RADIX;
30 /// ```
31 #[stable(feature = "rust1", since = "1.0.0")]
32 pub const RADIX: u32 = f64::RADIX;
33
34 /// Number of significant digits in base 2.
35 /// Use [`f64::MANTISSA_DIGITS`](../../std/primitive.f64.html#associatedconstant.MANTISSA_DIGITS) instead.
36 ///
37 /// # Examples
38 ///
39 /// ```rust
40 /// // deprecated way
41 /// let d = std::f64::MANTISSA_DIGITS;
42 ///
43 /// // intended way
44 /// let d = f64::MANTISSA_DIGITS;
45 /// ```
46 #[stable(feature = "rust1", since = "1.0.0")]
47 pub const MANTISSA_DIGITS: u32 = f64::MANTISSA_DIGITS;
48
49 /// Approximate number of significant digits in base 10.
50 /// Use [`f64::DIGITS`](../../std/primitive.f64.html#associatedconstant.DIGITS) instead.
51 ///
52 /// # Examples
53 ///
54 /// ```rust
55 /// // deprecated way
56 /// let d = std::f64::DIGITS;
57 ///
58 /// // intended way
59 /// let d = f64::DIGITS;
60 /// ```
61 #[stable(feature = "rust1", since = "1.0.0")]
62 pub const DIGITS: u32 = f64::DIGITS;
63
64 /// [Machine epsilon] value for `f64`.
65 /// Use [`f64::EPSILON`](../../std/primitive.f64.html#associatedconstant.EPSILON) instead.
66 ///
67 /// This is the difference between `1.0` and the next larger representable number.
68 ///
69 /// [Machine epsilon]: https://en.wikipedia.org/wiki/Machine_epsilon
70 ///
71 /// # Examples
72 ///
73 /// ```rust
74 /// // deprecated way
75 /// let e = std::f64::EPSILON;
76 ///
77 /// // intended way
78 /// let e = f64::EPSILON;
79 /// ```
80 #[stable(feature = "rust1", since = "1.0.0")]
81 pub const EPSILON: f64 = f64::EPSILON;
82
83 /// Smallest finite `f64` value.
84 /// Use [`f64::MIN`](../../std/primitive.f64.html#associatedconstant.MIN) instead.
85 ///
86 /// # Examples
87 ///
88 /// ```rust
89 /// // deprecated way
90 /// let min = std::f64::MIN;
91 ///
92 /// // intended way
93 /// let min = f64::MIN;
94 /// ```
95 #[stable(feature = "rust1", since = "1.0.0")]
96 pub const MIN: f64 = f64::MIN;
97
98 /// Smallest positive normal `f64` value.
99 /// Use [`f64::MIN_POSITIVE`](../../std/primitive.f64.html#associatedconstant.MIN_POSITIVE) instead.
100 ///
101 /// # Examples
102 ///
103 /// ```rust
104 /// // deprecated way
105 /// let min = std::f64::MIN_POSITIVE;
106 ///
107 /// // intended way
108 /// let min = f64::MIN_POSITIVE;
109 /// ```
110 #[stable(feature = "rust1", since = "1.0.0")]
111 pub const MIN_POSITIVE: f64 = f64::MIN_POSITIVE;
112
113 /// Largest finite `f64` value.
114 /// Use [`f64::MAX`](../../std/primitive.f64.html#associatedconstant.MAX) instead.
115 ///
116 /// # Examples
117 ///
118 /// ```rust
119 /// // deprecated way
120 /// let max = std::f64::MAX;
121 ///
122 /// // intended way
123 /// let max = f64::MAX;
124 /// ```
125 #[stable(feature = "rust1", since = "1.0.0")]
126 pub const MAX: f64 = f64::MAX;
127
128 /// One greater than the minimum possible normal power of 2 exponent.
129 /// Use [`f64::MIN_EXP`](../../std/primitive.f64.html#associatedconstant.MIN_EXP) instead.
130 ///
131 /// # Examples
132 ///
133 /// ```rust
134 /// // deprecated way
135 /// let min = std::f64::MIN_EXP;
136 ///
137 /// // intended way
138 /// let min = f64::MIN_EXP;
139 /// ```
140 #[stable(feature = "rust1", since = "1.0.0")]
141 pub const MIN_EXP: i32 = f64::MIN_EXP;
142
143 /// Maximum possible power of 2 exponent.
144 /// Use [`f64::MAX_EXP`](../../std/primitive.f64.html#associatedconstant.MAX_EXP) instead.
145 ///
146 /// # Examples
147 ///
148 /// ```rust
149 /// // deprecated way
150 /// let max = std::f64::MAX_EXP;
151 ///
152 /// // intended way
153 /// let max = f64::MAX_EXP;
154 /// ```
155 #[stable(feature = "rust1", since = "1.0.0")]
156 pub const MAX_EXP: i32 = f64::MAX_EXP;
157
158 /// Minimum possible normal power of 10 exponent.
159 /// Use [`f64::MIN_10_EXP`](../../std/primitive.f64.html#associatedconstant.MIN_10_EXP) instead.
160 ///
161 /// # Examples
162 ///
163 /// ```rust
164 /// // deprecated way
165 /// let min = std::f64::MIN_10_EXP;
166 ///
167 /// // intended way
168 /// let min = f64::MIN_10_EXP;
169 /// ```
170 #[stable(feature = "rust1", since = "1.0.0")]
171 pub const MIN_10_EXP: i32 = f64::MIN_10_EXP;
172
173 /// Maximum possible power of 10 exponent.
174 /// Use [`f64::MAX_10_EXP`](../../std/primitive.f64.html#associatedconstant.MAX_10_EXP) instead.
175 ///
176 /// # Examples
177 ///
178 /// ```rust
179 /// // deprecated way
180 /// let max = std::f64::MAX_10_EXP;
181 ///
182 /// // intended way
183 /// let max = f64::MAX_10_EXP;
184 /// ```
185 #[stable(feature = "rust1", since = "1.0.0")]
186 pub const MAX_10_EXP: i32 = f64::MAX_10_EXP;
187
188 /// Not a Number (NaN).
189 /// Use [`f64::NAN`](../../std/primitive.f64.html#associatedconstant.NAN) instead.
190 ///
191 /// # Examples
192 ///
193 /// ```rust
194 /// // deprecated way
195 /// let nan = std::f64::NAN;
196 ///
197 /// // intended way
198 /// let nan = f64::NAN;
199 /// ```
200 #[stable(feature = "rust1", since = "1.0.0")]
201 pub const NAN: f64 = f64::NAN;
202
203 /// Infinity (∞).
204 /// Use [`f64::INFINITY`](../../std/primitive.f64.html#associatedconstant.INFINITY) instead.
205 ///
206 /// # Examples
207 ///
208 /// ```rust
209 /// // deprecated way
210 /// let inf = std::f64::INFINITY;
211 ///
212 /// // intended way
213 /// let inf = f64::INFINITY;
214 /// ```
215 #[stable(feature = "rust1", since = "1.0.0")]
216 pub const INFINITY: f64 = f64::INFINITY;
217
218 /// Negative infinity (−∞).
219 /// Use [`f64::NEG_INFINITY`](../../std/primitive.f64.html#associatedconstant.NEG_INFINITY) instead.
220 ///
221 /// # Examples
222 ///
223 /// ```rust
224 /// // deprecated way
225 /// let ninf = std::f64::NEG_INFINITY;
226 ///
227 /// // intended way
228 /// let ninf = f64::NEG_INFINITY;
229 /// ```
230 #[stable(feature = "rust1", since = "1.0.0")]
231 pub const NEG_INFINITY: f64 = f64::NEG_INFINITY;
232
233 /// Basic mathematical constants.
234 #[stable(feature = "rust1", since = "1.0.0")]
235 pub mod consts {
236     // FIXME: replace with mathematical constants from cmath.
237
238     /// Archimedes' constant (π)
239     #[stable(feature = "rust1", since = "1.0.0")]
240     pub const PI: f64 = 3.14159265358979323846264338327950288_f64;
241
242     /// The full circle constant (τ)
243     ///
244     /// Equal to 2π.
245     #[stable(feature = "tau_constant", since = "1.47.0")]
246     pub const TAU: f64 = 6.28318530717958647692528676655900577_f64;
247
248     /// π/2
249     #[stable(feature = "rust1", since = "1.0.0")]
250     pub const FRAC_PI_2: f64 = 1.57079632679489661923132169163975144_f64;
251
252     /// π/3
253     #[stable(feature = "rust1", since = "1.0.0")]
254     pub const FRAC_PI_3: f64 = 1.04719755119659774615421446109316763_f64;
255
256     /// π/4
257     #[stable(feature = "rust1", since = "1.0.0")]
258     pub const FRAC_PI_4: f64 = 0.785398163397448309615660845819875721_f64;
259
260     /// π/6
261     #[stable(feature = "rust1", since = "1.0.0")]
262     pub const FRAC_PI_6: f64 = 0.52359877559829887307710723054658381_f64;
263
264     /// π/8
265     #[stable(feature = "rust1", since = "1.0.0")]
266     pub const FRAC_PI_8: f64 = 0.39269908169872415480783042290993786_f64;
267
268     /// 1/π
269     #[stable(feature = "rust1", since = "1.0.0")]
270     pub const FRAC_1_PI: f64 = 0.318309886183790671537767526745028724_f64;
271
272     /// 2/π
273     #[stable(feature = "rust1", since = "1.0.0")]
274     pub const FRAC_2_PI: f64 = 0.636619772367581343075535053490057448_f64;
275
276     /// 2/sqrt(π)
277     #[stable(feature = "rust1", since = "1.0.0")]
278     pub const FRAC_2_SQRT_PI: f64 = 1.12837916709551257389615890312154517_f64;
279
280     /// sqrt(2)
281     #[stable(feature = "rust1", since = "1.0.0")]
282     pub const SQRT_2: f64 = 1.41421356237309504880168872420969808_f64;
283
284     /// 1/sqrt(2)
285     #[stable(feature = "rust1", since = "1.0.0")]
286     pub const FRAC_1_SQRT_2: f64 = 0.707106781186547524400844362104849039_f64;
287
288     /// Euler's number (e)
289     #[stable(feature = "rust1", since = "1.0.0")]
290     pub const E: f64 = 2.71828182845904523536028747135266250_f64;
291
292     /// log<sub>2</sub>(10)
293     #[stable(feature = "extra_log_consts", since = "1.43.0")]
294     pub const LOG2_10: f64 = 3.32192809488736234787031942948939018_f64;
295
296     /// log<sub>2</sub>(e)
297     #[stable(feature = "rust1", since = "1.0.0")]
298     pub const LOG2_E: f64 = 1.44269504088896340735992468100189214_f64;
299
300     /// log<sub>10</sub>(2)
301     #[stable(feature = "extra_log_consts", since = "1.43.0")]
302     pub const LOG10_2: f64 = 0.301029995663981195213738894724493027_f64;
303
304     /// log<sub>10</sub>(e)
305     #[stable(feature = "rust1", since = "1.0.0")]
306     pub const LOG10_E: f64 = 0.434294481903251827651128918916605082_f64;
307
308     /// ln(2)
309     #[stable(feature = "rust1", since = "1.0.0")]
310     pub const LN_2: f64 = 0.693147180559945309417232121458176568_f64;
311
312     /// ln(10)
313     #[stable(feature = "rust1", since = "1.0.0")]
314     pub const LN_10: f64 = 2.30258509299404568401799145468436421_f64;
315 }
316
317 #[lang = "f64"]
318 #[cfg(not(test))]
319 impl f64 {
320     /// The radix or base of the internal representation of `f64`.
321     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
322     pub const RADIX: u32 = 2;
323
324     /// Number of significant digits in base 2.
325     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
326     pub const MANTISSA_DIGITS: u32 = 53;
327     /// Approximate number of significant digits in base 10.
328     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
329     pub const DIGITS: u32 = 15;
330
331     /// [Machine epsilon] value for `f64`.
332     ///
333     /// This is the difference between `1.0` and the next larger representable number.
334     ///
335     /// [Machine epsilon]: https://en.wikipedia.org/wiki/Machine_epsilon
336     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
337     pub const EPSILON: f64 = 2.2204460492503131e-16_f64;
338
339     /// Smallest finite `f64` value.
340     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
341     pub const MIN: f64 = -1.7976931348623157e+308_f64;
342     /// Smallest positive normal `f64` value.
343     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
344     pub const MIN_POSITIVE: f64 = 2.2250738585072014e-308_f64;
345     /// Largest finite `f64` value.
346     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
347     pub const MAX: f64 = 1.7976931348623157e+308_f64;
348
349     /// One greater than the minimum possible normal power of 2 exponent.
350     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
351     pub const MIN_EXP: i32 = -1021;
352     /// Maximum possible power of 2 exponent.
353     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
354     pub const MAX_EXP: i32 = 1024;
355
356     /// Minimum possible normal power of 10 exponent.
357     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
358     pub const MIN_10_EXP: i32 = -307;
359     /// Maximum possible power of 10 exponent.
360     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
361     pub const MAX_10_EXP: i32 = 308;
362
363     /// Not a Number (NaN).
364     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
365     pub const NAN: f64 = 0.0_f64 / 0.0_f64;
366     /// Infinity (∞).
367     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
368     pub const INFINITY: f64 = 1.0_f64 / 0.0_f64;
369     /// Negative infinity (−∞).
370     #[stable(feature = "assoc_int_consts", since = "1.43.0")]
371     pub const NEG_INFINITY: f64 = -1.0_f64 / 0.0_f64;
372
373     /// Returns `true` if this value is `NaN`.
374     ///
375     /// ```
376     /// let nan = f64::NAN;
377     /// let f = 7.0_f64;
378     ///
379     /// assert!(nan.is_nan());
380     /// assert!(!f.is_nan());
381     /// ```
382     #[stable(feature = "rust1", since = "1.0.0")]
383     #[rustc_const_unstable(feature = "const_float_classify", issue = "72505")]
384     #[inline]
385     pub const fn is_nan(self) -> bool {
386         self != self
387     }
388
389     // FIXME(#50145): `abs` is publicly unavailable in libcore due to
390     // concerns about portability, so this implementation is for
391     // private use internally.
392     #[inline]
393     #[rustc_const_unstable(feature = "const_float_classify", issue = "72505")]
394     const fn abs_private(self) -> f64 {
395         f64::from_bits(self.to_bits() & 0x7fff_ffff_ffff_ffff)
396     }
397
398     /// Returns `true` if this value is positive infinity or negative infinity, and
399     /// `false` otherwise.
400     ///
401     /// ```
402     /// let f = 7.0f64;
403     /// let inf = f64::INFINITY;
404     /// let neg_inf = f64::NEG_INFINITY;
405     /// let nan = f64::NAN;
406     ///
407     /// assert!(!f.is_infinite());
408     /// assert!(!nan.is_infinite());
409     ///
410     /// assert!(inf.is_infinite());
411     /// assert!(neg_inf.is_infinite());
412     /// ```
413     #[stable(feature = "rust1", since = "1.0.0")]
414     #[rustc_const_unstable(feature = "const_float_classify", issue = "72505")]
415     #[inline]
416     pub const fn is_infinite(self) -> bool {
417         self.abs_private() == Self::INFINITY
418     }
419
420     /// Returns `true` if this number is neither infinite nor `NaN`.
421     ///
422     /// ```
423     /// let f = 7.0f64;
424     /// let inf: f64 = f64::INFINITY;
425     /// let neg_inf: f64 = f64::NEG_INFINITY;
426     /// let nan: f64 = f64::NAN;
427     ///
428     /// assert!(f.is_finite());
429     ///
430     /// assert!(!nan.is_finite());
431     /// assert!(!inf.is_finite());
432     /// assert!(!neg_inf.is_finite());
433     /// ```
434     #[stable(feature = "rust1", since = "1.0.0")]
435     #[rustc_const_unstable(feature = "const_float_classify", issue = "72505")]
436     #[inline]
437     pub const fn is_finite(self) -> bool {
438         // There's no need to handle NaN separately: if self is NaN,
439         // the comparison is not true, exactly as desired.
440         self.abs_private() < Self::INFINITY
441     }
442
443     /// Returns `true` if the number is neither zero, infinite,
444     /// [subnormal], or `NaN`.
445     ///
446     /// ```
447     /// let min = f64::MIN_POSITIVE; // 2.2250738585072014e-308f64
448     /// let max = f64::MAX;
449     /// let lower_than_min = 1.0e-308_f64;
450     /// let zero = 0.0f64;
451     ///
452     /// assert!(min.is_normal());
453     /// assert!(max.is_normal());
454     ///
455     /// assert!(!zero.is_normal());
456     /// assert!(!f64::NAN.is_normal());
457     /// assert!(!f64::INFINITY.is_normal());
458     /// // Values between `0` and `min` are Subnormal.
459     /// assert!(!lower_than_min.is_normal());
460     /// ```
461     /// [subnormal]: https://en.wikipedia.org/wiki/Denormal_number
462     #[stable(feature = "rust1", since = "1.0.0")]
463     #[rustc_const_unstable(feature = "const_float_classify", issue = "72505")]
464     #[inline]
465     pub const fn is_normal(self) -> bool {
466         matches!(self.classify(), FpCategory::Normal)
467     }
468
469     /// Returns the floating point category of the number. If only one property
470     /// is going to be tested, it is generally faster to use the specific
471     /// predicate instead.
472     ///
473     /// ```
474     /// use std::num::FpCategory;
475     ///
476     /// let num = 12.4_f64;
477     /// let inf = f64::INFINITY;
478     ///
479     /// assert_eq!(num.classify(), FpCategory::Normal);
480     /// assert_eq!(inf.classify(), FpCategory::Infinite);
481     /// ```
482     #[stable(feature = "rust1", since = "1.0.0")]
483     #[rustc_const_unstable(feature = "const_float_classify", issue = "72505")]
484     pub const fn classify(self) -> FpCategory {
485         const EXP_MASK: u64 = 0x7ff0000000000000;
486         const MAN_MASK: u64 = 0x000fffffffffffff;
487
488         let bits = self.to_bits();
489         match (bits & MAN_MASK, bits & EXP_MASK) {
490             (0, 0) => FpCategory::Zero,
491             (_, 0) => FpCategory::Subnormal,
492             (0, EXP_MASK) => FpCategory::Infinite,
493             (_, EXP_MASK) => FpCategory::Nan,
494             _ => FpCategory::Normal,
495         }
496     }
497
498     /// Returns `true` if `self` has a positive sign, including `+0.0`, `NaN`s with
499     /// positive sign bit and positive infinity.
500     ///
501     /// ```
502     /// let f = 7.0_f64;
503     /// let g = -7.0_f64;
504     ///
505     /// assert!(f.is_sign_positive());
506     /// assert!(!g.is_sign_positive());
507     /// ```
508     #[stable(feature = "rust1", since = "1.0.0")]
509     #[rustc_const_unstable(feature = "const_float_classify", issue = "72505")]
510     #[inline]
511     pub const fn is_sign_positive(self) -> bool {
512         !self.is_sign_negative()
513     }
514
515     #[stable(feature = "rust1", since = "1.0.0")]
516     #[rustc_deprecated(since = "1.0.0", reason = "renamed to is_sign_positive")]
517     #[inline]
518     #[doc(hidden)]
519     pub fn is_positive(self) -> bool {
520         self.is_sign_positive()
521     }
522
523     /// Returns `true` if `self` has a negative sign, including `-0.0`, `NaN`s with
524     /// negative sign bit and negative infinity.
525     ///
526     /// ```
527     /// let f = 7.0_f64;
528     /// let g = -7.0_f64;
529     ///
530     /// assert!(!f.is_sign_negative());
531     /// assert!(g.is_sign_negative());
532     /// ```
533     #[stable(feature = "rust1", since = "1.0.0")]
534     #[rustc_const_unstable(feature = "const_float_classify", issue = "72505")]
535     #[inline]
536     pub const fn is_sign_negative(self) -> bool {
537         self.to_bits() & 0x8000_0000_0000_0000 != 0
538     }
539
540     #[stable(feature = "rust1", since = "1.0.0")]
541     #[rustc_deprecated(since = "1.0.0", reason = "renamed to is_sign_negative")]
542     #[inline]
543     #[doc(hidden)]
544     pub fn is_negative(self) -> bool {
545         self.is_sign_negative()
546     }
547
548     /// Takes the reciprocal (inverse) of a number, `1/x`.
549     ///
550     /// ```
551     /// let x = 2.0_f64;
552     /// let abs_difference = (x.recip() - (1.0 / x)).abs();
553     ///
554     /// assert!(abs_difference < 1e-10);
555     /// ```
556     #[stable(feature = "rust1", since = "1.0.0")]
557     #[inline]
558     pub fn recip(self) -> f64 {
559         1.0 / self
560     }
561
562     /// Converts radians to degrees.
563     ///
564     /// ```
565     /// let angle = std::f64::consts::PI;
566     ///
567     /// let abs_difference = (angle.to_degrees() - 180.0).abs();
568     ///
569     /// assert!(abs_difference < 1e-10);
570     /// ```
571     #[stable(feature = "rust1", since = "1.0.0")]
572     #[inline]
573     pub fn to_degrees(self) -> f64 {
574         // The division here is correctly rounded with respect to the true
575         // value of 180/π. (This differs from f32, where a constant must be
576         // used to ensure a correctly rounded result.)
577         self * (180.0f64 / consts::PI)
578     }
579
580     /// Converts degrees to radians.
581     ///
582     /// ```
583     /// let angle = 180.0_f64;
584     ///
585     /// let abs_difference = (angle.to_radians() - std::f64::consts::PI).abs();
586     ///
587     /// assert!(abs_difference < 1e-10);
588     /// ```
589     #[stable(feature = "rust1", since = "1.0.0")]
590     #[inline]
591     pub fn to_radians(self) -> f64 {
592         let value: f64 = consts::PI;
593         self * (value / 180.0)
594     }
595
596     /// Returns the maximum of the two numbers.
597     ///
598     /// ```
599     /// let x = 1.0_f64;
600     /// let y = 2.0_f64;
601     ///
602     /// assert_eq!(x.max(y), y);
603     /// ```
604     ///
605     /// If one of the arguments is NaN, then the other argument is returned.
606     #[stable(feature = "rust1", since = "1.0.0")]
607     #[inline]
608     pub fn max(self, other: f64) -> f64 {
609         intrinsics::maxnumf64(self, other)
610     }
611
612     /// Returns the minimum of the two numbers.
613     ///
614     /// ```
615     /// let x = 1.0_f64;
616     /// let y = 2.0_f64;
617     ///
618     /// assert_eq!(x.min(y), x);
619     /// ```
620     ///
621     /// If one of the arguments is NaN, then the other argument is returned.
622     #[stable(feature = "rust1", since = "1.0.0")]
623     #[inline]
624     pub fn min(self, other: f64) -> f64 {
625         intrinsics::minnumf64(self, other)
626     }
627
628     /// Rounds toward zero and converts to any primitive integer type,
629     /// assuming that the value is finite and fits in that type.
630     ///
631     /// ```
632     /// let value = 4.6_f64;
633     /// let rounded = unsafe { value.to_int_unchecked::<u16>() };
634     /// assert_eq!(rounded, 4);
635     ///
636     /// let value = -128.9_f64;
637     /// let rounded = unsafe { value.to_int_unchecked::<i8>() };
638     /// assert_eq!(rounded, i8::MIN);
639     /// ```
640     ///
641     /// # Safety
642     ///
643     /// The value must:
644     ///
645     /// * Not be `NaN`
646     /// * Not be infinite
647     /// * Be representable in the return type `Int`, after truncating off its fractional part
648     #[stable(feature = "float_approx_unchecked_to", since = "1.44.0")]
649     #[inline]
650     pub unsafe fn to_int_unchecked<Int>(self) -> Int
651     where
652         Self: FloatToInt<Int>,
653     {
654         // SAFETY: the caller must uphold the safety contract for
655         // `FloatToInt::to_int_unchecked`.
656         unsafe { FloatToInt::<Int>::to_int_unchecked(self) }
657     }
658
659     /// Raw transmutation to `u64`.
660     ///
661     /// This is currently identical to `transmute::<f64, u64>(self)` on all platforms.
662     ///
663     /// See `from_bits` for some discussion of the portability of this operation
664     /// (there are almost no issues).
665     ///
666     /// Note that this function is distinct from `as` casting, which attempts to
667     /// preserve the *numeric* value, and not the bitwise value.
668     ///
669     /// # Examples
670     ///
671     /// ```
672     /// assert!((1f64).to_bits() != 1f64 as u64); // to_bits() is not casting!
673     /// assert_eq!((12.5f64).to_bits(), 0x4029000000000000);
674     ///
675     /// ```
676     #[stable(feature = "float_bits_conv", since = "1.20.0")]
677     #[rustc_const_unstable(feature = "const_float_bits_conv", issue = "72447")]
678     #[inline]
679     pub const fn to_bits(self) -> u64 {
680         // SAFETY: `u64` is a plain old datatype so we can always transmute to it
681         unsafe { mem::transmute(self) }
682     }
683
684     /// Raw transmutation from `u64`.
685     ///
686     /// This is currently identical to `transmute::<u64, f64>(v)` on all platforms.
687     /// It turns out this is incredibly portable, for two reasons:
688     ///
689     /// * Floats and Ints have the same endianness on all supported platforms.
690     /// * IEEE-754 very precisely specifies the bit layout of floats.
691     ///
692     /// However there is one caveat: prior to the 2008 version of IEEE-754, how
693     /// to interpret the NaN signaling bit wasn't actually specified. Most platforms
694     /// (notably x86 and ARM) picked the interpretation that was ultimately
695     /// standardized in 2008, but some didn't (notably MIPS). As a result, all
696     /// signaling NaNs on MIPS are quiet NaNs on x86, and vice-versa.
697     ///
698     /// Rather than trying to preserve signaling-ness cross-platform, this
699     /// implementation favors preserving the exact bits. This means that
700     /// any payloads encoded in NaNs will be preserved even if the result of
701     /// this method is sent over the network from an x86 machine to a MIPS one.
702     ///
703     /// If the results of this method are only manipulated by the same
704     /// architecture that produced them, then there is no portability concern.
705     ///
706     /// If the input isn't NaN, then there is no portability concern.
707     ///
708     /// If you don't care about signaling-ness (very likely), then there is no
709     /// portability concern.
710     ///
711     /// Note that this function is distinct from `as` casting, which attempts to
712     /// preserve the *numeric* value, and not the bitwise value.
713     ///
714     /// # Examples
715     ///
716     /// ```
717     /// let v = f64::from_bits(0x4029000000000000);
718     /// assert_eq!(v, 12.5);
719     /// ```
720     #[stable(feature = "float_bits_conv", since = "1.20.0")]
721     #[rustc_const_unstable(feature = "const_float_bits_conv", issue = "72447")]
722     #[inline]
723     pub const fn from_bits(v: u64) -> Self {
724         // SAFETY: `u64` is a plain old datatype so we can always transmute from it
725         // It turns out the safety issues with sNaN were overblown! Hooray!
726         unsafe { mem::transmute(v) }
727     }
728
729     /// Return the memory representation of this floating point number as a byte array in
730     /// big-endian (network) byte order.
731     ///
732     /// # Examples
733     ///
734     /// ```
735     /// let bytes = 12.5f64.to_be_bytes();
736     /// assert_eq!(bytes, [0x40, 0x29, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00]);
737     /// ```
738     #[stable(feature = "float_to_from_bytes", since = "1.40.0")]
739     #[rustc_const_unstable(feature = "const_float_bits_conv", issue = "72447")]
740     #[inline]
741     pub const fn to_be_bytes(self) -> [u8; 8] {
742         self.to_bits().to_be_bytes()
743     }
744
745     /// Return the memory representation of this floating point number as a byte array in
746     /// little-endian byte order.
747     ///
748     /// # Examples
749     ///
750     /// ```
751     /// let bytes = 12.5f64.to_le_bytes();
752     /// assert_eq!(bytes, [0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x29, 0x40]);
753     /// ```
754     #[stable(feature = "float_to_from_bytes", since = "1.40.0")]
755     #[rustc_const_unstable(feature = "const_float_bits_conv", issue = "72447")]
756     #[inline]
757     pub const fn to_le_bytes(self) -> [u8; 8] {
758         self.to_bits().to_le_bytes()
759     }
760
761     /// Return the memory representation of this floating point number as a byte array in
762     /// native byte order.
763     ///
764     /// As the target platform's native endianness is used, portable code
765     /// should use [`to_be_bytes`] or [`to_le_bytes`], as appropriate, instead.
766     ///
767     /// [`to_be_bytes`]: #method.to_be_bytes
768     /// [`to_le_bytes`]: #method.to_le_bytes
769     ///
770     /// # Examples
771     ///
772     /// ```
773     /// let bytes = 12.5f64.to_ne_bytes();
774     /// assert_eq!(
775     ///     bytes,
776     ///     if cfg!(target_endian = "big") {
777     ///         [0x40, 0x29, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00]
778     ///     } else {
779     ///         [0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x29, 0x40]
780     ///     }
781     /// );
782     /// ```
783     #[stable(feature = "float_to_from_bytes", since = "1.40.0")]
784     #[rustc_const_unstable(feature = "const_float_bits_conv", issue = "72447")]
785     #[inline]
786     pub const fn to_ne_bytes(self) -> [u8; 8] {
787         self.to_bits().to_ne_bytes()
788     }
789
790     /// Return the memory representation of this floating point number as a byte array in
791     /// native byte order.
792     ///
793     /// [`to_ne_bytes`] should be preferred over this whenever possible.
794     ///
795     /// [`to_ne_bytes`]: #method.to_ne_bytes
796     ///
797     /// # Examples
798     ///
799     /// ```
800     /// #![feature(num_as_ne_bytes)]
801     /// let num = 12.5f64;
802     /// let bytes = num.as_ne_bytes();
803     /// assert_eq!(
804     ///     bytes,
805     ///     if cfg!(target_endian = "big") {
806     ///         &[0x40, 0x29, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00]
807     ///     } else {
808     ///         &[0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x29, 0x40]
809     ///     }
810     /// );
811     /// ```
812     #[unstable(feature = "num_as_ne_bytes", issue = "76976")]
813     #[inline]
814     pub fn as_ne_bytes(&self) -> &[u8; 8] {
815         // SAFETY: `f64` is a plain old datatype so we can always transmute to it
816         unsafe { &*(self as *const Self as *const _) }
817     }
818
819     /// Create a floating point value from its representation as a byte array in big endian.
820     ///
821     /// # Examples
822     ///
823     /// ```
824     /// let value = f64::from_be_bytes([0x40, 0x29, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00]);
825     /// assert_eq!(value, 12.5);
826     /// ```
827     #[stable(feature = "float_to_from_bytes", since = "1.40.0")]
828     #[rustc_const_unstable(feature = "const_float_bits_conv", issue = "72447")]
829     #[inline]
830     pub const fn from_be_bytes(bytes: [u8; 8]) -> Self {
831         Self::from_bits(u64::from_be_bytes(bytes))
832     }
833
834     /// Create a floating point value from its representation as a byte array in little endian.
835     ///
836     /// # Examples
837     ///
838     /// ```
839     /// let value = f64::from_le_bytes([0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x29, 0x40]);
840     /// assert_eq!(value, 12.5);
841     /// ```
842     #[stable(feature = "float_to_from_bytes", since = "1.40.0")]
843     #[rustc_const_unstable(feature = "const_float_bits_conv", issue = "72447")]
844     #[inline]
845     pub const fn from_le_bytes(bytes: [u8; 8]) -> Self {
846         Self::from_bits(u64::from_le_bytes(bytes))
847     }
848
849     /// Create a floating point value from its representation as a byte array in native endian.
850     ///
851     /// As the target platform's native endianness is used, portable code
852     /// likely wants to use [`from_be_bytes`] or [`from_le_bytes`], as
853     /// appropriate instead.
854     ///
855     /// [`from_be_bytes`]: #method.from_be_bytes
856     /// [`from_le_bytes`]: #method.from_le_bytes
857     ///
858     /// # Examples
859     ///
860     /// ```
861     /// let value = f64::from_ne_bytes(if cfg!(target_endian = "big") {
862     ///     [0x40, 0x29, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00]
863     /// } else {
864     ///     [0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x29, 0x40]
865     /// });
866     /// assert_eq!(value, 12.5);
867     /// ```
868     #[stable(feature = "float_to_from_bytes", since = "1.40.0")]
869     #[rustc_const_unstable(feature = "const_float_bits_conv", issue = "72447")]
870     #[inline]
871     pub const fn from_ne_bytes(bytes: [u8; 8]) -> Self {
872         Self::from_bits(u64::from_ne_bytes(bytes))
873     }
874
875     /// Returns an ordering between self and other values.
876     /// Unlike the standard partial comparison between floating point numbers,
877     /// this comparison always produces an ordering in accordance to
878     /// the totalOrder predicate as defined in IEEE 754 (2008 revision)
879     /// floating point standard. The values are ordered in following order:
880     /// - Negative quiet NaN
881     /// - Negative signaling NaN
882     /// - Negative infinity
883     /// - Negative numbers
884     /// - Negative subnormal numbers
885     /// - Negative zero
886     /// - Positive zero
887     /// - Positive subnormal numbers
888     /// - Positive numbers
889     /// - Positive infinity
890     /// - Positive signaling NaN
891     /// - Positive quiet NaN
892     ///
893     /// Note that this function does not always agree with the [`PartialOrd`]
894     /// and [`PartialEq`] implementations of `f64`. In particular, they regard
895     /// negative and positive zero as equal, while `total_cmp` doesn't.
896     ///
897     /// # Example
898     /// ```
899     /// #![feature(total_cmp)]
900     /// struct GoodBoy {
901     ///     name: String,
902     ///     weight: f64,
903     /// }
904     ///
905     /// let mut bois = vec![
906     ///     GoodBoy { name: "Pucci".to_owned(), weight: 0.1 },
907     ///     GoodBoy { name: "Woofer".to_owned(), weight: 99.0 },
908     ///     GoodBoy { name: "Yapper".to_owned(), weight: 10.0 },
909     ///     GoodBoy { name: "Chonk".to_owned(), weight: f64::INFINITY },
910     ///     GoodBoy { name: "Abs. Unit".to_owned(), weight: f64::NAN },
911     ///     GoodBoy { name: "Floaty".to_owned(), weight: -5.0 },
912     /// ];
913     ///
914     /// bois.sort_by(|a, b| a.weight.total_cmp(&b.weight));
915     /// # assert!(bois.into_iter().map(|b| b.weight)
916     /// #     .zip([-5.0, 0.1, 10.0, 99.0, f64::INFINITY, f64::NAN].iter())
917     /// #     .all(|(a, b)| a.to_bits() == b.to_bits()))
918     /// ```
919     #[unstable(feature = "total_cmp", issue = "72599")]
920     #[inline]
921     pub fn total_cmp(&self, other: &Self) -> crate::cmp::Ordering {
922         let mut left = self.to_bits() as i64;
923         let mut right = other.to_bits() as i64;
924
925         // In case of negatives, flip all the bits except the sign
926         // to achieve a similar layout as two's complement integers
927         //
928         // Why does this work? IEEE 754 floats consist of three fields:
929         // Sign bit, exponent and mantissa. The set of exponent and mantissa
930         // fields as a whole have the property that their bitwise order is
931         // equal to the numeric magnitude where the magnitude is defined.
932         // The magnitude is not normally defined on NaN values, but
933         // IEEE 754 totalOrder defines the NaN values also to follow the
934         // bitwise order. This leads to order explained in the doc comment.
935         // However, the representation of magnitude is the same for negative
936         // and positive numbers – only the sign bit is different.
937         // To easily compare the floats as signed integers, we need to
938         // flip the exponent and mantissa bits in case of negative numbers.
939         // We effectively convert the numbers to "two's complement" form.
940         //
941         // To do the flipping, we construct a mask and XOR against it.
942         // We branchlessly calculate an "all-ones except for the sign bit"
943         // mask from negative-signed values: right shifting sign-extends
944         // the integer, so we "fill" the mask with sign bits, and then
945         // convert to unsigned to push one more zero bit.
946         // On positive values, the mask is all zeros, so it's a no-op.
947         left ^= (((left >> 63) as u64) >> 1) as i64;
948         right ^= (((right >> 63) as u64) >> 1) as i64;
949
950         left.cmp(&right)
951     }
952 }