pub fn forward_jacobian_vec_prod_vec<F>( x: &Vec<F>, fs: &'_ dyn Fn(&Vec<F>) -> Result<Vec<F>, Error>, p: &[F], ) -> Result<Vec<F>, Error>where F: Float,