Lines Matching refs:nzs
3169 PetscInt nzs, j; in PCBDDCBenignDetectSaddlePoint() local
3171 PetscCall(ISGetLocalSize(pcbddc->benign_zerodiag_subs[i], &nzs)); in PCBDDCBenignDetectSaddlePoint()
3173 for (j = 0; j < nzs - 1; j++) nnz[idxs[j]] = 2; /* change on pressures */ in PCBDDCBenignDetectSaddlePoint()
3174 nnz[idxs[nzs - 1]] = nzs; /* last local pressure dof in subdomain */ in PCBDDCBenignDetectSaddlePoint()
3188 PetscInt nzs; in PCBDDCBenignDetectSaddlePoint() local
3190 PetscCall(ISGetLocalSize(pcbddc->benign_zerodiag_subs[s], &nzs)); in PCBDDCBenignDetectSaddlePoint()
3192 for (i = 0; i < nzs - 1; i++) { in PCBDDCBenignDetectSaddlePoint()
3197 cols[1] = idxs[nzs - 1]; in PCBDDCBenignDetectSaddlePoint()
3202 PetscCall(PetscMalloc1(nzs, &array)); in PCBDDCBenignDetectSaddlePoint()
3203 for (i = 0; i < nzs - 1; i++) array[i] = -1.; in PCBDDCBenignDetectSaddlePoint()
3204 array[nzs - 1] = 1.; in PCBDDCBenignDetectSaddlePoint()
3205 …PetscCall(MatSetValues(pcbddc->benign_change, 1, idxs + nzs - 1, nzs, idxs, array, INSERT_VALUES)); in PCBDDCBenignDetectSaddlePoint()
3207 pcbddc->benign_p0_lidx[s] = idxs[nzs - 1]; in PCBDDCBenignDetectSaddlePoint()