| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex66.c | 31 PetscScalar uu, tt; member
49 user.uu = 1.0; in main()
67 PetscScalar Hx, Hy, pi, uu, tt; in ComputeRHS() local
75 uu = user->uu; in ComputeRHS()
84 …for (i = xs; i < xs + xm; i++) array[j][i] = -PetscCosScalar(uu * pi * ((PetscReal)i + 0.5) * Hx) … in ComputeRHS()
|
| H A D | ex50.c | 35 PetscScalar uu, tt; member
53 user.uu = 1.0; in main()
71 PetscScalar Hx, Hy, pi, uu, tt; in ComputeRHS() local
79 uu = user->uu; in ComputeRHS()
88 …for (i = xs; i < xs + xm; i++) array[j][i] = -PetscCosScalar(uu * pi * ((PetscReal)i + 0.5) * Hx) … in ComputeRHS()
|
| H A D | ex2f.F90 | 114 Vec x, b, u, xx, bb, uu
340 PetscCallA(VecDuplicate(xx, uu, ierr))
341 PetscCallA(VecSet(uu, one, ierr))
342 PetscCallA(MatMult(AA, uu, bb, ierr))
385 PetscCallA(VecDestroy(uu, ierr))
|
| H A D | ex69.c | 25 PetscScalar *uu, *xx; in ComputeSolution() local
33 PetscCall(DMDAVecGetArray(da, u, &uu)); in ComputeSolution()
37 uu[j] = (xd * xd - 1.0) * PetscCosReal(5. * PETSC_PI * xd); in ComputeSolution()
40 PetscCall(DMDAVecRestoreArray(da, u, &uu)); in ComputeSolution()
|
| /petsc/src/ksp/ksp/impls/bcgs/qmrcgs/ |
| H A D | qmrcgs.c | 18 PetscScalar eta, rho1, rho2, alpha, eta2, omega, beta, cf, cf1, uu; in KSPSolve_QMRCGS() local
125 PetscCall(VecDotNorm2(S, T, &uu, &vv)); in KSPSolve_QMRCGS()
127 PetscCall(VecDot(S, S, &uu)); in KSPSolve_QMRCGS()
128 if (uu != 0.0) { in KSPSolve_QMRCGS()
151 if (uu == 0) { in KSPSolve_QMRCGS()
156 omega = uu / vv; /* omega <- uu/vv; */ in KSPSolve_QMRCGS()
|
| /petsc/src/ksp/pc/tests/ |
| H A D | ex5.c | 337 PetscScalar uu; in CalculateRhs() local
343 uu = 2.0 * h * h; in CalculateRhs()
344 PetscCall(VecSetValues(u, 1, &i, &uu, INSERT_VALUES)); in CalculateRhs()
353 PetscScalar uu; in CalculateSolution() local
360 uu = x * (1. - x); in CalculateSolution()
361 PetscCall(VecSetValues(*solution, 1, &i, &uu, INSERT_VALUES)); in CalculateSolution()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex11.c | 478 EulerNode *uu = (EulerNode *)u; in PhysicsSolution_Euler() local
484 for (i = 0; i < DIM; i++) uu->ru[i] = 0.0; /* zero out initial velocity */ in PhysicsSolution_Euler()
503 uu->r = rho * (p1 / press + gas1) / (gas1 * p1 / press + 1.0); in PhysicsSolution_Euler()
504 uu->ru[0] = ((uu->r - rho) * PetscSqrtReal(gamma * press / rho) * amach); in PhysicsSolution_Euler()
505 uu->E = p1 / (gamma - 1.0) + .5 / uu->r * uu->ru[0] * uu->ru[0]; in PhysicsSolution_Euler()
507 uu->r = 1.; /* rho = 1 */ in PhysicsSolution_Euler()
508 uu->E = p0 / (gamma - 1.0); in PhysicsSolution_Euler()
511 uu->r = eu->rhoR; in PhysicsSolution_Euler()
512 uu->E = p0 / (gamma - 1.0); in PhysicsSolution_Euler()
517 uu->r = 8.; in PhysicsSolution_Euler()
[all …]
|
| /petsc/src/dm/impls/swarm/tests/ |
| H A D | ex7.c | 16 Vec uu; member
91 PetscCall(MatCreateVecs(M_p, &matshellctx->uu, &matshellctx->ff)); in gridToParticles()
125 PetscCall(KSPSolve(ksp, rhs, matshellctx->uu)); in gridToParticles()
126 PetscCall(MatMult(M_p, matshellctx->uu, ff)); in gridToParticles()
129 PetscCall(VecDestroy(&matshellctx->uu)); in gridToParticles()
|
| /petsc/src/ksp/pc/impls/tfs/ |
| H A D | xxt.c | 170 PetscScalar *u, *uu, *v, *z, *w, alpha, alpha_w; in xxt_generate() local
231 uu = (PetscScalar *)malloc(m * sizeof(PetscScalar)); in xxt_generate()
286 PetscCall(PCTFS_rvec_zero(uu, m)); in xxt_generate()
293 PetscCallBLAS("BLASdot", uu[k] = BLASdot_(&dlen, u + off, &i1, x_ptr, &i1)); in xxt_generate()
298 PetscCall(PCTFS_ssgl_radd(uu, w, dim, stages)); in xxt_generate()
308 PetscCallBLAS("BLASaxpy", BLASaxpy_(&dlen, &uu[k], x_ptr, &i1, z + off, &i1)); in xxt_generate()
400 free(uu); in xxt_generate()
|
| H A D | xyt.c | 223 PetscScalar *u, *uu, *v, *z, *w, alpha, alpha_w; in xyt_generate() local
300 uu = (PetscScalar *)malloc(m * sizeof(PetscScalar)); in xyt_generate()
355 PetscCall(PCTFS_rvec_zero(uu, m)); in xyt_generate()
362 PetscCallBLAS("BLASdot", uu[k] = BLASdot_(&dlen, u + off, &i1, y_ptr, &i1)); in xyt_generate()
367 PetscCall(PCTFS_ssgl_radd(uu, w, dim, stages)); in xyt_generate()
377 PetscCallBLAS("BLASaxpy", BLASaxpy_(&dlen, &uu[k], x_ptr, &i1, z + off, &i1)); in xyt_generate()
509 free(uu); in xyt_generate()
|
| /petsc/src/mat/impls/h2opus/cuda/ |
| H A D | math2opus.cu | 1831 const PetscScalar *u, *v, *uu, *vv; in MatH2OpusLowRankUpdate() local
1871 uu = MatH2OpusGetThrustPointer(*a->xx); in MatH2OpusLowRankUpdate()
1872 PetscCall(PetscSFBcastBegin(usf, MPIU_SCALAR, u, (PetscScalar *)uu, MPI_REPLACE)); in MatH2OpusLowRankUpdate()
1873 PetscCall(PetscSFBcastEnd(usf, MPIU_SCALAR, u, (PetscScalar *)uu, MPI_REPLACE)); in MatH2OpusLowRankUpdate()
1874 } else uu = vv; in MatH2OpusLowRankUpdate()
1876 uu = u; in MatH2OpusLowRankUpdate()
1879 hlru_global(*a->hmatrix, uu, ldu, vv, ldv, U->cmap->N, s, handle); in MatH2OpusLowRankUpdate()
1898 uu = MatH2OpusGetThrustPointer(*a->xx_gpu); in MatH2OpusLowRankUpdate()
1899 PetscCall(PetscSFBcastBegin(usf, MPIU_SCALAR, u, (PetscScalar *)uu, MPI_REPLACE)); in MatH2OpusLowRankUpdate()
1900 PetscCall(PetscSFBcastEnd(usf, MPIU_SCALAR, u, (PetscScalar *)uu, MPI_REPLACE)); in MatH2OpusLowRankUpdate()
[all …]
|
| /petsc/src/ts/tests/ |
| H A D | ex30.c | 28 Vec uu; member
189 PetscCall(MatCreateVecs(M_p, &matshellctx->uu, &matshellctx->ff)); in gridToParticles()
217 PetscCall(VecDestroy(&matshellctx->uu)); in gridToParticles()
254 PetscCall(KSPSolve(ksp, work_ferhs, matshellctx->uu)); in gridToParticles()
256 PetscCall(MatMult(M_p, matshellctx->uu, ff)); in gridToParticles()
261 PetscCall(VecDestroy(&matshellctx->uu)); in gridToParticles()
|
| /petsc/src/ts/utils/dmplexlandau/ |
| H A D | plexland.c | 1014 static PetscErrorCode adapt(PetscInt grid, LandauCtx *ctx, Vec *uu) in adapt() argument
1023 PetscCall(adaptToleranceFEM(ctx->fe[0], *uu, type, grid, ctx, &newForest)); in adapt()
1026 PetscCall(VecDestroy(uu)); in adapt()
1027 PetscCall(DMCreateGlobalVector(newForest, uu)); in adapt()
1028 PetscCall(LandauSetInitialCondition(newForest, *uu, grid, 0, 1, ctx)); in adapt()
1031 PetscCall(PetscInfo(*uu, "No refinement\n")); in adapt()
1035 PetscCall(PetscObjectSetName((PetscObject)*uu, "uAMR")); in adapt()
|
| /petsc/src/ksp/ksp/impls/cg/pipecg2/ |
| H A D | pipecg2.c | 8 …_Private(Vec U, Vec W, Vec R, PetscInt normtype, PetscScalar *ru, PetscScalar *wu, PetscScalar *uu) in VecMergedDot_Private() argument
45 *uu = sumuu; in VecMergedDot_Private()
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | testcase3D.cas | 2205 (uu-stress/interpolate #t)
2787 (uu-stress/extrapolate/cells #t)
3551 (uu-stress/patch 1)
|