| /petsc/src/mat/impls/baij/seq/ |
| H A D | baij2.c | 285 const PetscScalar *x, *xb; in MatMult_SeqBAIJ_2() local 316 xb = x + 2 * (*idx++); in MatMult_SeqBAIJ_2() 317 x1 = xb[0]; in MatMult_SeqBAIJ_2() 318 x2 = xb[1]; in MatMult_SeqBAIJ_2() 338 const PetscScalar *x, *xb; in MatMult_SeqBAIJ_3() local 344 #pragma disjoint(*v, *z, *xb) in MatMult_SeqBAIJ_3() 373 xb = x + 3 * (*idx++); in MatMult_SeqBAIJ_3() 374 x1 = xb[0]; in MatMult_SeqBAIJ_3() 375 x2 = xb[1]; in MatMult_SeqBAIJ_3() 376 x3 = xb[2]; in MatMult_SeqBAIJ_3() [all …]
|
| H A D | baij.c | 222 const PetscScalar *b, *xb; in MatSOR_SeqBAIJ() local 546 xb = t; in MatSOR_SeqBAIJ() 547 } else xb = b; in MatSOR_SeqBAIJ() 553 s[0] = xb[i2]; in MatSOR_SeqBAIJ() 561 s[0] = xb[i2]; in MatSOR_SeqBAIJ() 573 s[0] = xb[i2]; in MatSOR_SeqBAIJ() 574 s[1] = xb[i2 + 1]; in MatSOR_SeqBAIJ() 584 s[0] = xb[i2]; in MatSOR_SeqBAIJ() 585 s[1] = xb[i2 + 1]; in MatSOR_SeqBAIJ() 601 s[0] = xb[i2]; in MatSOR_SeqBAIJ() [all …]
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij2.c | 291 const PetscScalar *x, *xb; in MatMult_SeqSBAIJ_2() local 304 xb = x; in MatMult_SeqSBAIJ_2() 308 x1 = xb[0]; in MatMult_SeqSBAIJ_2() 309 x2 = xb[1]; in MatMult_SeqSBAIJ_2() 331 xb += 2; in MatMult_SeqSBAIJ_2() 345 const PetscScalar *x, *xb; in MatMult_SeqSBAIJ_3() local 358 xb = x; in MatMult_SeqSBAIJ_3() 362 x1 = xb[0]; in MatMult_SeqSBAIJ_3() 363 x2 = xb[1]; in MatMult_SeqSBAIJ_3() 364 x3 = xb[2]; in MatMult_SeqSBAIJ_3() [all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex4.c | 119 PetscScalar f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb; in FormObjectiveLocal() local 130 xl = xr = xb = xt = xc; in FormObjectiveLocal() 137 xb = user->bottom[i + 1]; in FormObjectiveLocal() 138 } else xb = x[j - 1][i]; in FormObjectiveLocal() 151 d4 = (xc - xb); in FormObjectiveLocal() 217 PetscScalar f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormFunctionLocal() local 230 xlt = xrb = xl = xr = xb = xt = xc; in FormFunctionLocal() 238 xb = user->bottom[i + 1]; in FormFunctionLocal() 240 } else xb = x[j - 1][i]; in FormFunctionLocal() 258 d4 = (xc - xb); in FormFunctionLocal() [all …]
|
| H A D | ex58.c | 149 …PetscScalar f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormGradient() local 181 xlt = xrb = xl = xr = xb = xt = xc; in FormGradient() 189 xb = ctx->bottom[i + 1]; in FormGradient() 191 } else xb = x[j - 1][i]; in FormGradient() 209 d4 = (xc - xb); in FormGradient() 211 d6 = (xrb - xb); in FormGradient() 277 …PetscScalar f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormJacobian() local 312 xlt = xrb = xl = xr = xb = xt = xc; in FormJacobian() 322 xb = ctx->bottom[i + 1]; in FormJacobian() 324 } else xb = x[j - 1][i]; in FormJacobian() [all …]
|
| /petsc/src/tao/complementarity/tutorials/ |
| H A D | minsurf1.c | 142 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormConstraints() local 161 xlt = xrb = xl = xr = xb = xt = xc; in FormConstraints() 171 xb = user->bottom[i + 1]; in FormConstraints() 174 xb = x[row - mx]; in FormConstraints() 197 d4 = (xc - xb); in FormConstraints() 199 d6 = (xrb - xb); in FormConstraints() 264 …PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt… in FormJacobian() local 285 xlt = xrb = xl = xr = xb = xt = xc; in FormJacobian() 296 xb = user->bottom[i + 1]; in FormJacobian() 299 xb = x[row - mx]; in FormJacobian() [all …]
|
| /petsc/src/tao/unconstrained/tutorials/ |
| H A D | minsurf1.c | 131 …PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt… in FormFunctionGradient() local 148 xlt = xrb = xl = xr = xb = xt = xc; in FormFunctionGradient() 157 xb = user->bottom[i + 1]; in FormFunctionGradient() 160 xb = x[row - mx]; in FormFunctionGradient() 183 d4 = (xc - xb); in FormFunctionGradient() 185 d6 = (xrb - xb); in FormFunctionGradient() 310 …PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt… in QuadraticH() local 331 xlt = xrb = xl = xr = xb = xt = xc; in QuadraticH() 342 xb = user->bottom[i + 1]; in QuadraticH() 345 xb = x[row - mx]; in QuadraticH() [all …]
|
| H A D | minsurf2.c | 170 PetscReal f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb; in FormFunction() local 199 xb = user->bottom[i - xs + 1]; in FormFunction() 201 xb = x[j - 1][i]; in FormFunction() 219 d4 = (xc - xb); in FormFunction() 304 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormFunctionGradient() local 327 xlt = xrb = xl = xr = xb = xt = xc; in FormFunctionGradient() 337 xb = user->bottom[i - xs + 1]; in FormFunctionGradient() 340 xb = x[j - 1][i]; in FormFunctionGradient() 363 d4 = (xc - xb); in FormFunctionGradient() 365 d6 = (xrb - xb); in FormFunctionGradient() [all …]
|
| /petsc/src/tao/bound/tutorials/ |
| H A D | plate2f.F90 | 76 PetscReal xc, xl, xr, xt, xb, xlt, xrb 117 xb = xc 131 xb = bottom_v(1 + i - xs + 1) 134 xb = x_v(1 + row - gxm) 162 d4 = xc - xb 164 d6 = xrb - xb 308 PetscReal xc, xl, xr, xt, xb, xlt, xrb 356 xb = xc 370 xb = bottom_v(1 + i - xs + 1) 373 xb = x_v(1 + row - gxm) [all …]
|
| H A D | plate2.c | 207 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormFunctionGradient() local 238 xlt = xrb = xl = xr = xb = xt = xc; in FormFunctionGradient() 248 xb = bottom[i - xs + 1]; in FormFunctionGradient() 251 xb = x[row - gxm]; in FormFunctionGradient() 274 d4 = (xc - xb); in FormFunctionGradient() 276 d6 = (xrb - xb); in FormFunctionGradient() 423 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormHessian() local 460 xlt = xrb = xl = xr = xb = xt = xc; in FormHessian() 471 xb = bottom[i - xs + 1]; in FormHessian() 474 xb = x[row - gxm]; in FormHessian() [all …]
|
| /petsc/src/dm/impls/swarm/ |
| H A D | data_bucket.c | 408 PetscErrorCode DMSwarmDataBucketCopyPoint(const DMSwarmDataBucket xb, const PetscInt pid_x, const D… in DMSwarmDataBucketCopyPoint() argument 413 for (f = 0; f < xb->nfields; ++f) { in DMSwarmDataBucketCopyPoint() 417 PetscCall(DMSwarmDataFieldGetAccess(xb->field[f])); in DMSwarmDataBucketCopyPoint() 418 if (xb != yb) PetscCall(DMSwarmDataFieldGetAccess(yb->field[f])); in DMSwarmDataBucketCopyPoint() 419 PetscCall(DMSwarmDataFieldAccessPoint(xb->field[f], pid_x, &src)); in DMSwarmDataBucketCopyPoint() 421 PetscCall(PetscMemcpy(dest, src, xb->field[f]->atomic_size)); in DMSwarmDataBucketCopyPoint() 422 PetscCall(DMSwarmDataFieldRestoreAccess(xb->field[f])); in DMSwarmDataBucketCopyPoint() 423 if (xb != yb) PetscCall(DMSwarmDataFieldRestoreAccess(yb->field[f])); in DMSwarmDataBucketCopyPoint()
|
| /petsc/src/mat/impls/kaij/ |
| H A D | kaij.c | 861 const PetscScalar *b, *xb, *idiag; in MatSOR_SeqKAIJ() local 929 xb = t; in MatSOR_SeqKAIJ() 930 } else xb = b; in MatSOR_SeqKAIJ() 934 PetscCall(PetscArraycpy(w, xb + i2, bs)); in MatSOR_SeqKAIJ() 944 PetscCall(PetscArraycpy(w, xb + i2, bs)); in MatSOR_SeqKAIJ() 1039 xb = t; in MatSOR_SeqKAIJ() 1040 } else xb = b; in MatSOR_SeqKAIJ() 1044 if (xb == b) { in MatSOR_SeqKAIJ() 1103 PetscCall(PetscArraycpy(w, xb + i2, bs)); in MatSOR_SeqKAIJ()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | inode.c | 2500 const PetscScalar *xb, *b; in MatSOR_SeqAIJ_Inode() local 2704 xb = t; in MatSOR_SeqAIJ_Inode() 2706 } else xb = b; in MatSOR_SeqAIJ_Inode() 2720 sum1 = xb[row]; in MatSOR_SeqAIJ_Inode() 2740 sum1 = xb[row]; in MatSOR_SeqAIJ_Inode() 2741 sum2 = xb[row - 1]; in MatSOR_SeqAIJ_Inode() 2766 sum1 = xb[row]; in MatSOR_SeqAIJ_Inode() 2767 sum2 = xb[row - 1]; in MatSOR_SeqAIJ_Inode() 2768 sum3 = xb[row - 2]; in MatSOR_SeqAIJ_Inode() 2797 sum1 = xb[row]; in MatSOR_SeqAIJ_Inode() [all …]
|
| H A D | aij.c | 1859 const PetscScalar *b, *bs, *xb, *ts; in MatSOR_SeqAIJ() local 1953 xb = t; in MatSOR_SeqAIJ() 1955 } else xb = b; in MatSOR_SeqAIJ() 1961 sum = xb[i]; in MatSOR_SeqAIJ() 1963 if (xb == b) { in MatSOR_SeqAIJ() 1990 xb = t; in MatSOR_SeqAIJ() 1992 } else xb = b; in MatSOR_SeqAIJ() 1995 sum = xb[i]; in MatSOR_SeqAIJ() 1996 if (xb == b) { in MatSOR_SeqAIJ() 2011 if (xb == b) { in MatSOR_SeqAIJ()
|
| /petsc/src/sys/classes/draw/impls/x/ |
| H A D | xops.c | 19 int xa, ya, xb, yb, xmax = XiWin->w - 1, ymax = XiWin->h - 1; in PetscDrawSetViewport_X() local 25 xb = (int)(xr * (PetscReal)xmax); in PetscDrawSetViewport_X() 29 box.width = (unsigned short)(xb + 1 - xa); in PetscDrawSetViewport_X() 307 int xb = (int)(xr * xmax), yb = ymax - (int)(yl * ymax); in PetscDrawClear_X() local 308 unsigned int w = (unsigned int)(xb + 1 - xa); in PetscDrawClear_X()
|
| /petsc/src/mat/impls/sell/seq/ |
| H A D | sell.c | 1692 const PetscScalar *b, *xb; in MatSOR_SeqSELL() local 1722 xb = t; in MatSOR_SeqSELL() 1724 } else xb = b; in MatSOR_SeqSELL() 1728 sum = xb[i]; in MatSOR_SeqSELL() 1731 if (xb == b) { in MatSOR_SeqSELL() 1755 xb = t; in MatSOR_SeqSELL() 1757 } else xb = b; in MatSOR_SeqSELL() 1761 sum = xb[i]; in MatSOR_SeqSELL() 1762 if (xb == b) { in MatSOR_SeqSELL() 1773 if (xb == b) { in MatSOR_SeqSELL()
|
| /petsc/src/sys/classes/draw/impls/image/ |
| H A D | drawimage.c | 27 int xb = (int)(xr * xmax), yb = ymax - (int)(yl * ymax); in PetscDrawSetViewport_Image() local 28 PetscImageSetClip(img, xa, ya, xb + 1 - xa, yb + 1 - ya); in PetscDrawSetViewport_Image()
|
| /petsc/src/tao/constrained/impls/ipm/ |
| H A D | pdipm.c | 843 PetscInt offset, *xa, *xb, i, j, rstart, rend; in TaoSetup_PDIPM() local 1021 PetscCall(PetscMalloc2(pdipm->nx + pdipm->nce, &xa, 2 * pdipm->nci, &xb)); in TaoSetup_PDIPM() 1023 for (i = 0; i < 2 * pdipm->nci; i++) xb[i] = pdipm->off_lambdai + i; in TaoSetup_PDIPM() 1026 PetscCall(ISCreateGeneral(comm, 2 * pdipm->nci, xb, PETSC_OWN_POINTER, &pdipm->is2)); in TaoSetup_PDIPM()
|
| /petsc/src/mat/impls/sbaij/mpi/ |
| H A D | mpisbaij.c | 1657 Mat_SeqBAIJ *xb, *yb; in MatAXPY_MPISBAIJ() local 1666 xb = (Mat_SeqBAIJ *)xx->B->data; in MatAXPY_MPISBAIJ() 1668 PetscCall(PetscBLASIntCast(xb->nz, &bnz)); in MatAXPY_MPISBAIJ() 1669 PetscCallBLAS("BLASaxpy", BLASaxpy_(&bnz, &alpha, xb->a, &one, yb->a, &one)); in MatAXPY_MPISBAIJ()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plex.c | 338 PetscScalar *xa, *xb, *svals; in DMPlexVecView1D() local 343 PetscCall(DMPlexPointLocalRead(cdm, cone[1], coords, &xb)); in DMPlexVecView1D() 355 …PetscCall(PetscDrawLGAddCommonPoint(lg, 0.5 * (PetscRealPart(xa[0]) + PetscRealPart(xb[0])), vals)… in DMPlexVecView1D() 360 PetscCall(PetscDrawLGAddCommonPoint(lg, PetscRealPart(xb[0]), vals)); in DMPlexVecView1D() 366 PetscScalar *xa, *xb, *svals; in DMPlexVecView1D() local 371 PetscCall(DMPlexPointLocalRead(cdm, cone[1], coords, &xb)); in DMPlexVecView1D() 376 …PetscCall(PetscDrawLGAddCommonPoint(lg, 0.5 * (PetscRealPart(xa[0]) + PetscRealPart(xb[0])), vals)… in DMPlexVecView1D()
|