SF: add an MPI_Op argument to SFBcast

git: switch default branch from 'master' to 'main'

CHKERRQ() -> CHKERRMPI()

KSP: Reorganized monitors - Now the residual can be viewed using -ksp_monitor draw and a plot can be made using -ksp_monitor draw::draw_lg - This also works for the true residual and error - We no

KSP: Reorganized monitors - Now the residual can be viewed using -ksp_monitor draw and a plot can be made using -ksp_monitor draw::draw_lg - This also works for the true residual and error - We noe register monitors so that dispatch based on viewer and format works - Add KSPMonitorLGCreate(), and remove KSPMonitorLGResidualNorm*() and KSPMonitorLGTrueResidualNorm*() - Add KSPMonitorError(), used by -ksp_monitor_error - Add arguments to KSPMonitorSetFromOptions() to allow line graphs to be configured - Remove SNESMonitorLGCreate() and SNESMonitorLGResidualNorm() which are now handled by the default monitor - PetscViewerAndFormat now allows a payload - VecView_Plex() and VecView_Plex_Local() must propagate the flag indicating zero BC

show more ...

Merge branch 'jed/pcmg-transpose' into 'master'

PCMG: implement PCApplyTranspose_MG

See merge request petsc/petsc!3539

Add KSPSetUseExplicitTranspose

Added -ksp_use_explicittranspose to explicitly transpose the system in KSPSolveTranspose

PCMG: implement PCApplyTranspose_MG

A sample run:

$PETSC_ARCH/tests/ts/tutorials/advection-diffusion-reaction/ex5adj \
-pc_type mg -mg_levels_pc_type jacobi

SNES ex76: Add monitor for cell divergence

Merge branch 'hongzh/support-parameterized-mass-matrix' into 'master'

Add sensitivity suport for parameterized mass matrice

See merge request petsc/petsc!3500

Merge remote-tracking branch 'origin/release'

Add sensitivity suport for parameterized mass matrice

- Xdot now is calculated correctly before passing to TSComputeIJacobianP
- Always checkpoint th->X0 as the stage value for Theta methods. This a

Add sensitivity suport for parameterized mass matrice

- Xdot now is calculated correctly before passing to TSComputeIJacobianP
- Always checkpoint th->X0 as the stage value for Theta methods. This allows us to recover th->X (for endpoint!=1) from th->X0 and to calculate th->Xdot.

Commit-type: feature
/spend 4.0 hr
Thanks-to: salazardetro1@llnl.gov

show more ...

Merge branch 'emil/fix-description-ex25' into 'release'

modified description of ex25 to avoid IMEX confusion.

See merge request petsc/petsc!3109

Merge branch 'hzhang/dmnetwork-netcoupling' into 'master'

Hzhang/dmnetwork netcoupling

See merge request petsc/petsc!3051

New design for the coupling of sub-dmnetowrks.
See email discussions on Nov 23, 2019, Subject: Re: petsc | DMNetwork: Update manual. (!2297)
See dev.html for changes made
using sf, DMNetworkVariables

New design for the coupling of sub-dmnetowrks.
See email discussions on Nov 23, 2019, Subject: Re: petsc | DMNetwork: Update manual. (!2297)
See dev.html for changes made
using sf, DMNetworkVariablesSetUp() handles shared vertex easily

show more ...

TS ex77: Preliminary reactive flow simulation - The particles are passive right now, so we can check the accuracy of particle integration

Insure that PetscErrorPrintfHilight() does not use special charactors while running
test harness by skipping highlighting if PETSC_STDERR is PETSC_STDOUT.

Remove sed handling of special hightlightin

Insure that PetscErrorPrintfHilight() does not use special charactors while running
test harness by skipping highlighting if PETSC_STDERR is PETSC_STDOUT.

Remove sed handling of special hightlighting in output from PetscErrorPrintfHilight() because
it was not portable to all systems. And it would result say two different files where the same.

Remove detailed printing from ts/tutorials/advection-diffusion-reaction/ex3.c output because each
machine uses a slighly different number of function evaluations, etc.

Commit-type: portability-fix, testing-fix
/spend 3.5h

show more ...

Convert MPI error type to PETSc error with string message for all MPI calls

Now PETSc examples will ONLY return PETSc error codes and never MPI error codes directly so we can understand and post-pro

Convert MPI error type to PETSc error with string message for all MPI calls

Now PETSc examples will ONLY return PETSc error codes and never MPI error codes directly so we can understand and post-process their
errors better.

The test harness will now automatically retry tests that fail with MPI, this may help with Intel MPI that produces
seemingly random failures.

Commit-type: error-checking
/spend 30m

show more ...

Merge branch 'barry/2020-09-07/matstructure-unknown' into 'master'

Barry/2020 09 07/matstructure unknown

See merge request petsc/petsc!3182

TS ex76: Added Taylor-Green vortex test

Remove static on function prototypes that didn't belong.

Commit-type: bug-fix
/spend 5m
Reported-by: GitLab CI

Rewrite two TS examples in standard form u_t - F() instead of F() - u_t
Fix the two examples to use ARKIMEX methods with full implicit stages since the mass matrix for the two examples is not the ide

Rewrite two TS examples in standard form u_t - F() instead of F() - u_t
Fix the two examples to use ARKIMEX methods with full implicit stages since the mass matrix for the two examples is not the identity so the explicit stages will not work
Add equations for ts ex36 including mass matrix

Commit-type: bug-fix, testing-fix, style-fix, example
/spend 2h
Reported-by: Gitlab CI

Commit-type: example
/spend 5m
Thanks-to: Emil Constantinescu <emconsta@anl.gov>

show more ...

Fix code formating and VecGetArrayXXX()

Commit-type: testing-fix, style-fix
/spend 15m

Sets to ARKMEX fully implicit before setting the TSType to ARKIMEX hence flag is ignored.

Flag won't work anyway because the IFunction already contains the RHSFunction and fully implicit would make

Sets to ARKMEX fully implicit before setting the TSType to ARKIMEX hence flag is ignored.

Flag won't work anyway because the IFunction already contains the RHSFunction and fully implicit would make it contain it twice

Commit-type: bug-fix
/spend 10m

show more ...

TSView now shows the number of i and rhs function and Jacobian evaluations

Commit-type: feature
/spend 30m
Reported-by: Ed Bueler
Thanks-to: Jed Brown jed@jedbrown.org

DM: Added universal label - A universal label encodes the markers for a set of labels - We need this when using systems like TetGen that only take a single label - Added tests in Plex ex11