Searched refs:jac_pre (Results 1 – 1 of 1) sorted by relevance
160 static PetscErrorCode SNESVIComputeJacobian(Mat jac, Mat jac_pre, Vec Da, Vec Db) in SNESVIComputeJacobian() argument166 if (jac != jac_pre) { /* If jac and jac_pre are different */ in SNESVIComputeJacobian()167 PetscCall(MatDiagonalScale(jac_pre, Db, NULL)); in SNESVIComputeJacobian()168 PetscCall(MatDiagonalSet(jac_pre, Da, ADD_VALUES)); in SNESVIComputeJacobian()