Start PetscDTAltV

Added reference implementation of alternating multilinear forms

But I'm going to clean up the interface. I don't think an actual object is necessary
and we want to optimize the c

Start PetscDTAltV

Added reference implementation of alternating multilinear forms

But I'm going to clean up the interface. I don't think an actual object is necessary
and we want to optimize the cases N = 2, 3

PetscDTAltV: change representation order

so that when k != N / 2, there is no reordering in hodge star

PetscDTAltv: New basis for alterting multilinear forms.

For $k$-forms with $2 k \leq N$, the standard lexicographic
basis of $\{dx_{\sigma(1)} \wedge \dots \wedge dx_{\sigma(k)}\}$
is used.

For $2 k > N$, we modify this by first reversing the lexicographic
order, and then swapping $\sigma(1)$ and $\sigma(k)$ if $\sigma$ is odd.
While this involves more sign flips than the basic approach, it has two
good properties:

- it preserves the case as before that the dofs of a row vector /
covector are the dofs of the equivalent 1-form, but now
- the dofs of a contravariant vector (like an H(div)-conforming flux)
are the dofs of the equivalent (N-1)-form (importantly, this means
that the contravariant pullback is the same as the alternating form
pullback)

The only exception is 2D, for which the case is the same as before: the
1-form pullback is the covariant pullback, the contravariant pullback
can be obtained by Hodge star isometry before and afer the pullback.

Remove bad file from rename during merge

PetscDTAltV: interface simplifaction, undo specialized basis.

The fact that the N=2 case wasn't handle bothered me. Now I think the
thing to do is to let PetscDTAltVPullback() and PetscDTAltVInterior() take
negative form degrees, indicating that the pullback/interior product should be
conjugated with the hodge star.

Fix c++ PetscBool issues with PetscDT static inlines

PetscDTAltV: negative form index in pullback...

indicates $\star^{-1} \circ L* \circ$. This gives the contravariant
pullback by $ k=-(N-1)$ in all dimensions.

PetscDTAltV: stop using altv data structures

I'm going to remove the object

PetscDTAltV: Remove object

PetscDTAltV: avoid mallocs/frees for N <= 3

show more ...

Merge branch 'knepley/feature-move-fv' into 'master'

Knepley/feature move fv

See merge request petsc/petsc!1982

FV: Move to its own directory

Merge branch 'master' into jpusztay/feature-swarm-symplectic-example

Merge branch 'master' into jpustay/feature-swarm-example

Merged petsc/petsc into master

Merged petsc/petsc into master

Merge branch 'master' into arcowie-rem/feature-error-logging

* master: (345 commits)
Update dev notes for the coming release
Increase patchlevel to 3.9.4
doc fixes
Added some AIJSELL tests t

Merge branch 'master' into arcowie-rem/feature-error-logging

* master: (345 commits)
Update dev notes for the coming release
Increase patchlevel to 3.9.4
doc fixes
Added some AIJSELL tests that mirror the SELL tests.
Bib: More refs
Removed use of hardwired -1 to represent PETSC_NULL_XXX object, since the null value is now 0.
fix memleaks flagged by valgrind
ex62: disable long running test and workaround differences in iteration counts on various builds.
Add FAQ we can point to when users get error messages such as: "PETSc was configured with one OpenMPI mpi.h version but now appears to be compiling using a different OpenMPI mpi.h version"
swith test from mpich to openmpi as mpich ch3:sock is unuserable per PetscShmCommGet()
DM: Forgot to move include statement
fix nightlybuild issues - disable hypre [for openmp] and fix uninitialized variables
tests: enable openmp and hwloc with one of the mumps tests
Default BoomerAMG smoothers
Configure: Revert changes from debb90c, so that we are back to the original Libmesh
Configure: Reverted change for Libmesh from 83cd5d7
Split vecimpl.h into vecimpl.h and vecscatterimpl.h
SNES ex62: Turned off some tests for other precisions and shortened output of sensitive tests
Use mmap to allocate the shared pthread barrier if available
Add configure test for functional mmap
...

Conflicts:
include/petscconvest.h
src/snes/examples/tutorials/ex13.c
src/snes/interface/snes.c
src/snes/utils/convest.c

show more ...

Merge branch 'master' into knepley/fix-superdm-is

* master: (123 commits)
Fix a few more Fortran objects to have their NULL value match the C NULL value
remove dead code
Add variable point blo

Merge branch 'master' into knepley/fix-superdm-is

* master: (123 commits)
Fix a few more Fortran objects to have their NULL value match the C NULL value
remove dead code
Add variable point block Jacobi for MPIAIJ matrix.
should not pass user contexts through CHKFORTRANNULLOBJECT in Fortran stubs
fix mpiuni test
Remove CHKFORTRANNULLOBJECTDEREFERENCE since it is no longer needed with Fortran null equal to C null
PETSCFEOPENCL: declare PetscFEGetTabulation_Basic as external
PETSCFEOPENCL: declare PetscFESetUp_Basic as external
Change Fortran null pointer to match C null pointer
one more fix
jenkins: run 'cleantest' after examples are run - to recover diskspace
logical order bug in diagbrdn update fixed, tests updated again
replace new -petscspace_order instances with -petscspace_degree
PETSCSPACETENSOR: clear up naming confusion of subspaces
regression.py: remove configure code related to old test suite - now that SLEPc no longer needs it
PetscSpaceSetup_Poly(): make idempotent
Plex ex3: update tests
Lift changes to defunct dtfe.c from [PR #1074]
PetscSpace_Tensor: clarification of vector-valued
PetscSpace_Tensor: allow height subspaces, fixes tabulation initialization
...

show more ...

Merge branch 'master' into knepley/feature-pc-patch

* master: (131 commits)
fix mpiuni test
PETSCFEOPENCL: declare PetscFEGetTabulation_Basic as external
PETSCFEOPENCL: declare PetscFESetUp_Ba

Merge branch 'master' into knepley/feature-pc-patch

* master: (131 commits)
fix mpiuni test
PETSCFEOPENCL: declare PetscFEGetTabulation_Basic as external
PETSCFEOPENCL: declare PetscFESetUp_Basic as external
one more fix
jenkins: run 'cleantest' after examples are run - to recover diskspace
logical order bug in diagbrdn update fixed, tests updated again
replace new -petscspace_order instances with -petscspace_degree
PETSCSPACETENSOR: clear up naming confusion of subspaces
regression.py: remove configure code related to old test suite - now that SLEPc no longer needs it
PetscSpaceSetup_Poly(): make idempotent
Plex ex3: update tests
Lift changes to defunct dtfe.c from [PR #1074]
PetscSpace_Tensor: clarification of vector-valued
PetscSpace_Tensor: allow height subspaces, fixes tabulation initialization
PetscSpace: add a dimension member
PetscSpaceSetDegree(): now take maxDegree argument
-petscspacepoly_tensor -> -petscspace_type tensor
PetscSpaceEvaluate_Polynomial(): fix temp array dependencies
PETSCSPACETENSOR: A generic tensor product of spaces.
PetscSpace: -petscspace_order -> -petscspace_degree
...

Conflicts:
src/dm/dt/interface/dtfe.c
src/snes/examples/tutorials/ex13.c
src/snes/examples/tutorials/ex17.c
src/snes/examples/tutorials/ex62.c
src/snes/examples/tutorials/output/ex62_bddc_quad.out
src/snes/examples/tutorials/output/ex62_fetidp_2d_quad.out
src/snes/examples/tutorials/output/ex62_fetidp_2d_tri.out
src/snes/examples/tutorials/output/ex62_fetidp_3d_hex.out
src/snes/examples/tutorials/output/ex62_fetidp_3d_tet.out

show more ...

Merge remote-tracking branch 'origin/tisaac/feature-space-tensor'

PetscFE: create directory structure

It's becoming necessary, as we try to do more complex things,
to sharpen the interface between the things in dtfe, and make
extension more approachable. It's tim

PetscFE: create directory structure

It's becoming necessary, as we try to do more complex things,
to sharpen the interface between the things in dtfe, and make
extension more approachable. It's time.

show more ...

Merge commit 'd47bf9aaf1e5266cc3f0ff499f934c85788965a9' into jed/fix-matcreatempibaij/maint

Obtain changed (v3.5) handling of MatMPIBAIJSetPreallocationCSR to match
documentation.

Merge branch 'master' into stefano_zampini/feature-gmsh-periodic

* master: (1016 commits)
metis: setself.downloaddirnames for change in 6d90219d8fa4f8b2fd242cc6ff9008358d87c26f
fix compile warni

Merge branch 'master' into stefano_zampini/feature-gmsh-periodic

* master: (1016 commits)
metis: setself.downloaddirnames for change in 6d90219d8fa4f8b2fd242cc6ff9008358d87c26f
fix compile warning: /sandbox/petsc/petsc.clone/src/snes/examples/tests/ex2.c:250:8: warning: ‘pointsAllProcs’ may be used uninitialized in this function [-Wmaybe-uninitialized] ierr = DMInterpolationSetUp(interpolator, dm, pointsAllProcs);CHKERRQ(ierr);
migrate 'build requires: !mpiuni' to new format
migrate 'build requires: !complex' to new format
workarround 'division by zero' compiler warning on windows C:\cygwin64\home\petsc\PETSC~3.CLO\src\ts\examples\TUTORI~1\ex11.c(416): warning #39: division by zero if (uL->h < 0 || uR->h < 0) {for (i=0; i<1+dim; i++) flux[i] = 0./0.; return;} /* SETERRQ(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"Reconstructed thickness is negative"); */
dos2unix to prevent warnings on opensolaris "/export/home/petsc/petsc.clone/src/ksp/ksp/examples/tests/ex53.c", line 6: warning: invalid white space character in directive
fix output file - corresponding to code changes in efd4aadf157bf1ba2d80c2be092fcf4247860003
disable test on windows
update log from opensolaris so that the diff is clean on linux aswell
update grep in sync with change in efd4aadf157bf1ba2d80c2be092fcf4247860003
Fixed code that did not correctly set current color for block insertion
Update output files that did not get fixed after previous merges
Change MPIUni to not immediately abort on operations it cannot perform, instead return error code
Add filter for example that can produce 1 different iteration count
Example add system calloc, bad formating of %d in print, commented out code
Always use (void) for functions that take no arguments, never ()
always need to caste PetscReal variables to (double) in print statements
updates to manual page for clarity
update help menu for MatRARt() -- MatRARt_mpiaij_mpiaij is supported now, but not efficient.
fix nightlybuild error with c1cbb1def3f8ddca04c3998b9c06618fd1ffef35 on arch-linux-IPL64 src/ksp/pc/examples/tests ex5_1
...

Conflicts:
include/petscdm.h
src/dm/impls/plex/plexcreate.c
src/dm/impls/plex/plexgmsh.c
src/dm/impls/plex/plexhdf5.c

show more ...

Merge branch 'master' into knepley/feature-dm-bc-components

* master: (320 commits)
fix compile warning: /sandbox/petsc/petsc.clone/src/snes/examples/tests/ex2.c:250:8: warning: ‘pointsAllProcs’ m

Merge branch 'master' into knepley/feature-dm-bc-components

* master: (320 commits)
fix compile warning: /sandbox/petsc/petsc.clone/src/snes/examples/tests/ex2.c:250:8: warning: ‘pointsAllProcs’ may be used uninitialized in this function [-Wmaybe-uninitialized] ierr = DMInterpolationSetUp(interpolator, dm, pointsAllProcs);CHKERRQ(ierr);
migrate 'build requires: !mpiuni' to new format
migrate 'build requires: !complex' to new format
workarround 'division by zero' compiler warning on windows C:\cygwin64\home\petsc\PETSC~3.CLO\src\ts\examples\TUTORI~1\ex11.c(416): warning #39: division by zero if (uL->h < 0 || uR->h < 0) {for (i=0; i<1+dim; i++) flux[i] = 0./0.; return;} /* SETERRQ(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"Reconstructed thickness is negative"); */
dos2unix to prevent warnings on opensolaris "/export/home/petsc/petsc.clone/src/ksp/ksp/examples/tests/ex53.c", line 6: warning: invalid white space character in directive
fix output file - corresponding to code changes in efd4aadf157bf1ba2d80c2be092fcf4247860003
disable test on windows
update log from opensolaris so that the diff is clean on linux aswell
update grep in sync with change in efd4aadf157bf1ba2d80c2be092fcf4247860003
Fixed code that did not correctly set current color for block insertion
Update output files that did not get fixed after previous merges
Change MPIUni to not immediately abort on operations it cannot perform, instead return error code
Add filter for example that can produce 1 different iteration count
Example add system calloc, bad formating of %d in print, commented out code
Always use (void) for functions that take no arguments, never ()
always need to caste PetscReal variables to (double) in print statements
updates to manual page for clarity
update help menu for MatRARt() -- MatRARt_mpiaij_mpiaij is supported now, but not efficient.
fix nightlybuild error with c1cbb1def3f8ddca04c3998b9c06618fd1ffef35 on arch-linux-IPL64 src/ksp/pc/examples/tests ex5_1
output has changed with e6d0a238963c2a97dd04845ea512b529992c7cdb so update the output file
...

show more ...

Merge branch 'knepley/test-process-arg' into knepley/feature-adapt-interface

* knepley/test-process-arg: (49 commits)
Test: Can now use make test NP=<num procs>
Testing: Fix error reporting
I

Merge branch 'knepley/test-process-arg' into knepley/feature-adapt-interface

* knepley/test-process-arg: (49 commits)
Test: Can now use make test NP=<num procs>
Testing: Fix error reporting
I didn't realize there was a hacky special case for dense B matrix
Fix example to support -snes_mf_operator
cleanup src/ksp/ksp/examples/tutorials/network/makefile
add missing input in the help menu of PetscOptionsView() Reported-by: Richard Katz <richard.katz@earth.ox.ac.uk>
Fix docs for MatMatMult() that C matrix can only be created by user for sequential case.
protect '-' in macros
fix compile error with --download-petsc4py
Plex: Removed useless CHKERRQ
Configure: Replace use of PETSc modules
Plex ex7: Damn, quad precision is named __float128
TSSetFromOptions: Give subtypes a chance to set the default adapt type
TS ex11.c: fix corrupt filename initialization
small fix
Plex ex1: Update output
DT: Fixed Fortran module - Add DT module to DM and Petsc modules - Fix makefile to generate DT F90 module - Fix PetscQuadratureGet/Restore/SetData() F90 defs
Plex ex19: Updated output
Plex ex7: Update output
p4est: Update to new API
...

show more ...

Merge branch 'knepley/fix-fe-fortran'

* knepley/fix-fe-fortran:
DT: Fixed Fortran module - Add DT module to DM and Petsc modules - Fix makefile to generate DT F90 module - Fix PetscQuadratureGet/R

Merge branch 'knepley/fix-fe-fortran'

* knepley/fix-fe-fortran:
DT: Fixed Fortran module - Add DT module to DM and Petsc modules - Fix makefile to generate DT F90 module - Fix PetscQuadratureGet/Restore/SetData() F90 defs
PetscFE: Added Fortran bindings

show more ...

DT: Fixed Fortran module
- Add DT module to DM and Petsc modules
- Fix makefile to generate DT F90 module
- Fix PetscQuadratureGet/Restore/SetData() F90 defs

Merge branch 'master' into pr263/PierreBdR/ts-domain-error-in-stage/master

Merge branch 'barry/fix-nonew-notcollective/maint' into jed/mat-assembly-perf

Jed this fucking 6+ month running thing in next that is not in master is a royal fucking pain in the ass and
an abuse of

Merge branch 'barry/fix-nonew-notcollective/maint' into jed/mat-assembly-perf

Jed this fucking 6+ month running thing in next that is not in master is a royal fucking pain in the ass and
an abuse of git. I think it is not justified to ever have anything in next for more than a few weeks at most.
Either take it out of next if it is broken or put it in master if it is not broken. Hanging around in next but
not master for ever is not a good policy.

show more ...

Merge branch 'master' into knepley/solkx

* master: (6933 commits)
Bib: Added reference
Includegraph: Updated for compatibility with new folder layout.
Add SNESLineSearchReason to fortran inclu

Merge branch 'master' into knepley/solkx

* master: (6933 commits)
Bib: Added reference
Includegraph: Updated for compatibility with new folder layout.
Add SNESLineSearchReason to fortran includes Also truncate names to fit 32 char fortran limit
fix compile warnings /usr/home/balay/petsc.clone-2/src/vec/vec/interface/vector.c:1944: warning: division by zero in '1.0e+0 / 0.' /usr/home/balay/petsc.clone-2/src/ksp/ksp/interface/itfunc.c:508: warning: 'vec_rhs' may be used uninitialized in this function
updated output for new trust region initial size
initial size of the trust region is set as a percent of the norm of the initial guess, it should not be related to the initial norm of the function (which is kind of nuts).
The default bounds for SNESVISetVariableBounds() in the manual page were reversed from correct values
Bib: Added refs
Bib: Update MPICH webpage
Plex: Forgot to close file
removed nonexistent ex46f from makefile test rule
Revert "fixed bad merge into master"
fixed bad merge into master
fix memory leak in mkl_pardiso fix some formatting in mkl_pardiso code
mv runex111 from TESTEXAMPLES_C to TESTEXAMPLES_DATAFILESPATH
Mat+Doc: More Fortran docs
Plex ex3: Run tests with the correct number of field components
Bib: Added TetGen cite
parmetis: make sure 'ldd libparmetis.so' points to the correct libmetis.so that its linked with.
small fix as reported by the nightly test
...

Conflicts:
config/builder.py

show more ...

Merge branch 'master' into tisaac/dmforest

directory layout changed again: merging to switching branches back and forth
doesn't cause reconfigure to have a hard time

* master: (357 commits)
remov

Merge branch 'master' into tisaac/dmforest

directory layout changed again: merging to switching branches back and forth
doesn't cause reconfigure to have a hard time

* master: (357 commits)
removed nonexistent ex46f from makefile test rule
Revert "fixed bad merge into master"
fixed bad merge into master
fix memory leak in mkl_pardiso fix some formatting in mkl_pardiso code
mv runex111 from TESTEXAMPLES_C to TESTEXAMPLES_DATAFILESPATH
Mat+Doc: More Fortran docs
Bib: Added TetGen cite
parmetis: make sure 'ldd libparmetis.so' points to the correct libmetis.so that its linked with.
small fix as reported by the nightly test
add a test for MatPtAP()
small cleanup; add tests
cleanup MatPtAPSymbolic_MPIAIJ_MPIAIJ()
Fix nonzerostate tracking in all MATMPI types.
update location of test matrices in faq.html file
unused variable still declared
Fix for MatLoad_MPIDense() from sparsely stored matrix also fix for MatLoad_MPIDense_DenseInFile() when matrix has not had its local sizes preset by user
Preserve the B nonzerostate during disassembly.
TSEvent: Fix arch-mswin nightly errors due to declaring variables in the middle of a function.
Fix error in parallel loading of dense format matrix
removed KSPSetNullSpace() since mat already carries that information and having two locations for it is confusing, unnecessary and prone to error added MatSetTransposeNullSpace() for providing the null space of the transpose of the matrix and added support in KSPSolve() to remove this null space from the right handside of the linear system thus making the system consistent.
...

show more ...

Merge branch 'master' into knepley/feature-composable-callbacks

* master: (151 commits)
PetscMalloc: Fix handling of zero-sized allocations
add note to DMGetLocal/GlobalVector() about intended s

Merge branch 'master' into knepley/feature-composable-callbacks

* master: (151 commits)
PetscMalloc: Fix handling of zero-sized allocations
add note to DMGetLocal/GlobalVector() about intended short-time usage of the vectors.
sf: use patched openmpi tarball that has a potential fix for sf example hang
SNES problems that set VI bounds must use a VI solver since the other solvers ignore the bounds
pounders viewer did not cast PetscReal to double in printf nor use %D instead of %d
fix ts ex20adj
fix ts theta
fixed another variable declaration in the middle of a code block
example main() had PetscInt return type and argument that should be int
DMDAVTKWriteAll_VTS: Fix handling of higher-dimension coordinates
fixed example that declared variables in the middle of a function
C programs do not return PetscInt nor do they take an argument that is PetscInt, fixed example that did
You cannot declare variables in the middle of functions in PETSc.
add ierr checking to PetscOptionsBegin/End() in example
configure: add error check if user specifies PETSC_ARCH=-dbg [i.e starts with a '-']
scale norm of vector difference by norm of vector to insure it is small enough to consider the two vectors equal
missing PetscViewerDestroy() added
MatNestGetISs: use custom fortranstub to handle PETSC_NULL_OBJECT
update location of open-mx website on FAQ page
If the user sets bounds and calls a solver that does not support them, throw an error.
...

show more ...

Merge branch 'pr192/Fande-Kong/scalable-matincreaseoverlap/master' of https://bitbucket.org/petsc/petsc

Merged petsc/petsc into master