| /petsc/src/vec/vec/impls/seq/ |
| H A D | bvec1.c | 125 PetscScalar *zz; in VecAXPBYPCZ_Seq() local 130 PetscCall(VecGetArray(zin, &zz)); in VecAXPBYPCZ_Seq() 132 for (PetscInt i = 0; i < n; ++i) zz[i] = xx[i] + beta * yy[i] + gamma * zz[i]; in VecAXPBYPCZ_Seq() 134 for (PetscInt i = 0; i < n; ++i) zz[i] = alpha * xx[i] + beta * yy[i] + zz[i]; in VecAXPBYPCZ_Seq() 136 for (PetscInt i = 0; i < n; ++i) zz[i] = alpha * xx[i] + beta * yy[i]; in VecAXPBYPCZ_Seq() 139 for (PetscInt i = 0; i < n; ++i) zz[i] = alpha * xx[i] + beta * yy[i] + gamma * zz[i]; in VecAXPBYPCZ_Seq() 144 PetscCall(VecRestoreArray(zin, &zz)); in VecAXPBYPCZ_Seq()
|
| /petsc/src/dm/tests/ |
| H A D | ex27.c | 26 Vec coords, x, y, z, xx, yy, zz; in main() local 85 PetscCall(DMGetGlobalVector(da, &zz)); in main() 86 PetscCall(PetscObjectSetName((PetscObject)zz, "FFTW reconstructed vector")); in main() 127 PetscCall(MatMult(AA, xx, zz)); in main() 133 PetscCall(MatMultTranspose(AA, yy, zz)); in main() 147 PetscCall(VecView(zz, PETSC_VIEWER_STDOUT_WORLD)); in main() 159 PetscCall(VecScale(zz, s)); in main() 160 PetscCall(VecAXPY(xx, -1.0, zz)); in main() 173 PetscCall(VecAXPY(z, -1.0, zz)); in main() 184 PetscCall(DMRestoreGlobalVector(da, &zz)); in main()
|
| H A D | ex28.c | 28 Vec coords, x, y, z, xx, yy, zz; in main() local 94 PetscCall(DMGetGlobalVector(da, &zz)); in main() 95 PetscCall(PetscObjectSetName((PetscObject)zz, "FFTW reconstructed vector")); in main() 158 …PetscCall(VecStrideScatterAll(zzsplit, zz, INSERT_VALUES)); /*YES! 'Scatter' means 'collect' (or m… in main() 171 PetscCall(VecView(zz, PETSC_VIEWER_STDOUT_WORLD)); in main() 183 PetscCall(VecScale(zz, s)); in main() 184 PetscCall(VecAXPY(xx, -1.0, zz)); in main() 197 PetscCall(VecAXPY(z, -1.0, zz)); in main() 208 PetscCall(DMRestoreGlobalVector(da, &zz)); in main()
|
| H A D | ex36.c | 61 PetscScalar xx, yy, zz; in DAApplyConformalMapping() local 105 xx = yy = zz = 0.0; in DAApplyConformalMapping() 109 zz += Ni[p] * zn[p]; in DAApplyConformalMapping() 113 XX[dim * i + 2] = zz; in DAApplyConformalMapping() 201 PetscScalar xx, yy, zz; in DAApplyTrilinearMapping() local 234 xx = yy = zz = 0.0; in DAApplyTrilinearMapping() 238 zz += Ni[p] * zn[p]; in DAApplyTrilinearMapping() 242 XX[k][j][i].z = zz; in DAApplyTrilinearMapping()
|
| /petsc/src/mat/impls/maij/ |
| H A D | maij.c | 194 …tic inline PetscErrorCode MatMult_MatMultAdd_SeqMAIJ_Template(Mat A, Vec xx, Vec yy, Vec zz, int N) in MatMult_MatMultAdd_SeqMAIJ_Template() argument 211 if (mult_add && yy != zz) PetscCall(VecCopy(yy, zz)); in MatMult_MatMultAdd_SeqMAIJ_Template() 214 PetscCall(VecGetArray(zz, &z)); in MatMult_MatMultAdd_SeqMAIJ_Template() 216 PetscCall(VecGetArrayWrite(zz, &z)); in MatMult_MatMultAdd_SeqMAIJ_Template() 248 PetscCall(VecRestoreArray(zz, &z)); in MatMult_MatMultAdd_SeqMAIJ_Template() 250 PetscCall(VecRestoreArrayWrite(zz, &z)); in MatMult_MatMultAdd_SeqMAIJ_Template() 255 …rorCode MatMultTranspose_MatMultTransposeAdd_SeqMAIJ_Template(Mat A, Vec xx, Vec yy, Vec zz, int N) in MatMultTranspose_MatMultTransposeAdd_SeqMAIJ_Template() argument 272 if (yy != zz) PetscCall(VecCopy(yy, zz)); in MatMultTranspose_MatMultTransposeAdd_SeqMAIJ_Template() 274 PetscCall(VecSet(zz, 0.0)); in MatMultTranspose_MatMultTransposeAdd_SeqMAIJ_Template() 277 PetscCall(VecGetArray(zz, &z)); in MatMultTranspose_MatMultTransposeAdd_SeqMAIJ_Template() [all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baij2.c | 236 PetscErrorCode MatMult_SeqBAIJ_1(Mat A, Vec xx, Vec zz) in MatMult_SeqBAIJ_1() argument 248 PetscCall(VecGetArrayWrite(zz, &z)); in MatMult_SeqBAIJ_1() 276 PetscCall(VecRestoreArrayWrite(zz, &z)); in MatMult_SeqBAIJ_1() 281 PetscErrorCode MatMult_SeqBAIJ_2(Mat A, Vec xx, Vec zz) in MatMult_SeqBAIJ_2() argument 293 PetscCall(VecGetArrayWrite(zz, &zarray)); in MatMult_SeqBAIJ_2() 329 PetscCall(VecRestoreArrayWrite(zz, &zarray)); in MatMult_SeqBAIJ_2() 334 PetscErrorCode MatMult_SeqBAIJ_3(Mat A, Vec xx, Vec zz) in MatMult_SeqBAIJ_3() argument 349 PetscCall(VecGetArrayWrite(zz, &zarray)); in MatMult_SeqBAIJ_3() 390 PetscCall(VecRestoreArrayWrite(zz, &zarray)); in MatMult_SeqBAIJ_3() 395 PetscErrorCode MatMult_SeqBAIJ_4(Mat A, Vec xx, Vec zz) in MatMult_SeqBAIJ_4() argument [all …]
|
| /petsc/src/mat/impls/baij/seq/baijmkl/ |
| H A D | baijmkl.c | 271 static PetscErrorCode MatMultAdd_SeqBAIJMKL_SpMV2(Mat A, Vec xx, Vec yy, Vec zz) in MatMultAdd_SeqBAIJMKL_SpMV2() argument 283 PetscCall(VecCopy(yy, zz)); in MatMultAdd_SeqBAIJMKL_SpMV2() 288 PetscCall(VecGetArrayPair(yy, zz, &y, &z)); in MatMultAdd_SeqBAIJMKL_SpMV2() 296 if (zz == yy) { in MatMultAdd_SeqBAIJMKL_SpMV2() 309 PetscCall(VecRestoreArrayPair(yy, zz, &y, &z)); in MatMultAdd_SeqBAIJMKL_SpMV2() 313 static PetscErrorCode MatMultTransposeAdd_SeqBAIJMKL_SpMV2(Mat A, Vec xx, Vec yy, Vec zz) in MatMultTransposeAdd_SeqBAIJMKL_SpMV2() argument 326 PetscCall(VecCopy(yy, zz)); in MatMultTransposeAdd_SeqBAIJMKL_SpMV2() 331 PetscCall(VecGetArrayPair(yy, zz, &y, &z)); in MatMultTransposeAdd_SeqBAIJMKL_SpMV2() 339 if (zz == yy) { in MatMultTransposeAdd_SeqBAIJMKL_SpMV2() 352 PetscCall(VecRestoreArrayPair(yy, zz, &y, &z)); in MatMultTransposeAdd_SeqBAIJMKL_SpMV2()
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij2.c | 287 PetscErrorCode MatMult_SeqSBAIJ_2(Mat A, Vec xx, Vec zz) in MatMult_SeqSBAIJ_2() argument 298 PetscCall(VecSet(zz, zero)); in MatMult_SeqSBAIJ_2() 301 PetscCall(VecGetArray(zz, &z)); in MatMult_SeqSBAIJ_2() 336 PetscCall(VecRestoreArray(zz, &z)); in MatMult_SeqSBAIJ_2() 341 PetscErrorCode MatMult_SeqSBAIJ_3(Mat A, Vec xx, Vec zz) in MatMult_SeqSBAIJ_3() argument 352 PetscCall(VecSet(zz, zero)); in MatMult_SeqSBAIJ_3() 355 PetscCall(VecGetArray(zz, &z)); in MatMult_SeqSBAIJ_3() 394 PetscCall(VecRestoreArray(zz, &z)); in MatMult_SeqSBAIJ_3() 399 PetscErrorCode MatMult_SeqSBAIJ_4(Mat A, Vec xx, Vec zz) in MatMult_SeqSBAIJ_4() argument 410 PetscCall(VecSet(zz, zero)); in MatMult_SeqSBAIJ_4() [all …]
|
| H A D | relax.h | 54 PetscErrorCode MatMult_SeqSBAIJ_1_ushort(Mat A, Vec xx, Vec zz) in MatMult_SeqSBAIJ_1_ushort() argument 56 PetscErrorCode MatMult_SeqSBAIJ_1(Mat A, Vec xx, Vec zz) in MatMult_SeqSBAIJ_1_ushort() 77 PetscCall(VecSet(zz, 0.0)); in MatMult_SeqSBAIJ_1_ushort() 79 PetscCall(VecGetArray(zz, &z)); in MatMult_SeqSBAIJ_1_ushort() 116 PetscCall(VecRestoreArray(zz, &z)); in MatMult_SeqSBAIJ_1_ushort()
|
| /petsc/src/mat/impls/aij/seq/aijmkl/ |
| H A D | aijmkl.c | 483 static PetscErrorCode MatMultAdd_SeqAIJMKL(Mat A, Vec xx, Vec yy, Vec zz) in MatMultAdd_SeqAIJMKL() argument 505 PetscCall(VecGetArrayPair(yy, zz, &y, &z)); in MatMultAdd_SeqAIJMKL() 511 if (zz == yy) { in MatMultAdd_SeqAIJMKL() 525 PetscCall(VecRestoreArrayPair(yy, zz, &y, &z)); in MatMultAdd_SeqAIJMKL() 531 PetscErrorCode MatMultAdd_SeqAIJMKL_SpMV2(Mat A, Vec xx, Vec yy, Vec zz) in MatMultAdd_SeqAIJMKL_SpMV2() argument 546 PetscCall(VecGetArrayPair(yy, zz, &y, &z)); in MatMultAdd_SeqAIJMKL_SpMV2() 548 PetscCall(VecRestoreArrayPair(yy, zz, &y, &z)); in MatMultAdd_SeqAIJMKL_SpMV2() 553 PetscCall(VecGetArrayPair(yy, zz, &y, &z)); in MatMultAdd_SeqAIJMKL_SpMV2() 562 if (zz == yy) { in MatMultAdd_SeqAIJMKL_SpMV2() 575 PetscCall(VecRestoreArrayPair(yy, zz, &y, &z)); in MatMultAdd_SeqAIJMKL_SpMV2() [all …]
|
| /petsc/src/ts/tests/ |
| H A D | ex3.c | 341 PetscScalar dd, dl, zip, zipq, zz, b_z, bb_z, bij; in femBg() local 385 zz = zquad[il][iquad]; in femBg() 389 b_z = bspl(z, zz, il, iq, nli, 2); in femBg() 395 bb_z = bspl(z, zz, il, iqq, nli, 2); in femBg() 421 PetscScalar dd, dl, zip, zipq, zz, bb, bbb, aij; in femA() local 466 zz = rquad[il][iquad]; in femA() 470 bb = bspl(z, zz, il, iq, nli, 1); in femA() 475 bbb = bspl(z, zz, il, iqq, nli, 1); in femA()
|
| /petsc/src/mat/impls/aij/seq/aijsell/ |
| H A D | aijsell.c | 167 static PetscErrorCode MatMultAdd_SeqAIJSELL(Mat A, Vec xx, Vec yy, Vec zz) in MatMultAdd_SeqAIJSELL() argument 173 PetscCall(MatMultAdd_SeqSELL(aijsell->S, xx, yy, zz)); in MatMultAdd_SeqAIJSELL() 177 static PetscErrorCode MatMultTransposeAdd_SeqAIJSELL(Mat A, Vec xx, Vec yy, Vec zz) in MatMultTransposeAdd_SeqAIJSELL() argument 183 PetscCall(MatMultTransposeAdd_SeqSELL(aijsell->S, xx, yy, zz)); in MatMultTransposeAdd_SeqAIJSELL()
|
| /petsc/src/mat/impls/aij/seq/seqviennacl/ |
| H A D | aijviennacl.cxx | 207 static PetscErrorCode MatMultAdd_SeqAIJViennaCL(Mat A, Vec xx, Vec yy, Vec zz) in MatMultAdd_SeqAIJViennaCL() argument 221 PetscCall(VecViennaCLGetArrayWrite(zz, &zgpu)); in MatMultAdd_SeqAIJViennaCL() 225 if (zz != yy) *zgpu = *ygpu + *viennaclstruct->tempvec; in MatMultAdd_SeqAIJViennaCL() 231 PetscCall(VecViennaCLRestoreArrayWrite(zz, &zgpu)); in MatMultAdd_SeqAIJViennaCL() 238 PetscCall(VecCopy_SeqViennaCL(yy, zz)); in MatMultAdd_SeqAIJViennaCL()
|
| /petsc/src/mat/impls/kaij/ |
| H A D | kaij.c | 652 static PetscErrorCode MatMultAdd_SeqKAIJ(Mat A, Vec xx, Vec yy, Vec zz) in MatMultAdd_SeqKAIJ() argument 664 PetscCall(VecSet(zz, 0.0)); in MatMultAdd_SeqKAIJ() 666 PetscCall(VecCopy(yy, zz)); in MatMultAdd_SeqKAIJ() 671 PetscCall(VecGetArray(zz, &y)); in MatMultAdd_SeqKAIJ() 717 PetscCall(VecRestoreArray(zz, &y)); in MatMultAdd_SeqKAIJ() 1143 static PetscErrorCode MatMultAdd_MPIKAIJ(Mat A, Vec xx, Vec yy, Vec zz) in MatMultAdd_MPIKAIJ() argument 1149 PetscCall(VecSet(zz, 0.0)); in MatMultAdd_MPIKAIJ() 1151 PetscCall(VecCopy(yy, zz)); in MatMultAdd_MPIKAIJ() 1156 PetscCall((*b->AIJ->ops->multadd)(b->AIJ, xx, zz, zz)); in MatMultAdd_MPIKAIJ() 1158 PetscCall((*b->OAIJ->ops->multadd)(b->OAIJ, b->w, zz, zz)); in MatMultAdd_MPIKAIJ()
|
| /petsc/src/mat/impls/aij/seq/kokkos/ |
| H A D | aijkok.kokkos.cxx | 468 static PetscErrorCode MatMultAdd_SeqAIJKokkos(Mat A, Vec xx, Vec yy, Vec zz) in MatMultAdd_SeqAIJKokkos() argument 477 …if (zz != yy) PetscCall(VecCopy(yy, zz)); // depending on yy's sync flags, zz might get its latest… in MatMultAdd_SeqAIJKokkos() 479 PetscCall(VecGetKokkosView(zz, &zv)); // do after VecCopy(yy, zz) to get the latest data on device in MatMultAdd_SeqAIJKokkos() 483 PetscCall(VecRestoreKokkosView(zz, &zv)); in MatMultAdd_SeqAIJKokkos() 490 static PetscErrorCode MatMultTransposeAdd_SeqAIJKokkos(Mat A, Vec xx, Vec yy, Vec zz) in MatMultTransposeAdd_SeqAIJKokkos() argument 501 if (zz != yy) PetscCall(VecCopy(yy, zz)); in MatMultTransposeAdd_SeqAIJKokkos() 503 PetscCall(VecGetKokkosView(zz, &zv)); in MatMultTransposeAdd_SeqAIJKokkos() 514 PetscCall(VecRestoreKokkosView(zz, &zv)); in MatMultTransposeAdd_SeqAIJKokkos() 521 static PetscErrorCode MatMultHermitianTransposeAdd_SeqAIJKokkos(Mat A, Vec xx, Vec yy, Vec zz) in MatMultHermitianTransposeAdd_SeqAIJKokkos() argument 532 if (zz != yy) PetscCall(VecCopy(yy, zz)); in MatMultHermitianTransposeAdd_SeqAIJKokkos() [all …]
|
| /petsc/src/dm/impls/stag/tutorials/ |
| H A D | ex4.c | 362 const PetscReal zz = PetscRealPart(z) / d - 0.5; in GetRho_sinker_box3() local 364 …return (PetscAbsReal(xx) > half_width || PetscAbsReal(yy) > half_width || PetscAbsReal(zz) > half_… in GetRho_sinker_box3() 372 const PetscReal zz = PetscRealPart(z) / d - 0.5; in GetEta_sinker_box3() local 374 …return (PetscAbsReal(xx) > half_width || PetscAbsReal(yy) > half_width || PetscAbsReal(zz) > half_… in GetEta_sinker_box3() 382 const PetscReal zz = PetscRealPart(z) / d - 0.5; in GetRho_sinker_sphere3() local 384 return (xx * xx + yy * yy + zz * zz > half_width * half_width) ? ctx->rho1 : ctx->rho2; in GetRho_sinker_sphere3() 392 const PetscReal zz = PetscRealPart(z) / d - 0.5; in GetEta_sinker_sphere3() local 394 return (xx * xx + yy * yy + zz * zz > half_width * half_width) ? ctx->eta1 : ctx->eta2; in GetEta_sinker_sphere3() 402 const PetscReal zz = PetscRealPart(z) / d - 0.5; in GetEta_blob3() local 403 return ctx->eta1 + ctx->eta2 * PetscExpScalar(-20.0 * (xx * xx + yy * yy + zz * zz)); in GetEta_blob3() [all …]
|
| /petsc/src/mat/impls/aij/mpi/mpihipsparse/ |
| H A D | mpiaijhipsparse.hip.cxx | 337 static PetscErrorCode MatMultAdd_MPIAIJHIPSPARSE(Mat A, Vec xx, Vec yy, Vec zz) in MatMultAdd_MPIAIJHIPSPARSE() argument 343 PetscCall((*a->A->ops->multadd)(a->A, xx, yy, zz)); in MatMultAdd_MPIAIJHIPSPARSE() 345 PetscCall((*a->B->ops->multadd)(a->B, a->lvec, zz, zz)); in MatMultAdd_MPIAIJHIPSPARSE()
|
| /petsc/src/mat/impls/dense/seq/cupm/ |
| H A D | matseqdensecupm.hpp | 805 …se_Seq_CUPM<T>::MatMultAddColumnRange_Dispatch_(Mat A, Vec xx, Vec yy, Vec zz, PetscInt c_start, P… in MatMultAddColumnRange_Dispatch_() argument 812 Vec x = xx, y = yy, z = zz; in MatMultAddColumnRange_Dispatch_() 834 if (zz != yy) { in MatMultAddColumnRange_Dispatch_() 835 …PetscCall(PetscObjectTypeCompareAny(PetscObjectCast(zz), &ziscupm, VecSeq_CUPM::VECSEQCUPM(), VecS… in MatMultAddColumnRange_Dispatch_() 836 if (!ziscupm || zz->boundtocpu) { in MatMultAddColumnRange_Dispatch_() 837 PetscCall(VecCreate(PetscObjectComm(PetscObjectCast(zz)), &z)); in MatMultAddColumnRange_Dispatch_() 838 PetscCall(VecSetLayout(z, zz->map)); in MatMultAddColumnRange_Dispatch_() 866 if (z != zz) { in MatMultAddColumnRange_Dispatch_() 867 PetscCall(VecCopy(z, zz)); in MatMultAddColumnRange_Dispatch_() 886 inline PetscErrorCode MatDense_Seq_CUPM<T>::MatMultAdd_Dispatch_(Mat A, Vec xx, Vec yy, Vec zz) noe… in MatMultAdd_Dispatch_() argument [all …]
|
| /petsc/src/mat/impls/aij/mpi/mpicusparse/ |
| H A D | mpiaijcusparse.cu | 337 static PetscErrorCode MatMultAdd_MPIAIJCUSPARSE(Mat A, Vec xx, Vec yy, Vec zz) in MatMultAdd_MPIAIJCUSPARSE() argument 343 PetscCall((*a->A->ops->multadd)(a->A, xx, yy, zz)); in MatMultAdd_MPIAIJCUSPARSE() 345 PetscCall((*a->B->ops->multadd)(a->B, a->lvec, zz, zz)); in MatMultAdd_MPIAIJCUSPARSE()
|
| /petsc/src/ts/impls/implicit/sundials/ |
| H A D | sundials.c | 54 Vec rr = cvode->w1, zz = cvode->w2; in TSPSolve_Sundials_Petsc() local 62 PetscCall(VecPlaceArray(zz, z_data)); in TSPSolve_Sundials_Petsc() 66 PetscCall(PCApply(pc, rr, zz)); in TSPSolve_Sundials_Petsc() 68 PetscCall(VecResetArray(zz)); in TSPSolve_Sundials_Petsc()
|
| /petsc/src/mat/impls/dense/mpi/ |
| H A D | mpidense.c | 511 static PetscErrorCode MatMultAddColumnRange_MPIDense(Mat mat, Vec xx, Vec yy, Vec zz, PetscInt c_st… in MatMultAddColumnRange_MPIDense() argument 526 …mnRange_C", (Mat, Vec, Vec, Vec, PetscInt, PetscInt), (mdn->A, mdn->lvec, yy, zz, c_start, c_end)); in MatMultAddColumnRange_MPIDense() 530 static PetscErrorCode MatMultAdd_MPIDense(Mat mat, Vec xx, Vec yy, Vec zz) in MatMultAdd_MPIDense() argument 545 PetscCall((*mdn->A->ops->multadd)(mdn->A, mdn->lvec, yy, zz)); in MatMultAdd_MPIDense() 603 …atMultHermitianTransposeAddColumnRange_MPIDense(Mat A, Vec xx, Vec yy, Vec zz, PetscInt c_start, P… in MatMultHermitianTransposeAddColumnRange_MPIDense() argument 612 PetscCall(VecCopy(yy, zz)); in MatMultHermitianTransposeAddColumnRange_MPIDense() 616 PetscCall(VecGetArrayAndMemType(zz, &ay, &aymtype)); in MatMultHermitianTransposeAddColumnRange_MPIDense() 620 PetscCall(VecRestoreArrayAndMemType(zz, &ay)); in MatMultHermitianTransposeAddColumnRange_MPIDense() 624 static PetscErrorCode MatMultTransposeAddKernel_MPIDense(Mat A, Vec xx, Vec yy, Vec zz, PetscBool h… in MatMultTransposeAddKernel_MPIDense() argument 633 PetscCall(VecCopy(yy, zz)); in MatMultTransposeAddKernel_MPIDense() [all …]
|
| /petsc/src/mat/impls/sell/mpi/ |
| H A D | mpisell.c | 390 static PetscErrorCode MatMultAdd_MPISELL(Mat A, Vec xx, Vec yy, Vec zz) in MatMultAdd_MPISELL() argument 396 PetscCall((*a->A->ops->multadd)(a->A, xx, yy, zz)); in MatMultAdd_MPISELL() 398 PetscCall((*a->B->ops->multadd)(a->B, a->lvec, zz, zz)); in MatMultAdd_MPISELL() 457 static PetscErrorCode MatMultTransposeAdd_MPISELL(Mat A, Vec xx, Vec yy, Vec zz) in MatMultTransposeAdd_MPISELL() argument 465 PetscCall((*a->A->ops->multtransposeadd)(a->A, xx, yy, zz)); in MatMultTransposeAdd_MPISELL() 467 PetscCall(VecScatterBegin(a->Mvctx, a->lvec, zz, ADD_VALUES, SCATTER_REVERSE)); in MatMultTransposeAdd_MPISELL() 468 PetscCall(VecScatterEnd(a->Mvctx, a->lvec, zz, ADD_VALUES, SCATTER_REVERSE)); in MatMultTransposeAdd_MPISELL()
|
| /petsc/src/ksp/ksp/utils/lmvm/blas_cyclic/ |
| H A D | blas_cyclic.c | 556 PETSC_INTERN PetscErrorCode MatMultAddColumnRange(Mat A, Vec xx, Vec zz, Vec yy, PetscInt c_start, … in MatMultAddColumnRange() argument 560 …PetscUseMethod(A, "MatMultAddColumnRange_C", (Mat, Vec, Vec, Vec, PetscInt, PetscInt), (A, xx, zz,… in MatMultAddColumnRange() 574 PETSC_INTERN PetscErrorCode MatMultHermitianTransposeAddColumnRange(Mat A, Vec xx, Vec zz, Vec yy, … in MatMultHermitianTransposeAddColumnRange() argument 578 …nsposeAddColumnRange_C", (Mat, Vec, Vec, Vec, PetscInt, PetscInt), (A, xx, zz, yy, c_start, c_end)… in MatMultHermitianTransposeAddColumnRange()
|
| /petsc/src/ts/utils/dmplexlandau/ |
| H A D | plexland.c | 211 …vJ, *xx = (PetscReal *)ctx->SData_d.x, *yy = (PetscReal *)ctx->SData_d.y, *zz = (PetscReal *)ctx->… in LandauFormJacobian_Internal() local 341 …const PetscReal vj[3] = {xx[jpidx_glb], yy[jpidx_glb], zz ? zz[jpidx_glb] : 0}, wj = ww[jpi… in LandauFormJacobian_Internal() 359 …PetscReal U[3][3], z = zz[ipidx], mask = (PetscAbs(vj[0] - x) < 100 * PETSC_SQRT_MACHINE_EPSILON &… in LandauFormJacobian_Internal() 1699 PetscReal *invJe, *ww, *xx, *yy, *zz = NULL, *invJ_a; in CreateStaticData() local 1716 if (dim == 3) PetscCall(PetscMalloc1(nip_glb, &zz)); in CreateStaticData() 1783 if (dim == 3) zz[gidx] = vj[qj * dim + 2]; in CreateStaticData() 1784 … xx[gidx] + yy[gidx] * yy[gidx] + zz[gidx] * zz[gidx]), ej, qj, dg2_c2, (double)xx[gidx], (double)… in CreateStaticData() 1797 if (dim == 3) zz[gidx] = eGradPhi[2]; in CreateStaticData() 1802 if (dim == 3) zz[gidx] = vj[qj * dim + 2]; in CreateStaticData() 1816 …et, nu_alpha, nu_beta, invMass, (PetscReal *)ctx->lambdas, invJ_a, xx, yy, zz, ww, &ctx->SData_d)); in CreateStaticData() [all …]
|
| /petsc/src/dm/impls/swarm/tests/ |
| H A D | ex7.c | 31 PetscErrorCode MatMultAddMtM_SeqAIJ(Mat MtM, Vec xx, Vec yy, Vec zz) in MatMultAddMtM_SeqAIJ() argument 39 PetscCall(MatMultAdd(matshellctx->MpTrans, matshellctx->ff, yy, zz)); in MatMultAddMtM_SeqAIJ()
|