Use regularised Jacobian when solving pressure system.

This commit is contained in:
Bård Skaflestad 2013-05-03 12:26:02 +02:00
parent cd9c9d5cde
commit dddb30b5a5

View File

@ -260,7 +260,7 @@ int main()
M matr = residual.derivative()[0];
matr.coeffRef(0,0) *= 2.0;
matr.makeCompressed();
solver.compute(residual.derivative()[0]);
solver.compute(matr);
// if (solver.info() != Eigen::Succeeded) {
// std::cerr << "Decomposition error!\n";
// return 1;