rfa 0.5.9

A port ERFA to Rust.
Documentation 
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

///  Multiply two r-matrices.
///
///  Given:
///   * a first r-matrix
///   * b second r-matrix
///
///  Returned:
///   * atb a * b
///
///  # Note:
///     It is permissible to re-use the same array for any of the
///     arguments.
///  This revision:  2021 May 11

pub fn rxr(a: &[[f64; 3]; 3], b: &[[f64; 3]; 3], atb: &mut[[f64; 3]; 3])
{
    for i in 0..3 {
       for j in 0..3 {
          let mut w = 0.0;
          for k in  0..3 {
             w +=  a[i][k] * b[k][j];
          }
          atb[i][j] = w;
       }
    }
}