Sys: add PetscMalloc1 macro, array allocation without redundant types

The type is inferred from the pointer return type. This patch is
automated via the following script:

git grep -l 'PetscMalloc(

Sys: add PetscMalloc1 macro, array allocation without redundant types

The type is inferred from the pointer return type. This patch is
automated via the following script:

git grep -l 'PetscMalloc(.*sizeof' src | xargs perl -pi -e 's@PetscMalloc\(([^,;]*[^,; ]) *\* *sizeof\([^,;()]+\),@PetscMalloc1($1,@'

This commit contains an additional bug-fix in csrperm.c, fixing pointer
arity. The code was introduced in 2006, but the allocation could not
have been correct at any time. This probably means that
MatDuplicate_SeqAIJPERM has never been tested.

a54129beb540034ba105796c682d589e7e1111f2
Richard Tran Mills <rmills@ornl.gov>

Added MATSEQCSRPERM support for MatDuplicate() and conversion to/from
MATSEQAIJ. Note that these changes are not quite debugged.

show more ...

Merge branch 'barry/minor-draw-improvements'

added -ts_monitor_envelope and TSMonitorEnvelope() to track extremes of solution vector
added support to extchem.c to print maximums of species
updated sample run for extchem.c for displaying only "i

added -ts_monitor_envelope and TSMonitorEnvelope() to track extremes of solution vector
added support to extchem.c to print maximums of species
updated sample run for extchem.c for displaying only "important species"

show more ...

added TSMonitorLGSetTransform() and added transform to mole fraction in extchem.c

added -ts_lg_monitor_solution_variables and TSMonitorLGSetDisplayVariables()

changed TSMonitorLGSolutionSetLegend() to TSMonitorLGSetVariableNames()

add support for providing legends for TSMonitorLGSolution()

Merge branch 'barry/minor-draw-improvements' into jed/chem

added -draw_save_final_image and fixed -lg_indicate_data_points
added PetscDrawLGSetFromOptions()

introduce XXXViewFromOptions() and use consistently when possible

Merge branch 'knepley/example-ts-friction'

* knepley/example-ts-friction:
TS ex34: Now reproducing plots from Berkeley paper - Fixed coordinate DM screwups
DMDA: Fixed bug with 1D DMDAGetRay()

Merge branch 'knepley/example-ts-friction'

* knepley/example-ts-friction:
TS ex34: Now reproducing plots from Berkeley paper - Fixed coordinate DM screwups
DMDA: Fixed bug with 1D DMDAGetRay()
TS: Allow the DMDARay viewer to operator in LG mode - Added TSMonitorLGDMDARay() activated by -ts_monitor_lg_dmda_ray - Added TSMonitorLGCtx to TSMonitorDMDARayCtx
DMDA: Allow DMDAGetRay() to work in 1D
TS: Make TSMonitorLGCtx visible to other code
TS: Added an example of DR friction

show more ...

Merge branch 'master' into knepley/feature-dmda-section

* master: (528 commits)
Doc: Fixed configure doc
Nightly summary: Added configure section.
TSSetErrorIfStepFails: fix __FUNCT__
TS: fi

Merge branch 'master' into knepley/feature-dmda-section

* master: (528 commits)
Doc: Fixed configure doc
Nightly summary: Added configure section.
TSSetErrorIfStepFails: fix __FUNCT__
TS: fix for TSGetRHSMats_Private() in the creation of the Arhs and Brhs matrices
protect logging of objects (which uses global variables) with PETSC_USE_LOG
error handler API had been previously changed but this example was not updated to new API.
DMPlex: Fixed bug in coordinate permutation
PetscSection: Doc fix for permutation
Bib: Added Top500 and fixed entry
fun3d: update PetscMallocValidate() usage
bugfix of MatCreateSeqAIJFromTriple() contributed by Nicolas Marsic <nicolas.marsic@montefiore.ulg.ac.be>
bib: rename ref with duplicate key
Webpage: Corrected '-dm_mat_type cusp' to 'dm_mat_type aijcusp'.
completed removing server code that was not removed in previous commit
removed PETSc webserver and configure package files that supported it
Allow calling MatGetBlockSize[s]() before matrix preallocation
SNESLINESEARCHBT: Set the norms when exiting early due to negligible step.
fix outdated PetscOptionsList
SNESQN: only monitor real part of dot product (fails with C++ complex)
Fix for multiple definitions in mat/color/interface/ custom fortran bindings
...

Conflicts:
config/builder.py
include/petscdmplex.h
src/dm/impls/da/dageometry.c
src/dm/impls/plex/examples/tests/ex4.c
src/dm/impls/plex/plex.c
src/dm/impls/plex/plexinterpolate.c

show more ...

Merge branch 'maint'

* maint:
TSSetErrorIfStepFails: fix __FUNCT__

TSSetErrorIfStepFails: fix __FUNCT__

Reported-by: Richard Katz <richard.katz@earth.ox.ac.uk>

Merge branch 'jed/fix-ts-rhs-mats'

* jed/fix-ts-rhs-mats:
TS: fix for TSGetRHSMats_Private() in the creation of the Arhs and Brhs matrices

TS: fix for TSGetRHSMats_Private() in the creation of the Arhs and Brhs matrices

When an RHS Jacobian is provided, matrices must be copied or duplicated
for Arhs and Brhs. Commit e1244c694be870577d

TS: fix for TSGetRHSMats_Private() in the creation of the Arhs and Brhs matrices

When an RHS Jacobian is provided, matrices must be copied or duplicated
for Arhs and Brhs. Commit e1244c694be870577ddfeb27d69a78d09ad24070

TS: add TSRHSJacobianSetReuse, fixes TSComputeRHSJacobianConstant

changed the logic for how to recompute RHS Jacobians, which had the
side-effect of causing TSGetRHSMats_Private to duplicate Brhs even if A
was being reused and A == B. This commit fixes the logic to reuse Arhs
as Brhs when they should match.

This commit cherry-picks the change to src/ts/interface/ts.c from
3162be4ab4fd0bc00dad2f51f4d64a4f9f8a384e.

Reported-by: Jay Jay Billings <jayjaybillings@gmail.com>

show more ...

Merge branch 'barry/saws'

Merge branch 'master' into knepley/feature-fem-dgspace

* master: (593 commits)
Bib: Added Top500 and fixed entry
fun3d: update PetscMallocValidate() usage
bib: rename ref with duplicate key

Merge branch 'master' into knepley/feature-fem-dgspace

* master: (593 commits)
Bib: Added Top500 and fixed entry
fun3d: update PetscMallocValidate() usage
bib: rename ref with duplicate key
Webpage: Corrected '-dm_mat_type cusp' to 'dm_mat_type aijcusp'.
Allow calling MatGetBlockSize[s]() before matrix preallocation
SNESLINESEARCHBT: Set the norms when exiting early due to negligible step.
SNESQN: only monitor real part of dot product (fails with C++ complex)
Corrected the #include statement in the man pages for the PetscObjectComposedDataSet/Getwhatever routines.
DMPlex: Added doc for DMPlexGetHybridBounds()
winzip: detect if winzip is used to extract petsc.tar.gz and error out.
configure: remove self.archIndependent as its not being used. [it is used to create externalpackages/package/arch dir which is unused]
configure: fix packages using package-dir/PETSC_ARCH as out-of-source-build location Perhaps this should be changed to package-dir/build?
configure: remove dead code previously used to download BuildSystem
configure: check if compilerDefines [and compilerFixes] exist in framework before using them
configure: save/restore reconfigure.py when --with-clean is used
configure: add --with-clean option to delete buildfiles/externalpackages in PETSC_DIR/PETSC_ARCH and potential externalpackages in --with-external-package-dir/PETSC_ARCH
configure: With --with-externalpackages-dir=dir store/build packages in dir/PETSC_ARCH
Add MatSeqSBAIJSetPreallocationCSR()
configure: switch to using PETSC_DIR/PETSC_ARCH/externalpackages by updating package.py to use externalpackages.py as the externalPackagesDirProvider
configure: move dead code configureExternalPackagesDir() from petscdir.py to externalpackagesdir.py Also set PETSC_DIR/PETSC_ARCH/externalpackages as the externalpackagesdir
...

show more ...

Merge branch 'master' into barry/saws

Conflicts:
src/ksp/pc/impls/gamg/gamg.c
src/sys/classes/viewer/impls/ams/ams.c
src/sys/objects/pinit.c

TS: Allow the DMDARay viewer to operator in LG mode
- Added TSMonitorLGDMDARay() activated by -ts_monitor_lg_dmda_ray
- Added TSMonitorLGCtx to TSMonitorDMDARayCtx

TS: Make TSMonitorLGCtx visible to other code

Merge branch 'emil/ts-poststage'

* emil/ts-poststage:
Added a post-stage function for time-integrators

Added a post-stage function for time-integrators

Merge branch 'master' into knepley/feature-dmda-section

* master: (397 commits)
PetscSynchronizedFGets: fix deadlock at EOF
Compiler: Fix warnings from MPI impls which do not initialize outputs

Merge branch 'master' into knepley/feature-dmda-section

* master: (397 commits)
PetscSynchronizedFGets: fix deadlock at EOF
Compiler: Fix warnings from MPI impls which do not initialize outputs
Increase patchlevel to 3.4.3
SNES: Now work vectors come from the DM
SNES ex62: Remove code generation
SNES: Move setup involving snes->vec_sol from SNESSetUp() to SNESolve()
Compiler: Fix warnings from MPI impls which do not initialize outputs
DMPlex ex7: Added missing test output
SNES ex52: Removed old Jacobian stuff, and fixed call to DMPlexProjectFunction() - Fixed calls for new element handling
SNES ex12: Added a performance profiling mode
DMPlex ex8: Fixed leak
PetscSection: Must reset the section when changing the number of fields
DMPlex ex7: Fixed test output - Corrected orientations in interpolation
DMPlex ex1: Fixed test output - Uniform refinement for quads changed
DMPlex: Fix overagressive checks
PC: Removed support graph PC
configure: both downloadonWindows and worksonWindows refer to MS compilers. Also fix isWindows() -> isWindows(CC).
MatXAIJSetPreallocation: use array[] notation to help out Fortran
__float128: when blas/lapack is not found instruct using f2cblaslapack
Configure: Package.downloadonWindows is supposed to mean that Windows compilers work. not Cygwin - Also changed doc for Package.worksonWindows, which does refer to Cygwin
...

Conflicts:
src/dm/impls/da/dalocal.c

show more ...

Merge branch 'master' into knepley/feature-plex-refine-3d

* master: (149 commits)
Bib: Added papers referencing PETSc from ResearchGate
Bib: Fixed urls for some refs Suggested-by: Michele De Ste

Merge branch 'master' into knepley/feature-plex-refine-3d

* master: (149 commits)
Bib: Added papers referencing PETSc from ResearchGate
Bib: Fixed urls for some refs Suggested-by: Michele De Stefano
fixed bug that has appeared recently of false positive errors when there are zero local equations on coarse grids in gamg-agg.
MatSchurComplement: Must initialize KSP package
bib: fix inbook citation
Bib: Updated some refs
SNES: fix copy/paste in man page
fortran: fix KSPCONVERGEDDEFAULT and SNESCONVERGEDSKIP in fortran includes
Mat SBAIJ: avoid invalid memory read in MatSOR_SBAIJ
Revert "MatMPIBAIJSetPreallocationCSR() handled the input matrix entries incorrectly, it treated them as if each row of values"
PCGAMG: fix corruption for multiple solves when MatSetNearNullSpace is used
PCGAMG: remove distracting comment
KSP ex49: restrict near-null space for elasticity example with sym BC
MatMPIBAIJSetPreallocationCSR() handled the input matrix entries incorrectly, it treated them as if each row of values was provided ordered by row within the block rows instead of by block as it was originally designed to be. Also added more discussion to the manual pages on the expected order of the values
Fix accidental creation of the matrices in SNESNCG used for linear analogue LS.
PCGAMG: add man page for PCGAMGType
fortran: fix snesconvergedskip_() kspconvergedskip_() kspconvergeddefault*_()
rm //
Sys: new_ir changed from PetscInt to PetscObjectState - so fix mallocs to use the new type
SNES: Fixed package init in SNESLineSearchCreate()
...

Conflicts:
include/petscdmplex.h
src/dm/impls/plex/examples/tests/ex5.c
src/dm/impls/plex/plex.c
src/dm/interface/dm.c

show more ...