Perform iterative refinement on the polished solution:
(repeat)
1. (K + dK) * dz = b - K*z
2. z <- z + dz
@param work Solver workspace
@param p Private variable for solving linear system
@param z Initial z value
@param b RHS of the linear system
@return Exitflag
Perform iterative refinement on the polished solution: (repeat) 1. (K + dK) * dz = b - K*z 2. z <- z + dz @param work Solver workspace @param p Private variable for solving linear system @param z Initial z value @param b RHS of the linear system @return Exitflag