Lines Matching refs:pcis
2136 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCBenignRemoveInterior() local
2149 PetscCall(VecScatterBegin(pcis->global_to_D, r, pcis->vec1_D, INSERT_VALUES, SCATTER_FORWARD)); in PCBDDCBenignRemoveInterior()
2150 PetscCall(VecScatterEnd(pcis->global_to_D, r, pcis->vec1_D, INSERT_VALUES, SCATTER_FORWARD)); in PCBDDCBenignRemoveInterior()
2152 PetscCall(KSPSolve(pcbddc->ksp_D, pcis->vec1_D, pcis->vec2_D)); in PCBDDCBenignRemoveInterior()
2154 PetscCall(KSPCheckSolve(pcbddc->ksp_D, pc, pcis->vec2_D)); in PCBDDCBenignRemoveInterior()
2156 PetscCall(VecScatterBegin(pcis->global_to_D, pcis->vec2_D, z, INSERT_VALUES, SCATTER_REVERSE)); in PCBDDCBenignRemoveInterior()
2157 PetscCall(VecScatterEnd(pcis->global_to_D, pcis->vec2_D, z, INSERT_VALUES, SCATTER_REVERSE)); in PCBDDCBenignRemoveInterior()
2263 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCBenignShellMat() local
2275 PetscCall(PetscMalloc1(pcis->n, &work)); in PCBDDCBenignShellMat()
2277 PetscCall(MatSetSizes(A_IB, pcis->n - pcis->n_B, pcis->n_B, PETSC_DECIDE, PETSC_DECIDE)); in PCBDDCBenignShellMat()
2294 PetscCall(ISLocalToGlobalMappingCreateIS(pcis->is_I_local, &N_to_D)); in PCBDDCBenignShellMat()
2300 ctx->A = pcis->A_IB; in PCBDDCBenignShellMat()
2305 pcis->A_IB = A_IB; in PCBDDCBenignShellMat()
2309 pcbddc->benign_original_mat = pcis->A_BI; in PCBDDCBenignShellMat()
2310 pcis->A_BI = A_BI; in PCBDDCBenignShellMat()
2313 PetscCall(MatShellGetContext(pcis->A_IB, &ctx)); in PCBDDCBenignShellMat()
2314 PetscCall(MatDestroy(&pcis->A_IB)); in PCBDDCBenignShellMat()
2315 pcis->A_IB = ctx->A; in PCBDDCBenignShellMat()
2317 PetscCall(MatDestroy(&pcis->A_BI)); in PCBDDCBenignShellMat()
2318 pcis->A_BI = pcbddc->benign_original_mat; in PCBDDCBenignShellMat()
2714 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCBenignCheck() local
2728 PetscCall(VecSet(pcis->vec1_N, 0.)); in PCBDDCBenignCheck()
2729 PetscCall(PetscMalloc1(pcis->n, &vals)); in PCBDDCBenignCheck()
2733 PetscCall(VecSetValues(pcis->vec1_N, nz, idxs, vals, INSERT_VALUES)); in PCBDDCBenignCheck()
2734 PetscCall(VecAssemblyBegin(pcis->vec1_N)); in PCBDDCBenignCheck()
2735 PetscCall(VecAssemblyEnd(pcis->vec1_N)); in PCBDDCBenignCheck()
2737 PetscCall(VecSetRandom(pcis->vec2_N, NULL)); in PCBDDCBenignCheck()
2739 PetscCall(VecSetValues(pcis->vec2_N, nz, idxs, vals, INSERT_VALUES)); in PCBDDCBenignCheck()
2741 PetscCall(ISGetIndices(pcis->is_B_local, &idxs)); in PCBDDCBenignCheck()
2742 for (i = 0; i < pcis->n_B; i++) vals[i] = 0.; in PCBDDCBenignCheck()
2743 PetscCall(VecSetValues(pcis->vec2_N, pcis->n_B, idxs, vals, INSERT_VALUES)); in PCBDDCBenignCheck()
2744 PetscCall(ISRestoreIndices(pcis->is_B_local, &idxs)); in PCBDDCBenignCheck()
2751 PetscCall(VecSetValues(pcis->vec2_N, n, idxs, vals, INSERT_VALUES)); in PCBDDCBenignCheck()
2754 PetscCall(VecAssemblyBegin(pcis->vec2_N)); in PCBDDCBenignCheck()
2755 PetscCall(VecAssemblyEnd(pcis->vec2_N)); in PCBDDCBenignCheck()
2756 PetscCall(VecDuplicate(pcis->vec1_N, &vec3_N)); in PCBDDCBenignCheck()
2759 PetscCall(MatMult(A, pcis->vec1_N, vec3_N)); in PCBDDCBenignCheck()
2760 PetscCall(VecDot(vec3_N, pcis->vec2_N, &vals[0])); in PCBDDCBenignCheck()
2766 PetscCall(PetscCalloc1(pcis->n, &count)); in PCBDDCBenignCheck()
2767 PetscCall(ISGetIndices(pcis->is_B_local, &idxs)); in PCBDDCBenignCheck()
2768 for (i = 0; i < pcis->n_B; i++) count[idxs[i]]++; in PCBDDCBenignCheck()
2769 PetscCall(ISRestoreIndices(pcis->is_B_local, &idxs)); in PCBDDCBenignCheck()
2778 PetscCall(VecSetRandom(pcis->vec1_global, NULL)); in PCBDDCBenignCheck()
2780 PetscCall(PCBDDCBenignGetOrSetP0(pc, pcis->vec1_global, PETSC_FALSE)); in PCBDDCBenignCheck()
2782 PetscCall(PCBDDCBenignGetOrSetP0(pc, pcis->vec1_global, PETSC_TRUE)); in PCBDDCBenignCheck()
4068 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCSetUpLocalWorkVectors() local
4074 n_R = pcis->n - pcbddc->n_vertices; in PCBDDCSetUpLocalWorkVectors()
4075 PetscCall(VecGetType(pcis->vec1_N, &impVecType)); in PCBDDCSetUpLocalWorkVectors()
4083 PetscCall(VecCreate(PetscObjectComm((PetscObject)pcis->vec1_N), &pcbddc->vec1_R)); in PCBDDCSetUpLocalWorkVectors()
4093 PetscCall(VecCreate(PetscObjectComm((PetscObject)pcis->vec1_N), &pcbddc->vec1_P)); in PCBDDCSetUpLocalWorkVectors()
4102 PetscCall(VecCreate(PetscObjectComm((PetscObject)pcis->vec1_N), &pcbddc->vec1_C)); in PCBDDCSetUpLocalWorkVectors()
4278 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCSetUpCorrection() local
4317 n_B = pcis->n_B; in PCBDDCSetUpCorrection()
4318 n_D = pcis->n - n_B; in PCBDDCSetUpCorrection()
4319 n_R = pcis->n - n_vertices; in PCBDDCSetUpCorrection()
4323 …PetscCall(ISGlobalToLocalMappingApply(pcis->BtoNmap, IS_GTOLM_DROP, n_vertices, pcbddc->local_prim… in PCBDDCSetUpCorrection()
4422 PetscCall(ISGetIndices(pcis->is_B_local, &idx)); in PCBDDCSetUpCorrection()
4431 PetscCall(ISRestoreIndices(pcis->is_B_local, &idx)); in PCBDDCSetUpCorrection()
4492 PetscCall(VecCreate(PetscObjectComm((PetscObject)pcis->vec1_N), &dummy_vec)); in PCBDDCSetUpCorrection()
4494 PetscCall(VecSetType(dummy_vec, ((PetscObject)pcis->vec1_N)->type_name)); in PCBDDCSetUpCorrection()
4512 …PetscCall(MatCreateSubMatrix(pcbddc->ConstraintMatrix, is_C, pcis->is_B_local, MAT_INITIAL_MATRIX,… in PCBDDCSetUpCorrection()
4674 PetscCall(ISComplement(tis, 0, pcis->n, &is_aux)); in PCBDDCSetUpCorrection()
4677 PetscCall(ISComplement(pcbddc->is_R_local, 0, pcis->n, &is_aux)); in PCBDDCSetUpCorrection()
4700 PetscCall(ISGetIndices(pcis->is_I_local, &idxs)); in PCBDDCSetUpCorrection()
4702 …dc->benign_n; i++) PetscCall(PetscFindInt(pcbddc->benign_p0_lidx[i], pcis->n - pcis->n_B, idxs, &p… in PCBDDCSetUpCorrection()
4703 PetscCall(ISRestoreIndices(pcis->is_I_local, &idxs)); in PCBDDCSetUpCorrection()
5088 …PetscCall(MatCreateSubMatrix(T, pcis->is_B_local, NULL, MAT_INITIAL_MATRIX, &pcbddc->coarse_phi_B)… in PCBDDCSetUpCorrection()
5089 …if (pcbddc->switch_static || pcbddc->dbg_flag) PetscCall(MatCreateSubMatrix(T, pcis->is_I_local, N… in PCBDDCSetUpCorrection()
5103 …PetscCall(MatCreateSubMatrix(pcbddc->benign_B0, is_dummy, pcis->is_B_local, MAT_INITIAL_MATRIX, &B… in PCBDDCSetUpCorrection()
5268 … PetscCall(MatCreateSubMatrix(A, pcis->is_I_local, pcis->is_I_local, MAT_INITIAL_MATRIX, &A_II)); in PCBDDCSetUpCorrection()
5269 … PetscCall(MatCreateSubMatrix(A, pcis->is_I_local, pcis->is_B_local, MAT_INITIAL_MATRIX, &A_IB)); in PCBDDCSetUpCorrection()
5270 … PetscCall(MatCreateSubMatrix(A, pcis->is_B_local, pcis->is_I_local, MAT_INITIAL_MATRIX, &A_BI)); in PCBDDCSetUpCorrection()
5271 … PetscCall(MatCreateSubMatrix(A, pcis->is_B_local, pcis->is_B_local, MAT_INITIAL_MATRIX, &A_BB)); in PCBDDCSetUpCorrection()
5274 PetscCall(MatConvert(pcis->A_II, checkmattype, MAT_INITIAL_MATRIX, &A_II)); in PCBDDCSetUpCorrection()
5275 PetscCall(MatConvert(pcis->A_IB, checkmattype, MAT_INITIAL_MATRIX, &A_IB)); in PCBDDCSetUpCorrection()
5276 PetscCall(MatConvert(pcis->A_BI, checkmattype, MAT_INITIAL_MATRIX, &A_BI)); in PCBDDCSetUpCorrection()
5277 PetscCall(MatConvert(pcis->A_BB, checkmattype, MAT_INITIAL_MATRIX, &A_BB)); in PCBDDCSetUpCorrection()
5322 …PetscCall(MatCreateSubMatrix(pcbddc->benign_B0, is_dummy, pcis->is_B_local, MAT_INITIAL_MATRIX, &B… in PCBDDCSetUpCorrection()
5347 …PetscCall(MatCreateSubMatrix(pcbddc->ConstraintMatrix, is_dummy, pcis->is_B_local, MAT_INITIAL_MAT… in PCBDDCSetUpCorrection()
5354 …PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, pcis->n_B, pcbddc->local_primal_size - pcbddc->benign… in PCBDDCSetUpCorrection()
5422 PetscCall(PetscObjectSetName((PetscObject)pcis->is_I_local,"I")); in PCBDDCSetUpCorrection()
5423 PetscCall(ISView(pcis->is_I_local,viewer)); in PCBDDCSetUpCorrection()
5424 PetscCall(PetscObjectSetName((PetscObject)pcis->is_B_local,"B")); in PCBDDCSetUpCorrection()
5425 PetscCall(ISView(pcis->is_B_local,viewer)); in PCBDDCSetUpCorrection()
5437 …PetscCall(PetscObjectTypeCompareAny((PetscObject)pcis->vec1_N, &iscuda, VECCUDA, VECMPICUDA, VECSE… in PCBDDCSetUpCorrection()
5438 …PetscCall(PetscObjectTypeCompareAny((PetscObject)pcis->vec1_N, &iship, VECHIP, VECMPIHIP, VECSEQHI… in PCBDDCSetUpCorrection()
5439 …PetscCall(PetscObjectTypeCompareAny((PetscObject)pcis->vec1_N, &iskokkos, VECKOKKOS, VECMPIKOKKOS,… in PCBDDCSetUpCorrection()
5710 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCSetUpLocalScatters() local
5733 n_B = pcis->n_B; in PCBDDCSetUpLocalScatters()
5734 n_D = pcis->n - n_B; in PCBDDCSetUpLocalScatters()
5741 PetscCall(PetscMalloc1(pcis->n - n_vertices, &idx_R_local)); in PCBDDCSetUpLocalScatters()
5742 PetscCall(PetscBTCreate(pcis->n, &bitmask)); in PCBDDCSetUpLocalScatters()
5745 for (i = 0, n_R = 0; i < pcis->n; i++) { in PCBDDCSetUpLocalScatters()
5762 PetscCall(PetscMalloc1(pcis->n / bs, &vary)); in PCBDDCSetUpLocalScatters()
5763 PetscCall(PetscArrayzero(vary, pcis->n / bs)); in PCBDDCSetUpLocalScatters()
5767 for (i = 0; i < pcis->n / bs; i++) { in PCBDDCSetUpLocalScatters()
5809 …T ", dirichlet_size = %" PetscInt_FMT ", boundary_size = %" PetscInt_FMT "\n", pcis->n, n_D, n_B)); in PCBDDCSetUpLocalScatters()
5821 PetscCall(PetscMalloc1(pcis->n_B - n_vertices, &aux_array1)); in PCBDDCSetUpLocalScatters()
5822 PetscCall(PetscMalloc1(pcis->n_B - n_vertices, &aux_array2)); in PCBDDCSetUpLocalScatters()
5823 PetscCall(ISGetIndices(pcis->is_I_local, (const PetscInt **)&is_indices)); in PCBDDCSetUpLocalScatters()
5825 PetscCall(ISRestoreIndices(pcis->is_I_local, (const PetscInt **)&is_indices)); in PCBDDCSetUpLocalScatters()
5830 PetscCall(ISGetIndices(pcis->is_B_local, (const PetscInt **)&is_indices)); in PCBDDCSetUpLocalScatters()
5834 PetscCall(ISRestoreIndices(pcis->is_B_local, (const PetscInt **)&is_indices)); in PCBDDCSetUpLocalScatters()
5836 PetscCall(VecScatterCreate(pcbddc->vec1_R, is_aux1, pcis->vec1_B, is_aux2, &pcbddc->R_to_B)); in PCBDDCSetUpLocalScatters()
5846 PetscCall(VecScatterCreate(pcbddc->vec1_R, is_aux1, pcis->vec1_D, (IS)0, &pcbddc->R_to_D)); in PCBDDCSetUpLocalScatters()
5858 …PetscCall(VecScatterCreate(pcbddc->vec1_R, tis, pcis->vec1_B, reuse_solver->is_B, &pcbddc->R_to_B)… in PCBDDCSetUpLocalScatters()
5862 PetscCall(VecScatterCreate(pcbddc->vec1_R, tis, pcis->vec1_D, (IS)0, &pcbddc->R_to_D)); in PCBDDCSetUpLocalScatters()
5936 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCSetUpLocalSolvers() local
5995 PetscCall(PCBDDCBenignProject(pc, pcis->is_I_local, pcis->is_I_local, &A_IIn)); in PCBDDCSetUpLocalSolvers()
5996 PetscCall(MatDestroy(&pcis->A_II)); in PCBDDCSetUpLocalSolvers()
5997 pcis->A_II = A_IIn; in PCBDDCSetUpLocalSolvers()
6001 if (isset) PetscCall(MatSetOption(pcis->A_II, MAT_SYMMETRIC, issym)); in PCBDDCSetUpLocalSolvers()
6004 n_D = pcis->n - pcis->n_B; in PCBDDCSetUpLocalSolvers()
6014 PetscCall(PetscObjectTypeCompare((PetscObject)pcis->pA_II, MATSEQSBAIJ, &issbaij)); in PCBDDCSetUpLocalSolvers()
6023 PetscCall(MatSetOptionsPrefix(pcis->pA_II, ((PetscObject)pcbddc->ksp_D)->prefix)); in PCBDDCSetUpLocalSolvers()
6024 PetscCall(MatViewFromOptions(pcis->pA_II, NULL, "-mat_view")); in PCBDDCSetUpLocalSolvers()
6025 PetscCall(KSPSetOperators(pcbddc->ksp_D, pcis->A_II, pcis->pA_II)); in PCBDDCSetUpLocalSolvers()
6028 PetscCall(MatGetNearNullSpace(pcis->pA_II, &nnsp)); in PCBDDCSetUpLocalSolvers()
6030 PetscCall(MatNullSpacePropagateAny_Private(pcbddc->local_mat, pcis->is_I_local, pcis->pA_II)); in PCBDDCSetUpLocalSolvers()
6032 PetscCall(MatGetNearNullSpace(pcis->pA_II, &nnsp)); in PCBDDCSetUpLocalSolvers()
6040 PetscCall(ISGetLocalSize(pcis->is_I_local, &nl)); in PCBDDCSetUpLocalSolvers()
6041 PetscCall(ISGetIndices(pcis->is_I_local, &idxs)); in PCBDDCSetUpLocalSolvers()
6046 PetscCall(ISRestoreIndices(pcis->is_I_local, &idxs)); in PCBDDCSetUpLocalSolvers()
6064 PetscCall(KSPDestroy(&pcis->ksp_D)); in PCBDDCSetUpLocalSolvers()
6065 pcis->ksp_D = pcbddc->ksp_D; in PCBDDCSetUpLocalSolvers()
6227 PetscCall(VecSetRandom(pcis->vec1_D, NULL)); in PCBDDCSetUpLocalSolvers()
6228 PetscCall(MatMult(pcis->A_II, pcis->vec1_D, pcis->vec2_D)); in PCBDDCSetUpLocalSolvers()
6229 PetscCall(KSPSolve(pcbddc->ksp_D, pcis->vec2_D, pcis->vec2_D)); in PCBDDCSetUpLocalSolvers()
6230 PetscCall(KSPCheckSolve(pcbddc->ksp_D, pc, pcis->vec2_D)); in PCBDDCSetUpLocalSolvers()
6231 PetscCall(VecAXPY(pcis->vec1_D, m_one, pcis->vec2_D)); in PCBDDCSetUpLocalSolvers()
6232 PetscCall(VecNorm(pcis->vec1_D, NORM_INFINITY, &value)); in PCBDDCSetUpLocalSolvers()
6340 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCApplyInterfacePreconditioner() local
6347 PetscCall(MatMultTranspose(pcbddc->coarse_phi_B, pcis->vec1_B, pcbddc->vec1_P)); in PCBDDCApplyInterfacePreconditioner()
6348 …if (pcbddc->switch_static) PetscCall(MatMultTransposeAdd(pcbddc->coarse_phi_D, pcis->vec1_D, pcbdd… in PCBDDCApplyInterfacePreconditioner()
6350 PetscCall(MatMultTranspose(pcbddc->coarse_psi_B, pcis->vec1_B, pcbddc->vec1_P)); in PCBDDCApplyInterfacePreconditioner()
6351 …if (pcbddc->switch_static) PetscCall(MatMultTransposeAdd(pcbddc->coarse_psi_D, pcis->vec1_D, pcbdd… in PCBDDCApplyInterfacePreconditioner()
6431 …apply_coarse_only) PetscCall(PCBDDCSolveSubstructureCorrection(pc, pcis->vec1_B, pcis->vec1_D, app… in PCBDDCApplyInterfacePreconditioner()
6439 PetscCall(MatMultAdd(pcbddc->coarse_psi_B, pcbddc->vec1_P, pcis->vec1_B, pcis->vec1_B)); in PCBDDCApplyInterfacePreconditioner()
6440 …ch_static) PetscCall(MatMultAdd(pcbddc->coarse_psi_D, pcbddc->vec1_P, pcis->vec1_D, pcis->vec1_D)); in PCBDDCApplyInterfacePreconditioner()
6442 PetscCall(MatMultAdd(pcbddc->coarse_phi_B, pcbddc->vec1_P, pcis->vec1_B, pcis->vec1_B)); in PCBDDCApplyInterfacePreconditioner()
6443 …ch_static) PetscCall(MatMultAdd(pcbddc->coarse_phi_D, pcbddc->vec1_P, pcis->vec1_D, pcis->vec1_D)); in PCBDDCApplyInterfacePreconditioner()
6456 PetscCall(MatMult(pcbddc->coarse_psi_B, pcbddc->vec1_P, pcis->vec1_B)); in PCBDDCApplyInterfacePreconditioner()
6458 PetscCall(MatMult(pcbddc->coarse_phi_B, pcbddc->vec1_P, pcis->vec1_B)); in PCBDDCApplyInterfacePreconditioner()
6526 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCConstraintsSetUp() local
6668 PetscCall(VecDuplicate(pcis->vec1_N, &localnearnullsp[k])); in PCBDDCConstraintsSetUp()
6958 …PetscCall(ISGlobalToLocalMappingApply(pcis->BtoNmap, IS_GTOLM_DROP, constraints_idxs_ptr[total_cou… in PCBDDCConstraintsSetUp()
6965 …es(pcbddc->ConstraintMatrix, pcbddc->local_primal_size, pcis->n, pcbddc->local_primal_size, pcis->… in PCBDDCConstraintsSetUp()
7075 PetscCall(MatSetSizes(localChangeOfBasisMatrix, pcis->n, pcis->n, pcis->n, pcis->n)); in PCBDDCConstraintsSetUp()
7077 PetscCall(PetscMalloc1(pcis->n, &nnz)); in PCBDDCConstraintsSetUp()
7079 for (i = 0; i < pcis->n; i++) nnz[i] = 1; in PCBDDCConstraintsSetUp()
7104 …for (i = 0; i < pcis->n; i++) PetscCall(MatSetValue(localChangeOfBasisMatrix, i, i, 1.0, INSERT_VA… in PCBDDCConstraintsSetUp()
7173 PetscCall(PetscBTCreate(pcis->n_B, &is_primal)); in PCBDDCConstraintsSetUp()
7175 …PetscCall(ISGlobalToLocalMappingApply(pcis->BtoNmap, IS_GTOLM_DROP, total_primal_vertices, pcbddc-… in PCBDDCConstraintsSetUp()
7344 PetscCall(VecGetSize(pcis->vec1_global, &global_size)); in PCBDDCConstraintsSetUp()
7345 PetscCall(VecGetLocalSize(pcis->vec1_global, &local_size)); in PCBDDCConstraintsSetUp()
7352 PetscCall(VecSet(pcis->vec1_global, 0.0)); in PCBDDCConstraintsSetUp()
7353 PetscCall(VecSet(pcis->vec1_N, 1.0)); in PCBDDCConstraintsSetUp()
7354 …PetscCall(VecScatterBegin(matis->rctx, pcis->vec1_N, pcis->vec1_global, ADD_VALUES, SCATTER_REVERS… in PCBDDCConstraintsSetUp()
7355 …PetscCall(VecScatterEnd(matis->rctx, pcis->vec1_N, pcis->vec1_global, ADD_VALUES, SCATTER_REVERSE)… in PCBDDCConstraintsSetUp()
7356 PetscCall(VecReciprocal(pcis->vec1_global)); in PCBDDCConstraintsSetUp()
7357 PetscCall(MatDiagonalScale(pcbddc->ChangeOfBasisMatrix, pcis->vec1_global, NULL)); in PCBDDCConstraintsSetUp()
7364 PetscCall(VecDuplicate(pcis->vec1_global, &x)); in PCBDDCConstraintsSetUp()
7365 PetscCall(VecDuplicate(pcis->vec1_global, &x_change)); in PCBDDCConstraintsSetUp()
7367 PetscCall(VecCopy(x, pcis->vec1_global)); in PCBDDCConstraintsSetUp()
7368 PetscCall(VecScatterBegin(matis->rctx, x, pcis->vec1_N, INSERT_VALUES, SCATTER_FORWARD)); in PCBDDCConstraintsSetUp()
7369 PetscCall(VecScatterEnd(matis->rctx, x, pcis->vec1_N, INSERT_VALUES, SCATTER_FORWARD)); in PCBDDCConstraintsSetUp()
7370 PetscCall(MatMult(localChangeOfBasisMatrix, pcis->vec1_N, pcis->vec2_N)); in PCBDDCConstraintsSetUp()
7371 PetscCall(VecScatterBegin(matis->rctx, pcis->vec2_N, x, INSERT_VALUES, SCATTER_REVERSE)); in PCBDDCConstraintsSetUp()
7372 PetscCall(VecScatterEnd(matis->rctx, pcis->vec2_N, x, INSERT_VALUES, SCATTER_REVERSE)); in PCBDDCConstraintsSetUp()
7373 PetscCall(MatMult(pcbddc->ChangeOfBasisMatrix, pcis->vec1_global, x_change)); in PCBDDCConstraintsSetUp()
7389 PetscCall(ISLocalToGlobalMappingApplyIS(pcis->BtoNmap, sub_schurs->is_Ej_all, &is_all_N)); in PCBDDCConstraintsSetUp()
7413 PetscCall(VecGetArrayRead(pcis->D, &array)); in PCBDDCConstraintsSetUp()
7424 PetscCall(VecRestoreArrayRead(pcis->D, &array)); in PCBDDCConstraintsSetUp()
7471 PetscCall(VecCopy(matis->counter, pcis->vec1_N)); in PCBDDCConstraintsSetUp()
7472 PetscCall(VecReciprocal(pcis->vec1_N)); in PCBDDCConstraintsSetUp()
7476 PetscCall(MatDiagonalScale(M, pcis->vec1_N, NULL)); in PCBDDCConstraintsSetUp()
7478 PetscCall(MatCreateSeqAIJ(PETSC_COMM_SELF, pcis->n, pcis->n, 1, NULL, &M)); in PCBDDCConstraintsSetUp()
7479 PetscCall(MatDiagonalSet(M, pcis->vec1_N, INSERT_VALUES)); in PCBDDCConstraintsSetUp()
7498 …PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)pc), pcis->n, gidxs, PETSC_COPY_VALUES, &is… in PCBDDCConstraintsSetUp()
7504 …bddc->ChangeOfBasisMatrix && !pcbddc->work_change) PetscCall(VecDuplicate(pcis->vec1_global, &pcbd… in PCBDDCConstraintsSetUp()
8515 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCSetUpCoarseSolver() local
8591 im_active = !!pcis->n; in PCBDDCSetUpCoarseSolver()
8754 … PetscCall(MatCreateSubMatrix(loc_divudotp, dummy, pcis->is_B_local, MAT_INITIAL_MATRIX, &B)); in PCBDDCSetUpCoarseSolver()
9286 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCComputePrimalNumbering() local
9296 PetscCall(ISLocalToGlobalMappingApplyIS(pcis->mapping, subset_n, &subset)); in PCBDDCComputePrimalNumbering()
9365 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCComputeFakeChange() local
9379 pcisf->is_B_local = pcis->is_B_local; in PCBDDCComputeFakeChange()
9380 pcisf->vec1_N = pcis->vec1_N; in PCBDDCComputeFakeChange()
9381 pcisf->BtoNmap = pcis->BtoNmap; in PCBDDCComputeFakeChange()
9382 pcisf->n = pcis->n; in PCBDDCComputeFakeChange()
9383 pcisf->n_B = pcis->n_B; in PCBDDCComputeFakeChange()
9421 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCSetUpSubSchurs() local
9463 …PetscCall(MatCreateSchurComplement(pcis->A_II, pcis->pA_II, pcis->A_IB, pcis->A_BI, pcis->A_BB, &S… in PCBDDCSetUpSubSchurs()
9501 if (!pcbddc->use_deluxe_scaling) scaling = pcis->D; in PCBDDCSetUpSubSchurs()
9534 PC_IS *pcis = (PC_IS *)pc->data; in PCBDDCInitSubSchurs() local
9577 …sInit(pcbddc->sub_schurs, ((PetscObject)pc)->prefix, pcis->is_I_local, pcis->is_B_local, graph, pc… in PCBDDCInitSubSchurs()