| /petsc/src/ksp/ksp/impls/gmres/agmres/ |
| H A D | agmresleja.c | 25 PetscScalar rd, id, pd, max; in KSPAGMRESLejaCfpdMax() local
29 pd = 1.0; in KSPAGMRESLejaCfpdMax()
36 pd = pd * PetscSqrtReal(rd * rd + id * id); in KSPAGMRESLejaCfpdMax()
38 if (max < pd) { in KSPAGMRESLejaCfpdMax()
40 max = pd; in KSPAGMRESLejaCfpdMax()
42 pd = 1.0; in KSPAGMRESLejaCfpdMax()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexfvm.c | 91 PetscInt ghost, c, pd, d, numChildren, numCells; in DMPlexReconstructGradients_Internal() local
109 for (pd = 0; pd < dof; ++pd) { in DMPlexReconstructGradients_Internal()
110 PetscScalar delta = cx[1][pd] - cx[0][pd]; in DMPlexReconstructGradients_Internal()
113 if (cgrad[0]) cgrad[0][pd * dim + d] += fg->grad[0][d] * delta; in DMPlexReconstructGradients_Internal()
114 if (cgrad[1]) cgrad[1][pd * dim + d] -= fg->grad[1][d] * delta; in DMPlexReconstructGradients_Internal()
127 PetscInt coneSize, f, pd, d; in DMPlexReconstructGradients_Internal() local
147 for (pd = 0; pd < dof; ++pd) /* Scalar limiter applied to each component separately */ in DMPlexReconstructGradients_Internal()
148 for (d = 0; d < dim; ++d) cgrad[pd * dim + d] *= cellPhi[pd]; in DMPlexReconstructGradients_Internal()
|
| /petsc/src/ksp/ksp/impls/cg/pipecg2/ |
| H A D | pipecg2.c | 86 …PetscScalar *PETSC_RESTRICT pc, *PETSC_RESTRICT pd, *PETSC_RESTRICT pg0, *PETSC_RESTRICT ph0, *PET… in VecMergedOpsShort_Private() local
97 PetscCall(VecGetArray(vd, (PetscScalar **)&pd)); in VecMergedOpsShort_Private()
122 pd[j] = ph0[j]; in VecMergedOpsShort_Private()
131 pn[j] = pn[j] - alpha0 * pd[j]; in VecMergedOpsShort_Private()
143 pd[j] = ph0[j] + beta1 * pd[j]; in VecMergedOpsShort_Private()
150 pn[j] = pn[j] - alpha1 * pd[j]; in VecMergedOpsShort_Private()
177 pd[j] = ph0[j]; in VecMergedOpsShort_Private()
186 pn[j] = pn[j] - alpha0 * pd[j]; in VecMergedOpsShort_Private()
198 pd[j] = ph0[j] + beta1 * pd[j]; in VecMergedOpsShort_Private()
205 pn[j] = pn[j] - alpha1 * pd[j]; in VecMergedOpsShort_Private()
[all …]
|
| /petsc/src/binding/petsc4py/ |
| H A D | setup.py | 176 pd = os.environ[pkg.upper() + '_DIR']
178 depends += glob_join(pd, 'include', '*.h')
179 depends += glob_join(pd, 'include', pkg, 'private', '*.h')
180 depends += glob_join(pd, pa, 'include', '%sconf.h' % pkg)
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpiptap.c | 576 …data, *ao = (Mat_SeqAIJ *)a->B->data, *p_oth = (Mat_SeqAIJ *)P_oth->data, *pd = (Mat_SeqAIJ *)p->A… in MatPtAPSymbolicComputeOneRowOfAP_private() local
593 nzpi = pd->i[row + 1] - pd->i[row]; in MatPtAPSymbolicComputeOneRowOfAP_private()
594 pj = pd->j + pd->i[row]; in MatPtAPSymbolicComputeOneRowOfAP_private()
635 …data, *ao = (Mat_SeqAIJ *)a->B->data, *p_oth = (Mat_SeqAIJ *)P_oth->data, *pd = (Mat_SeqAIJ *)p->A… in MatPtAPNumericComputeOneRowOfAP_private() local
653 nzpi = pd->i[row + 1] - pd->i[row]; in MatPtAPNumericComputeOneRowOfAP_private()
654 pj = pd->j + pd->i[row]; in MatPtAPNumericComputeOneRowOfAP_private()
655 pa = pd->a + pd->i[row]; in MatPtAPNumericComputeOneRowOfAP_private()
698 Mat_SeqAIJ *cd, *co, *po = (Mat_SeqAIJ *)p->B->data, *pd = (Mat_SeqAIJ *)p->A->data; in MatPtAPNumeric_MPIAIJ_MPIXAIJ_allatonce() local
801 nzi = pd->i[ii + 1] - pd->i[ii]; in MatPtAPNumeric_MPIAIJ_MPIXAIJ_allatonce()
808 pdj = pd->j + pd->i[ii]; in MatPtAPNumeric_MPIAIJ_MPIXAIJ_allatonce()
[all …]
|
| H A D | mpimatmatmult.c | 1593 Mat_SeqAIJ *pd = (Mat_SeqAIJ *)p->A->data, *po = (Mat_SeqAIJ *)p->B->data; in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ() local
1672 pnz = pd->i[i + 1] - pd->i[i]; in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
1673 pdJ = pd->j + pd->i[i]; in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
1674 pA = pd->a + pd->i[i]; in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
|
| H A D | mpiaij.c | 5423 Mat_SeqAIJ *pd = (Mat_SeqAIJ *)p->A->data, *po = (Mat_SeqAIJ *)p->B->data, *p_oth; in MatCreateSeqSubMatrixWithRows_Private() local
5470 nrcols[i * 2 + 0] = pd->i[i + 1] - pd->i[i]; in MatCreateSeqSubMatrixWithRows_Private()
5540 …PetscCall(PetscSFSetGraph(sf, pd->i[plocalsize], dntotalcols, ilocal, PETSC_OWN_POINTER, iremote, … in MatCreateSeqSubMatrixWithRows_Private()
5556 for (i = 0; i < pd->i[plocalsize]; i++) pd->j[i] += pcstart; in MatCreateSeqSubMatrixWithRows_Private()
5557 PetscCall(PetscSFBcastBegin(sf, MPIU_INT, pd->j, p_oth->j, MPI_REPLACE)); in MatCreateSeqSubMatrixWithRows_Private()
5564 PetscCall(PetscSFBcastEnd(sf, MPIU_INT, pd->j, p_oth->j, MPI_REPLACE)); in MatCreateSeqSubMatrixWithRows_Private()
5566 for (i = 0; i < pd->i[plocalsize]; i++) pd->j[i] -= pcstart; in MatCreateSeqSubMatrixWithRows_Private()
|
| /petsc/lib/petsc/bin/ |
| H A D | petsc_tas_analysis.py | 23 import pandas as pd namespace
327 df = pd.read_csv(fileName)
774 if not pd.isna(field.fieldData['Errors'][0]):
|
| /petsc/share/petsc/bin/ |
| H A D | dmnetwork_view.py | 381 data = pd.read_csv(filename, skipinitialspace=True, dtype=str)
498 import pandas as pd namespace
|
| /petsc/src/ts/tutorials/ |
| H A D | ex18.c | 80 PetscInt vd, pd, d; in ProcessOptions() local
98 pd = options->porosityDist; in ProcessOptions()
99 …ty distribution type", "ex18.c", porosityDist, 5, porosityDist[options->porosityDist], &pd, NULL)); in ProcessOptions()
100 options->porosityDist = (PorosityDistribution)pd; in ProcessOptions()
|
| /petsc/src/dm/impls/plex/cgns/ |
| H A D | plexcgns2.c | 419 static PetscErrorCode PetscCGNSDataType(PetscDataType pd, CGNS_ENUMT(DataType_t) * cd) in PetscCGNSDataType() argument
422 switch (pd) { in PetscCGNSDataType()
433 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "Data type %s", PetscDataTypes[pd]); in PetscCGNSDataType()
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | petscmat.pxi | 1115 if st == "pd": return MAT_SHIFT_POSITIVE_DEFINITE
|