9 #ifndef __IP_AUGSYSTEMSOLVER_HPP__
10 #define __IP_AUGSYSTEMSOLVER_HPP__
52 const std::string& prefix) = 0;
83 Index numberOfNegEVals)
85 std::vector<SmartPtr<const Vector> > rhs_xV(1);
87 std::vector<SmartPtr<const Vector> > rhs_sV(1);
89 std::vector<SmartPtr<const Vector> > rhs_cV(1);
91 std::vector<SmartPtr<const Vector> > rhs_dV(1);
93 std::vector<SmartPtr<Vector> > sol_xV(1);
95 std::vector<SmartPtr<Vector> > sol_sV(1);
97 std::vector<SmartPtr<Vector> > sol_cV(1);
99 std::vector<SmartPtr<Vector> > sol_dV(1);
101 return MultiSolve(W, W_factor, D_x, delta_x, D_s, delta_s, J_c, D_c, delta_c,
102 J_d, D_d, delta_d, rhs_xV, rhs_sV, rhs_cV, rhs_dV,
103 sol_xV, sol_sV, sol_cV, sol_dV, check_NegEVals,
132 Index numberOfNegEVals)
146 for (
Index i=0; i<nrhs; i++) {
147 retval =
Solve(W, W_factor, D_x, delta_x, D_s, delta_s, J_c, D_c, delta_c,
149 *rhs_xV[i], *rhs_sV[i], *rhs_cV[i], *rhs_dV[i],
150 *sol_xV[i], *sol_sV[i], *sol_cV[i], *sol_dV[i],
151 check_NegEVals, numberOfNegEVals);