Lines Matching refs:idxs
166 const PetscInt *idxs, *ii, *jj, *iit, *jjt; in PCBDDCNedelecSupport() local
257 PetscCall(ISGetIndices(nedfieldlocal, &idxs)); in PCBDDCNedelecSupport()
259 if (PetscRealPart(vals[idxs[i]]) > 2.) eidxs[cum++] = idxs[i]; in PCBDDCNedelecSupport()
261 PetscCall(ISRestoreIndices(nedfieldlocal, &idxs)); in PCBDDCNedelecSupport()
309 PetscCall(ISGetIndices(nedfieldlocal, &idxs)); in PCBDDCNedelecSupport()
310 for (i = 0; i < ne; i++) matis->sf_leafdata[idxs[i]] = 1; in PCBDDCNedelecSupport()
311 PetscCall(ISRestoreIndices(nedfieldlocal, &idxs)); in PCBDDCNedelecSupport()
341 PetscCall(ISGetIndices(nedfieldlocal, &idxs)); in PCBDDCNedelecSupport()
343 if (matis->sf_leafdata[idxs[i]] == 1) tbz[cum++] = i; in PCBDDCNedelecSupport()
344 PetscCall(ISRestoreIndices(nedfieldlocal, &idxs)); in PCBDDCNedelecSupport()
351 PetscCall(ISLocalToGlobalMappingGetIndices(el2g, &idxs)); in PCBDDCNedelecSupport()
352 PetscCall(ISCreateGeneral(comm, ne, idxs, PETSC_USE_POINTER, &lned)); in PCBDDCNedelecSupport()
355 PetscCall(ISLocalToGlobalMappingRestoreIndices(el2g, &idxs)); in PCBDDCNedelecSupport()
429 PetscCall(ISLocalToGlobalMappingGetIndices(vl2g, &idxs)); in PCBDDCNedelecSupport()
430 PetscCall(PetscSFSetGraphLayout(sfv, lGis->cmap, nv, NULL, PETSC_OWN_POINTER, idxs)); in PCBDDCNedelecSupport()
431 PetscCall(ISLocalToGlobalMappingRestoreIndices(vl2g, &idxs)); in PCBDDCNedelecSupport()
478 PetscCall(ISGetIndices(is, &idxs)); in PCBDDCNedelecSupport()
480 if (idxs[i] >= 0 && idxs[i] < ne) { in PCBDDCNedelecSupport()
481 PetscCall(PetscBTSet(btb, idxs[i])); in PCBDDCNedelecSupport()
482 PetscCall(PetscBTSet(btbd, idxs[i])); in PCBDDCNedelecSupport()
485 PetscCall(ISRestoreIndices(is, &idxs)); in PCBDDCNedelecSupport()
497 PetscCall(ISGetIndices(is, &idxs)); in PCBDDCNedelecSupport()
499 if (idxs[i] >= 0 && idxs[i] < ne) PetscCall(PetscBTSet(btb, idxs[i])); in PCBDDCNedelecSupport()
501 PetscCall(ISRestoreIndices(is, &idxs)); in PCBDDCNedelecSupport()
579 PetscCall(ISGetIndices(elements_corners, &idxs)); in PCBDDCNedelecSupport()
580 for (i = 0; i < cum; i++) PetscCall(PetscBTSet(btvcand, idxs[i])); in PCBDDCNedelecSupport()
581 PetscCall(ISRestoreIndices(elements_corners, &idxs)); in PCBDDCNedelecSupport()
617 const PetscInt *idxs; in PCBDDCNedelecSupport() local
622 PetscCall(ISGetIndices(pcbddc->local_subs[i], &idxs)); in PCBDDCNedelecSupport()
624 const PetscInt e = idxs[j]; in PCBDDCNedelecSupport()
629 PetscCall(ISRestoreIndices(pcbddc->local_subs[i], &idxs)); in PCBDDCNedelecSupport()
881 PetscCall(ISGetIndices(nedfieldlocal, &idxs)); in PCBDDCNedelecSupport()
883 PetscCall(PetscBTSet(btf, idxs[i])); in PCBDDCNedelecSupport()
884 iia[idxs[i] + 1] = ii[i + 1] - ii[i]; in PCBDDCNedelecSupport()
900 PetscInt e = idxs[i]; in PCBDDCNedelecSupport()
904 PetscCall(ISRestoreIndices(nedfieldlocal, &idxs)); in PCBDDCNedelecSupport()
935 PetscCall(ISGetIndices(primals, &idxs)); in PCBDDCNedelecSupport()
936 for (i = 0; i < cum; i++) marks[idxs[i]] = nee + 1; in PCBDDCNedelecSupport()
937 PetscCall(ISRestoreIndices(primals, &idxs)); in PCBDDCNedelecSupport()
948 PetscCall(ISGetIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
949 for (j = 0; j < size; j++) marks[idxs[j]] = mark; in PCBDDCNedelecSupport()
950 PetscCall(ISRestoreIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
963 PetscCall(ISGetIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
969 PetscInt k, ee = idxs[j]; in PCBDDCNedelecSupport()
996 PetscCall(ISRestoreIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1004 PetscCall(ISGetIndices(primals, &idxs)); in PCBDDCNedelecSupport()
1005 for (i = 0; i < cum; i++) marks[idxs[i]] = 0; in PCBDDCNedelecSupport()
1006 PetscCall(ISRestoreIndices(primals, &idxs)); in PCBDDCNedelecSupport()
1030 PetscCall(ISGetIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1033 PetscInt k, ee = idxs[j]; in PCBDDCNedelecSupport()
1040 PetscCall(ISRestoreIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1065 PetscCall(ISGetIndices(primals, &idxs)); in PCBDDCNedelecSupport()
1066 PetscCall(PetscArraycpy(newprimals, idxs, cum)); in PCBDDCNedelecSupport()
1067 PetscCall(ISRestoreIndices(primals, &idxs)); in PCBDDCNedelecSupport()
1076 PetscCall(ISGetIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1079 PetscInt k, ee = idxs[j]; in PCBDDCNedelecSupport()
1107 PetscInt k, ee = idxs[0], *tmarks; in PCBDDCNedelecSupport()
1118 if (tmarks[idxs[j]] > 1) { in PCBDDCNedelecSupport()
1119 …) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Edge dof set to primal %" PetscInt_FMT "\n", idxs[j])); in PCBDDCNedelecSupport()
1120 newprimals[cum++] = idxs[j]; in PCBDDCNedelecSupport()
1125 PetscCall(ISRestoreIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1162 PetscCall(ISGetIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1163 for (j = 0; j < size; j++) marks[idxs[j]] = mark; in PCBDDCNedelecSupport()
1164 PetscCall(ISRestoreIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1180 PetscCall(ISGetIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1182 PetscInt k, ee = idxs[j]; in PCBDDCNedelecSupport()
1186 PetscCall(ISRestoreIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1221 PetscCall(ISGetIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1224 PetscInt k, ee = idxs[j]; in PCBDDCNedelecSupport()
1236 PetscCall(ISLocalToGlobalMappingApply(fl2g, 1, idxs, &e)); in PCBDDCNedelecSupport()
1238 e = idxs[0]; in PCBDDCNedelecSupport()
1240 …dge %" PetscInt_FMT " (astart %" PetscInt_FMT ", estart %" PetscInt_FMT ")", found, i, e, idxs[0]); in PCBDDCNedelecSupport()
1250 cedges[i] = idxs[size - 1]; in PCBDDCNedelecSupport()
1251 PetscCall(ISRestoreIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1349 const PetscInt *idxs; in PCBDDCNedelecSupport() local
1353 PetscCall(ISGetIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1354 PetscCall(VecSetValuesLocal(w, nl, idxs, NULL, INSERT_VALUES)); in PCBDDCNedelecSupport()
1355 PetscCall(ISRestoreIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1433 const PetscInt *idxs; in PCBDDCNedelecSupport() local
1437 PetscCall(ISGetIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1438 PetscCall(PetscArraycpy(eedgesidxs + cum, idxs, j - 1)); /* last on the edge is primal */ in PCBDDCNedelecSupport()
1439 PetscCall(ISRestoreIndices(eedges[i], &idxs)); in PCBDDCNedelecSupport()
1601 PetscInt ncc, *idxs; in PCBDDCNedelecSupport() local
1605 PetscCall(ISGetIndices(pcbddc->nedclocal, (const PetscInt **)&idxs)); in PCBDDCNedelecSupport()
1608 idxs = cedges; in PCBDDCNedelecSupport()
1631 if (pcbddc->nedclocal) PetscCall(ISRestoreIndices(pcbddc->nedclocal, (const PetscInt **)&idxs)); in PCBDDCNedelecSupport()
1760 PetscInt n, *mmask, *mask, *idxs, nmr, nr; in PCBDDCComputeNoNetFlux() local
1768 PetscCall(PetscCalloc3(nmr, &mmask, n, &mask, n, &idxs)); in PCBDDCComputeNoNetFlux()
1780 idxs[i] = i; in PCBDDCComputeNoNetFlux()
1782 PetscCall(VecSetValuesLocal(quad_vec, n, idxs, array, ADD_VALUES)); in PCBDDCComputeNoNetFlux()
1784 PetscCall(PetscFree3(mmask, mask, idxs)); in PCBDDCComputeNoNetFlux()
2100 const PetscInt *idxs; in PCBDDCConsistencyCheckIS() local
2113 PetscCall(ISGetIndices(*is, &idxs)); in PCBDDCConsistencyCheckIS()
2115 if (-1 < idxs[i] && idxs[i] < n) matis->sf_leafdata[idxs[i]] = 1; in PCBDDCConsistencyCheckIS()
2116 PetscCall(ISRestoreIndices(*is, &idxs)); in PCBDDCConsistencyCheckIS()
2191 const PetscInt *idxs; in PCBDDCBenignMatMult_Private_Private() local
2194 PetscCall(ISGetIndices(ctx->benign_zerodiag_subs[i], &idxs)); in PCBDDCBenignMatMult_Private_Private()
2197 val = ctx->work[idxs[nz - 1]]; in PCBDDCBenignMatMult_Private_Private()
2199 sum += ctx->work[idxs[j]]; in PCBDDCBenignMatMult_Private_Private()
2200 ctx->work[idxs[j]] += val; in PCBDDCBenignMatMult_Private_Private()
2203 for (j = 0; j < nz - 1; j++) sum += ctx->work[idxs[j]]; in PCBDDCBenignMatMult_Private_Private()
2205 ctx->work[idxs[nz - 1]] -= sum; in PCBDDCBenignMatMult_Private_Private()
2206 PetscCall(ISRestoreIndices(ctx->benign_zerodiag_subs[i], &idxs)); in PCBDDCBenignMatMult_Private_Private()
2224 const PetscInt *idxs; in PCBDDCBenignMatMult_Private_Private() local
2227 PetscCall(ISGetIndices(ctx->benign_zerodiag_subs[i], &idxs)); in PCBDDCBenignMatMult_Private_Private()
2228 val = -ay[idxs[nz - 1]]; in PCBDDCBenignMatMult_Private_Private()
2232 sum += ay[idxs[j]]; in PCBDDCBenignMatMult_Private_Private()
2233 ay[idxs[j]] += val; in PCBDDCBenignMatMult_Private_Private()
2235 ay[idxs[nz - 1]] += sum; in PCBDDCBenignMatMult_Private_Private()
2237 for (j = 0; j < nz - 1; j++) ay[idxs[j]] += val; in PCBDDCBenignMatMult_Private_Private()
2238 ay[idxs[nz - 1]] = 0.; in PCBDDCBenignMatMult_Private_Private()
2240 PetscCall(ISRestoreIndices(ctx->benign_zerodiag_subs[i], &idxs)); in PCBDDCBenignMatMult_Private_Private()
2724 const PetscInt *idxs; in PCBDDCBenignCheck() local
2731 PetscCall(ISGetIndices(zerodiag, &idxs)); in PCBDDCBenignCheck()
2733 PetscCall(VecSetValues(pcis->vec1_N, nz, idxs, vals, INSERT_VALUES)); in PCBDDCBenignCheck()
2739 PetscCall(VecSetValues(pcis->vec2_N, nz, idxs, vals, INSERT_VALUES)); in PCBDDCBenignCheck()
2740 PetscCall(ISRestoreIndices(zerodiag, &idxs)); in PCBDDCBenignCheck()
2741 PetscCall(ISGetIndices(pcis->is_B_local, &idxs)); 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()
2749 PetscCall(ISGetIndices(dirIS, &idxs)); in PCBDDCBenignCheck()
2751 PetscCall(VecSetValues(pcis->vec2_N, n, idxs, vals, INSERT_VALUES)); in PCBDDCBenignCheck()
2752 PetscCall(ISRestoreIndices(dirIS, &idxs)); 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()
2770 PetscCall(ISGetIndices(zerodiag, &idxs)); in PCBDDCBenignCheck()
2771 …ck(!count[idxs[i]], PETSC_COMM_SELF, PETSC_ERR_SUP, "Benign trick can not be applied! pressure dof… in PCBDDCBenignCheck()
2772 PetscCall(ISRestoreIndices(zerodiag, &idxs)); in PCBDDCBenignCheck()
2840 const PetscInt *idxs; in PCBDDCBenignDetectSaddlePoint() local
2845 PetscCall(ISGetIndices(pcbddc->ISForDofsLocal[pp[p]], &idxs)); in PCBDDCBenignDetectSaddlePoint()
2850 for (i = 0; i < npl / bs; i++) bidxs[i] = idxs[bs * i + b]; in PCBDDCBenignDetectSaddlePoint()
2855 PetscCall(ISRestoreIndices(pcbddc->ISForDofsLocal[pp[p]], &idxs)); in PCBDDCBenignDetectSaddlePoint()
2912 const PetscInt *idxs; in PCBDDCBenignDetectSaddlePoint() local
2922 PetscCall(ISGetIndices(zerodiag, &idxs)); in PCBDDCBenignDetectSaddlePoint()
2926 for (j = 0; j < nl; j++) array[idxs[j]] = 1.; in PCBDDCBenignDetectSaddlePoint()
2931 for (j = 0; j < nl; j++) array[idxs[j]] = 0.; in PCBDDCBenignDetectSaddlePoint()
2933 PetscCall(ISRestoreIndices(zerodiag, &idxs)); in PCBDDCBenignDetectSaddlePoint()
2968 PetscCall(ISGetIndices(subs[i], &idxs)); in PCBDDCBenignDetectSaddlePoint()
2970 for (j = 0; j < nl; j++) array[idxs[j]] = 1.; in PCBDDCBenignDetectSaddlePoint()
2972 PetscCall(ISRestoreIndices(subs[i], &idxs)); in PCBDDCBenignDetectSaddlePoint()
2986 const PetscInt *idxs; in PCBDDCBenignDetectSaddlePoint() local
2989 PetscCall(ISGetIndices(pcbddc->NeumannBoundariesLocal, &idxs)); in PCBDDCBenignDetectSaddlePoint()
2990 PetscCall(ISGlobalToLocalMappingApply(l2g, IS_GTOLM_DROP, nneu, idxs, &nzb, NULL)); in PCBDDCBenignDetectSaddlePoint()
2991 PetscCall(ISRestoreIndices(pcbddc->NeumannBoundariesLocal, &idxs)); in PCBDDCBenignDetectSaddlePoint()
3072 PetscInt ns, *idxs; in PCBDDCBenignDetectSaddlePoint() local
3074 PetscCall(ISGetIndices(zerodiag_subs[i], (const PetscInt **)&idxs)); in PCBDDCBenignDetectSaddlePoint()
3075 PetscCall(PetscArraycpy(new_idxs + nzn, idxs, ns)); in PCBDDCBenignDetectSaddlePoint()
3076 PetscCall(ISRestoreIndices(zerodiag_subs[i], (const PetscInt **)&idxs)); in PCBDDCBenignDetectSaddlePoint()
3168 const PetscInt *idxs; in PCBDDCBenignDetectSaddlePoint() local
3172 PetscCall(ISGetIndices(pcbddc->benign_zerodiag_subs[i], &idxs)); 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()
3175 PetscCall(ISRestoreIndices(pcbddc->benign_zerodiag_subs[i], &idxs)); in PCBDDCBenignDetectSaddlePoint()
3187 const PetscInt *idxs; in PCBDDCBenignDetectSaddlePoint() local
3191 PetscCall(ISGetIndices(pcbddc->benign_zerodiag_subs[s], &idxs)); in PCBDDCBenignDetectSaddlePoint()
3196 cols[0] = idxs[i]; in PCBDDCBenignDetectSaddlePoint()
3197 cols[1] = idxs[nzs - 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()
3208 PetscCall(ISRestoreIndices(pcbddc->benign_zerodiag_subs[s], &idxs)); in PCBDDCBenignDetectSaddlePoint()
3304 PetscInt *idxs, j, nz, cum; in PCBDDCBenignPopOrPushB0() local
3308 PetscCall(ISGetIndices(pcbddc->benign_zerodiag_subs[i], (const PetscInt **)&idxs)); in PCBDDCBenignPopOrPushB0()
3310 PetscCall(VecSetValues(matis->x, nz, idxs, vals, INSERT_VALUES)); in PCBDDCBenignPopOrPushB0()
3326 PetscCall(ISRestoreIndices(pcbddc->benign_zerodiag_subs[i], (const PetscInt **)&idxs)); in PCBDDCBenignPopOrPushB0()
3466 const PetscInt *idxs; in PCBDDCAdaptiveSelection() local
3468 PetscCall(ISGetIndices(sub_schurs->is_vertices, &idxs)); in PCBDDCAdaptiveSelection()
3471 pcbddc->adaptive_constraints_idxs[cum] = idxs[cum]; in PCBDDCAdaptiveSelection()
3476 PetscCall(ISRestoreIndices(sub_schurs->is_vertices, &idxs)); in PCBDDCAdaptiveSelection()
3493 const PetscInt *idxs; in PCBDDCAdaptiveSelection() local
3514 PetscCall(ISGetIndices(sub_schurs->is_subs[i], &idxs)); in PCBDDCAdaptiveSelection()
3531 const PetscInt *idxs; in PCBDDCAdaptiveSelection() local
3534 PetscCall(ISGetIndices(sub_schurs->change_primal_sub[i], &idxs)); in PCBDDCAdaptiveSelection()
3536 *(Sarray + cumarray + idxs[k] * (subset_size + 1)) *= -1.0; in PCBDDCAdaptiveSelection()
3537 *(Starray + cumarray + idxs[k] * (subset_size + 1)) = 0.0; in PCBDDCAdaptiveSelection()
3539 PetscCall(ISRestoreIndices(sub_schurs->change_primal_sub[i], &idxs)); in PCBDDCAdaptiveSelection()
3585 …PetscInt nc = 0, c = pcbddc->mat_graph->nodes[idxs[0]].count, w = pcbddc->mat_graph->nodes[idxs[0]… in PCBDDCAdaptiveSelection()
3916 …cbddc->adaptive_constraints_idxs + pcbddc->adaptive_constraints_idxs_ptr[cum], idxs, subset_size)); in PCBDDCAdaptiveSelection()
3921 PetscCall(ISRestoreIndices(sub_schurs->is_subs[i], &idxs)); in PCBDDCAdaptiveSelection()
4698 const PetscInt *idxs; in PCBDDCSetUpCorrection() local
4700 PetscCall(ISGetIndices(pcis->is_I_local, &idxs)); in PCBDDCSetUpCorrection()
4702 … i++) PetscCall(PetscFindInt(pcbddc->benign_p0_lidx[i], pcis->n - pcis->n_B, idxs, &p0_lidx_I[i])); in PCBDDCSetUpCorrection()
4703 PetscCall(ISRestoreIndices(pcis->is_I_local, &idxs)); in PCBDDCSetUpCorrection()
4793 const PetscInt *idxs, *idxs_zero; in PCBDDCSetUpCorrection() local
4798 PetscCall(MatGetRow(A_RV_bcorr, i, &n, &idxs, &vals)); in PCBDDCSetUpCorrection()
4801 PetscInt k, col = idxs[j]; in PCBDDCSetUpCorrection()
4804 PetscCall(MatRestoreRow(A_RV_bcorr, i, &n, &idxs, &vals)); in PCBDDCSetUpCorrection()
4974 const PetscInt *idxs, *idxs_zero; in PCBDDCSetUpCorrection() local
4983 PetscCall(MatGetRow(A_RV_bcorr, i, &n, &idxs, &vals)); in PCBDDCSetUpCorrection()
4986 …t k = 0; k < n_vertices; k++) PetscCall(MatSetValue(*coarse_submat, idx_V[idxs[j]], idx_V[k], val … in PCBDDCSetUpCorrection()
4988 PetscCall(MatRestoreRow(A_RV_bcorr, i, &n, &idxs, &vals)); in PCBDDCSetUpCorrection()
5473 const PetscInt *idxs; in MatCreateSubMatrixUnsorted() local
5479 PetscCall(ISGetIndices(isrow, &idxs)); in MatCreateSubMatrixUnsorted()
5480 PetscCall(PetscSortIntWithPermutation(rsize, idxs, idxs_perm_r)); in MatCreateSubMatrixUnsorted()
5481 for (i = 0; i < rsize; i++) idxs_sorted[i] = idxs[idxs_perm_r[i]]; in MatCreateSubMatrixUnsorted()
5482 PetscCall(ISRestoreIndices(isrow, &idxs)); in MatCreateSubMatrixUnsorted()
5494 const PetscInt *idxs; in MatCreateSubMatrixUnsorted() local
5500 PetscCall(ISGetIndices(iscol, &idxs)); in MatCreateSubMatrixUnsorted()
5501 PetscCall(PetscSortIntWithPermutation(csize, idxs, idxs_perm_c)); in MatCreateSubMatrixUnsorted()
5502 for (i = 0; i < csize; i++) idxs_sorted[i] = idxs[idxs_perm_c[i]]; in MatCreateSubMatrixUnsorted()
5503 PetscCall(ISRestoreIndices(iscol, &idxs)); in MatCreateSubMatrixUnsorted()
6037 const PetscInt *idxs; in PCBDDCSetUpLocalSolvers() local
6041 PetscCall(ISGetIndices(pcis->is_I_local, &idxs)); in PCBDDCSetUpLocalSolvers()
6044 for (d = 0; d < cdim; d++) scoords[i * cdim + d] = coords[idxs[i] * cdim + d]; in PCBDDCSetUpLocalSolvers()
6046 PetscCall(ISRestoreIndices(pcis->is_I_local, &idxs)); in PCBDDCSetUpLocalSolvers()
6185 const PetscInt *idxs; in PCBDDCSetUpLocalSolvers() local
6189 PetscCall(ISGetIndices(pcbddc->is_R_local, &idxs)); in PCBDDCSetUpLocalSolvers()
6192 for (d = 0; d < cdim; d++) scoords[i * cdim + d] = coords[idxs[i] * cdim + d]; in PCBDDCSetUpLocalSolvers()
6194 PetscCall(ISRestoreIndices(pcbddc->is_R_local, &idxs)); in PCBDDCSetUpLocalSolvers()
7639 const PetscInt *idxs; in PCBDDCAnalyzeInterface() local
7643 PetscCall(ISGetIndices(pcbddc->local_subs[i], &idxs)); in PCBDDCAnalyzeInterface()
7644 for (j = 0; j < nl; j++) local_subs[idxs[j]] = i; in PCBDDCAnalyzeInterface()
7645 PetscCall(ISRestoreIndices(pcbddc->local_subs[i], &idxs)); in PCBDDCAnalyzeInterface()
8653 const PetscInt *idxs; in PCBDDCSetUpCoarseSolver() local
8681 PetscCall(ISGetIndices(pcbddc->ISForDofsLocal[i], &idxs)); in PCBDDCSetUpCoarseSolver()
8682 PetscCall(ISGlobalToLocalMappingApply(tmap, IS_GTOLM_DROP, tsize, idxs, &nout, tidxs)); in PCBDDCSetUpCoarseSolver()
8683 PetscCall(ISRestoreIndices(pcbddc->ISForDofsLocal[i], &idxs)); in PCBDDCSetUpCoarseSolver()
8686 PetscCall(ISGetIndices(pcbddc->nedclocal, &idxs)); in PCBDDCSetUpCoarseSolver()
8687 PetscCall(ISGlobalToLocalMappingApply(tmap, IS_GTOLM_DROP, tsize, idxs, &nout, tidxs)); in PCBDDCSetUpCoarseSolver()
8689 PetscCall(ISRestoreIndices(pcbddc->nedclocal, &idxs)); in PCBDDCSetUpCoarseSolver()
8699 PetscCall(ISGetIndices(pcbddc->NeumannBoundariesLocal, &idxs)); in PCBDDCSetUpCoarseSolver()
8700 PetscCall(ISGlobalToLocalMappingApply(tmap, IS_GTOLM_DROP, tsize, idxs, &nout, tidxs)); in PCBDDCSetUpCoarseSolver()
8701 PetscCall(ISRestoreIndices(pcbddc->NeumannBoundariesLocal, &idxs)); in PCBDDCSetUpCoarseSolver()
8710 PetscCall(ISGetIndices(corners, &idxs)); in PCBDDCSetUpCoarseSolver()
8711 PetscCall(ISGlobalToLocalMappingApply(tmap, IS_GTOLM_DROP, tsize, idxs, &nout, tidxs)); in PCBDDCSetUpCoarseSolver()
8713 PetscCall(ISRestoreIndices(corners, &idxs)); in PCBDDCSetUpCoarseSolver()
9057 const PetscInt *idxs; in PCBDDCSetUpCoarseSolver() local
9073 PetscCall(ISGetIndices(is, &idxs)); in PCBDDCSetUpCoarseSolver()
9076 … for (d = 0; d < cdim; d++) coords[cdim * i + d] = pcbddc->mat_graph->coords[cdim * idxs[i] + d]; in PCBDDCSetUpCoarseSolver()
9078 PetscCall(ISRestoreIndices(is, &idxs)); in PCBDDCSetUpCoarseSolver()
9082 PetscCall(ISGetIndices(corners, &idxs)); in PCBDDCSetUpCoarseSolver()
9083 PetscCall(VecSetValuesBlocked(gv, n, idxs, coords, INSERT_VALUES)); in PCBDDCSetUpCoarseSolver()
9084 PetscCall(ISRestoreIndices(corners, &idxs)); in PCBDDCSetUpCoarseSolver()
9326 PetscInt i, lsize, *idxs, n; in PCBDDCGlobalToLocal() local
9334 PetscCall(ISGetIndices(globalis, (const PetscInt **)&idxs)); in PCBDDCGlobalToLocal()
9337 if (idxs) { /* multilevel guard */ in PCBDDCGlobalToLocal()
9339 PetscCall(VecSetValues(gwork, lsize, idxs, vals, INSERT_VALUES)); in PCBDDCGlobalToLocal()
9342 PetscCall(ISRestoreIndices(globalis, (const PetscInt **)&idxs)); in PCBDDCGlobalToLocal()
9353 PetscCall(PetscMalloc1(lsize, &idxs)); in PCBDDCGlobalToLocal()
9355 if (PetscRealPart(vals[i]) > 0.5) idxs[lsize++] = i; in PCBDDCGlobalToLocal()
9358 …PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)gwork), lsize, idxs, PETSC_OWN_POINTER, &lo… in PCBDDCGlobalToLocal()
9542 PetscInt vsize, *idxs; in PCBDDCInitSubSchurs() local
9546 PetscCall(ISGetIndices(verticesIS, (const PetscInt **)&idxs)); in PCBDDCInitSubSchurs()
9547 …PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)pc), vsize, idxs, PETSC_COPY_VALUES, &verti… in PCBDDCInitSubSchurs()
9548 PetscCall(ISRestoreIndices(verticesIS, (const PetscInt **)&idxs)); in PCBDDCInitSubSchurs()
9588 const PetscInt *idxs; in PCBDDCViewGlobalIS() local
9598 PetscCall(ISGetIndices(is, &idxs)); in PCBDDCViewGlobalIS()
9600 if (idxs[i] < 0 || idxs[i] >= ln) continue; in PCBDDCViewGlobalIS()
9601 matis->sf_leafdata[idxs[i]] = 1; in PCBDDCViewGlobalIS()
9603 PetscCall(ISRestoreIndices(is, &idxs)); in PCBDDCViewGlobalIS()
9814 const PetscInt *idxs; in MatAIJExtractRows() local
9821 PetscCall(ISGetIndices(rows, &idxs)); in MatAIJExtractRows()
9822 PetscCall(PetscSFSetGraphLayout(sf, A->rmap, ni, NULL, PETSC_USE_POINTER, idxs)); in MatAIJExtractRows()
9823 PetscCall(ISRestoreIndices(rows, &idxs)); in MatAIJExtractRows()