| /petsc/src/dm/tests/ |
| H A D | ex36.c | 60 PetscScalar *XX; in DAApplyConformalMapping() local
78 PetscCall(VecGetArray(Gcoords, &XX)); in DAApplyConformalMapping()
87 PetscScalar xi = XX[dim * i]; in DAApplyConformalMapping()
88 PetscScalar eta = XX[dim * i + 1]; in DAApplyConformalMapping()
89 PetscScalar zeta = XX[dim * i + 2]; in DAApplyConformalMapping()
111 XX[dim * i] = xx; in DAApplyConformalMapping()
112 XX[dim * i + 1] = yy; in DAApplyConformalMapping()
113 XX[dim * i + 2] = zz; in DAApplyConformalMapping()
119 xx = XX[dim * i] - 0.8; in DAApplyConformalMapping()
120 yy = XX[dim * i + 1] + 1.5; in DAApplyConformalMapping()
[all …]
|
| /petsc/src/ts/trajectory/impls/visualization/ |
| H A D | trajvisualization.c | 40 Vec XX; in TSTrajectorySet_Visualization() local
41 PetscCall((*tj->transform)(tj->transformctx, X, &XX)); in TSTrajectorySet_Visualization()
42 PetscCall(VecView(XX, viewer)); in TSTrajectorySet_Visualization()
43 PetscCall(VecDestroy(&XX)); in TSTrajectorySet_Visualization()
54 Vec XX; in TSTrajectorySet_Visualization() local
55 PetscCall((*tj->transform)(tj->transformctx, X, &XX)); in TSTrajectorySet_Visualization()
56 PetscCall(VecView(XX, viewer)); in TSTrajectorySet_Visualization()
57 PetscCall(VecDestroy(&XX)); in TSTrajectorySet_Visualization()
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex26.c | 118 const PetscScalar *xx, *XX; in main() local
121 PetscCall(MatDenseGetArrayRead(X, &XX)); in main()
124 if (PetscAbsScalar(xx[i] - XX[nlocal * n + i]) > PETSC_SMALL) { in main()
125 … (double)PetscImaginaryPart(xx[i]), (double)PetscRealPart(XX[i]), (double)PetscImaginaryPart(XX[i]… in main()
129 PetscCall(MatDenseRestoreArrayRead(X, &XX)); in main()
|
| /petsc/src/ksp/ksp/impls/gmres/dgmres/ |
| H A D | dgmres.c | 71 XX = NULL; in KSPSetUp_DGMRES()
271 if (XX) { in KSPDestroy_DGMRES()
272 PetscCall(VecDestroyVecs(neig1, &XX)); in KSPDestroy_DGMRES()
563 if (!XX) PetscCall(VecDuplicateVecs(VEC_VV(0), neig1, &XX)); in KSPDGMRESComputeDeflationData_DGMRES()
564 for (j = 0; j < neig; j++) PetscCall(VecMAXPBY(XX[j], n, &SR[j * N], 0, &VEC_VV(0))); in KSPDGMRESComputeDeflationData_DGMRES()
573 PetscCall(VecDot(XX[j], UU[i], &alpha)); in KSPDGMRESComputeDeflationData_DGMRES()
575 PetscCall(VecAXPY(XX[j], -alpha, UU[i])); in KSPDGMRESComputeDeflationData_DGMRES()
581 for (j = 0; j < neig; j++) PetscCall(KSP_PCApplyBAorAB(ksp, XX[j], MX[j], VEC_TEMP_MATOP)); in KSPDGMRESComputeDeflationData_DGMRES()
592 for (j = 0; j < neig; j++) PetscCall(VecMDot(MX[j], neig, XX, &XMX[j * neig1])); in KSPDGMRESComputeDeflationData_DGMRES()
600 for (j = 0; j < r; j++) PetscCall(VecMDot(MU[j], neig, XX, &XMU[j * neig1])); in KSPDGMRESComputeDeflationData_DGMRES()
[all …]
|
| H A D | dgmresimpl.h | 72 #define XX dgmres->X macro
|
| /petsc/src/tao/pde_constrained/tutorials/ |
| H A D | elliptic.c | 540 Vec XX, YY, ZZ, XXwork, YYwork, ZZwork, UTwork; in EllipticInitialize() local
666 PetscCall(VecCreate(PETSC_COMM_WORLD, &XX)); in EllipticInitialize()
668 PetscCall(VecSetSizes(XX, ysubnlocal, n)); in EllipticInitialize()
670 PetscCall(VecSetFromOptions(XX)); in EllipticInitialize()
673 PetscCall(VecDuplicate(XX, &YY)); in EllipticInitialize()
674 PetscCall(VecDuplicate(XX, &ZZ)); in EllipticInitialize()
675 PetscCall(VecDuplicate(XX, &XXwork)); in EllipticInitialize()
676 PetscCall(VecDuplicate(XX, &YYwork)); in EllipticInitialize()
677 PetscCall(VecDuplicate(XX, &ZZwork)); in EllipticInitialize()
678 PetscCall(VecDuplicate(XX, &UTwork)); in EllipticInitialize()
[all …]
|
| H A D | hyperbolic.c | 676 Vec XX, YY, XXwork, YYwork, yi, uxi, ui, bc; in HyperbolicInitialize() local
790 PetscCall(VecCreate(PETSC_COMM_WORLD, &XX)); in HyperbolicInitialize()
792 PetscCall(VecSetSizes(XX, PETSC_DECIDE, n)); in HyperbolicInitialize()
794 PetscCall(VecSetFromOptions(XX)); in HyperbolicInitialize()
797 PetscCall(VecDuplicate(XX, &YY)); in HyperbolicInitialize()
798 PetscCall(VecDuplicate(XX, &XXwork)); in HyperbolicInitialize()
799 PetscCall(VecDuplicate(XX, &YYwork)); in HyperbolicInitialize()
800 PetscCall(VecDuplicate(XX, &user->d)); in HyperbolicInitialize()
801 PetscCall(VecDuplicate(XX, &user->dwork)); in HyperbolicInitialize()
803 PetscCall(VecGetOwnershipRange(XX, &istart, &iend)); in HyperbolicInitialize()
[all …]
|
| H A D | parabolic.c | 659 Vec XX, YY, ZZ, XXwork, YYwork, ZZwork, UTwork, yi, di, bc; in ParabolicInitialize() local
698 PetscCall(VecCreate(PETSC_COMM_WORLD, &XX)); in ParabolicInitialize()
700 PetscCall(VecSetSizes(XX, PETSC_DECIDE, n)); in ParabolicInitialize()
702 PetscCall(VecSetFromOptions(XX)); in ParabolicInitialize()
705 PetscCall(VecDuplicate(XX, &YY)); in ParabolicInitialize()
706 PetscCall(VecDuplicate(XX, &ZZ)); in ParabolicInitialize()
707 PetscCall(VecDuplicate(XX, &XXwork)); in ParabolicInitialize()
708 PetscCall(VecDuplicate(XX, &YYwork)); in ParabolicInitialize()
709 PetscCall(VecDuplicate(XX, &ZZwork)); in ParabolicInitialize()
710 PetscCall(VecDuplicate(XX, &UTwork)); in ParabolicInitialize()
[all …]
|
| /petsc/src/ksp/pc/impls/bjacobi/bjkokkos/ |
| H A D | bjkokkos.kokkos.cxx | 101 PetscScalar *XX = ptr; in BJSolve_TFQMR() local
161 XX[rowb - start] = 0; in BJSolve_TFQMR()
246 …parallel_for(Kokkos::TeamVectorRange(team, Nblk), [=](int idx) { XX[idx] = XX[idx] + eta * D[idx];… in BJSolve_TFQMR()
311 …parallel_for(Kokkos::TeamVectorRange(team, Nblk), [=](int idx) { XX[idx] = Diag[idx] * XX[idx]; }); in BJSolve_TFQMR()
316 glb_x[rowa] = XX[rowb - start]; in BJSolve_TFQMR()
343 PetscScalar *XX = ptr; in BJSolve_BICG() local
386 XX[rowb - start] = 0; in BJSolve_BICG()
460 XX[idx] = XX[idx] + a * Pr[idx]; in BJSolve_BICG()
512 glb_x[rowa] = XX[rowb - start]; in BJSolve_BICG()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex49.c | 985 Vec ff, XX; in solve_elasticity_2d() local
990 PetscCall(VecDuplicate(ff, &XX)); in solve_elasticity_2d()
995 PetscCall(KSPSolve(ksp_E, ff, XX)); in solve_elasticity_2d()
1000 PetscCall(VecScatterCreate(XX, NULL, X, is, &scat)); in solve_elasticity_2d()
1001 PetscCall(VecScatterBegin(scat, XX, X, INSERT_VALUES, SCATTER_FORWARD)); in solve_elasticity_2d()
1002 PetscCall(VecScatterEnd(scat, XX, X, INSERT_VALUES, SCATTER_FORWARD)); in solve_elasticity_2d()
1007 PetscCall(VecDestroy(&XX)); in solve_elasticity_2d()
|
| /petsc/doc/manual/ |
| H A D | vec.md | 570 DMDAVecGetKokkosOffsetView(DM dm, Vec vec, Kokkos::View<const PetscScalar*XX*, MemorySpace> *ov)
573 where `*XX*` can contain any number of `*`. This allows one to write very natural Kokkos multi-dime…
|