sofars 0.6.0

Pure Rust implementation of the IAU SOFA library
Documentation
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
use super::cr;

/// Multiply two r-matrices.
pub fn rxr(a: &[[f64; 3]; 3], b: &[[f64; 3]; 3], atb: &mut [[f64; 3]; 3]) {
    let wm = &mut [[0.0; 3]; 3];

    for i in 0..3 {
        for j in 0..3 {
            wm[i][j] = 0.0;
            for k in 0..3 {
                wm[i][j] += a[i][k] * b[k][j];
            }
        }
    }

    cr(wm, atb)
}