still debugging, it seems the two rows of Jacobian and nfix are correctly added. Also, when I output J, v, and f in phase3 to matlab, it no longer says matrix singular. However, I'm still getting that Pastix arc error even though the matrix is no longer singular!!!