Lines Matching refs:p
25 …cErrorCode PetscPDFMaxwellBoltzmann1D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscPDFMaxwellBoltzmann1D() argument
27 p[0] = PetscSqrtReal(2. / PETSC_PI) * PetscExpReal(-0.5 * PetscSqr(x[0])); in PetscPDFMaxwellBoltzmann1D()
47 …cErrorCode PetscCDFMaxwellBoltzmann1D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscCDFMaxwellBoltzmann1D() argument
49 p[0] = PetscErfReal(x[0] / PETSC_SQRT2); in PetscCDFMaxwellBoltzmann1D()
69 …cErrorCode PetscPDFMaxwellBoltzmann2D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscPDFMaxwellBoltzmann2D() argument
71 p[0] = x[0] * PetscExpReal(-0.5 * PetscSqr(x[0])); in PetscPDFMaxwellBoltzmann2D()
91 …cErrorCode PetscCDFMaxwellBoltzmann2D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscCDFMaxwellBoltzmann2D() argument
93 p[0] = 1. - PetscExpReal(-0.5 * PetscSqr(x[0])); in PetscCDFMaxwellBoltzmann2D()
113 …cErrorCode PetscPDFMaxwellBoltzmann3D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscPDFMaxwellBoltzmann3D() argument
115 p[0] = PetscSqrtReal(2. / PETSC_PI) * PetscSqr(x[0]) * PetscExpReal(-0.5 * PetscSqr(x[0])); in PetscPDFMaxwellBoltzmann3D()
135 …cErrorCode PetscCDFMaxwellBoltzmann3D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscCDFMaxwellBoltzmann3D() argument
137 …p[0] = PetscErfReal(x[0] / PETSC_SQRT2) - PetscSqrtReal(2. / PETSC_PI) * x[0] * PetscExpReal(-0.5 … in PetscCDFMaxwellBoltzmann3D()
157 PetscErrorCode PetscPDFGaussian1D(const PetscReal x[], const PetscReal scale[], PetscReal p[]) in PetscPDFGaussian1D() argument
160 …p[0] = PetscSqrtReal(1. / (2. * PETSC_PI)) * PetscExpReal(-0.5 * PetscSqr(x[0] / … in PetscPDFGaussian1D()
164 PetscErrorCode PetscCDFGaussian1D(const PetscReal x[], const PetscReal scale[], PetscReal p[]) in PetscCDFGaussian1D() argument
167 p[0] = 0.5 * (1. + PetscErfReal(x[0] / PETSC_SQRT2 / sigma)); in PetscCDFGaussian1D()
191 PetscErrorCode PetscPDFSampleGaussian1D(const PetscReal p[], const PetscReal unused[], PetscReal x[… in PetscPDFSampleGaussian1D() argument
193 const PetscReal q = 2 * p[0] - 1.; in PetscPDFSampleGaussian1D()
234 PetscErrorCode PetscPDFGaussian2D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscPDFGaussian2D() argument
236 p[0] = (1. / PETSC_PI) * PetscExpReal(-0.5 * (PetscSqr(x[0]) + PetscSqr(x[1]))); in PetscPDFGaussian2D()
257 PetscErrorCode PetscPDFSampleGaussian2D(const PetscReal p[], const PetscReal unused[], PetscReal x[… in PetscPDFSampleGaussian2D() argument
259 const PetscReal mag = PetscSqrtReal(-2.0 * PetscLogReal(p[0])); in PetscPDFSampleGaussian2D()
260 x[0] = mag * PetscCosReal(2.0 * PETSC_PI * p[1]); in PetscPDFSampleGaussian2D()
261 x[1] = mag * PetscSinReal(2.0 * PETSC_PI * p[1]); in PetscPDFSampleGaussian2D()
281 PetscErrorCode PetscPDFGaussian3D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscPDFGaussian3D() argument
283 …p[0] = (1. / PETSC_PI * PetscSqrtReal(PETSC_PI)) * PetscExpReal(-0.5 * (PetscSqr(x[0]) + PetscSqr(… in PetscPDFGaussian3D()
304 PetscErrorCode PetscPDFSampleGaussian3D(const PetscReal p[], const PetscReal unused[], PetscReal x[… in PetscPDFSampleGaussian3D() argument
306 PetscCall(PetscPDFSampleGaussian1D(p, unused, x)); in PetscPDFSampleGaussian3D()
307 PetscCall(PetscPDFSampleGaussian2D(&p[1], unused, &x[1])); in PetscPDFSampleGaussian3D()
327 PetscErrorCode PetscPDFConstant1D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscPDFConstant1D() argument
329 p[0] = x[0] > -1. && x[0] <= 1. ? 0.5 : 0.; in PetscPDFConstant1D()
349 PetscErrorCode PetscCDFConstant1D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscCDFConstant1D() argument
351 p[0] = x[0] <= -1. ? 0. : (x[0] > 1. ? 1. : 0.5 * (x[0] + 1.)); in PetscCDFConstant1D()
371 PetscErrorCode PetscPDFSampleConstant1D(const PetscReal p[], const PetscReal unused[], PetscReal x[… in PetscPDFSampleConstant1D() argument
373 x[0] = 2. * p[0] - 1.; in PetscPDFSampleConstant1D()
393 PetscErrorCode PetscPDFConstant2D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscPDFConstant2D() argument
395 p[0] = x[0] > -1. && x[0] <= 1. && x[1] > -1. && x[1] <= 1. ? 0.25 : 0.; in PetscPDFConstant2D()
415 PetscErrorCode PetscCDFConstant2D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscCDFConstant2D() argument
417 …p[0] = x[0] <= -1. || x[1] <= -1. ? 0. : (x[0] > 1. ? 1. : 0.5 * (x[0] + 1.)) * (x[1] > 1. ? 1. : … in PetscCDFConstant2D()
437 PetscErrorCode PetscPDFSampleConstant2D(const PetscReal p[], const PetscReal unused[], PetscReal x[… in PetscPDFSampleConstant2D() argument
439 x[0] = 2. * p[0] - 1.; in PetscPDFSampleConstant2D()
440 x[1] = 2. * p[1] - 1.; in PetscPDFSampleConstant2D()
460 PetscErrorCode PetscPDFConstant3D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscPDFConstant3D() argument
462 …p[0] = x[0] > -1. && x[0] <= 1. && x[1] > -1. && x[1] <= 1. && x[2] > -1. && x[2] <= 1. ? 0.125 : … in PetscPDFConstant3D()
482 PetscErrorCode PetscCDFConstant3D(const PetscReal x[], const PetscReal unused[], PetscReal p[]) in PetscCDFConstant3D() argument
484 …p[0] = x[0] <= -1. || x[1] <= -1. || x[2] <= -1. ? 0. : (x[0] > 1. ? 1. : 0.5 * (x[0] + 1.)) * (x[… in PetscCDFConstant3D()
504 PetscErrorCode PetscPDFSampleConstant3D(const PetscReal p[], const PetscReal unused[], PetscReal x[… in PetscPDFSampleConstant3D() argument
506 x[0] = 2. * p[0] - 1.; in PetscPDFSampleConstant3D()
507 x[1] = 2. * p[1] - 1.; in PetscPDFSampleConstant3D()
508 x[2] = 2. * p[2] - 1.; in PetscPDFSampleConstant3D()
703 for (PetscInt p = 0; p < n; ++p) { in PetscProbComputeKSStatistic_Internal()
704 const PetscReal x = val[p]; in PetscProbComputeKSStatistic_Internal()
705 const PetscReal w = wgt ? wgt[p] : 1. / n; in PetscProbComputeKSStatistic_Internal()
939 for (PetscInt p = 0; p < n; ++p) { in PetscProbComputeKSStatisticMagnitude() local
942 for (PetscInt d = 0; d < dim; ++d) mag += PetscSqr(PetscRealPart(a[p * dim + d])); in PetscProbComputeKSStatisticMagnitude()
943 speed[p] = PetscSqrtReal(mag); in PetscProbComputeKSStatisticMagnitude()