Home
last modified time | relevance | path

Searched refs:phi (Results 1 – 25 of 43) sorted by relevance

12

/petsc/src/dm/impls/moab/
H A Ddmmbfem.cxx110 …npts, const PetscReal *quad, PetscReal *phypts, PetscReal *jxw, PetscReal *phi, PetscReal *dphidx,… in Compute_Lagrange_Basis_1D_Internal() argument
128 phi[0 + offset] = (1.0 - r); in Compute_Lagrange_Basis_1D_Internal()
129 phi[1 + offset] = (r); in Compute_Lagrange_Basis_1D_Internal()
137 if (phypts) phypts[3 * j + 0] += phi[i + offset] * vertices[0]; in Compute_Lagrange_Basis_1D_Internal()
156 phi[0 + offset] = 1.0 + r * (2.0 * r - 3.0); in Compute_Lagrange_Basis_1D_Internal()
157 phi[1 + offset] = 4.0 * r * (1.0 - r); in Compute_Lagrange_Basis_1D_Internal()
158 phi[2 + offset] = r * (2.0 * r - 1.0); in Compute_Lagrange_Basis_1D_Internal()
166 if (phypts) phypts[3 * j + 0] += phi[i + offset] * vertices[0]; in Compute_Lagrange_Basis_1D_Internal()
220 …npts, const PetscReal *quad, PetscReal *phypts, PetscReal *jxw, PetscReal *phi, PetscReal *dphidx,… in Compute_Lagrange_Basis_2D_Internal() argument
228 PetscCall(PetscArrayzero(phi, npts)); in Compute_Lagrange_Basis_2D_Internal()
[all …]
/petsc/src/ksp/ksp/tutorials/amrex/
H A DMyTest.cxx55 for (int ilev = 0; ilev <= max_level; ++ilev) mleb.setLevelBC(ilev, &phi[ilev]); in solve()
65 …for (int ilev = 0; ilev <= max_level; ++ilev) mleb.setEBDirichlet(ilev, phi[ilev], bcoef_eb[ilev]); in solve()
79 mlmg.solve(amrex::GetVecOfPtrs(phi), amrex::GetVecOfConstPtrs(rhs), tol_rel, tol_abs); in solve()
88 MultiFab::Copy(plotmf[ilev], phi[ilev], 0, 0, 1, 0); in writePlotfile()
162 phi.resize(nlevels); in initData()
174 phi[ilev].define(grids[ilev], dmap[ilev], 1, 1, MFInfo(), *factory[ilev]); in initData()
183 phi[ilev].setVal(0.0); in initData()
203 phi[ilev].setVal(10.0); in initData()
204 phi[ilev].setVal(0.0, 0, 1, 0); // set interior in initData()
207 for (MFIter mfi(phi[ilev]); mfi.isValid(); ++mfi) { in initData()
[all …]
H A DMyTest.H51 amrex::Vector<amrex::MultiFab> phi; variable
/petsc/src/ksp/ksp/tests/output/
H A Dex5_mat_type-lmvmsymbroyden.out19 Convex factor phi = 0.125
37 Convex factor phi = 0.125
55 Convex factor phi = 0.125
73 Convex factor phi = 0.125
91 Convex factor phi = 0.125
109 Convex factor phi = 0.125
127 Convex factor phi = 0.125
145 Convex factor phi = 0.125
163 Convex factor phi = 0.125
181 Convex factor phi = 0.125
[all …]
/petsc/src/snes/impls/vi/ss/
H A Dviss.c17 PetscErrorCode SNESVIComputeMeritFunction(Vec phi, PetscReal *merit, PetscReal *phinorm) in SNESVIComputeMeritFunction() argument
20 PetscCall(VecNormBegin(phi, NORM_2, phinorm)); in SNESVIComputeMeritFunction()
21 PetscCall(VecNormEnd(phi, NORM_2, phinorm)); in SNESVIComputeMeritFunction()
53 PetscErrorCode SNESVIComputeFunction(SNES snes, Vec X, Vec phi, void *functx) in SNESVIComputeFunction() argument
68 PetscCall(VecGetArray(phi, &phi_arr)); in SNESVIComputeFunction()
88 PetscCall(VecRestoreArray(phi, &phi_arr)); in SNESVIComputeFunction()
188 static PetscErrorCode SNESVIComputeMeritFunctionGradient(Mat H, Vec phi, Vec dpsi) in SNESVIComputeMeritFunctionGradient() argument
191 PetscCall(MatMultTranspose(H, phi, dpsi)); in SNESVIComputeMeritFunctionGradient()
228 PetscCall(SNESComputeFunction(snes, X, vi->phi)); in SNESSolve_VINEWTONSSLS()
236 PetscCall(SNESVIComputeMeritFunction(vi->phi, &vi->merit, &vi->phinorm)); in SNESSolve_VINEWTONSSLS()
[all …]
H A Dvissimpl.h13 Vec phi; /* pointer to semismooth function */ member
/petsc/src/dm/dt/space/impls/subspace/
H A Dspacesubspace.c127 PetscReal *phi, *psi; in PetscSpaceEvaluate_Subspace() local
129 PetscCall(DMGetWorkArray(sp->dm, origNc * origDim * origDim, MPIU_REAL, &phi)); in PetscSpaceEvaluate_Subspace()
136 for (k = 0; k < origNc * origDim; k++) phi[k] = 0.; in PetscSpaceEvaluate_Subspace()
139 …for (l = 0; l < origNc * origDim * origDim; l++) phi[l] += inH[(j * origNb + k) * origNc * origDim… in PetscSpaceEvaluate_Subspace()
146 …) * subDim + n] += subsp->Jx[m * subDim + l] * subsp->Jx[o * subDim + n] * phi[(k * origDim + m) *… in PetscSpaceEvaluate_Subspace()
154 … m < PetscMin(subDim, origDim); m++) psi[(k * subDim + l) * subDim + m] += phi[(k * origDim + l) *… in PetscSpaceEvaluate_Subspace()
172 PetscCall(DMRestoreWorkArray(sp->dm, origNc * origDim, MPIU_REAL, &phi)); in PetscSpaceEvaluate_Subspace()
176 PetscReal *phi, *psi; in PetscSpaceEvaluate_Subspace() local
178 PetscCall(DMGetWorkArray(sp->dm, origNc * origDim, MPIU_REAL, &phi)); in PetscSpaceEvaluate_Subspace()
185 for (k = 0; k < origNc * origDim; k++) phi[k] = 0.; in PetscSpaceEvaluate_Subspace()
[all …]
/petsc/src/snes/tutorials/
H A Dex33.c10 PetscReal phi; member
67 PetscReal phi = user->phi; in FormFunctionLocal() local
96 …f[i].s = phi * (u[i].s - uold[i].s) + (dt / dx) * ((lambdaWet / lambda) * u[i].v - (lambdaWetL / l… in FormFunctionLocal()
146 user.phi = 1.0; in main()
H A Dex48.c204 … PetscReal hx, PetscReal hy, const PetscReal dz[PETSC_RESTRICT], PetscReal phi[PETSC_RESTRICT], Pe… in HexComputeGeometry()
221 phi[i] = HexQInterp[q][i]; in HexComputeGeometry()
710 static void PointwiseNonlinearity(THI thi, const Node n[PETSC_RESTRICT], const PetscReal phi[PETSC_… in PointwiseNonlinearity()
720 *u += phi[l] * n[l].u; in PointwiseNonlinearity()
721 *v += phi[l] * n[l].v; in PointwiseNonlinearity()
731 static void PointwiseNonlinearity2D(THI thi, Node n[], PetscReal phi[], PetscReal dphi[4][2], Petsc… in PointwiseNonlinearity2D() argument
741 *u += phi[l] * n[l].u; in PointwiseNonlinearity2D()
742 *v += phi[l] * n[l].v; in PointwiseNonlinearity2D()
791 PetscReal dz[3], phi[8], dphi[8][3], jw, eta, deta; in THIFunctionLocal() local
794 HexComputeGeometry(q, hx, hy, dz, phi, dphi, &jw); in THIFunctionLocal()
[all …]
/petsc/src/ksp/ksp/utils/lmvm/symbrdn/
H A Dsymbrdn.c9 PetscScalar phi = numer / (numer + psi * stBs * ytHy); in PhiFromPsi() local
10 return phi; in PhiFromPsi()
65 PetscScalar phi = Phi_t == SYMBROYDEN_PRODUCTS_PHI ? lsb->phi_scalar : lsb->psi_scalar; in MatLMVMSymBroydenGetConvexFactor() local
70 if (phi != PETSC_DETERMINE) { in MatLMVMSymBroydenGetConvexFactor()
76 …for (PetscInt i = start; i < next; i++) PetscCall(LMProductsInsertNextDiagonalValue(*Phi, i, phi)); in MatLMVMSymBroydenGetConvexFactor()
261 PetscScalar phi; in SymBroydenCompactDenseUpdateArrays() local
279 phi = PhiFromPsi(Psi[i], yts, sBis, YtHkY[i]); in SymBroydenCompactDenseUpdateArrays()
281 phi = Phi[i]; in SymBroydenCompactDenseUpdateArrays()
283 alpha = PetscRealPart(-(1.0 - phi) / sBis); in SymBroydenCompactDenseUpdateArrays()
284 beta = -phi / yts; in SymBroydenCompactDenseUpdateArrays()
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex36.cxx265 PetscReal *phi, *phypts, *jxw; in ComputeRHS_MOAB() local
277 PetscCall(PetscMalloc3(user->VPERE * npoints, &phi, npoints * 3, &phypts, npoints, &jxw)); in ComputeRHS_MOAB()
299 PetscCall(DMMoabFEMComputeBasis(user->dim, nconn, vpos, quadratureObj, phypts, jxw, phi, NULL)); in ComputeRHS_MOAB()
307 for (i = 0; i < nconn; ++i) localv[i] += jxw[q] * phi[offset + i] * ff; in ComputeRHS_MOAB()
343 PetscCall(PetscFree3(phi, phypts, jxw)); in ComputeRHS_MOAB()
361 PetscReal *phi, *dphi[3], *phypts, *jxw; in ComputeMatrix_MOAB() local
376 …PetscCall(PetscMalloc6(user->VPERE * npoints, &phi, user->VPERE * npoints, &dphi[0], user->VPERE *… in ComputeMatrix_MOAB()
394 PetscCall(DMMoabFEMComputeBasis(user->dim, nconn, vpos, quadratureObj, phypts, jxw, phi, dphi)); in ComputeMatrix_MOAB()
409 …fset + j] + dphi[2][offset + i] * dphi[2][offset + j]) + alpha * (phi[offset + i] * phi[offset + j… in ComputeMatrix_MOAB()
448 PetscCall(PetscFree6(phi, dphi[0], dphi[1], dphi[2], phypts, jxw)); in ComputeMatrix_MOAB()
H A Dex35.cxx294 PetscReal *phi, *phypts, *jxw; in ComputeRHS() local
306 PetscCall(PetscMalloc3(user->VPERE * npoints, &phi, npoints * 3, &phypts, npoints, &jxw)); in ComputeRHS()
334 PetscCall(DMMoabFEMComputeBasis(2, nconn, vpos, quadratureObj, phypts, jxw, phi, NULL)); in ComputeRHS()
340 for (i = 0; i < nconn; ++i) localv[i] += jxw[q] * phi[q * nconn + i] * ff; in ComputeRHS()
379 PetscCall(PetscFree3(phi, phypts, jxw)); in ComputeRHS()
397 PetscReal *phi, *dphi[2], *phypts, *jxw; in ComputeMatrix() local
412 …PetscCall(PetscMalloc5(user->VPERE * npoints, &phi, user->VPERE * npoints, &dphi[0], user->VPERE *… in ComputeMatrix()
434 PetscCall(DMMoabFEMComputeBasis(2, nconn, vpos, quadratureObj, phypts, jxw, phi, dphi)); in ComputeMatrix()
488 PetscCall(PetscFree5(phi, dphi[0], dphi[1], phypts, jxw)); in ComputeMatrix()
/petsc/share/petsc/datafiles/meshes/
H A Dgmsh-B2tri.geo22 phi = Pi/180 * GetValue("Rotation (degrees)", 0);
23 Rotate {{0, 0, 1}, {Cx, Cy, 0}, phi} { Point{1:9}; }
/petsc/src/ksp/ksp/utils/lmvm/tests/
H A Dex2.c11 PetscReal phi = 0.618; in main() local
36 if (is_sym) PetscCall(MatLMVMSymBroydenSetPhi(sym, phi)); in main()
37 if (is_symbad) PetscCall(MatLMVMSymBadBroydenSetPsi(sym, phi)); in main()
H A Dex1.c231 static PetscErrorCode SymmetricBroydenUpdate_Explicit(Mat B, PetscReal phi, Vec s, Vec y) in SymmetricBroydenUpdate_Explicit() argument
241 PetscCall(RankOneAXPY(B, (yts + phi * stBs) / (yts * yts), y, y)); in SymmetricBroydenUpdate_Explicit()
242 PetscCall(RankOneAXPY(B, -phi / yts, y, Bs)); in SymmetricBroydenUpdate_Explicit()
243 PetscCall(RankOneAXPY(B, -phi / yts, Bs, y)); in SymmetricBroydenUpdate_Explicit()
244 PetscCall(RankOneAXPY(B, (phi - 1.0) / stBs, Bs, Bs)); in SymmetricBroydenUpdate_Explicit()
315 …te)(Mat, PetscReal, Vec, Vec), PetscErrorCode (*H_update)(Mat, PetscReal, Vec, Vec), PetscReal phi) in TestUpdate() argument
405 for (PetscInt i = oldest; i < next; i++) PetscCall((*B_update)(B_k_exp, phi, dxs[i], dfs[i])); in TestUpdate()
426 for (PetscInt i = oldest; i < next; i++) PetscCall((*H_update)(H_k_exp, phi, dfs[i], dxs[i])); in TestUpdate()
489 PetscReal phi = 0.0; in main() local
589 PetscCall(MatLMVMSymBroydenGetPhi(B, &phi)); in main()
[all …]
/petsc/src/dm/impls/swarm/tests/
H A Dex6.c269 Vec phi, locPhi, rho, f; in ComputeFieldAtParticles_Primal() local
316 PetscCall(DMGetGlobalVector(dm, &phi)); in ComputeFieldAtParticles_Primal()
317 PetscCall(PetscObjectSetName((PetscObject)phi, "potential")); in ComputeFieldAtParticles_Primal()
318 PetscCall(VecSet(phi, 0.0)); in ComputeFieldAtParticles_Primal()
319 PetscCall(SNESSolve(snes, rho, phi)); in ComputeFieldAtParticles_Primal()
321 PetscCall(VecViewFromOptions(phi, NULL, "-phi_view")); in ComputeFieldAtParticles_Primal()
324 PetscCall(DMGlobalToLocalBegin(dm, phi, INSERT_VALUES, locPhi)); in ComputeFieldAtParticles_Primal()
325 PetscCall(DMGlobalToLocalEnd(dm, phi, INSERT_VALUES, locPhi)); in ComputeFieldAtParticles_Primal()
326 PetscCall(DMRestoreGlobalVector(dm, &phi)); in ComputeFieldAtParticles_Primal()
378 Vec phi, locPhi, rho, f, temp_rho; in ComputeFieldAtParticles_Mixed() local
[all …]
/petsc/src/ts/tutorials/power_grid/output/
H A Dex3sa_2.out2 trajectory sensitivity: d[phi(tf)]/d[pm] d[omega(tf)]/d[pm]
/petsc/src/dm/dt/fv/interface/
H A Dfv.c333 PetscErrorCode PetscLimiterLimit(PetscLimiter lim, PetscReal flim, PetscReal *phi) in PetscLimiterLimit() argument
337 PetscAssertPointer(phi, 3); in PetscLimiterLimit()
338 PetscUseTypeMethod(lim, limit, flim, phi); in PetscLimiterLimit()
373 static PetscErrorCode PetscLimiterLimit_Sin(PetscLimiter lim, PetscReal f, PetscReal *phi) in PetscLimiterLimit_Sin() argument
376 *phi = PetscSinReal(PETSC_PI * PetscMax(0, PetscMin(f, 1))); in PetscLimiterLimit_Sin()
441 static PetscErrorCode PetscLimiterLimit_Zero(PetscLimiter lim, PetscReal f, PetscReal *phi) in PetscLimiterLimit_Zero() argument
444 *phi = 0.0; in PetscLimiterLimit_Zero()
509 static PetscErrorCode PetscLimiterLimit_None(PetscLimiter lim, PetscReal f, PetscReal *phi) in PetscLimiterLimit_None() argument
512 *phi = 1.0; in PetscLimiterLimit_None()
577 static PetscErrorCode PetscLimiterLimit_Minmod(PetscLimiter lim, PetscReal f, PetscReal *phi) in PetscLimiterLimit_Minmod() argument
[all …]
/petsc/src/ts/tutorials/hamiltonian/
H A Dex4.c424 Vec rho, rhohat, phi; in SetupContext() local
473 PetscCall(DMGetNamedGlobalVector(ctx->dmPot, "phi", &phi)); in SetupContext()
474 PetscCall(PetscObjectSetName((PetscObject)phi, "potential")); in SetupContext()
475 PetscCall(DMRestoreNamedGlobalVector(ctx->dmPot, "phi", &phi)); in SetupContext()
646 PetscReal C, gamma, omega, phi; in ComputeEmaxResidual() local
654 phi = PetscRealPart(a[3]); in ComputeEmaxResidual()
656 …x[i]) - C * PetscExpReal(-gamma * ctx->t[i]) * PetscAbsReal(PetscCosReal(omega * ctx->t[i] - phi)); in ComputeEmaxResidual()
667 PetscReal C, gamma, omega, phi; in ComputeEmaxJacobian() local
675 phi = PetscRealPart(a[3]); in ComputeEmaxJacobian()
681 …* 4 + 0] = -PetscExpReal(-gamma * ctx->t[k]) * PetscAbsReal(PetscCosReal(omega * ctx->t[k] - phi)); in ComputeEmaxJacobian()
[all …]
H A Dex2.c339 Vec rho, rhohat, phi; in SetupContext() local
388 PetscCall(DMGetNamedGlobalVector(ctx->dmPot, "phi", &phi)); in SetupContext()
389 PetscCall(PetscObjectSetName((PetscObject)phi, "potential")); in SetupContext()
390 PetscCall(DMRestoreNamedGlobalVector(ctx->dmPot, "phi", &phi)); in SetupContext()
550 Vec phi; in MonitorEField() local
552 PetscCall(DMGetNamedGlobalVector(ctx->dmPot, "phi", &phi)); in MonitorEField()
555 PetscCall(DMPlexComputeIntegralFEM(ctx->dmPot, phi, &intESq, ctx)); in MonitorEField()
556 PetscCall(DMRestoreNamedGlobalVector(ctx->dmPot, "phi", &phi)); in MonitorEField()
843 Vec rho, rhohat, phi; in MonitorPoisson() local
859 PetscCall(DMGetNamedGlobalVector(ctx->dmPot, "phi", &phi)); in MonitorPoisson()
[all …]
/petsc/src/dm/impls/plex/tutorials/
H A Dex19.c247 PetscReal xold, yold, zold, rold, phi, theta; in main() local
254 phi = PetscAtan2Real(yold, xold); in main()
280 …cpCoordData[faceCPStartRow + (3 * ii) + 0] += -0.05 * deltaR * PetscCosReal(phi) * PetscSinReal(th… in main()
281 …cpCoordData[faceCPStartRow + (3 * ii) + 1] += -0.05 * deltaR * PetscSinReal(phi) * PetscSinReal(th… in main()
/petsc/src/tao/unconstrained/impls/bmrm/
H A Dbmrm.c6 static PetscReal phi(PetscReal *, PetscInt, PetscReal, PetscReal *, PetscReal, PetscReal *, Pe…
671 static PetscReal phi(PetscReal *x, PetscInt n, PetscReal lambda, PetscReal *a, PetscReal b, PetscRe… in phi() function
715 r = phi(x, n, lambda, a, b, c, l, u); in project()
729 r = phi(x, n, lambda, a, b, c, l, u); in project()
737 r = phi(x, n, lambda, a, b, c, l, u); in project()
745 r = phi(x, n, lambda, a, b, c, l, u); in project()
753 r = phi(x, n, lambda, a, b, c, l, u); in project()
767 r = phi(x, n, lambda, a, b, c, l, u); in project()
808 r = phi(x, n, lambda, a, b, c, l, u); in project()
/petsc/src/tao/constrained/impls/ipm/
H A Dipm.h26 PetscScalar phi; member
H A Dipm.c49 PetscCall(TaoLogConvergenceHistory(tao, ipmP->kkt_f, ipmP->phi, 0.0, tao->ksp_its)); in TaoSolve_IPM()
50 PetscCall(TaoMonitor(tao, tao->niter, ipmP->kkt_f, ipmP->phi, 0.0, 1.0)); in TaoSolve_IPM()
165 phi_target = ipmP->dec * ipmP->phi; in TaoSolve_IPM()
179 if (ipmP->phi <= phi_target) break; in TaoSolve_IPM()
183 PetscCall(TaoLogConvergenceHistory(tao, ipmP->kkt_f, ipmP->phi, 0.0, tao->ksp_its)); in TaoSolve_IPM()
184 PetscCall(TaoMonitor(tao, tao->niter, ipmP->kkt_f, ipmP->phi, 0.0, stepsize)); in TaoSolve_IPM()
594 ipmP->phi = norm; in IPMComputeKKT()
597 ipmP->phi += norm; in IPMComputeKKT()
601 ipmP->phi += norm; in IPMComputeKKT()
603 ipmP->phi += norm; in IPMComputeKKT()
[all …]
/petsc/src/ts/impls/implicit/glle/
H A Dglle.h28 PetscScalar *phi; /* dim=[3][s] for estimating higher moments, see B,J,W 2007 */ member

12