Function integer_matrix_power 

pub fn integer_matrix_power<F>(
    a: &ArrayView2<'_, F>,
    p: i32,
) -> LinalgResult<Array2<F>>
where
    F: Float + NumAssign + Sum + One + Send + Sync + ScalarOperand + 'static,

Compute matrix power for small integer powers using repeated squaring