| /petsc/src/snes/tutorials/ |
| H A D | ex19.c | 239 PetscScalar u, uxx, uyy, vx, vy, avx, avy, vxp, vxm, vyp, vym; in FormFunctionLocal() local 334 uxx = (2.0 * u - x[j][i - 1].u - x[j][i + 1].u) * hydhx; in FormFunctionLocal() 336 f[j][i].u = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx; in FormFunctionLocal() 340 uxx = (2.0 * u - x[j][i - 1].v - x[j][i + 1].v) * hydhx; in FormFunctionLocal() 342 f[j][i].v = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy; in FormFunctionLocal() 346 uxx = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx; in FormFunctionLocal() 348 …f[j][i].omega = uxx + uyy + (vxp * (u - x[j][i - 1].omega) + vxm * (x[j][i + 1].omega - u)) * hy +… in FormFunctionLocal() 352 uxx = (2.0 * u - x[j][i - 1].temp - x[j][i + 1].temp) * hydhx; in FormFunctionLocal() 354 …f[j][i].temp = uxx + uyy + prandtl * ((vxp * (u - x[j][i - 1].temp) + vxm * (x[j][i + 1].temp - u)… in FormFunctionLocal() 380 PetscScalar u, uxx, uyy, vx, vy, avx, avy, vxp, vxm, vyp, vym; in NonlinearGS() local [all …]
|
| H A D | ex5m.m | 17 uxx = (2.0*u - localX(1:m-2,2:n-1) - localX(3:m,2:n-1))*hydhx; variable 22 localF(2:m-1,2:n-1) = uxx + uyy - sc*exp(u);
|
| H A D | ex35.c | 249 PetscScalar **x, F, J, u, uxx, uyy; in NonlinearGS() local 293 uxx = (2.0 * u - x[j][i - 1] - x[j][i + 1]) * hydhx; in NonlinearGS() 295 F = uxx + uyy; in NonlinearGS()
|
| H A D | ex46.c | 124 PetscScalar u, ux, uy, uxx, uyy; in FormFunctionLocal() local 148 uxx = (2.0 * u - x[j][i - 1] - x[j][i + 1]) * hydhx; in FormFunctionLocal() 150 …f[j][i] = D * (uxx + uyy) - (K * funcA(x[j][i], user) * PetscSqrtScalar(ux * ux + uy * uy) + funcU… in FormFunctionLocal()
|
| H A D | ex5f.F90 | 171 PetscScalar u, uxx, uyy 197 uxx = hydhx*(two*u - x(i - 1, j) - x(i + 1, j)) 199 f(i, j) = uxx + uyy - sc*exp(u)
|
| H A D | ex5f90t.F90 | 230 PetscScalar u, uxx, uyy 249 uxx = hydhx*(two*u - x(i - 1, j) - x(i + 1, j)) 251 f(i, j) = uxx + uyy - sc*exp(u)
|
| H A D | ex5f90.F90 | 240 PetscScalar u, uxx, uyy 259 uxx = hydhx*(two*u - x(i - 1, j) - x(i + 1, j)) 261 f(i, j) = uxx + uyy - sc*exp(u)
|
| /petsc/src/ts/tutorials/ |
| H A D | ex26.c | 263 PetscScalar u, udot, uxx, uyy, vx, vy, avx, avy, vxp, vxm, vyp, vym; in FormIFunctionLocal() local 359 uxx = (2.0 * u - x[j][i - 1].u - x[j][i + 1].u) * hydhx; in FormIFunctionLocal() 361 f[j][i].u = udot + uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx; in FormIFunctionLocal() 366 uxx = (2.0 * u - x[j][i - 1].v - x[j][i + 1].v) * hydhx; in FormIFunctionLocal() 368 f[j][i].v = udot + uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy; in FormIFunctionLocal() 372 uxx = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx; in FormIFunctionLocal() 374 …f[j][i].omega = (xdot[j][i].omega + uxx + uyy + (vxp * (u - x[j][i - 1].omega) + vxm * (x[j][i + 1… in FormIFunctionLocal() 378 uxx = (2.0 * u - x[j][i - 1].temp - x[j][i + 1].temp) * hydhx; in FormIFunctionLocal() 380 …f[j][i].temp = (xdot[j][i].temp + uxx + uyy + prandtl * ((vxp * (u - x[j][i - 1].temp) + vxm * (x[… in FormIFunctionLocal()
|
| H A D | ex7.c | 129 PetscScalar u, uxx, uyy, **x, **f; in FormFunction() local 172 uxx = (two * u - x[j][i - 1] - x[j][i + 1]) * sx; in FormFunction() 175 …f[j][i] = -u * (uxx + uyy) - (4.0 - 1.0) * ((x[j][i + 1] - x[j][i - 1]) * (x[j][i + 1] - x[j][i - … in FormFunction()
|
| H A D | ex15.c | 161 PetscScalar u, uxx, uyy, **uarray, **f, **udot; in FormIFunction() local 220 uxx = (-2.0 * u + uarray[j][i - 1] + uarray[j][i + 1]); in FormIFunction() 224 …uxx = 2.0 * uxx / 3.0 + (0.5 * (uarray[j - 1][i - 1] + uarray[j - 1][i + 1] + uarray[j + 1][i - 1]… in FormIFunction() 227 f[j][i] = udot[j][i] - (uxx * sx + uyy * sy); in FormIFunction()
|
| H A D | ex12.c | 136 PetscScalar u, uxx, uyy, v, ***x, ***f; in FormFunction() local 182 uxx = (-2.0 * u + x[j][i - 1][0] + x[j][i + 1][0]) * sx; in FormFunction() 185 f[j][i][1] = uxx + uyy; in FormFunction()
|
| H A D | ex1.c | 218 PetscScalar ut, ub, ul, ur, u, uxx, uyy, sc, *f; in FormFunction() local 246 uxx = (-ur + two * u - ul) * hydhx; in FormFunction() 248 f[row] = -uxx + -uyy + sc * lambda * PetscExpScalar(u); in FormFunction()
|
| H A D | ex1f.F90 | 38 PetscScalar uxx, uyy, sc 67 uxx = (-ur + two*u - ul)*hydhx 69 ff(row) = -uxx - uyy + sc*lambda*exp(u)
|
| H A D | ex13.c | 122 PetscScalar u, uxx, uyy, **uarray, **f; in RHSFunction() local 159 uxx = (-two * u + uarray[j][i - 1] + uarray[j][i + 1]) * sx; in RHSFunction() 161 f[j][i] = uxx + uyy; in RHSFunction()
|
| H A D | ex27.c | 191 PetscScalar u, uxx, uyy; in FormIFunctionLocal() local 251 if (i != info->mx - 1) uxx = (2.0 * u - x[j][i - 1].sp[l] - x[j][i + 1].sp[l]) * hydhx; in FormIFunctionLocal() 252 else uxx = u - x[j][i - 1].sp[l]; in FormIFunctionLocal() 258 f[j][i].sp[l] += f_advect + ctx->dispersivity * (uxx + uyy); in FormIFunctionLocal()
|
| /petsc/src/snes/tests/ |
| H A D | ex69.c | 214 PetscScalar u, uxx, uyy, vx, vy, avx, avy, vxp, vxm, vyp, vym; in FormFunctionLocal() local 315 uxx = (2.0 * u - x[j][i - 1].u - x[j][i + 1].u) * hydhx; in FormFunctionLocal() 317 f[j][i].u = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx; in FormFunctionLocal() 321 uxx = (2.0 * u - x[j][i - 1].v - x[j][i + 1].v) * hydhx; in FormFunctionLocal() 323 f[j][i].v = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy; in FormFunctionLocal() 327 uxx = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx; in FormFunctionLocal() 329 …f[j][i].omega = uxx + uyy + (vxp * (u - x[j][i - 1].omega) + vxm * (x[j][i + 1].omega - u)) * hy +… in FormFunctionLocal() 333 uxx = (2.0 * u - x[j][i - 1].temp - x[j][i + 1].temp) * hydhx; in FormFunctionLocal() 335 …f[j][i].temp = uxx + uyy + prandtl * ((vxp * (u - x[j][i - 1].temp) + vxm * (x[j][i + 1].temp - u)… in FormFunctionLocal()
|
| H A D | ex1f.F90 | 255 PetscScalar u, uxx, uyy 275 uxx = hydhx*(two*u - x(i - 1, j) - x(i + 1, j)) 277 f(i, j) = uxx + uyy - sc*exp(u)
|
| /petsc/src/ts/tutorials/advection-diffusion-reaction/ |
| H A D | ex5adj_mf.c | 65 PetscScalar uc, uxx, uyy, vc, vxx, vyy, ucb, vcb; in MyRHSMatMultTranspose() local 100 uxx = (-2.0 * ucb + x[j][i - 1].u + x[j][i + 1].u) * sx; in MyRHSMatMultTranspose() 106 y[j][i].u = appctx->D1 * (uxx + uyy) - ucb * (vc * vc + appctx->gamma) + vc * vc * vcb; in MyRHSMatMultTranspose() 124 PetscScalar uc, uxx, uyy, vc, vxx, vyy, ucb, vcb; in MyIMatMultTranspose() local 159 uxx = (-2.0 * ucb + x[j][i - 1].u + x[j][i + 1].u) * sx; in MyIMatMultTranspose() 165 y[j][i].u = appctx->D1 * (uxx + uyy) - ucb * (vc * vc + appctx->gamma) + vc * vc * vcb; in MyIMatMultTranspose() 185 PetscScalar uc, uxx, uyy, vc, vxx, vyy, ucb, vcb; in MyIMatMult() local 220 uxx = (-2.0 * ucb + x[j][i - 1].u + x[j][i + 1].u) * sx; in MyIMatMult() 226 y[j][i].u = appctx->D1 * (uxx + uyy) - (vc * vc + appctx->gamma) * ucb - 2.0 * uc * vc * vcb; in MyIMatMult()
|
| H A D | reaction_diffusion.c | 33 PetscScalar uc, uxx, uyy, vc, vxx, vyy; in RHSFunction() local 72 uxx = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx; in RHSFunction() 77 f[j][i].u = appctx->D1 * (uxx + uyy) - uc * vc * vc + appctx->gamma * (1.0 - uc); in RHSFunction() 241 PetscScalar uc, uxx, uyy, vc, vxx, vyy; in IFunction() local 281 uxx = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx; in IFunction() 286 … f[j][i].u = udot[j][i].u - (appctx->D1 * (uxx + uyy) - uc * vc * vc + appctx->gamma * (1.0 - uc)); in IFunction()
|
| /petsc/src/mat/impls/h2opus/cuda/ |
| H A D | math2opus.cu | 342 PetscScalar *xx, *yy, *uxx, *uyy; in MatMultNKernel_H2OPUS() local 371 uxx = MatH2OpusGetThrustPointer(*h2opus->xx); in MatMultNKernel_H2OPUS() 373 PetscCall(PetscSFBcastBegin(bsf, MPIU_SCALAR, xx, uxx, MPI_REPLACE)); in MatMultNKernel_H2OPUS() 374 PetscCall(PetscSFBcastEnd(bsf, MPIU_SCALAR, xx, uxx, MPI_REPLACE)); in MatMultNKernel_H2OPUS() 376 uxx = xx; in MatMultNKernel_H2OPUS() 383 …ansA ? H2Opus_Trans : H2Opus_NoTrans, */ h2opus->s, *h2opus->dist_hmatrix, uxx, blda, 0.0, uyy, cl… in MatMultNKernel_H2OPUS() 387 …hgemv(transA ? H2Opus_Trans : H2Opus_NoTrans, h2opus->s, *h2opus->hmatrix, uxx, blda, 0.0, uyy, cl… in MatMultNKernel_H2OPUS() 410 uxx = MatH2OpusGetThrustPointer(*h2opus->xx_gpu); in MatMultNKernel_H2OPUS() 412 PetscCall(PetscSFBcastBegin(bsf, MPIU_SCALAR, xx, uxx, MPI_REPLACE)); in MatMultNKernel_H2OPUS() 413 PetscCall(PetscSFBcastEnd(bsf, MPIU_SCALAR, xx, uxx, MPI_REPLACE)); in MatMultNKernel_H2OPUS() [all …]
|
| /petsc/src/ts/tests/ |
| H A D | ex21.c | 38 PetscScalar ut, u, ue, uw, un, us, uxx, uyy; in FormIFunctionLocal() local 62 uxx = (uw - 2.0 * u + ue) / (hx * hx); in FormIFunctionLocal() 64 f[j][i] = ut - uxx - uyy - lambda * PetscExpScalar(u); in FormIFunctionLocal()
|
| /petsc/src/ts/tutorials/autodiff/ |
| H A D | adr_ex5adj_mf.cxx | 253 PetscScalar uc, uxx, uyy, vc, vxx, vyy; in IFunctionLocalPassive() local 271 uxx = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx; in IFunctionLocalPassive() 276 … f[j][i].u = udot[j][i].u - appctx->D1 * (uxx + uyy) + uc * vc * vc - appctx->gamma * (1.0 - uc); in IFunctionLocalPassive() 293 adouble uc, uxx, uyy, vc, vxx, vyy; in IFunctionActive() local 367 uxx = (-2.0 * uc + u_a[j][i - 1].u + u_a[j][i + 1].u) * sx; in IFunctionActive() 372 … f_a[j][i].u = udot[j][i].u - appctx->D1 * (uxx + uyy) + uc * vc * vc - appctx->gamma * (1.0 - uc); in IFunctionActive()
|
| H A D | adr_ex5adj.cxx | 365 PetscScalar uc, uxx, uyy, vc, vxx, vyy; in IFunctionLocalPassive() local 383 uxx = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx; in IFunctionLocalPassive() 388 … f[j][i].u = udot[j][i].u - appctx->D1 * (uxx + uyy) + uc * vc * vc - appctx->gamma * (1.0 - uc); in IFunctionLocalPassive() 405 adouble uc, uxx, uyy, vc, vxx, vyy; in IFunctionActive() local 479 uxx = (-2.0 * uc + u_a[j][i - 1].u + u_a[j][i + 1].u) * sx; in IFunctionActive() 484 … f_a[j][i].u = udot[j][i].u - appctx->D1 * (uxx + uyy) + uc * vc * vc - appctx->gamma * (1.0 - uc); in IFunctionActive() 532 PetscScalar uc, uxx, uyy, vc, vxx, vyy; in RHSFunctionPassive() local 575 uxx = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx; in RHSFunctionPassive() 580 f[j][i].u = appctx->D1 * (uxx + uyy) - uc * vc * vc + appctx->gamma * (1.0 - uc); in RHSFunctionPassive() 613 adouble uc, uxx, uyy, vc, vxx, vyy; in RHSFunctionActive() local [all …]
|
| /petsc/src/binding/petsc4py/demo/legacy/wrap-f2py/ |
| H A D | Bratu2D.F90 | 97 PetscScalar u, uxx, uyy 115 uxx = (two*u - x(i - 1, j) - x(i + 1, j))*hydhx 117 f(i, j) = uxx + uyy - lambda*exp(u)*sc
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex14f.F90 | 153 PetscScalar u, uxx, uyy 198 uxx = (two*u - xx(row - 1) - xx(row + 1))*hydhx 200 ff(rowf) = uxx + uyy - sc*exp(u)
|