A 30 old_code/2015linfty/dune_bdf.hxx Eigen::Vector3d operator() (const Eigen::Matrix3d& A, const Eigen::Vector3d& b) const { A 31 old_code/2015linfty/dune_bdf.hxx return A.fullPivLu().solve( (-1) * b);