| /petsc/src/mat/impls/sell/seq/ |
| H A D | sell.h | 125 #define MatSetValue_SeqSELL_Private(A, row, col, value, addv, orow, ocol, cp, vp, lastcol, low, hig… argument
140 if (addv == ADD_VALUES) *(vp + a->sliceheight * _i) += value; \
141 else *(vp + a->sliceheight * _i) = value; \
164 vp = new_val + a->sliidx[row / a->sliceheight] + (row % a->sliceheight); \
185 *(vp + a->sliceheight * (ii + 1)) = *(vp + a->sliceheight * ii); \
188 *(vp + a->sliceheight * _i) = value; \
|
| H A D | sell.c | 1048 MatScalar *vp; in MatGetValues_SeqSELL() local
1057 vp = a->val + shift; /* pointer to the row */ in MatGetValues_SeqSELL()
1072 *v++ = *(vp + a->sliceheight * i); in MatGetValues_SeqSELL()
1440 MatScalar *vp; in MatAssemblyEnd_SeqSELL() local
1456 …vp = PetscSafePointerPlusOffset(a->val, shift); /* pointer to the nonzero … in MatAssemblyEnd_SeqSELL()
1481 vp[a->sliceheight * k + row_in_slice] = (MatScalar)0; in MatAssemblyEnd_SeqSELL()
1538 MatScalar *vp, value; in MatSetValues_SeqSELL() local
1551 vp = a->val + shift; /* pointer to the row */ in MatSetValues_SeqSELL()
1579 if (is == ADD_VALUES) *(vp + a->sliceheight * i) += value; in MatSetValues_SeqSELL()
1580 else *(vp + a->sliceheight * i) = value; in MatSetValues_SeqSELL()
[all …]
|
| /petsc/src/ts/tests/ |
| H A D | ex27.c | 287 static PetscErrorCode ComputeGradS(PetscInt dim, PetscInt Np, const PetscReal vp[], const PetscReal… in ComputeGradS() argument
306 … < dim; ++d) integral[d] += (-1. / epsilon) * PetscAbsReal(vp[d] - vc_l[d]) * Gaussian(dim, vp, ep… in ComputeGradS()
314 static PetscErrorCode QCompute(PetscInt dim, const PetscReal vp[], const PetscReal vq[], PetscReal … in QCompute() argument
320 DMPlex_WaxpyD_Internal(dim, -1.0, vq, vp, xi); in QCompute()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexglvis.c | 1011 PetscInt vp, gvp; in DMPlexView_GLVis_ASCII() local
1013 for (vp = 0, p = vStart; p < vEnd; p++) { in DMPlexView_GLVis_ASCII()
1021 if (parent != p) vp++; in DMPlexView_GLVis_ASCII()
1023 PetscCallMPI(MPIU_Allreduce(&vp, &gvp, 1, MPIU_INT, MPI_SUM, PetscObjectComm((PetscObject)dm))); in DMPlexView_GLVis_ASCII()
1029 PetscCall(PetscViewerASCIIPrintf(viewer, "%" PetscInt_FMT "\n", vp)); in DMPlexView_GLVis_ASCII()
1052 vp--; in DMPlexView_GLVis_ASCII()
1104 vp--; in DMPlexView_GLVis_ASCII()
1115 …PetscCheck(!vp, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Unexpected %" PetscInt_FMT " hanging vertices", … in DMPlexView_GLVis_ASCII()
|
| /petsc/src/ts/interface/sensitivity/ |
| H A D | tssen.c | 1733 PetscErrorCode TSForwardSetIntegralGradients(TS ts, PetscInt numfwdint, Vec vp[]) in TSForwardSetIntegralGradients() argument
1740 ts->vecs_integral_sensip = vp; in TSForwardSetIntegralGradients()
1758 PetscErrorCode TSForwardGetIntegralGradients(TS ts, PetscInt *numfwdint, Vec *vp[]) in TSForwardGetIntegralGradients() argument
1762 PetscAssertPointer(vp, 3); in TSForwardGetIntegralGradients()
1764 if (vp) *vp = ts->vecs_integral_sensip; in TSForwardGetIntegralGradients()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex11.c | 1627 PetscReal vp[3]; in initLinearWave() local
1646 projecttoprim(vp, wcp, rv); in initLinearWave()
1647 ux->r = vp[0]; /* density */ in initLinearWave()
1648 ux->ru[0] = vp[0] * vp[1]; /* x momentum */ in initLinearWave()
1651 ux->E = vp[2] / (gamma - 1.) + 0.5 * vp[0] * vp[1] * vp[1]; in initLinearWave()
|
| /petsc/src/ksp/ksp/impls/cg/pipecg2/ |
| H A D | pipecg2.c | 82 static PetscErrorCode VecMergedOpsShort_Private(Vec vx, Vec vr, Vec vz, Vec vw, Vec vp, Vec vq, Vec… in VecMergedOpsShort_Private() argument
94 PetscCall(VecGetArray(vp, (PetscScalar **)&pp)); in VecMergedOpsShort_Private()
284 PetscCall(VecRestoreArray(vp, (PetscScalar **)&pp)); in VecMergedOpsShort_Private()
304 static PetscErrorCode VecMergedOps_Private(Vec vx, Vec vr, Vec vz, Vec vw, Vec vp, Vec vq, Vec vc, … in VecMergedOps_Private() argument
316 PetscCall(VecGetArray(vp, (PetscScalar **)&pp)); in VecMergedOps_Private()
513 PetscCall(VecRestoreArray(vp, (PetscScalar **)&pp)); in VecMergedOps_Private()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex70.c | 701 PetscErrorCode MaterialPoint_AdvectRK1(DM dm_vp, Vec vp, PetscReal dt, DM dm_mpoint) in MaterialPoint_AdvectRK1() argument
721 PetscCall(DMGlobalToLocalBegin(dm_vp, vp, INSERT_VALUES, vp_l)); in MaterialPoint_AdvectRK1()
722 PetscCall(DMGlobalToLocalEnd(dm_vp, vp, INSERT_VALUES, vp_l)); in MaterialPoint_AdvectRK1()
|
| /petsc/src/ksp/ksp/utils/dm/ |
| H A D | dmproject.c | 1088 const PetscReal vp = v[p * dim + d]; in DMSwarmRemap_Colella_Internal() local
1091 vs[d] = PetscFloorReal((vp - vmin[d]) / hv[d]); in DMSwarmRemap_Colella_Internal()
|
| /petsc/src/ksp/pc/impls/bddc/ |
| H A D | bddcprivate.c | 8733 Vec vp[1]; in PCBDDCSetUpCoarseSolver() local
8737 vp[0] = NULL; in PCBDDCSetUpCoarseSolver()
8740 PetscCall(VecCreate(PetscObjectComm((PetscObject)pc), &vp[0])); in PCBDDCSetUpCoarseSolver()
8741 PetscCall(VecSetSizes(vp[0], pcbddc->local_primal_size, PETSC_DECIDE)); in PCBDDCSetUpCoarseSolver()
8742 PetscCall(VecSetType(vp[0], VECSTANDARD)); in PCBDDCSetUpCoarseSolver()
8760 PetscCall(VecGetArray(vp[0], &array)); in PCBDDCSetUpCoarseSolver()
8764 PetscCall(VecRestoreArray(vp[0], &array)); in PCBDDCSetUpCoarseSolver()
8777 …dc->coarse_subassembling, 0, restr, full_restr, PETSC_TRUE, &coarse_mat, nis, isarray, nvecs, vp)); in PCBDDCSetUpCoarseSolver()
8779 …coarse_subassembling, 0, restr, full_restr, PETSC_FALSE, &coarse_mat_is, nis, isarray, nvecs, vp)); in PCBDDCSetUpCoarseSolver()
8781 if (vp[0]) { /* vp[0] could have been placed on a different set of processes */ in PCBDDCSetUpCoarseSolver()
[all …]
|