Function quad_risk

crates/openquant/src/portfolio_optimization.rs:195–198

(cov: &DMatrix<f64>, w: &[f64])

}

fn quad_risk(cov: &DMatrix<f64>, w: &[f64]) -> f64 {
    let wv = DVector::from_vec(w.to_vec());
    (wv.transpose() * cov * wv)[(0, 0)]
}

fn ones(n: usize) -> DVector<f64> {
    DVector::from_element(n, 1.0)