Function compute_optimized_constants 

pub fn compute_optimized_constants<F: Field, const N: usize>(
    mds: &[[F; N]; N],
    rounds_p: usize,
    partial_rc: &[[F; N]],
) -> ([F; N], Vec<F>, [[F; N]; N], Vec<[F; N]>, Vec<[F; N]>)

Compute all optimized partial round constants from raw parameters.

Combines the round constant compression and sparse matrix factorization into a single entry point, keeping the individual helpers private.

Returns

A tuple of:

• The compressed first-round constant vector (WIDTH elements).
• The optimized scalar round constants (RP-1 entries).
• The dense transition matrix m_i.
• The per-round sparse v vectors.
• The per-round sparse ŵ vectors.