Lines Matching refs:cum
168 PetscInt i, j, extmem, cum, maxsize, nee; in PCBDDCNedelecSupport() local
258 for (i = 0, cum = 0; i < ne; i++) { in PCBDDCNedelecSupport()
259 if (PetscRealPart(vals[idxs[i]]) > 2.) eidxs[cum++] = idxs[i]; in PCBDDCNedelecSupport()
263 for (i = 0, cum = 0; i < ne; i++) { in PCBDDCNedelecSupport()
264 if (PetscRealPart(vals[i]) > 2.) eidxs[cum++] = i; in PCBDDCNedelecSupport()
268 PetscCall(ISCreateGeneral(comm, cum, eidxs, PETSC_COPY_VALUES, &enedfieldlocal)); in PCBDDCNedelecSupport()
330 for (i = 0, cum = 0; i < pc->pmat->rmap->n; i++) { in PCBDDCNedelecSupport()
331 if (matis->sf_rootdata[i] < 2) matis->sf_rootdata[cum++] = i + rst; in PCBDDCNedelecSupport()
334 PetscCall(MatZeroRows(G, cum, matis->sf_rootdata, 0., NULL, NULL)); in PCBDDCNedelecSupport()
342 for (i = 0, cum = 0; i < ne; i++) in PCBDDCNedelecSupport()
343 if (matis->sf_leafdata[idxs[i]] == 1) tbz[cum++] = i; in PCBDDCNedelecSupport()
345 PetscCall(ISLocalToGlobalMappingApply(el2g, cum, tbz, tbz)); in PCBDDCNedelecSupport()
346 PetscCall(MatZeroRows(G, cum, tbz, 0., NULL, NULL)); in PCBDDCNedelecSupport()
477 PetscCall(ISGetLocalSize(is, &cum)); in PCBDDCNedelecSupport()
479 for (i = 0; i < cum; i++) { in PCBDDCNedelecSupport()
496 PetscCall(ISGetLocalSize(is, &cum)); in PCBDDCNedelecSupport()
498 for (i = 0; i < cum; i++) { in PCBDDCNedelecSupport()
522 cum = 0; in PCBDDCNedelecSupport()
526 marks[cum++] = i; in PCBDDCNedelecSupport()
532 marks[cum++] = i; in PCBDDCNedelecSupport()
553 marks[cum++] = i; in PCBDDCNedelecSupport()
572 PetscCall(MatZeroRows(lGe, cum, marks, 0., NULL, NULL)); in PCBDDCNedelecSupport()
578 PetscCall(ISGetLocalSize(elements_corners, &cum)); in PCBDDCNedelecSupport()
580 for (i = 0; i < cum; i++) PetscCall(PetscBTSet(btvcand, idxs[i])); in PCBDDCNedelecSupport()
774 cum = 0; in PCBDDCNedelecSupport()
776 if (sfvleaves[i]) vmarks[cum++] = i; in PCBDDCNedelecSupport()
778 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, cum, vmarks, PETSC_COPY_VALUES, &tbz)); in PCBDDCNedelecSupport()
791 cum = 0; in PCBDDCNedelecSupport()
794 vmarks[cum++] = i; in PCBDDCNedelecSupport()
796 } else if (!PetscBTLookup(btvi, i)) vmarks[cum++] = i; in PCBDDCNedelecSupport()
802 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, cum, vmarks, PETSC_COPY_VALUES, &tbz)); in PCBDDCNedelecSupport()
807 PetscCall(MatZeroRows(lGt, cum, vmarks, 0., NULL, NULL)); in PCBDDCNedelecSupport()
823 cum = 0; in PCBDDCNedelecSupport()
825 if (PetscUnlikely(PetscBTLookup(bte, i))) marks[cum++] = i; in PCBDDCNedelecSupport()
827 if (fl2g) PetscCall(ISLocalToGlobalMappingApply(fl2g, cum, marks, marks)); in PCBDDCNedelecSupport()
828 PetscCall(ISCreateGeneral(comm, cum, marks, PETSC_COPY_VALUES, &primals)); in PCBDDCNedelecSupport()
934 PetscCall(ISGetLocalSize(primals, &cum)); in PCBDDCNedelecSupport()
936 for (i = 0; i < cum; i++) marks[idxs[i]] = nee + 1; in PCBDDCNedelecSupport()
1003 PetscCall(ISGetLocalSize(primals, &cum)); in PCBDDCNedelecSupport()
1005 for (i = 0; i < cum; i++) marks[idxs[i]] = 0; in PCBDDCNedelecSupport()
1026 cum = 0; in PCBDDCNedelecSupport()
1036 if (!PetscBTLookup(btv, vv)) extrow[cum++] = vv; in PCBDDCNedelecSupport()
1041 PetscCall(PetscSortRemoveDupsInt(&cum, extrow)); in PCBDDCNedelecSupport()
1042 PetscCall(ISLocalToGlobalMappingApply(vl2g, cum, extrow, gidxs)); in PCBDDCNedelecSupport()
1043 PetscCall(PetscSortIntWithArray(cum, gidxs, extrow)); in PCBDDCNedelecSupport()
1044 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, cum, extrow, PETSC_COPY_VALUES, &extcols[i])); in PCBDDCNedelecSupport()
1047 if (cum != size - 1 || found != 2) { in PCBDDCNedelecSupport()
1064 PetscCall(ISGetLocalSize(primals, &cum)); in PCBDDCNedelecSupport()
1066 PetscCall(PetscArraycpy(newprimals, idxs, cum)); in PCBDDCNedelecSupport()
1093 newprimals[cum++] = ee2; in PCBDDCNedelecSupport()
1120 newprimals[cum++] = idxs[j]; in PCBDDCNedelecSupport()
1131 PetscCall(PetscSortRemoveDupsInt(&cum, newprimals)); in PCBDDCNedelecSupport()
1133 PetscCall(ISLocalToGlobalMappingApply(fl2g, cum, newprimals, newprimals)); in PCBDDCNedelecSupport()
1139 PetscCall(ISCreateGeneral(comm, cum, newprimals, PETSC_COPY_VALUES, &primals)); in PCBDDCNedelecSupport()
1176 cum = 0; in PCBDDCNedelecSupport()
1184 if (!PetscBTLookup(btv, jj[k])) extrow[cum++] = jj[k]; in PCBDDCNedelecSupport()
1187 PetscCall(PetscSortRemoveDupsInt(&cum, extrow)); in PCBDDCNedelecSupport()
1188 PetscCall(ISLocalToGlobalMappingApply(vl2g, cum, extrow, gidxs)); in PCBDDCNedelecSupport()
1189 PetscCall(PetscSortIntWithArray(cum, gidxs, extrow)); in PCBDDCNedelecSupport()
1190 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, cum, extrow, PETSC_COPY_VALUES, &extcols[i])); in PCBDDCNedelecSupport()
1191 if (cum != size - 1) { in PCBDDCNedelecSupport()
1315 cum = 0; in PCBDDCNedelecSupport()
1320 cum = PetscMax(cum, size); in PCBDDCNedelecSupport()
1327 PetscCall(PetscMalloc2((5 + cum + maxsize) * maxsize, &work, maxsize, &rwork)); in PCBDDCNedelecSupport()
1422 cum = 0; in PCBDDCNedelecSupport()
1428 cum += j - 1; in PCBDDCNedelecSupport()
1430 PetscCall(PetscMalloc1(PetscMax(cum, pc->pmat->rmap->n), &eedgesidxs)); in PCBDDCNedelecSupport()
1431 cum = 0; in PCBDDCNedelecSupport()
1438 PetscCall(PetscArraycpy(eedgesidxs + cum, idxs, j - 1)); /* last on the edge is primal */ in PCBDDCNedelecSupport()
1440 cum += j - 1; in PCBDDCNedelecSupport()
1442 PetscCall(ISLocalToGlobalMappingApply(al2g, cum, eedgesidxs, eedgesidxs)); in PCBDDCNedelecSupport()
1443 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, cum, eedgesidxs, PETSC_USE_POINTER, &galleedges)); in PCBDDCNedelecSupport()
1452 PetscCall(VecSetValues(E, cum, eedgesidxs, NULL, INSERT_VALUES)); in PCBDDCNedelecSupport()
1456 for (i = 0, cum = 0; i < pc->pmat->rmap->n; i++) in PCBDDCNedelecSupport()
1457 if (evals[i] == 0.0) eedgesidxs[cum++] = i + pc->pmat->rmap->rstart; in PCBDDCNedelecSupport()
1459 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, cum, eedgesidxs, PETSC_COPY_VALUES, &is_E_to_zero)); in PCBDDCNedelecSupport()
1548 cum = 0; in PCBDDCNedelecSupport()
1553 cum = 1; in PCBDDCNedelecSupport()
1556 …PetscCall(VecScatterBegin(nnsp_vscat, nnsp_vecs[i], quadvecs[i + cum], INSERT_VALUES, SCATTER_FORW… in PCBDDCNedelecSupport()
1557 …PetscCall(VecScatterEnd(nnsp_vscat, nnsp_vecs[i], quadvecs[i + cum], INSERT_VALUES, SCATTER_FORWAR… in PCBDDCNedelecSupport()
1561 cum = 0; in PCBDDCNedelecSupport()
1565 cum = 1; in PCBDDCNedelecSupport()
1568 PetscCall(VecCopy(onnsp_vecs[i], quadvecs[i + cum + new_nnsp_nvecs])); in PCBDDCNedelecSupport()
1569 PetscCall(VecISSet(quadvecs[i + cum + new_nnsp_nvecs], is_E_to_zero, 0)); in PCBDDCNedelecSupport()
1610 cum = 0; in PCBDDCNedelecSupport()
1611 while (cum < nee && cedges[cum] < 0) cum++; in PCBDDCNedelecSupport()
1618 if (cum != nee && oqueue[ocptr[i + 1] - 1] == cedges[cum]) { /* this cc has a primal dof */ in PCBDDCNedelecSupport()
1620 graph->queue[graph->cptr[ncc]] = cedges[cum]; in PCBDDCNedelecSupport()
1623 cum++; in PCBDDCNedelecSupport()
1624 while (cum < nee && cedges[cum] < 0) cum++; in PCBDDCNedelecSupport()
2522 PetscInt j, cum; in PCBDDCDetectDisconnectedComponents() local
2526 cum = 0; in PCBDDCDetectDisconnectedComponents()
2532 adjncy_filtered[cum + xadj_filtered[i]++] = adjncy[j]; in PCBDDCDetectDisconnectedComponents()
2535 xadj_filtered[i] = cum; in PCBDDCDetectDisconnectedComponents()
2536 cum += t; in PCBDDCDetectDisconnectedComponents()
2575 PetscInt *ids, cum, cump, *cids, *pids; in PCBDDCDetectDisconnectedComponents() local
2645 for (i = 0, cump = 0, cum = 0; i < graph->ncc; i++) { in PCBDDCDetectDisconnectedComponents()
2662 else if (!PetscBTLookup(btv, off + s)) ids[cum++] = off + s; in PCBDDCDetectDisconnectedComponents()
2673 else if (!PetscBTLookup(btv, off + s)) ids[cum++] = off + s; in PCBDDCDetectDisconnectedComponents()
2680 cids[i + 1] = cum; in PCBDDCDetectDisconnectedComponents()
3304 PetscInt *idxs, j, nz, cum; in PCBDDCBenignPopOrPushB0() local
3316 cum = 0; in PCBDDCBenignPopOrPushB0()
3319 vals[cum] = array[j]; in PCBDDCBenignPopOrPushB0()
3320 idxs_ins[cum] = j; in PCBDDCBenignPopOrPushB0()
3321 cum++; in PCBDDCBenignPopOrPushB0()
3324 PetscCall(MatSetValues(pcbddc->benign_B0, 1, &i, cum, idxs_ins, vals, INSERT_VALUES)); in PCBDDCBenignPopOrPushB0()
3361 PetscInt i, nmax, nmin, nv, cum, mss, cum2, cumarray, maxneigs; in PCBDDCAdaptiveSelection() local
3404 cum = cum2 = 0; in PCBDDCAdaptiveSelection()
3411 cum += subset_size; in PCBDDCAdaptiveSelection()
3457 … nv + sub_schurs->n_subs + 1, &pcbddc->adaptive_constraints_data_ptr, nv + cum, &pcbddc->adaptive_… in PCBDDCAdaptiveSelection()
3462 cum = cumarray = 0; in PCBDDCAdaptiveSelection()
3469 for (cum = 0; cum < nv; cum++) { in PCBDDCAdaptiveSelection()
3470 pcbddc->adaptive_constraints_n[cum] = 1; in PCBDDCAdaptiveSelection()
3471 pcbddc->adaptive_constraints_idxs[cum] = idxs[cum]; in PCBDDCAdaptiveSelection()
3472 pcbddc->adaptive_constraints_data[cum] = 1.0; in PCBDDCAdaptiveSelection()
3473 … pcbddc->adaptive_constraints_idxs_ptr[cum + 1] = pcbddc->adaptive_constraints_idxs_ptr[cum] + 1; in PCBDDCAdaptiveSelection()
3474 … pcbddc->adaptive_constraints_data_ptr[cum + 1] = pcbddc->adaptive_constraints_data_ptr[cum] + 1; in PCBDDCAdaptiveSelection()
3899 …(pcbddc->adaptive_constraints_data + pcbddc->adaptive_constraints_data_ptr[cum], eigv + eigs_start… in PCBDDCAdaptiveSelection()
3907 …ddc->adaptive_constraints_data[ii * subset_size + j + pcbddc->adaptive_constraints_data_ptr[cum]]); in PCBDDCAdaptiveSelection()
3908 …ddc->adaptive_constraints_data[ii * subset_size + j + pcbddc->adaptive_constraints_data_ptr[cum]]); in PCBDDCAdaptiveSelection()
3911 …c->adaptive_constraints_data[ii * subset_size + j + pcbddc->adaptive_constraints_data_ptr[cum]]))); in PCBDDCAdaptiveSelection()
3916 …(pcbddc->adaptive_constraints_idxs + pcbddc->adaptive_constraints_idxs_ptr[cum], idxs, subset_size… in PCBDDCAdaptiveSelection()
3917 …pcbddc->adaptive_constraints_idxs_ptr[cum + 1] = pcbddc->adaptive_constraints_idxs_ptr[cum] + subs… in PCBDDCAdaptiveSelection()
3918 …pcbddc->adaptive_constraints_data_ptr[cum + 1] = pcbddc->adaptive_constraints_data_ptr[cum] + subs… in PCBDDCAdaptiveSelection()
3919 cum++; in PCBDDCAdaptiveSelection()
6537 PetscInt i, j, k, total_counts, total_counts_cc, cum; in PCBDDCConstraintsSetUp() local
7004 cum = total_primal_vertices; in PCBDDCConstraintsSetUp()
7007 pcbddc->local_primal_ref_node[cum] = constraints_idxs[constraints_idxs_ptr[i]]; in PCBDDCConstraintsSetUp()
7008 pcbddc->local_primal_ref_mult[cum] = constraints_n[i]; in PCBDDCConstraintsSetUp()
7009 cum++; in PCBDDCConstraintsSetUp()