| /petsc/src/dm/tests/ |
| H A D | ex10k.kokkos.cxx | 10 using PetscScalarKokkosOffsetView2D = Kokkos::Experimental::OffsetView<PetscScalar **, Kokkos:… 11 using ConstPetscScalarKokkosOffsetView2D = Kokkos::Experimental::OffsetView<const PetscScalar **, K… 13 using PetscScalarKokkosOffsetView3D = Kokkos::Experimental::OffsetView<PetscScalar ***, Kokkos… 14 using ConstPetscScalarKokkosOffsetView3D = Kokkos::Experimental::OffsetView<const PetscScalar ***, … 24 using NodeKokkosOffsetView2D = Kokkos::Experimental::OffsetView<Node **, Kokkos::LayoutRight, … 25 using ConstNodeKokkosOffsetView2D = Kokkos::Experimental::OffsetView<const Node **, Kokkos::LayoutR… 41 Show how to manage a multi-component DMDA with DMDAVecGetKokkosOffsetViewDOF in main() 84 …nt c = 0; c < dof; c++) garray[j][i][c] = (larray[j][i - 1][c] + larray[j][i + 1][c] + larray[j - … in main() 95 … PetscInt i, PetscInt c) { gview(j, i, c) = (lview(j, i - 1, c) + lview(j, i + 1, c) + lview(j - 1… in main() 100 PetscCall(VecAXPY(g, -1.0, gg)); in main() [all …]
|
| H A D | ex2k.kokkos.cxx | 4 …On a machine with AMD EPYC-7452 CPUs, we got this data using one MPI rank and a serial-only Kokkos: 6 ------------------------------------------ 8 ------------------------------------------ 9 32 4.6464E-05 4.7451E-05 1.6880E-04 10 64 2.5654E-04 2.5164E-04 5.6780E-04 11 128 1.9383E-03 1.8878E-03 4.7938E-03 12 256 1.4450E-02 1.3619E-02 3.7778E-02 13 512 1.1580E-01 1.1551E-01 2.8428E-01 16 Overall, C is -2% ~ 5% faster than PETSc. But Kokkos is 1.6~3.6x slower than PETSc 28 using PetscScalarKokkosOffsetView3D = Kokkos::Experimental::OffsetView<PetscScalar ***, Kokkos… [all …]
|
| /petsc/doc/changes/ |
| H A D | 315.md | 6 - Add `PetscKokkosInitializeCheck()`, which initializes Kokkos if it 8 - Add support for `-debug_terminal` Terminal to use Apple's Terminal 9 instead of xterm, allowing use of cut-paste 10 - Make Terminal the default device to display the debugger on Apple 12 - Add PetscHasExternalPackage() determining whether PETSc has been 18 - On macOS, `MACOS_FIREWALL=1` or `MACOS_FIREWALL_REFRESH=1` can 21 `make -f gmakefile.test help` for details 23 - `./configure --with-macos-firewall-rules` makes 26 - Change `--download-petsc4py` to `--with-petsc4py` to have PETSc build 29 - Add `--download-mmg` and `--download-parmmg`, 3D unstructured mesh [all …]
|
| H A D | 317.md | 13 - PETSc now requires a C99 compliant C compiler in all cases. Previously C99 was only required when… 14 - PETSc now requires a C++11 compliant C++ compiler. Note this requirement is only enforced if C++ … 15 - PETSc now requires at least Microsoft Visual Studio 2015 when using the Microsoft Visual C/C++ Co… 16 - Change `MPIU_Allreduce()` to always returns a PETSc error code that should be checked with `CHKER… 21 - Change minimum value of `--with-cxx-dialect` argument from "03" to "11" 22 - C++ dialect will now also be inferred from compiler flags, although users will be warned that the… 23 - Change C++ dialect flag option to be consistent with compiler flags; `--with-cxx-dialect=gnu++14`… 24 - Fix for requesting no C++ dialect flag via `--with-cxx-dialect=0`. Previously `configure` would b… 25 - Fix a number of corner-cases when handling C++ dialect detection 26 - Remove deprecated `PETSC_VERSION_PATCH` so as to not have confusion with patch releases where the… [all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex3k.kokkos.cxx | 7 User-defined application context 11 Vec F; /* right-hand side of PDE */ 15 /* ------------------------------------------------------------------- */ 17 FormInitialGuess - Computes initial guess. 20 . x - the solution vector 31 /* ------------------------------------------------------------------- */ 33 CpuFunction - Evaluates nonlinear function, F(x) on CPU 36 . snes - the SNES context 37 . x - input vector 38 . ctx - optional user-defined context, as set by SNESSetFunction() [all …]
|
| H A D | ex55k.kokkos.cxx | 14 using ConstPetscScalarKokkosOffsetView2D = Kokkos::Experimental::OffsetView<const PetscScalar **, K… 15 using PetscScalarKokkosOffsetView2D = Kokkos::Experimental::OffsetView<PetscScalar **, Kokkos:… 26 PetscReal x = PetscRealPart(c->x), y = PetscRealPart(c->y); in MMSSolution1() 27 u[0] = x * (1 - x) * y * (1 - y); in MMSSolution1() 33 PetscReal x = PetscRealPart(c->x), y = PetscRealPart(c->y); in MMSForcing1() 34 f[0] = 2 * x * (1 - x) + 2 * y * (1 - y) - user_param * PetscExpReal(x * (1 - x) * y * (1 - y)); in MMSForcing1() 41 …PetscInt xs = info->xs, ys = info->ys, xm = info->xm, ym = info->ym, mx = info->mx, my = info->my; in FormFunctionLocalVec() 42 PetscReal user_param = user->param; in FormFunctionLocalVec() 49 lambda = user->param; in FormFunctionLocalVec() 50 hx = 1.0 / (PetscReal)(info->mx - 1); in FormFunctionLocalVec() [all …]
|
| H A D | ex30.c | 1 static const char help[] = "Steady-state 2D subduction flow, pressure and temperature solver.\n\ 3 ---------------------------------ex30 help---------------------------------\n\ 4 -OPTION <DEFAULT> = (UNITS) DESCRIPTION.\n\n\ 5 -width <320> = (km) width of domain.\n\ 6 -depth <300> = (km) depth of domain.\n\ 7 -slab_dip <45> = (degrees) dip angle of the slab (determines the grid aspect ratio).\n\ 8 -lid_depth <35> = (km) depth of the static conductive lid.\n\ 9 -fault_depth <35> = (km) depth of slab-wedge mechanical coupling\n\ 12 -ni <82> = grid cells in x-direction. (nj adjusts to accommodate\n\ 13 the slab dip & depth). DO NOT USE -da_grid_x option!!!\n\ [all …]
|
| /petsc/ |
| H A D | .gitlab-ci.yml | 2 # stage-1 take only a few minutes; they do not run the full test suite or external packages. 4 # stage-2 runs on MCS systems and may take 10 to 15 minutes. They run the full test suite but with … 6 # stage-3 runs on MCS systems and may take an hour or more. They run the full test suite and heavil… 8 # The stage-(n) tests are only started if all of the stage-(n-1) tests run without error 13 - stage-1 14 - stage-2 15 - stage-3 16 - stage-4 19 GIT_CLEAN_FLAGS: -ffdxq 21 BASE_EXTRA_OPTIONS: -nox -nox_warning -malloc_dump [all …]
|
| /petsc/doc/community/meetings/2025/ |
| H A D | index.md | 1 --- 3 --- 11 May 20-21, 2025, 101 Davis Hall, University of Buffalo, NY, USA ([105 White Rd, Amherst, NY 14260](… 15 - Monday, May 19 - Tutorial (tutorials begin at 9am) 16 - Tuesday, May 20 - Meeting (begin at 9am) 17 - Wednesday, May 21 - Meeting (ends around 5pm) 28 …-------- | ---------------------------------------------------------------------------------------… 39 …-------- | ---------------------------------------------------------------------------------------… 42 | 9:30 am | Dense Broyden-Fletcher-Goldfarb-Shanno (BFGS) … 43 …-Boundary Adaptive Mesh Refinement … [all …]
|
| /petsc/src/ksp/ksp/impls/gmres/dgmres/ |
| H A D | dgmres.c | 53 KSP_DGMRES *dgmres = (KSP_DGMRES *)ksp->data; in KSPSetUp_DGMRES() 54 PetscInt neig = dgmres->neig + EIG_OFFSET; in KSPSetUp_DGMRES() 55 PetscInt max_k = dgmres->max_k + 1; in KSPSetUp_DGMRES() 59 if (!dgmres->neig) PetscFunctionReturn(PETSC_SUCCESS); in KSPSetUp_DGMRES() 63 dgmres->wr = NULL; in KSPSetUp_DGMRES() 64 dgmres->wi = NULL; in KSPSetUp_DGMRES() 65 dgmres->perm = NULL; in KSPSetUp_DGMRES() 66 dgmres->modul = NULL; in KSPSetUp_DGMRES() 67 dgmres->Q = NULL; in KSPSetUp_DGMRES() 68 dgmres->Z = NULL; in KSPSetUp_DGMRES() [all …]
|
| /petsc/src/ksp/ksp/impls/gmres/agmres/ |
| H A D | agmres.c | 27 KSP_AGMRES *agmres = (KSP_AGMRES *)ksp->data; in KSPSetUp_AGMRES() 28 PetscInt neig = agmres->neig; in KSPSetUp_AGMRES() 29 const PetscInt max_k = agmres->max_k; in KSPSetUp_AGMRES() 31 PetscInt lwork = PetscMax(8 * N + 16, 4 * neig * (N - neig)); in KSPSetUp_AGMRES() 36 agmres->q_preallocate = PETSC_TRUE; /* No allocation on the fly */ in KSPSetUp_AGMRES() 38 ksp->calc_sings = PETSC_TRUE; in KSPSetUp_AGMRES() 39 agmres->max_k = N; /* Set the augmented size to be allocated in KSPSetup_GMRES */ in KSPSetUp_AGMRES() 41 agmres->max_k = max_k; in KSPSetUp_AGMRES() 45 …PetscCall(PetscCalloc4(max_k, &agmres->Rshift, max_k, &agmres->Ishift, hes, &agmres->Rloc, (N + 1)… in KSPSetUp_AGMRES() 46 PetscCall(PetscMalloc3(N + 1, &agmres->tau, lwork, &agmres->work, N + 1, &agmres->nrs)); in KSPSetUp_AGMRES() [all …]
|
| /petsc/src/ksp/ksp/impls/minres/ |
| H A D | minres.c | 8 … " title={MINRES-QLP: A Krylov subspace method for indefinite or singular symmetric systems},\n" 9 … " author={Choi, Sou-Cheng T and Paige, Christopher C and Saunders, Michael A},\n" 13 " pages={1810--1836},\n" 25 // The following arrays are of size ksp->maxit 38 if (ksp->calc_sings) { in KSPSetUp_MINRES() 39 KSP_MINRES *minres = (KSP_MINRES *)ksp->data; in KSPSetUp_MINRES() 40 PetscInt maxit = ksp->max_it; in KSPSetUp_MINRES() 41 PetscCall(PetscFree4(minres->e, minres->d, minres->ee, minres->dd)); in KSPSetUp_MINRES() 42 …PetscCall(PetscMalloc4(maxit + 1, &minres->e, maxit, &minres->d, maxit, &minres->ee, maxit, &minre… in KSPSetUp_MINRES() 44 ksp->ops->computeextremesingularvalues = KSPComputeExtremeSingularValues_MINRES; in KSPSetUp_MINRES() [all …]
|
| /petsc/config/BuildSystem/config/packages/ |
| H A D | MPI.py | 45 # defaults to --with-mpi=yes 54 # support MPI-3 process shared memory 68 …nt('MPI', '-with-mpiexec=<prog>', nargs.Arg(None, None, 'The utility … 69 …t('MPI', '-with-mpiexec-tail=<prog>', nargs.Arg(None, None, 'The utility y… 70 …help.addArgument('MPI', '-with-mpi-compilers=<bool>', nargs.ArgBool(None, 1… 71 …t('MPI', '-known-mpi-shared-libraries=<bool>', nargs.ArgBool(None, None, 'Indicates… 72 …help.addArgument('MPI', '-with-mpi-ftn-module=<mpi or mpi_f08>', nargs.ArgStr… 106 raise RuntimeError('--download-mpich libraries cannot be used') 109 raise RuntimeError('--download-openmpi libraries cannot be used') 113 dirExp = re.compile(r'mpi(ch)?(-.*)?') [all …]
|
| /petsc/doc/community/meetings/2023/ |
| H A D | index.md | 1 --- 3 --- 9 ```{image} https://petsc.gitlab.io/annual-meetings/2023/GroupPhoto.jpg 14 June 5-7, 2023, at the [Hermann Hall Conference Center](https://www.iit.edu/event-services/meeting-… 19 [Parking use B5 (32nd & Federal St.)](https://www.iit.edu/cbsc/parking/visitor-and-event-parking). 21 Please test for Covid before attending the meeting and 24 … particular simulations. Finally, the meeting will be interspersed with mini-tutorials that will d… 28 - Monday, June 5: 1 pm to 5:30 pm 29 - Tuesday, June 6: 10:15 am to 5:30 pm 30 - Wednesday, June 7: 9 am to 3 pm [all …]
|
| /petsc/src/tao/leastsquares/impls/pounders/ |
| H A D | pounders.c | 16 PetscCall(MatMult(mfqP->subH, x, g)); in pounders_fg() 20 PetscCall(VecDot(mfqP->subb, x, &d2)); in pounders_fg() 24 PetscCall(VecAXPY(g, 1.0, mfqP->subb)); in pounders_fg() 30 TAO_POUNDERS *mfqP = (TAO_POUNDERS *)tao->data; in pounders_feval() 36 if (tao->res_weights_v) { in pounders_feval() 37 PetscCall(VecPointwiseMult(mfqP->workfvec, tao->res_weights_v, F)); in pounders_feval() 38 PetscCall(VecDot(mfqP->workfvec, mfqP->workfvec, fsum)); in pounders_feval() 39 } else if (tao->res_weights_w) { in pounders_feval() 41 for (i = 0; i < tao->res_weights_n; i++) { in pounders_feval() 42 row = tao->res_weights_rows[i]; in pounders_feval() [all …]
|
| /petsc/doc/ |
| H A D | petsc.bib | 3 % bibtool petsc.bib -- expand.macros=on -- print.line.length=100 -- pass.comments=on 17 % LiteralHTML: <a name="nano"><H3><center>Nano-simulations</center></H3> 19 @Misc{ semver-webpage, 25 title = {{D}-stability and {K}aps-{R}entrop-methods}, 30 pages = {229--237}, 41 pages = {93--113}, 55 …title = {Generalized {R}unge--{K}utta methods of order four with stepsize control for sti… 59 pages = {55--68}, 73 …gorithm for Structural Dynamics With Improved Numerical Dissipation: The Generalized-alpha Method}, 77 pages = {371-375}, [all …]
|
| /petsc/src/ksp/pc/impls/bddc/ |
| H A D | bddcprivate.c | 68 PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, nr, nr - i, NULL, B)); in MatDenseOrthogonalRangeOrComplement() 70 PetscCall(PetscArraycpy(data, U + nr * i, (nr - i) * nr)); in MatDenseOrthogonalRangeOrComplement() 155 PC_BDDC *pcbddc = (PC_BDDC *)pc->data; in PCBDDCNedelecSupport() 156 Mat_IS *matis = (Mat_IS *)pc->pmat->data; in PCBDDCNedelecSupport() 180 order = pcbddc->nedorder; in PCBDDCNedelecSupport() 181 conforming = pcbddc->conforming; in PCBDDCNedelecSupport() 182 field = pcbddc->nedfield; in PCBDDCNedelecSupport() 183 global = pcbddc->nedglobal; in PCBDDCNedelecSupport() 188 …PetscOptionsBegin(PetscObjectComm((PetscObject)pc), ((PetscObject)pc)->prefix, "BDDC Nedelec optio… in PCBDDCNedelecSupport() 189 …PetscCall(PetscOptionsBool("-pc_bddc_nedelec_field_primal", "All edge dofs set as primals: Toselli… in PCBDDCNedelecSupport() [all …]
|