Rust における倍精度浮動小数点数の基礎

Rust 開発において、実数値を扱う際に最も頻繁に利用されるのが f64 型です。これは IEEE 754 規格に準拠した 64 ビットの倍精度浮動小数点数であり、約 15 から 17 桁の十進有効数字を表現できます。科学技術計算、物理シミュレーション、あるいは金融分野における高精度な数値処理が必要な场合に不可欠な型となります。

メモリ構造と表現範囲

f64 は 8 バイト（64 ビット）のメモリ領域を占有します。内部構造は、符号ビット 1 ビット、指数部 11 ビット、仮数部 52 ビットで構成されています。この構成により、極端に大きな値（例：1.7e308）や極めて小さな値（例：5e-324）まで広範な範囲をカバーすることが可能です。

変数の宣言と型推論

変数を定義する際、明示的に型を指定するか、接尾辞を使用して型を確定させることができます。文脈によっては型推論が働きますが、数値リテラルの場合は注意が必要です。

let circumference_ratio: f64 = 3.141592653589793;
let natural_base = 2.71828_f64; // 接尾辞で型を固定
let default_zero = 0.0;         // 文脈により f64 と推論される

// 整数との違いに注意
let integer_val = 5;   // i32 型として推論
let float_val = 5.0;   // f64 型として推論

Rust は型安全性を重視するため、整数から浮動小数点数への暗黙的な変換は行われません。必要に応じて明示的なキャストが必要です。

標準ライブラリによる数学演算

f64 型には、標準ライブラリを通じて豊富な数学関数がメソッドとして用意されています。これらを利用することで、複雑な計算も簡潔に記述可能です。

let base = 4.0_f64;
let exponent = 0.5_f64;

// べき乗計算（4 の 0.5 乗 = 平方根）
let result_pow = base.powf(exponent); // 2.0

// 平方根メソッド
let sqrt_val = 25.0_f64.sqrt(); // 5.0

// 三角関数（ラジアン単位）
let angle_rad = std::f64::consts::PI / 4.0;
let cos_val = angle_rad.cos(); // √2/2 ≈ 0.707

// 対数と指数関数
let log_val = base.ln();      // ln(4)
let exp_val = natural_base.exp(); // e^e

また、std::f64::consts モジュールを利用すれば、円周率 PI やネイピア数 E などの数学定数を高精度で呼び出すことができます。

特殊な数値状態の扱い

浮動小数点演算では、計算結果が定義できない場合や範囲を超えた場合に、特殊な値が生成されることがあります。代表的なものに NaN（Not a Number）および無限大（INFINITY, NEG_INFINITY）があります。

let invalid_op = 0.0_f64 / 0.0; // NaN
let overflow_val = f64::INFINITY;

// 特殊値の検出メソッド
println!("{} is NaN? {}", invalid_op, invalid_op.is_nan());       // true
println!("{} is finite? {}", overflow_val, overflow_val.is_finite()); // false
println!("{} is infinite? {}", overflow_val, overflow_val.is_infinite()); // true

重要な注意点として、NaN は自身を含めあらゆる値と等しくないと評価されます。したがって、等値演算子 == での判定は不可であり、必ず is_nan() メソッドを使用する必要があります。

精度誤差と比較処理のベストプラクティス

二进制浮動小数点数の仕組み上、十進法の特定の小数値は正確に表現できません。例えば、0.1 と 0.2 の和は厳密には 0.3 になりません。

let sum = 0.1_f64 + 0.2_f64;
println!("{} == 0.3? {}", sum, sum == 0.3); // false が返る

この問題に対処するには、許容誤差（イプシロン）を設定した比較関数を作成するのが一般的です。

fn check_approximation(left: f64, right: f64, tolerance: f64) -> bool {
    (left - right).abs() < tolerance
}

let calculated = 0.1 + 0.2;
let is_equal = check_approximation(calculated, 0.3, 1e-10);
println!("Approx equal? {}", is_equal); // true

数値計算を行う際は、浮動小数点数の本質的な限界を理解し、厳密な等値比較を避ける設計が求められます。

出力フォーマットと文字列化

値を文字列として出力する際、フォーマット指定子を用いて表示形式を制御できます。桁数や指数表記の切り替えが可能です。

let ratio = 100.0 / 3.0;

println!("{:.2}", ratio);   // 33.33（小数点以下 2 桁）
println!("{:.5}", ratio);   // 33.33333（小数点以下 5 桁）
println!("{:e}", ratio);    // 3.333333e1（指数表記）
println!("{:.2e}", ratio);  // 3.33e1（指数表記・精度指定）