pub fn linear_lcorner(rho: &VectorF64, eta: &VectorF64) -> (Value, usize)
Returns (Value, idx).
(Value, idx)