Function circular_dichroism 

fn circular_dichroism(
    matrix: &DMatrix<Complex>,
    pol: &[u8],
    ks: [f64; 2],
) -> Result<(f64, MetricResidual)>

Absorption circular dichroism with its matrix and wavenumber gradients.