Home
last modified time | relevance | path

Searched refs:rho (Results 1 – 25 of 105) sorted by relevance

12345

/petsc/src/ts/tutorials/advection-diffusion-reaction/
H A Dex4.c29 PetscScalar rho, c; member
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()
203 …rho = udot[i].rho - appctx->epsilon * rhoxx + kcxrhox - appctx->mu * PetscAbsScalar(rho) * (1.0 - … in IFunction()
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex13.c33 PetscScalar *userx, *rho, *solution, *userb, hx, hy, x, y; in main() local
69 PetscCall(PetscMalloc1(N, &rho)); in main()
82 rho[Ii] = x; in main()
99 PetscCall(UserDoLinearSolver(rho, &userctx, userb, userx)); in main()
117 PetscCall(PetscFree(rho)); in main()
171 PetscErrorCode UserDoLinearSolver(PetscScalar *rho, UserCtx *userctx, PetscScalar *userb, PetscScal… in UserDoLinearSolver() argument
197 v = -.5 * (rho[Ii] + rho[J]) * hy2; in UserDoLinearSolver()
202 v = -.5 * (rho[Ii] + rho[J]) * hy2; in UserDoLinearSolver()
207 v = -.5 * (rho[Ii] + rho[J]) * hx2; in UserDoLinearSolver()
212 v = -.5 * (rho[Ii] + rho[J]) * hx2; in UserDoLinearSolver()
[all …]
H A Dex29.c36 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()
180 v[0] = 2.0 * rho * (HxdHy + HydHx); in ComputeMatrix()
[all …]
H A Dex13f90.F9072 subroutine UserDoLinearSolver(rho, userctx, userb, userx, ierr) argument
76 PetscScalar rho(*), userb(*), userx(*)
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
132 v = 2*rho(II + 1)*(hx2 + hy2)
247 PetscScalar, allocatable :: rho(:, :)
285 allocate (rho(m, n))
296 rho(i, j) = x
[all …]
H A Dex73.c39 PetscReal rho; member
55 user->rho = 1.0; in UserContextCreate()
56 PetscCall(PetscOptionsReal("-rho", "The conductivity", "ex29.c", user->rho, &user->rho, NULL)); in UserContextCreate()
974 …e ComputeRho(PetscInt i, PetscInt j, PetscInt mx, PetscInt my, PetscReal centerRho, PetscReal *rho) in ComputeRho() argument
978 *rho = centerRho; in ComputeRho()
980 *rho = 1.0; in ComputeRho()
991 PetscReal Hx, Hy, HydHx, HxdHy, rho; in ComputeMatrix_DMDA() local
999 centerRho = user->rho; in ComputeMatrix_DMDA()
1010 PetscCall(ComputeRho(i, j, mx, my, centerRho, &rho)); in ComputeMatrix_DMDA()
1013 v[0] = 2.0 * rho * (HxdHy + HydHx); in ComputeMatrix_DMDA()
[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-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 …]
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 …]
/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/ts/tutorials/
H A Dex9.c873 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
884 fl = (L.rho < rho) ? (rho - L.rho) / PetscSqrtScalar(L.rho * rho) /* shock */ in PhysicsRiemann_IsoGas_Exact()
885 : PetscLogScalar(rho) - PetscLogScalar(L.rho); /* rarefaction */ in PhysicsRiemann_IsoGas_Exact()
886 fr = (R.rho < rho) ? (rho - R.rho) / PetscSqrtScalar(R.rho * rho) /* shock */ in PhysicsRiemann_IsoGas_Exact()
887 : PetscLogScalar(rho) - PetscLogScalar(R.rho); /* rarefaction */ in PhysicsRiemann_IsoGas_Exact()
891 star.rho = rho; in PhysicsRiemann_IsoGas_Exact()
895 …dfl = (L.rho < rho) ? 1 / PetscSqrtScalar(L.rho * rho) * (1 - 0.5 * (rho - L.rho) / rho) : 1 / rho; in PhysicsRiemann_IsoGas_Exact()
896 …dfr = (R.rho < rho) ? 1 / PetscSqrtScalar(R.rho * rho) * (1 - 0.5 * (rho - R.rho) / rho) : 1 / rho; in PhysicsRiemann_IsoGas_Exact()
[all …]
/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/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/cgs/
H A Dcgs.c18 PetscScalar rho, rhoold, a, s, b; in KSPSolve_CGS() local
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/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()
/petsc/src/ksp/ksp/impls/bcgs/fbcgsr/
H A Dfbcgsr.c21 PetscReal rho; in KSPSolve_FBCGSR() local
83 PetscCall(VecNorm(R, NORM_2, &rho)); in KSPSolve_FBCGSR()
86 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = rho; in KSPSolve_FBCGSR()
106 tau = rho * rho; in KSPSolve_FBCGSR()
172 rho = PetscSqrtReal(PetscAbsScalar(xi1 - omega * xi2)); /* residual norm */ in KSPSolve_FBCGSR()
180 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = rho; in KSPSolve_FBCGSR()
/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/tao/complementarity/impls/ssls/
H A Dssils.c39 PetscReal delta, rho; in TaoSolve_SSILS() local
46 rho = ssls->rho; in TaoSolve_SSILS()
81 if (innerd <= delta * PetscPowReal(normd, rho)) { in TaoSolve_SSILS()
122 ssls->rho = 2.1; in TaoCreate_SSILS()
H A Dssfls.c25 PetscReal delta, rho; in TaoSolve_SSFLS() local
32 rho = ssls->rho; in TaoSolve_SSFLS()
73 if (innerd >= -delta * PetscPowReal(normd, rho)) { in TaoSolve_SSFLS()
132 ssls->rho = 2.1; in TaoCreate_SSFLS()
/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()
/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 …all(PetscPrintf(PETSC_COMM_WORLD, "The value of rho after loading is: %f\n", (double)params->rho)); 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 …le params: alpha=%g, beta=%g, rho=%g\n", (double)ldb->alpha, (double)ldb->beta, (double)ldb->rho)); in SymBroydenRescaleView()
542 (*ldb)->rho = 1.0; in SymBroydenRescaleCreate()
/petsc/src/tao/linesearch/impls/gpcglinesearch/
H A Dgpcglinesearch.c32 PetscReal d1, finit, actred, prered, rho, gdx; in TaoLineSearchApply_GPCG() local
77 PetscCall(VecStepBoundInfo(x, s, ls->lower, ls->upper, &rho, &actred, &d1)); in TaoLineSearchApply_GPCG()
80 rho = 0; in TaoLineSearchApply_GPCG()
127 rho = actred / prered; in TaoLineSearchApply_GPCG()
137 } else if (rho > ls->ftol) { in TaoLineSearchApply_GPCG()

12345