Enforcing Array Dimension Boundaries at Compile-Time via Const Generics

pub struct Matrix<const ROWS: usize, const COLS: usize> {
    pub data: [[f64; COLS]; ROWS],
}
 
impl<const ROWS: usize, const COLS: usize> Matrix<ROWS, COLS> {
    pub fn multiply<const OTHER_COLS: usize>(
        &self,
        other: &Matrix<COLS, OTHER_COLS>
    ) -> Matrix<ROWS, OTHER_COLS> {
        let mut result_data = [[0.0; OTHER_COLS]; ROWS];
        for r in 0..ROWS {
            for c in 0..OTHER_COLS {
                let mut sum = 0.0;
                for i in 0..COLS {
                    sum += self.data[r][i] * other.data[i][c];
                }
                result_data[r][c] = sum;
            }
        }
        Matrix { data: result_data }
    }
}