| /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 335 uyy = (2.0 * u - x[j - 1][i].u - x[j + 1][i].u) * hxdhy; in FormFunctionLocal() 336 f[j][i].u = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx; in FormFunctionLocal() 341 uyy = (2.0 * u - x[j - 1][i].v - x[j + 1][i].v) * hxdhy; in FormFunctionLocal() 342 f[j][i].v = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy; in FormFunctionLocal() 347 uyy = (2.0 * u - x[j - 1][i].omega - x[j + 1][i].omega) * hxdhy; 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() 353 uyy = (2.0 * u - x[j - 1][i].temp - x[j + 1][i].temp) * hxdhy; 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 | 18 uyy = (2.0*u - localX(2:m-1,1:n-2) - localX(2:m-1,3:n))*hxdhy; 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 294 uyy = (2.0 * u - x[j - 1][i] - x[j + 1][i]) * hxdhy; in NonlinearGS() 295 F = uxx + uyy; in NonlinearGS()
|
| H A D | ex46.c | 124 PetscScalar u, ux, uy, uxx, uyy; in FormFunctionLocal() local 149 uyy = (2.0 * u - x[j - 1][i] - x[j + 1][i]) * hxdhy; 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 198 uyy = hxdhy*(two*u - x(i, j - 1) - x(i, j + 1)) 199 f(i, j) = uxx + uyy - sc*exp(u)
|
| H A D | ex5f90t.F90 | 230 PetscScalar u, uxx, uyy 250 uyy = hxdhy*(two*u - x(i, j - 1) - x(i, j + 1)) 251 f(i, j) = uxx + uyy - sc*exp(u)
|
| H A D | ex5f90.F90 | 240 PetscScalar u, uxx, uyy 260 uyy = hxdhy*(two*u - x(i, j - 1) - x(i, j + 1)) 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 360 uyy = (2.0 * u - x[j - 1][i].u - x[j + 1][i].u) * hxdhy; in FormIFunctionLocal() 361 f[j][i].u = udot + uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx; in FormIFunctionLocal() 367 uyy = (2.0 * u - x[j - 1][i].v - x[j + 1][i].v) * hxdhy; in FormIFunctionLocal() 368 f[j][i].v = udot + uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy; in FormIFunctionLocal() 373 uyy = (2.0 * u - x[j - 1][i].omega - x[j + 1][i].omega) * hxdhy; 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() 379 uyy = (2.0 * u - x[j - 1][i].temp - x[j + 1][i].temp) * hxdhy; 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 173 uyy = (two * u - x[j - 1][i] - x[j + 1][i]) * sy; 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 221 uyy = (-2.0 * u + uarray[j - 1][i] + uarray[j + 1][i]); in FormIFunction() 225 …uyy = 2.0 * uyy / 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 183 uyy = (-2.0 * u + x[j - 1][i][0] + x[j + 1][i][0]) * sy; in FormFunction() 185 f[j][i][1] = uxx + uyy; in FormFunction()
|
| H A D | ex27.c | 191 PetscScalar u, uxx, uyy; in FormIFunctionLocal() local 247 if (j == 0) uyy = u - x[j + 1][i].sp[l]; in FormIFunctionLocal() 248 else if (j == info->my - 1) uyy = u - x[j - 1][i].sp[l]; in FormIFunctionLocal() 249 else uyy = (2.0 * u - x[j - 1][i].sp[l] - x[j + 1][i].sp[l]) * hxdhy; in FormIFunctionLocal() 258 f[j][i].sp[l] += f_advect + ctx->dispersivity * (uxx + uyy); in FormIFunctionLocal()
|
| H A D | ex1.c | 218 PetscScalar ut, ub, ul, ur, u, uxx, uyy, sc, *f; in FormFunction() local 247 uyy = (-ut + two * u - ub) * hxdhy; in FormFunction() 248 f[row] = -uxx + -uyy + sc * lambda * PetscExpScalar(u); in FormFunction()
|
| H A D | ex1f.F90 | 38 PetscScalar uxx, uyy, sc 68 uyy = (-ut + two*u - ub)*hxdhy 69 ff(row) = -uxx - uyy + sc*lambda*exp(u)
|
| H A D | ex13.c | 122 PetscScalar u, uxx, uyy, **uarray, **f; in RHSFunction() local 160 uyy = (-two * u + uarray[j - 1][i] + uarray[j + 1][i]) * sy; in RHSFunction() 161 f[j][i] = uxx + uyy; in RHSFunction()
|
| /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 316 uyy = (2.0 * u - x[j - 1][i].u - x[j + 1][i].u) * hxdhy; in FormFunctionLocal() 317 f[j][i].u = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx; in FormFunctionLocal() 322 uyy = (2.0 * u - x[j - 1][i].v - x[j + 1][i].v) * hxdhy; in FormFunctionLocal() 323 f[j][i].v = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy; in FormFunctionLocal() 328 uyy = (2.0 * u - x[j - 1][i].omega - x[j + 1][i].omega) * hxdhy; 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() 334 uyy = (2.0 * u - x[j - 1][i].temp - x[j + 1][i].temp) * hxdhy; 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 276 uyy = hxdhy*(two*u - x(i, j - 1) - x(i, j + 1)) 277 f(i, j) = uxx + uyy - sc*exp(u)
|
| /petsc/src/mat/impls/h2opus/cuda/ |
| H A D | math2opus.cu | 342 PetscScalar *xx, *yy, *uxx, *uyy; in MatMultNKernel_H2OPUS() local 372 uyy = MatH2OpusGetThrustPointer(*h2opus->yy); in MatMultNKernel_H2OPUS() 377 uyy = yy; in MatMultNKernel_H2OPUS() 383 …ans : H2Opus_NoTrans, */ h2opus->s, *h2opus->dist_hmatrix, uxx, blda, 0.0, uyy, clda, B->cmap->N, … in MatMultNKernel_H2OPUS() 387 …2Opus_Trans : H2Opus_NoTrans, h2opus->s, *h2opus->hmatrix, uxx, blda, 0.0, uyy, clda, B->cmap->N, … in MatMultNKernel_H2OPUS() 391 PetscCall(PetscSFReduceBegin(csf, MPIU_SCALAR, uyy, yy, MPI_REPLACE)); in MatMultNKernel_H2OPUS() 392 PetscCall(PetscSFReduceEnd(csf, MPIU_SCALAR, uyy, yy, MPI_REPLACE)); in MatMultNKernel_H2OPUS() 411 uyy = MatH2OpusGetThrustPointer(*h2opus->yy_gpu); in MatMultNKernel_H2OPUS() 416 uyy = yy; in MatMultNKernel_H2OPUS() 423 …: H2Opus_NoTrans, */ h2opus->s, *h2opus->dist_hmatrix_gpu, uxx, blda, 0.0, uyy, clda, B->cmap->N, … in MatMultNKernel_H2OPUS() [all …]
|
| /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 101 uyy = (-2.0 * ucb + x[j - 1][i].u + x[j + 1][i].u) * sy; 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 160 uyy = (-2.0 * ucb + x[j - 1][i].u + x[j + 1][i].u) * sy; 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 221 uyy = (-2.0 * ucb + x[j - 1][i].u + x[j + 1][i].u) * sy; 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 73 uyy = (-2.0 * uc + u[j - 1][i].u + u[j + 1][i].u) * sy; 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 282 uyy = (-2.0 * uc + u[j - 1][i].u + u[j + 1][i].u) * sy; 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/ts/tests/ |
| H A D | ex21.c | 38 PetscScalar ut, u, ue, uw, un, us, uxx, uyy; in FormIFunctionLocal() local 63 uyy = (un - 2.0 * u + us) / (hy * hy); 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 272 uyy = (-2.0 * uc + u[j - 1][i].u + u[j + 1][i].u) * sy; 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 368 uyy = (-2.0 * uc + u_a[j - 1][i].u + u_a[j + 1][i].u) * sy; 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 384 uyy = (-2.0 * uc + u[j - 1][i].u + u[j + 1][i].u) * sy; 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 480 uyy = (-2.0 * uc + u_a[j - 1][i].u + u_a[j + 1][i].u) * sy; 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 576 uyy = (-2.0 * uc + u[j - 1][i].u + u[j + 1][i].u) * sy; 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 116 uyy = (two*u - x(i, j - 1) - x(i, j + 1))*hxdhy 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 199 uyy = (two*u - xx(row - gxm) - xx(row + gxm))*hxdhy 200 ff(rowf) = uxx + uyy - sc*exp(u)
|