Home
last modified time | relevance | path

Searched refs:psi (Results 1 – 16 of 16) sorted by relevance

/petsc/src/ksp/ksp/tests/output/
H A Dex5_mat_type-lmvmsymbadbroyden.out19 Dual convex factor psi = 0.875
37 Dual convex factor psi = 0.875
55 Dual convex factor psi = 0.875
73 Dual convex factor psi = 0.875
91 Dual convex factor psi = 0.875
109 Dual convex factor psi = 0.875
127 Dual convex factor psi = 0.875
145 Dual convex factor psi = 0.875
163 Dual convex factor psi = 0.875
181 Dual convex factor psi = 0.875
[all …]
/petsc/src/ksp/ksp/impls/tfqmr/
H A Dtfqmr.c15 PetscReal dp, dpold, w, dpest, tau, psi, cm; in KSPSolve_TFQMR() local
80 psi = w / tau; in KSPSolve_TFQMR()
81 cm = 1.0 / PetscSqrtReal(1.0 + psi * psi); in KSPSolve_TFQMR()
82 tau = tau * psi * cm; in KSPSolve_TFQMR()
103 psiold = psi; in KSPSolve_TFQMR()
/petsc/src/dm/dt/space/impls/subspace/
H A Dspacesubspace.c127 PetscReal *phi, *psi; in PetscSpaceEvaluate_Subspace() local
130 PetscCall(DMGetWorkArray(sp->dm, origNc * subDim * subDim, MPIU_REAL, &psi)); in PetscSpaceEvaluate_Subspace()
137 for (k = 0; k < origNc * subDim; k++) psi[k] = 0.; in PetscSpaceEvaluate_Subspace()
146 …for (o = 0; o < origDim; o++) psi[(k * subDim + l) * subDim + n] += subsp->Jx[m * subDim + l] * su… in PetscSpaceEvaluate_Subspace()
154 …for (m = 0; m < PetscMin(subDim, origDim); m++) psi[(k * subDim + l) * subDim + m] += phi[(k * ori… in PetscSpaceEvaluate_Subspace()
161 …Nb + i) * subNc + k) * subDim * subDim + m] += subsp->Ju[k * origNc + l] * psi[l * subDim * subDim… in PetscSpaceEvaluate_Subspace()
166 …m * subDim; l++) H[((j * subNb + i) * subNc + k) * subDim * subDim + l] += psi[k * subDim * subDim… in PetscSpaceEvaluate_Subspace()
171 PetscCall(DMRestoreWorkArray(sp->dm, subNc * origDim, MPIU_REAL, &psi)); in PetscSpaceEvaluate_Subspace()
176 PetscReal *phi, *psi; in PetscSpaceEvaluate_Subspace() local
179 PetscCall(DMGetWorkArray(sp->dm, origNc * subDim, MPIU_REAL, &psi)); in PetscSpaceEvaluate_Subspace()
[all …]
/petsc/src/snes/tutorials/
H A Dex7.c193 …ion(PetscInt d, PetscBool forward, const PetscScalar U[], const PetscScalar psi[], PetscScalar f[]) in ComputeAction() argument
200 if (forward) DMPlex_Mult3D_Internal(U, 1, &psi[beta * 3], &tmp[beta * 3]); in ComputeAction()
201 else DMPlex_MultTranspose3D_Internal(U, 1, &psi[beta * 3], &tmp[beta * 3]); in ComputeAction()
400 Vec psi, psiHat; in TestFreeField() local
420 PetscCall(DMGetGlobalVector(dm, &psi)); in TestFreeField()
421 PetscCall(PetscObjectSetName((PetscObject)psi, "psi")); in TestFreeField()
432 if (constRhs) PetscCall(VecSet(psi, 1.)); in TestFreeField()
433 else PetscCall(VecSetRandom(psi, r)); in TestFreeField()
447 PetscCall(ComputeResidual(dm, psi, eta)); in TestFreeField()
448 PetscCall(VecViewFromOptions(psi, NULL, "-psi_view")); in TestFreeField()
[all …]
H A Dex9.c34 PetscReal psi(PetscReal x, PetscReal y) in psi() function
60 return (r <= afree) ? psi(x, y) /* active set; on the obstacle */ in u_exact()
157 aXl[j][i] = psi(x, y); in FormBounds()
/petsc/src/tao/complementarity/impls/ssls/
H A Dssils.c38 PetscReal psi, ndpsi, normd, innerd, t = 0; in TaoSolve_SSILS() local
55 …PetscCall(TaoLineSearchComputeObjectiveAndGradient(tao->linesearch, tao->solution, &psi, ssls->dps… in TaoSolve_SSILS()
91 …PetscCall(TaoLineSearchApply(tao->linesearch, tao->solution, &psi, ssls->dpsi, tao->stepdirection,… in TaoSolve_SSILS()
H A Dssfls.c24 PetscReal psi, ndpsi, normd, innerd, t = 0; in TaoSolve_SSFLS() local
42 …PetscCall(TaoLineSearchComputeObjectiveAndGradient(tao->linesearch, tao->solution, &psi, ssls->dps… in TaoSolve_SSFLS()
83 …PetscCall(TaoLineSearchApply(tao->linesearch, tao->solution, &psi, ssls->dpsi, tao->stepdirection,… in TaoSolve_SSFLS()
/petsc/src/ts/impls/implicit/glle/
H A Dglle.c128 … r, &scheme->beta, r, &scheme->gamma, 3 * s, &scheme->phi, 3 * r, &scheme->psi, r, &scheme->stage_… in TSGLLESchemeCreate()
253 scheme->psi[0 * r + 0] = 0.; in TSGLLESchemeCreate()
254 scheme->psi[1 * r + 0] = 0.; in TSGLLESchemeCreate()
255 scheme->psi[2 * r + 0] = 0.; in TSGLLESchemeCreate()
257 scheme->psi[0 * r + j] = 0.; in TSGLLESchemeCreate()
258 scheme->psi[1 * r + j] = 0.; in TSGLLESchemeCreate()
259 scheme->psi[2 * r + j] = 0.; in TSGLLESchemeCreate()
261 scheme->psi[0 * r + j] -= CPowF(c[k], j - 1) * scheme->phi[0 * s + k]; in TSGLLESchemeCreate()
262 scheme->psi[1 * r + j] -= CPowF(c[k], j - 1) * scheme->phi[1 * s + k]; in TSGLLESchemeCreate()
263 scheme->psi[2 * r + j] -= CPowF(c[k], j - 1) * scheme->phi[2 * s + k]; in TSGLLESchemeCreate()
[all …]
H A Dglle.h29 PetscScalar *psi; /* dim=[3][r-1], [0 psi^T] of B,J,W 2007 */ member
/petsc/src/ksp/ksp/tutorials/
H A Dex54f.F90263 PetscReal xx, yy, psi, cs, sn, c2, s2, dd(2, 2)
271 psi = dir(xx, yy)
273 cs = cos(psi)
274 sn = sin(psi)
/petsc/src/ksp/ksp/utils/lmvm/symbrdn/
H A Dsymbrdn.c6 static inline PetscScalar PhiFromPsi(PetscScalar psi, PetscScalar yts, PetscScalar stBs, PetscScala… in PhiFromPsi() argument
8 PetscScalar numer = (1.0 - psi) * PetscRealPart(PetscConj(yts) * yts); in PhiFromPsi()
9 PetscScalar phi = numer / (numer + psi * stBs * ytHy); in PhiFromPsi()
966 PetscErrorCode MatLMVMSymBadBroydenGetPsi(Mat B, PetscReal *psi) in MatLMVMSymBadBroydenGetPsi() argument
969 *psi = PETSC_DETERMINE; in MatLMVMSymBadBroydenGetPsi()
970 PetscTryMethod(B, "MatLMVMSymBadBroydenGetPsi_C", (Mat, PetscReal *), (B, psi)); in MatLMVMSymBadBroydenGetPsi()
974 static PetscErrorCode MatLMVMSymBadBroydenGetPsi_SymBrdn(Mat B, PetscReal *psi) in MatLMVMSymBadBroydenGetPsi_SymBrdn() argument
980 *psi = lsb->psi_scalar; in MatLMVMSymBadBroydenGetPsi_SymBrdn()
1005 PetscErrorCode MatLMVMSymBadBroydenSetPsi(Mat B, PetscReal psi) in MatLMVMSymBadBroydenSetPsi() argument
1008 PetscTryMethod(B, "MatLMVMSymBadBroydenSetPsi_C", (Mat, PetscReal), (B, psi)); in MatLMVMSymBadBroydenSetPsi()
[all …]
/petsc/src/tao/complementarity/impls/asls/
H A Dasils.c114 PetscReal psi, ndpsi, normd, innerd, t = 0; in TaoSolve_ASILS() local
128 …PetscCall(TaoLineSearchComputeObjectiveAndGradient(tao->linesearch, tao->solution, &psi, asls->dps… in TaoSolve_ASILS()
254 …PetscCall(TaoLineSearchApply(tao->linesearch, tao->solution, &psi, asls->dpsi, tao->stepdirection,… in TaoSolve_ASILS()
H A Dasfls.c114 PetscReal psi, ndpsi, normd, innerd, t = 0; in TaoSolve_ASFLS() local
131 …PetscCall(TaoLineSearchComputeObjectiveAndGradient(tao->linesearch, tao->solution, &psi, asls->dps… in TaoSolve_ASFLS()
259 …PetscCall(TaoLineSearchApply(tao->linesearch, tao->solution, &psi, asls->dpsi, tao->stepdirection,… in TaoSolve_ASFLS()
/petsc/src/ksp/pc/impls/bjacobi/bjkokkos/
H A Dbjkokkos.kokkos.cxx93 PetscReal dp, dpold, w, dpest, tau, psi, cm, r0; in BJSolve_TFQMR() local
233 psi = w / tau; in BJSolve_TFQMR()
234 cm = 1.0 / PetscSqrtReal(1.0 + psi * psi); in BJSolve_TFQMR()
235 tau = tau * psi * cm; in BJSolve_TFQMR()
280 psiold = psi; in BJSolve_TFQMR()
/petsc/doc/manual/
H A Dsnes.md533 % math::L^2 = \|\Delta x\|^2 + \psi^2 (\Delta\lambda)^2,
537 accumulated updates over the current load step, and $\psi^2$ is a
539 Generally, $\psi^2 > 0$ leads to a hyper-sphere constraint surface, while
540 $\psi^2 = 0$ leads to a hyper-cylinder constraint surface.
565 a &= \|\delta\mathbf x^Q\|^2 + \psi^2,\\
566 b &= 2\delta\mathbf x^Q\cdot (\Delta\mathbf x + \delta s\delta\mathbf x^F) + 2\psi^2 \Delta\lambda,…
567 c &= \|\Delta\mathbf x + \delta s\delta\mathbf x^F\|^2 + \psi^2 \Delta\lambda^2 - L^2.
576 \delta\lambda = \pm\frac{L}{\sqrt{\|\delta\mathbf x^Q\|^2 + \psi^2}},
583 \text{sign}(\delta\lambda) = \text{sign}\big(\delta\mathbf x^Q \cdot (\Delta\mathbf x)_{i-1} + \psi…
605 …mathbf x \cdot \delta \mathbf x^F}{\Delta\mathbf x \cdot \delta\mathbf x^Q + \psi^2 \Delta\lambda}.
/petsc/src/dm/impls/plex/tests/
H A Dex8.c271 PetscReal scale, phi, theta, psi = 0.0, norm; in CheckCell() local
289 const PetscReal cs = PetscCosReal(psi), ss = PetscSinReal(psi); in CheckCell()