revert BLASdotu_ change - and correctly fix the duplicate definition

fix -Wstrict-prototypes warnings

Reported-by: Till Martens <till.martens@dlr.de>

fix duplicate definitions warned by -Wredundant-decls

Also a bugfix? BLASdotu_ -> dotu

Reported-by: Till Martens <till.martens@dlr.de>

Merge branch 'barry/half-precision'

Support for half-precision using __fp16.

Requires repository trunk of gcc and ARM systems
Configure with --with-precision=__fp16 --download-f2cblaslapack --with-fc=0 --with-cxx=0
Many standard PETSc

Support for half-precision using __fp16.

Requires repository trunk of gcc and ARM systems
Configure with --with-precision=__fp16 --download-f2cblaslapack --with-fc=0 --with-cxx=0
Many standard PETSc tests fail due to lack of precision

Funded-by: ECP
Project: fp16
Time: 24 hours

show more ...

Merge branch 'stefano_zampini/feature-pcbddc-saddlepoint'

Conflicts:
src/mat/impls/aij/mpi/mpiaij.c
src/mat/impls/is/matis.c

Modified ex23.c for the API change of the function CheckMat

Merge branch 'master' into stefano_zampini/feature-pcbddc-saddlepoint

This is needed to deal with some changes in MatNullSpace from commit 9d2471e

Merge branch 'tisaac/dmplex-feature-real-to-ref'

This branch introduces DMPlexReferenceToCoordinates() and
DMPlexCoordinatesToReference() and push finite element reference element
points forward int

Merge branch 'tisaac/dmplex-feature-real-to-ref'

This branch introduces DMPlexReferenceToCoordinates() and
DMPlexCoordinatesToReference() and push finite element reference element
points forward into the mesh and inverting to pull them back.

DMPlexCoordinatesToReference() is currently implemented with a simple
Newton iteration. When the reference maps extend outside the element in
a non-unique way, a Newton iteration doesn't always do what we want.

What we really want are two properties:

- If there is a solution inside the reference element domain, converge
preferentially to that one.
- If there is no solution inside the reference element domain, diverge
outside of the domain.

Newton has no guarantee of which solution will be found, so
DMPlexCoordinatesToReference() can fail for very nonlinearly deformed
elements. Fixing this is a project for another time.

* tisaac/dmplex-feature-real-to-ref:
DMPlex tests ex22: larger tolerance, so that single precision passes
DMPlex tests ex22: use petsc portable implementation of snprintf
DMPlexComputeProject3Dto2D(): precision-aware bounds for orthogonalization success
DMPlex: fix bad declaration order in point location routine
DMPlex: complex and C++ fixes to recent point location routines
DMPlex: added point location to regression tests
DMPlex: reduce max iterations for point location.
DMPlex: indexing and work array size corrections for point location routines
DMPlex: add debugging info in point location routines
DMPlex tests ex22: fix memory leak
DMPlex: fix work array declarations in reference/coordinate transformations
DMPlex test ex22: added iso parametric test
DM: fix reference/reference order in DMSetCoordinateDM()
DMPlex tests: move my ex21 to ex22 to prepare for merge
DMPlex: implement point inversion for multilinear maps.
CoordinatesRealToRef(): correct stride
DMPlex: propose DMPlexCoordinatesToReference() interface

show more ...

DMPlex: implement point inversion for multilinear maps.

This is implemented by a fixed number of Newton iterations. In 2D, I
know that this is a contraction map for a non-degenerate element, and
wi

DMPlex: implement point inversion for multilinear maps.

This is implemented by a fixed number of Newton iterations. In 2D, I
know that this is a contraction map for a non-degenerate element, and
will always converge to the right solution. The same cannot be said in
3D. There are non-degenerate trilinearly mapped hexahedra for which any
Newton iteration, damped or undamped, can converge to a spurious
solution outside of the element, even if there is a solution inside the
element. One could go about a costly procedure of trying to determine
if the target coordinate is inside of the region defined by the mapped
boundary of the reference cell, but this is an expensive procedure. For
now, I am relying on the elements not being so deformed that Newton
doesn't converge. I spent much longer than I should have trying to come
up with a better solution to this problem. I should write the whole
thing up.

WIP: implementation of all point location routines

DMPlex: added complementary DMPlexReferenceToCoordinates() methods

For affine, this just wraps CoordinatesRefToReal.
For tensor, it uses the tensor structure.
In general, the basis space of the PetscDS is used to map points.

WIP: DMPlex: test point inversion.

It seems to work, but the test is flawed for simplicial elements

DMPlex: point location is working for implicit discretizations

Now to do the work for non-implicit (e.g. isoparametric) meshes

show more ...

BLASLAPACK: fix syr2k prototype

PETSCBLASLAPCK: added wrappers for symmetric rank updates

Merge branch 'master' into stefano_zampini/feature-pcbddc-saddlepoint

Conflicts:
src/ksp/pc/impls/bddc/bddc.c
src/ksp/pc/impls/bddc/bddcgraph.c
src/ksp/pc/impls/bddc/bddcprivate.c
src/ksp/pc/imp

Merge branch 'master' into stefano_zampini/feature-pcbddc-saddlepoint

Conflicts:
src/ksp/pc/impls/bddc/bddc.c
src/ksp/pc/impls/bddc/bddcgraph.c
src/ksp/pc/impls/bddc/bddcprivate.c
src/ksp/pc/impls/bddc/bddcscalingbasic.c
src/ksp/pc/impls/bddc/bddcschurs.c
src/mat/impls/is/matis.c

show more ...

Merge branch 'master' into barry/add-pcfieldsplit-logevent

Merge branch 'master' into barry/downloads

Merge branch 'tisaac/dmforest' of bitbucket.org:petsc/petsc into tisaac/dmforest

Merge branch 'master' into tisaac/dmforest

* master: (136 commits)
Make error message specific to exactly what matrix is supported
Prevent fblaslapack or f2cblaslapack from being compiled with t

Merge branch 'master' into tisaac/dmforest

* master: (136 commits)
Make error message specific to exactly what matrix is supported
Prevent fblaslapack or f2cblaslapack from being compiled with the -no-prec-dev Intel compiler option since it can produce incorrect results
Fix for PCReset_HYPRE to make it reentrant.
replace default superlu_dist '-mat_superlu_dist_colperm METIS_AT_PLUS_A' with MMD_ATA to avoid different random numbers for nightly tests
fixes for merge from maint patch
Added PCReset_HYPRE() so it may be used with SNESVI solvers
remove Microsoft compiler message for paste macros used for handling Blas/Lapack macro name mangling for example
Increase the number of PetscViewerPushFormat() possible to 25 and use a macro PETSCVIEWERFORMATPUSHESMAX for the value
Configure: Remove redundant logging
Bib: Added ref
nightlybuids: fix permissions of script
MatRegisterDestroy() does not exist but it is listed in manual pages
externalpackages-git: save local changes if any - in git stash - before 'git checkout -f' Useful incase --with-package-commit=HEAD option was msissing [perhaps due to a typo]
nightlybuilds: add xsdk builds
Add architectures for testing xsdk package downloads (Trilinos, Hypre, SuperLU)
create another version of TurnForward() for adaptive time stepping specifically
Configure: Give proper default if HOME is missing from the environment
fix abbreviation
allow zeropivot for MatSOR_SeqAIJ_Inode, node size=1
add test case using adpative time stepping for checkpointing type 'NONE'
...

show more ...

Merge branch 'pr394/stevengcj/blassuffix/master'

Merge branch 'master' into barry/add-pcreset-hypre-for-master

Needed to move PCDestroy_HYPRE() added to master to PCReset_HYPRE() since that was introduced
as a bug fix for maint

add --blaslapack-suffix=xxx argument for user-defined symbol suffix, consolidate Unix mangling code into a single header file for BLAS/LAPACK

Merge branch 'master' into hongzh/checkpointing

Manual merges for a couple of PetscOptions changes in master

Merge branch 'tisaac/dm-label-promote' into tisaac/dmforest

* tisaac/dm-label-promote: (319 commits)
snes tutorial ex69: give mesh a name for regression purposes
DM: make DMLabel labels a part o

Merge branch 'tisaac/dm-label-promote' into tisaac/dmforest

* tisaac/dm-label-promote: (319 commits)
snes tutorial ex69: give mesh a name for regression purposes
DM: make DMLabel labels a part of all DMs
builder: update regression tests
builder: update some missing package requirements
DM: set coordinate vector name
ts tutorials ex11_sa: make it compile
mark variable as unused so compiler doesn't warn about unused value
configure: use 'git cat-file -e' to check for existance of a commit in localrepo. ['git rev-parse' does not give error for a hash thats not in repo]
petsc4py: revert to manually syncing petsc4py with petsc master changes
configure: revert gitcommit default value to Null. If a pacakge using giturl does not change this value - raise an error
fix: revert change in a44bcdf that breaks with pgf90
update test output for fixed to hypre interface that now actually pass the interp vectors to BoomerAMG
Error message should say use negative threshold value to prevent filtering out zeros
Calls to HYPRE_BoomerAMGSetInterpVectors() did not actually copy the vector data pointer from the PETSc vector to the hypre vector thus resulting in the hypre interp vectors being identically zero and hence not affecting the quality of the AMG hierarchy produced/
DMPlex: Silence overzealous compiler warnings
add row/col orientation to MatConvert_Elemental_Dense
fix formatting error in previous commit [broke nightlybuilds]
fix formating of TSTHETA manual page and add missing command line option
run hwloc with MPIEXEC and binding so it runs on the correct machine
fix: memory leaks [in valgrind test]
...

show more ...

Merge branch 'master' into barry-allow-vecset-nan-3.6

Merge branch 'stefano_zampini/feature-pardiso-schur' of https://bitbucket.org/petsc/petsc into stefano_zampini/feature-pardiso-schur

Merge branch 'master' into knepley/fix-examples-plex

* master: (112 commits)
correct expression for implict ODE/DAE in table 11 of users manual
remove warning about deprecated MPI call in exampl

Merge branch 'master' into knepley/fix-examples-plex

* master: (112 commits)
correct expression for implict ODE/DAE in table 11 of users manual
remove warning about deprecated MPI call in examples
rm variables declared in middle of subroutine
tolerance for equality in PETSc test should depend on machine epsilon
tolerance for acccuracy in test example should depend on machine epsilon
increase convergence tolerance to get same result in single precision as double for convergence history
update no change output for gasm test which had not been properly updated
remove C++ // comment
example mistakenly used MPI_COMM_WORLD as if it was rank in print statement thus producing output only with MPI Uni
grep out name of DM object since on some systems the address printed has extra zeros thus producing a false difference in output
increase tolerance for test example
build: revert el6-4 build to cg - as el6 compilers give link errors
fix error type
insure similar convergence history for single, double, quad
fix ex104.c based on Barry's suggestion
increase example convergence tolerance to get matching convergence history in single and double precision
increase convergence tolerance so that test example has similar convergence history in single and double precisions
changed from cg to gmres for eigenestimator since in single precision cg produced a possible indefinite matrix report
change test example tolerance to depend on machine epsilon so as to give correct results in various precisions
loosen convergence tolerance to have similar convergence history with single, double, quad precisions
...

show more ...

Merge branch 'master' into shri/add-tsevent-to-powergrid-examples