Lines Matching refs:benign_n
2189 for (i = 0; i < ctx->benign_n; i++) { in PCBDDCBenignMatMult_Private_Private()
2222 for (i = 0; i < ctx->benign_n; i++) { in PCBDDCBenignMatMult_Private_Private()
2274 …if (!pcbddc->benign_change || !pcbddc->benign_n || pcbddc->benign_change_explicit) PetscFunctionRe… in PCBDDCBenignShellMat()
2286 ctx->benign_n = pcbddc->benign_n; in PCBDDCBenignShellMat()
2295 PetscCall(PetscMalloc1(pcbddc->benign_n, &ctx->benign_zerodiag_subs)); in PCBDDCBenignShellMat()
2296 …for (i = 0; i < pcbddc->benign_n; i++) PetscCall(ISGlobalToLocalMappingApplyIS(N_to_D, IS_GTOLM_DR… in PCBDDCBenignShellMat()
2322 for (i = 0; i < ctx->benign_n; i++) PetscCall(ISDestroy(&ctx->benign_zerodiag_subs[i])); in PCBDDCBenignShellMat()
2340 PetscCall(MatZeroRowsColumns(An, pcbddc->benign_n, pcbddc->benign_p0_lidx, 1.0, NULL, NULL)); in PCBDDCBenignProject()
2779 for (i = 0; i < pcbddc->benign_n; i++) pcbddc->benign_p0[i] = -PetscGlobalRank - i; in PCBDDCBenignCheck()
2781 for (i = 0; i < pcbddc->benign_n; i++) pcbddc->benign_p0[i] = 1; in PCBDDCBenignCheck()
2783 for (i = 0; i < pcbddc->benign_n; i++) { in PCBDDCBenignCheck()
2795 PetscInt nz, n, benign_n, bsp = 1; in PCBDDCBenignDetectSaddlePoint() local
2803 for (n = 0; n < pcbddc->benign_n; n++) PetscCall(ISDestroy(&pcbddc->benign_zerodiag_subs[n])); in PCBDDCBenignDetectSaddlePoint()
2891 benign_n = 0; in PCBDDCBenignDetectSaddlePoint()
3007 … PetscCall(ISLocalToGlobalMappingApplyIS(l2g, t_zerodiag_subs, &zerodiag_subs[benign_n])); in PCBDDCBenignDetectSaddlePoint()
3008 benign_n++; in PCBDDCBenignDetectSaddlePoint()
3033 benign_n = 1; in PCBDDCBenignDetectSaddlePoint()
3034 PetscCall(PetscMalloc1(benign_n, &zerodiag_subs)); in PCBDDCBenignDetectSaddlePoint()
3043 if (!benign_n) { in PCBDDCBenignDetectSaddlePoint()
3057 if (benign_n == 1) { in PCBDDCBenignDetectSaddlePoint()
3064 for (i = 0; i < benign_n; i++) { in PCBDDCBenignDetectSaddlePoint()
3071 for (i = 0; i < benign_n; i++) { in PCBDDCBenignDetectSaddlePoint()
3147 pcbddc->benign_n = benign_n; in PCBDDCBenignDetectSaddlePoint()
3151 have_null = (PetscBool)(!!pcbddc->benign_n); in PCBDDCBenignDetectSaddlePoint()
3157 if (pcbddc->benign_n) { in PCBDDCBenignDetectSaddlePoint()
3167 for (i = 0; i < pcbddc->benign_n; i++) { in PCBDDCBenignDetectSaddlePoint()
3183 …tscCall(PetscMalloc3(pcbddc->benign_n, &pcbddc->benign_p0_lidx, pcbddc->benign_n, &pcbddc->benign_… in PCBDDCBenignDetectSaddlePoint()
3185 for (s = 0; s < pcbddc->benign_n; s++) { in PCBDDCBenignDetectSaddlePoint()
3224 …PetscCall(ISLocalToGlobalMappingApply(matis->rmapping, pcbddc->benign_n, pcbddc->benign_p0_lidx, p… in PCBDDCBenignDetectSaddlePoint()
3238 …PetscCall(PetscSFSetGraphLayout(pcbddc->benign_sf, pc->pmat->rmap, pcbddc->benign_n, NULL, PETSC_O… in PCBDDCBenignGetOrSetP0()
3264 if (!pcbddc->benign_n) PetscFunctionReturn(PETSC_SUCCESS); in PCBDDCBenignPopOrPushB0()
3273 …PetscCall(ISCreateGeneral(PETSC_COMM_SELF, pcbddc->benign_n, pcbddc->benign_p0_lidx, PETSC_COPY_VA… in PCBDDCBenignPopOrPushB0()
3291 PetscCall(MatSetSizes(pcbddc->benign_B0, pcbddc->benign_n, n, PETSC_DECIDE, PETSC_DECIDE)); in PCBDDCBenignPopOrPushB0()
3292 PetscCall(PetscMalloc1(pcbddc->benign_n, &nnz)); in PCBDDCBenignPopOrPushB0()
3293 for (i = 0; i < pcbddc->benign_n; i++) { in PCBDDCBenignPopOrPushB0()
3302 for (i = 0; i < pcbddc->benign_n; i++) { in PCBDDCBenignPopOrPushB0()
3335 for (PetscInt i = 0; i < pcbddc->benign_n; i++) { in PCBDDCBenignPopOrPushB0()
4058 for (i = 0; i < pcbddc->benign_n; i++) PetscCall(ISDestroy(&pcbddc->benign_zerodiag_subs[i])); in PCBDDCResetSolvers()
4073 n_constraints = pcbddc->local_primal_size - pcbddc->benign_n - pcbddc->n_vertices; in PCBDDCSetUpLocalWorkVectors()
4311 …PetscCheck(pcbddc->symmetric_primal || !pcbddc->benign_n, PETSC_COMM_SELF, PETSC_ERR_SUP, "Non-sym… in PCBDDCSetUpCorrection()
4316 n_constraints = pcbddc->local_primal_size - pcbddc->benign_n - n_vertices; in PCBDDCSetUpCorrection()
4331 PetscCheck(!pcbddc->benign_n, PETSC_COMM_SELF, PETSC_ERR_SUP, "Not yet implemented"); in PCBDDCSetUpCorrection()
4477 need_benign_correction = (PetscBool)(!!reuse_solver->benign_n); in PCBDDCSetUpCorrection()
4559 PetscCall(PetscArrayzero(reuse_solver->benign_save_vals, pcbddc->benign_n)); in PCBDDCSetUpCorrection()
4697 if (pcbddc->benign_n && (pcbddc->switch_static || pcbddc->dbg_flag)) { in PCBDDCSetUpCorrection()
4701 PetscCall(PetscMalloc1(pcbddc->benign_n, &p0_lidx_I)); in PCBDDCSetUpCorrection()
4702 …for (i = 0; i < pcbddc->benign_n; i++) PetscCall(PetscFindInt(pcbddc->benign_p0_lidx[i], pcis->n -… in PCBDDCSetUpCorrection()
4743 PetscCall(PetscMalloc1(pcbddc->benign_n, &idxs_p0)); in PCBDDCSetUpCorrection()
4745 …PetscCall(ISGlobalToLocalMappingApply(RtoN, IS_GTOLM_DROP, pcbddc->benign_n, pcbddc->benign_p0_lid… in PCBDDCSetUpCorrection()
4746 …pcbddc->benign_n, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Error in R numbering for benign p0! %" PetscIn… in PCBDDCSetUpCorrection()
4791 for (i = 0; i < reuse_solver->benign_n; i++) { in PCBDDCSetUpCorrection()
4972 for (i = 0; i < reuse_solver->benign_n; i++) { in PCBDDCSetUpCorrection()
5021 if (pcbddc->benign_n) { in PCBDDCSetUpCorrection()
5022 …< n_vertices; i++) PetscCall(MatSetValues(pcbddc->coarse_phi_D, pcbddc->benign_n, p0_lidx_I, 1, &i… in PCBDDCSetUpCorrection()
5073 if (pcbddc->benign_n) { in PCBDDCSetUpCorrection()
5074 …for (i = 0; i < n_constraints; i++) PetscCall(MatSetValues(B2, pcbddc->benign_n, p0_lidx_I, 1, &i,… in PCBDDCSetUpCorrection()
5096 if (pcbddc->benign_n) { in PCBDDCSetUpCorrection()
5102 PetscCall(ISCreateStride(PETSC_COMM_SELF, pcbddc->benign_n, 0, 1, &is_dummy)); in PCBDDCSetUpCorrection()
5108 for (j = 0; j < pcbddc->benign_n; j++) { in PCBDDCSetUpCorrection()
5109 PetscInt primal_idx = pcbddc->local_primal_size - pcbddc->benign_n + j; in PCBDDCSetUpCorrection()
5111 …PetscCall(MatSetValue(*coarse_submat, primal_idx, i, data[i * pcbddc->benign_n + j], INSERT_VALUES… in PCBDDCSetUpCorrection()
5112 …PetscCall(MatSetValue(*coarse_submat, i, primal_idx, data[i * pcbddc->benign_n + j], INSERT_VALUES… in PCBDDCSetUpCorrection()
5265 if (pcbddc->benign_n && !pcbddc->benign_change_explicit) { in PCBDDCSetUpCorrection()
5315 if (pcbddc->benign_n) { in PCBDDCSetUpCorrection()
5321 PetscCall(ISCreateStride(PETSC_COMM_SELF, pcbddc->benign_n, 0, 1, &is_dummy)); in PCBDDCSetUpCorrection()
5327 for (j = 0; j < pcbddc->benign_n; j++) { in PCBDDCSetUpCorrection()
5328 PetscInt primal_idx = pcbddc->local_primal_size - pcbddc->benign_n + j; in PCBDDCSetUpCorrection()
5330 data[primal_idx * pcbddc->local_primal_size + i] += data2[i * pcbddc->benign_n + j]; in PCBDDCSetUpCorrection()
5331 data[i * pcbddc->local_primal_size + primal_idx] += data2[i * pcbddc->benign_n + j]; in PCBDDCSetUpCorrection()
5346 …PetscCall(ISCreateStride(PETSC_COMM_SELF, pcbddc->local_primal_size - pcbddc->benign_n, 0, 1, &is_… in PCBDDCSetUpCorrection()
5348 if (!pcbddc->benign_n) { /* TODO: add benign case */ in PCBDDCSetUpCorrection()
5354 …teSeqDense(PETSC_COMM_SELF, pcis->n_B, pcbddc->local_primal_size - pcbddc->benign_n, data, &tmat)); in PCBDDCSetUpCorrection()
5811 … pcbddc->local_primal_size - n_vertices - pcbddc->benign_n, pcbddc->local_primal_size)); in PCBDDCSetUpLocalScatters()
5990 if (pcbddc->benign_n && !pcbddc->benign_change_explicit) { in PCBDDCSetUpLocalSolvers()
6137 … if (reuse_solver->benign_n) { /* we are not using the explicit change of basis on the pressures */ in PCBDDCSetUpLocalSolvers()
6358 if (pcbddc->benign_n) { in PCBDDCApplyInterfacePreconditioner()
6363 …for (j = 0; j < pcbddc->benign_n; j++) array[pcbddc->local_primal_size - pcbddc->benign_n + j] += … in PCBDDCApplyInterfacePreconditioner()
6446 if (pcbddc->benign_n) { in PCBDDCApplyInterfacePreconditioner()
6451 …for (j = 0; j < pcbddc->benign_n; j++) pcbddc->benign_p0[j] = array[pcbddc->local_primal_size - pc… in PCBDDCApplyInterfacePreconditioner()
6954 …PetscCall(PetscMalloc1(pcbddc->local_primal_size + pcbddc->benign_n, &pcbddc->primal_indices_local… in PCBDDCConstraintsSetUp()
6993 …>local_primal_size_cc + pcbddc->benign_n, &pcbddc->local_primal_ref_node, pcbddc->local_primal_siz… in PCBDDCConstraintsSetUp()
7508 for (i = 0; i < pcbddc->benign_n; i++) { in PCBDDCConstraintsSetUp()
9129 …for (i = 0; i < pcbddc->benign_n; i++) PetscCall(VecSetValue(pcbddc->vec1_P, pcbddc->local_primal_… in PCBDDCSetUpCoarseSolver()
9472 PetscInt benign_n; in PCBDDCSetUpSubSchurs() local
9484 benign_n = pcbddc->benign_n; in PCBDDCSetUpSubSchurs()
9486 benign_n = 0; in PCBDDCSetUpSubSchurs()
9514 …pcbddc->adaptive_selection, reuse_solvers, pcbddc->benign_saddle_point, benign_n, pcbddc->benign_p… in PCBDDCSetUpSubSchurs()
9518 …_layers, scaling, pcbddc->adaptive_selection, reuse_solvers, pcbddc->benign_saddle_point, benign_n, in PCBDDCSetUpSubSchurs()