Home
last modified time | relevance | path

Searched full:rho (Results 1 – 25 of 120) sorted by relevance

12345

/petsc/src/ts/tutorials/advection-diffusion-reaction/
H A Dex4.c29 PetscScalar rho, c; member
75 PetscCall(DMDASetFieldName(da, 0, "rho")); in main()
141 PetscScalar rho, c, rhoxx, cxx, cx, rhox, kcxrhox; in IFunction() local
175 f[0].rho = udot[0].rho; /* u[0].rho - 0.0; */ in IFunction()
181 f[Mx - 1].rho = udot[Mx - 1].rho; /* u[Mx-1].rho - 1.0; */ in IFunction()
190 rho = u[i].rho; in IFunction()
191 rhoxx = (-2.0 * rho + u[i - 1].rho + u[i + 1].rho) * sx; in IFunction()
196 rhox = .5 * (u[i + 1].rho - u[i - 1].rho) / hx; in IFunction()
198 kcxrhox = appctx->kappa * (cxx * rho + cx * rhox); in IFunction()
200 …rhox = appctx->kappa * ((u[i + 1].c - u[i].c) * u[i + 1].rho - (u[i].c - u[i - 1].c) * u[i].rho) *… in IFunction()
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex29.c6 -div \rho grad u = f, 0 < x,y < 1,
36 PetscReal rho; member
61 user.rho = 1.0; in main()
62 PetscCall(PetscOptionsReal("-rho", "The conductivity", "ex29.c", user.rho, &user.rho, NULL)); in main()
142 …e ComputeRho(PetscInt i, PetscInt j, PetscInt mx, PetscInt my, PetscReal centerRho, PetscReal *rho) in ComputeRho() argument
146 *rho = centerRho; in ComputeRho()
148 *rho = 1.0; in ComputeRho()
159 PetscReal Hx, Hy, HydHx, HxdHy, rho; in ComputeMatrix() local
166 centerRho = user->rho; in ComputeMatrix()
177 PetscCall(ComputeRho(i, j, mx, my, centerRho, &rho)); in ComputeMatrix()
[all …]
H A Dex13f90.F9068 ! Solves -div (rho grad psi) = F using finite differences.
69 ! rho is a 2-dimensional array of size m by n, stored in Fortran
72 subroutine UserDoLinearSolver(rho, userctx, userb, userx, ierr) argument
76 PetscScalar rho(*), userb(*), userx(*)
103 ! -div rho grad
104 ! on an m by n grid with zero Dirichlet boundary conditions. The rho
114 v = -0.5*(rho(II + 1) + rho(JJ + 1))*hy2
119 v = -0.5*(rho(II + 1) + rho(JJ + 1))*hy2
124 v = -0.5*(rho(II + 1) + rho(JJ + 1))*hx2
129 v = -0.5*(rho(II + 1) + rho(JJ + 1))*hx2
[all …]
H A Dex13.c33 PetscScalar *userx, *rho, *solution, *userb, hx, hy, x, y; in main() local
69 PetscCall(PetscMalloc1(N, &rho)); in main()
72 Fill up the array rho[] with the function rho(x,y) = x; fill the in main()
82 rho[Ii] = x; in main()
99 PetscCall(UserDoLinearSolver(rho, &userctx, userb, userx)); in main()
117 PetscCall(PetscFree(rho)); in main()
166 Solves -div (rho grad psi) = F using finite differences.
167 rho is a 2-dimensional array of size m by n, stored in Fortran
171 PetscErrorCode UserDoLinearSolver(PetscScalar *rho, UserCtx *userctx, PetscScalar *userb, PetscScal… in UserDoLinearSolver() argument
186 -div rho grad in UserDoLinearSolver()
[all …]
H A Dex35.cxx4 -div \rho grad u = f, 0 < x,y < 1,
43 …Run with different values of \rho and \nu (problem 1) to control diffusion and gaussian source spr…
45 mpiexec -n $NP ./ex35 -problem 2 -n 20 -nu 0.02 -rho 0.01
46 mpiexec -n $NP ./ex35 -problem 2 -n 40 -nu 0.01 -rho 0.005 -io -ksp_monitor
47 mpiexec -n $NP ./ex35 -problem 2 -n 80 -nu 0.01 -rho 0.005 -io -ksp_monitor -pc_type hypre
48 mpiexec -n $NP ./ex35 -problem 2 -n 160 -bc neumann -nu 0.005 -rho 0.01 -io
49 mpiexec -n $NP ./ex35 -problem 2 -n 320 -bc neumann -nu 0.001 -rho 1 -io
79 … Usage: ./ex35 -problem 2 -n 80 -nu 0.01 -rho 0.005 -io -ksp_monitor -pc_type gamg\n \
95 PetscReal rho; member
227 return user->rho; in ComputeDiffusionCoefficient()
[all …]
H A Dex36.cxx4 -div \rho grad u + \alpha u = f, 0 < x,y,z < 1,
16 and \rho = 1.0, \alpha = 10.0 uniformly in the domain.
49 …Run with different values of \rho and \nu (problem 1) to control diffusion and gaussian source spr…
51 mpiexec -n $NP ./ex36 -problem 2 -n 20 -nu 0.02 -rho 0.01
52 mpiexec -n $NP ./ex36 -problem 2 -n 40 -nu 0.01 -rho 0.005 -io -ksp_monitor
53 mpiexec -n $NP ./ex36 -problem 2 -n 80 -nu 0.01 -rho 0.005 -io -ksp_monitor -pc_type gamg
54 mpiexec -n $NP ./ex36 -problem 2 -n 160 -bc neumann -nu 0.005 -rho 0.01 -io
55 mpiexec -n $NP ./ex36 -problem 2 -n 320 -bc neumann -nu 0.001 -rho 1 -io
78 PetscReal rho; member
212 …&& (coords[1] < 2.0 / 3.0) && (coords[2] > 1.0 / 3.0) && (coords[2] < 2.0 / 3.0)) return user->rho; in ComputeDiffusionCoefficient()
[all …]
/petsc/src/ksp/ksp/tests/output/
H A Dex5_mat_type-lmvmsymbroyden.out22 Rescale params: alpha=1., beta=0.5, rho=1.
40 Rescale params: alpha=1., beta=0.5, rho=1.
58 Rescale params: alpha=1., beta=0.5, rho=1.
76 Rescale params: alpha=1., beta=0.5, rho=1.
94 Rescale params: alpha=1., beta=0.5, rho=1.
112 Rescale params: alpha=1., beta=0.5, rho=1.
130 Rescale params: alpha=1., beta=0.5, rho=1.
148 Rescale params: alpha=1., beta=0.5, rho=1.
166 Rescale params: alpha=1., beta=0.5, rho=1.
184 Rescale params: alpha=1., beta=0.5, rho=1.
[all …]
H A Dex5_mat_type-lmvmdfp.out21 Rescale params: alpha=1., beta=0.5, rho=1.
38 Rescale params: alpha=1., beta=0.5, rho=1.
55 Rescale params: alpha=1., beta=0.5, rho=1.
72 Rescale params: alpha=1., beta=0.5, rho=1.
89 Rescale params: alpha=1., beta=0.5, rho=1.
106 Rescale params: alpha=1., beta=0.5, rho=1.
123 Rescale params: alpha=1., beta=0.5, rho=1.
140 Rescale params: alpha=1., beta=0.5, rho=1.
157 Rescale params: alpha=1., beta=0.5, rho=1.
174 Rescale params: alpha=1., beta=0.5, rho=1.
[all …]
H A Dex5_mat_type-lmvmsymbadbroyden.out22 Rescale params: alpha=1., beta=0.5, rho=1.
40 Rescale params: alpha=1., beta=0.5, rho=1.
58 Rescale params: alpha=1., beta=0.5, rho=1.
76 Rescale params: alpha=1., beta=0.5, rho=1.
94 Rescale params: alpha=1., beta=0.5, rho=1.
112 Rescale params: alpha=1., beta=0.5, rho=1.
130 Rescale params: alpha=1., beta=0.5, rho=1.
148 Rescale params: alpha=1., beta=0.5, rho=1.
166 Rescale params: alpha=1., beta=0.5, rho=1.
184 Rescale params: alpha=1., beta=0.5, rho=1.
[all …]
H A Dex5_mat_type-lmvmbfgs.out21 Rescale params: alpha=1., beta=0.5, rho=1.
38 Rescale params: alpha=1., beta=0.5, rho=1.
55 Rescale params: alpha=1., beta=0.5, rho=1.
72 Rescale params: alpha=1., beta=0.5, rho=1.
89 Rescale params: alpha=1., beta=0.5, rho=1.
106 Rescale params: alpha=1., beta=0.5, rho=1.
123 Rescale params: alpha=1., beta=0.5, rho=1.
140 Rescale params: alpha=1., beta=0.5, rho=1.
157 Rescale params: alpha=1., beta=0.5, rho=1.
174 Rescale params: alpha=1., beta=0.5, rho=1.
[all …]
H A Dex5_mat_type-lmvmdiagbroyden.out21 Rescale params: alpha=1., beta=0.5, rho=1.
38 Rescale params: alpha=1., beta=0.5, rho=1.
55 Rescale params: alpha=1., beta=0.5, rho=1.
72 Rescale params: alpha=1., beta=0.5, rho=1.
89 Rescale params: alpha=1., beta=0.5, rho=1.
106 Rescale params: alpha=1., beta=0.5, rho=1.
123 Rescale params: alpha=1., beta=0.5, rho=1.
140 Rescale params: alpha=1., beta=0.5, rho=1.
157 Rescale params: alpha=1., beta=0.5, rho=1.
174 Rescale params: alpha=1., beta=0.5, rho=1.
[all …]
/petsc/src/ts/tutorials/
H A Dex9.c13 " rho_t + (rho*u)_x = 0\n"
14 " (rho*u)_t + (rho*u^2 + c^2*rho)_x = 0\n"
647 * z = sqrt(rho*bulk), c = sqrt(rho/bulk)
658 PetscReal c; /* speed of sound: c = sqrt(bulk/rho) */
659 PetscReal z; /* impedance: z = sqrt(rho*bulk) */
777 …PetscCall(PetscOptionsReal("-physics_acoustics_c", "c = sqrt(bulk/rho)", "", user->c, &user->c, NU… in PhysicsCreate_Acoustics()
778 …PetscCall(PetscOptionsReal("-physics_acoustics_z", "z = sqrt(bulk*rho)", "", user->z, &user->z, NU… in PhysicsCreate_Acoustics()
873 PetscScalar rho, u; in PhysicsRiemann_IsoGas_Exact() member
878 …PetscCheck(L.rho > 0 && R.rho > 0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Reconstructed densi… in PhysicsRiemann_IsoGas_Exact()
881 PetscScalar res, tmp, rho = 0.5 * (L.rho + R.rho); /* initial guess */ in PhysicsRiemann_IsoGas_Exact() local
[all …]
/petsc/src/dm/impls/swarm/tutorials/
H A Dex1f90.F9015 Vec :: f, rho, rhs
64 PetscCallA(DMCreateGlobalVector(dm, rho, ierr))
66 PetscCallA(MatMultTranspose(M_p, f, rho, ierr))
70 PetscCallA(PetscObjectViewFromOptions(PetscObjectCast(rho), PETSC_NULL_OBJECT, '-rho_view', ierr))
77 PetscCallA(MatMult(M, rho, rhs, ierr)) ! this is what you would do for and FE solve
79 PetscCallA(VecCopy(rho, rhs, ierr)) ! Identity: M^1 M rho
98 PetscCallA(VecDestroy(rho, ierr))
/petsc/src/ksp/ksp/impls/bcgs/fbcgs/
H A Dfbcgs.c18 PetscScalar rho, rhoold, alpha, beta, omega, omegaold, d1; in KSPSolve_FBCGS() local
77 PetscCall(VecDot(R, RP, &rho)); /* rho <- (r,rp) */ in KSPSolve_FBCGS()
78 beta = (rho / rhoold) * (alpha / omegaold); in KSPSolve_FBCGS()
91 alpha = rho / d1; /* alpha <- rho / (v,rp) */ in KSPSolve_FBCGS()
123 rhoold = rho; in KSPSolve_FBCGS()
134 if (rho == 0.0) { in KSPSolve_FBCGS()
135 …mm((PetscObject)ksp), PETSC_ERR_NOT_CONVERGED, "KSPSolve breakdown due to zero rho inner product"); in KSPSolve_FBCGS()
137 PetscCall(PetscInfo(ksp, "Breakdown due to zero rho inner product\n")); in KSPSolve_FBCGS()
/petsc/src/ksp/ksp/impls/cgs/
H A Dcgs.c18 PetscScalar rho, rhoold, a, s, b; in KSPSolve_CGS() local
90 a = rhoold / s; /* a <- rho / s */ in KSPSolve_CGS()
96 PetscCall(VecDot(R, RP, &rho)); /* rho <- (r,rp) */ in KSPSolve_CGS()
97 KSPCheckDot(ksp, rho); in KSPSolve_CGS()
99 dp = PetscAbsScalar(rho); in KSPSolve_CGS()
114 b = rho / rhoold; /* b <- rho / rhoold */ in KSPSolve_CGS()
119 rhoold = rho; in KSPSolve_CGS()
/petsc/src/ksp/ksp/impls/bcgs/pipebcgs/
H A Dpipebcgs.c18 PetscScalar rho, rhoold, alpha, beta, omega = 0.0, d1, d2, d3; in KSPSolve_PIPEBCGS() local
75 PetscCall(VecDotBegin(R, RP, &rho)); /* rho <- (r,rp) */ in KSPSolve_PIPEBCGS()
79 PetscCall(VecDotEnd(R, RP, &rho)); in KSPSolve_PIPEBCGS()
87 alpha = rho / d2; in KSPSolve_PIPEBCGS()
142 rhoold = rho; in KSPSolve_PIPEBCGS()
145 …PetscCall(VecDotBegin(R, RP, &rho)); … in KSPSolve_PIPEBCGS()
155 PetscCall(VecDotEnd(R, RP, &rho)); in KSPSolve_PIPEBCGS()
162 beta = (rho / rhoold) * (alpha / omega); in KSPSolve_PIPEBCGS()
163 …alpha = rho / (d2 + beta * d1 - beta * omega * d3); /* alpha <- rho / (d2 + beta * d1 - beta * ome… in KSPSolve_PIPEBCGS()
189 if (rho == 0.0) { in KSPSolve_PIPEBCGS()
/petsc/src/ksp/ksp/impls/gmres/agmres/
H A Dagmresorthog.c126 PetscReal c, s, rho, Ajj, val, tt, old; in KSPAGMRESRoddec() local
152 PetscCallBLAS("BLASnrm2", rho = -PetscSign(Ajj) * BLASnrm2_(&len, &Qloc[j * nloc + j], &pas)); in KSPAGMRESRoddec()
153 if (rho == 0.0) tloc[j] = 0.0; in KSPAGMRESRoddec()
155 tloc[j] = (Ajj - rho) / rho; in KSPAGMRESRoddec()
157 val = 1.0 / (Ajj - rho); in KSPAGMRESRoddec()
165 Qloc[j * nloc + j] = rho; in KSPAGMRESRoddec()
179 PetscCall(KSPAGMRESRoddecGivens(&c, &s, &rho, 1)); in KSPAGMRESRoddec()
186 Qloc[d * nloc] = rho; in KSPAGMRESRoddec()
194 PetscCall(KSPAGMRESRoddecGivens(&c, &s, &rho, 1)); in KSPAGMRESRoddec()
200 Qloc[j * nloc + i] = rho; in KSPAGMRESRoddec()
[all …]
/petsc/src/snes/impls/tr/
H A Dtrimpl.h19 Given rho = (fk - fkp1) / (m(0) - m(pk))
22 rho < eta2 -> delta *= t1
23 rho > eta3 -> delta *= t2
26 The step is accepted if rho > eta1
/petsc/src/ksp/ksp/impls/bcgs/
H A Dbcgs.c21 PetscScalar rho, rhoold, alpha, beta, omega, omegaold, d1; in KSPSolve_BCGS() local
74 PetscCall(VecDot(R, RP, &rho)); /* rho <- (r,rp) */ in KSPSolve_BCGS()
75 beta = (rho / rhoold) * (alpha / omegaold); in KSPSolve_BCGS()
86 alpha = rho / d1; /* a <- rho / (v,rp) */ in KSPSolve_BCGS()
118 rhoold = rho; in KSPSolve_BCGS()
129 if (rho == 0.0) { in KSPSolve_BCGS()
132 PetscCall(PetscInfo(ksp, "Breakdown due to zero rho inner product\n")); in KSPSolve_BCGS()
/petsc/src/ksp/ksp/impls/tfqmr/
H A Dtfqmr.c14 PetscScalar rho, rhoold, a, s, b, eta, etaold, psiold, cf; in KSPSolve_TFQMR() local
70 a = rhoold / s; /* a <- rho / s */ in KSPSolve_TFQMR()
107 PetscCall(VecDot(R, RP, &rho)); /* rho <- (r,rp) */ in KSPSolve_TFQMR()
108 b = rho / rhoold; /* b <- rho / rhoold */ in KSPSolve_TFQMR()
114 rhoold = rho; in KSPSolve_TFQMR()
/petsc/src/sys/tutorials/output/
H A Dex5_1.out3 rho = 44. ; Density, kg/m^3
16 rho = 44. ; Density, kg/m^3
27 The value of rho after loading is: 44.000000
H A Dex5_yaml.out3 rho = 1. ; Density, kg/m^3
16 rho = 1. ; Density, kg/m^3
27 The value of rho after loading is: 1.000000
/petsc/src/sys/tutorials/
H A Dex5.c35 PetscReal rho; member
75 PetscCall(PetscBagRegisterReal(bag, &params->rho, 3.0, "rho", "Density, kg/m^3")); in main()
117 …tscCall(PetscPrintf(PETSC_COMM_WORLD, "The value of rho after loading is: %f\n", (double)params->r… in main()
/petsc/src/ksp/ksp/utils/lmvm/rescale/
H A Dsymbrdnrescale.c60 sigma = ldb->rho * signew + (1.0 - ldb->rho) * sigma; in SymBroydenRescaleUpdateScalar()
301 if (ldb->rho == 1.0) PetscCall(VecCopy(ldb->invDnew, invD)); in SymBroydenRescaleUpdateDiagonal()
302 else if (ldb->rho) PetscCall(VecAXPBY(invD, 1.0 - ldb->rho, ldb->rho, ldb->invDnew)); in SymBroydenRescaleUpdateDiagonal()
377 mctx->rho = bctx->rho; in SymBroydenRescaleCopy()
417 …l("-mat_lmvm_rho", "(developer) update limiter in the J0 scaling", "", ldb->rho, &ldb->rho, NULL)); in SymBroydenRescaleSetFromOptions()
426 …PetscCheck(!(ldb->rho < 0.0) && !(ldb->rho > 1.0), PetscObjectComm((PetscObject)B), PETSC_ERR_ARG_… in SymBroydenRescaleSetFromOptions()
499 …tf(pv, "Rescale params: alpha=%g, beta=%g, rho=%g\n", (double)ldb->alpha, (double)ldb->beta, (doub… in SymBroydenRescaleView()
542 (*ldb)->rho = 1.0; in SymBroydenRescaleCreate()
/petsc/src/dm/impls/swarm/tests/
H A Dex6.c269 Vec phi, locPhi, rho, f; in ComputeFieldAtParticles_Primal() local
290 /* Create the charges rho */ in ComputeFieldAtParticles_Primal()
291 PetscCall(DMGetGlobalVector(dm, &rho)); in ComputeFieldAtParticles_Primal()
292 PetscCall(PetscObjectSetName((PetscObject)rho, "rho")); in ComputeFieldAtParticles_Primal()
295 PetscCall(MatMultTranspose(M_p, f, rho)); in ComputeFieldAtParticles_Primal()
305 /* Remove constant from rho */ in ComputeFieldAtParticles_Primal()
306 PetscCall(VecGetSize(rho, &n)); in ComputeFieldAtParticles_Primal()
307 PetscCall(VecSum(rho, &sum)); in ComputeFieldAtParticles_Primal()
308 PetscCall(VecShift(rho, -sum / n)); in ComputeFieldAtParticles_Primal()
309 PetscCall(VecSum(rho, &sum)); in ComputeFieldAtParticles_Primal()
[all …]

12345