pub fn matrix_power<F>( a: &CsrMatrix<F>, power: i32, ) -> SparseResult<CsrMatrix<F>>where F: Float + NumAssign + Sum + SparseElement + 'static + Debug,
Matrix power