Fix errors revealed by gitlabci

Use PETSC_SMALL to determine when quadrature points are the same in PETSCDUALSPACELAGRANGE

Filter floating point output

Explicitly cast PetscBool to PetscInt before

Fix errors revealed by gitlabci

Use PETSC_SMALL to determine when quadrature points are the same in PETSCDUALSPACELAGRANGE

Filter floating point output

Explicitly cast PetscBool to PetscInt before printing

Do not construct vertex indices for simplex spaces in tensor cells

Right now the only way simplex spaces in tensor cells happens is discontinuous pressures
(typically in Q2P1). In this case there is no need to compute symmetries for the pressure
dofs.

show more ...

Add PetscDTIndexToBary, PetscDTBaryToIndex, and PetscDTNodeType

Move tests that missed the example excision

Add enum lists for PetscDTNodeTypes

Merge remote-tracking branch 'origin/psanan/docs-fixes-2020-03-25'

Doc fixes (prerequisites for Sphinx docs)

See merge request petsc/petsc!2633

docs: fix sowing list syntax/incomplete man page

doc fixes

Merge branch 'jed/promote-examples-tests-tutorials' [petsc/petsc!2610]

* jed/promote-examples-tests-tutorials:
Promote examples/{tests,tutorials}/ to {tests,tutorials}/

Promote examples/{tests,tutorials}/ to {tests,tutorials}/

This shortens paths and improves consistency between test target names
and paths to the source and output files. Most of the work was
compl

Promote examples/{tests,tutorials}/ to {tests,tutorials}/

This shortens paths and improves consistency between test target names
and paths to the source and output files. Most of the work was
completed by this script, followed by mild cleanup of nonconforming
cases.

for makefile in `git ls-files 'src/*makefile'`; do
if rg -q 'DIRS.*\bexamples\b' $makefile; then
base=$(dirname $makefile)
dirs=$(cd $base/examples && ls -d tests tutorials 2>/dev/null | xargs echo)
perl -pi -e "s#^(DIRS.*)\bexamples\b#\1${dirs}#" $makefile
git rm $base/examples/makefile
for t in $dirs; do
git mv $base/examples/$t $base/
perl -pi -e "s#^(LOCDIR[[:space:]]*=).*#\1 $base/$t/#" $base/$t
done
fi
done

git grep -l -E -z 'examples/(tutorials|tests)' | xargs -0 perl -pi -e 's#examples/(tutorials|tests)#\1#g'
git checkout @ \
src/docs/website/documentation/changes/ \
src/benchmarks/results/

show more ...

Merge branch 'jed/fix-unused-variable-warning' into 'master'

PetscDT: fix unused variable warning when not PETSC_HAVE_LGAMMA

See merge request petsc/petsc!2584

PetscDT: fix unused variable warning when not PETSC_HAVE_LGAMMA

Merge branch 'tisaac/feature-unify-quadrature' into 'master'

Unify Jacobi polynomial based computations in PetscDT

See merge request petsc/petsc!2574

Add PetscDTGaussJacobiQuadrature, PetscDTGaussLobattoJacobiQuadrature, PetscDTJacobiEval

Implement the Jacobi recurrence in exactly one place

Add PetscDTJacobiEval

squashme: Fixes to Jacobi polyno

Add PetscDTGaussJacobiQuadrature, PetscDTGaussLobattoJacobiQuadrature, PetscDTJacobiEval

Implement the Jacobi recurrence in exactly one place

Add PetscDTJacobiEval

squashme: Fixes to Jacobi polynomial changes

Use lgamma to determine the scaling coefficient in Newton's method for Gauss-Jacobi quadrature

Fix complex-to-real in postprocessing eigensolve in dt.c

Correctly choose newton for PetscGaussLobattoLegendreQuadrature

Use better initial guess for PetscDTGaussJacobiQuadrature_Newton_Internal

Work around Sun bug (nullus cimex sub sole novus)

show more ...

Rename PetscDTGaussJacobiQuadrature to PetscDTStroudConicalQuadrature and consolidate internal quadrature calculations

Revert incorrectly changed quadrature in dmmbfem.cxx

Merge branch 'balay/essl-reorg-PETSC_MISSING_LAPACK' into 'master'

Balay/essl reorg petsc missing lapack

See merge request petsc/petsc!2566

ESSL: reorganize code that handles PETSC_MISSING_LAPACK by using PetscMissingLapack()

And sync the missing function list with missing functions in current ESSL

If we find more LAPACK symbols missin

ESSL: reorganize code that handles PETSC_MISSING_LAPACK by using PetscMissingLapack()

And sync the missing function list with missing functions in current ESSL

If we find more LAPACK symbols missing in commonly used LAPACK impls,
we can add them back in the new format.

src/ksp/pc/impls/bddc/bddcprivate.c has alternate code for
PETSC_MISSING_LAPACK_GESVD. This code is retained for now [in case we
need to add this flag back in]

show more ...

Merge remote-tracking branch 'origin/maint'

Merge remote-tracking branch 'origin/balay/essl-fixes/maint' into maint

essl: fix build by protecting code with PETSC_MISSING_LAPACK_ORGQR etc flags

Reported-by: Victor Eijkhout <eijkhout@tacc.utexas.edu>

Merge remote-tracking branch 'origin/knepley/feature-dt-quad-classid'

Knepley/feature dt quad classid

See merge request petsc/petsc!2355

Merge branch 'tisaac/feature-alt-v' into 'master'

Tisaac/feature alt v

See merge request petsc/petsc!2357

Fix uninitialized, C++, and Fortran issues with PetscDTAltV, other small
fixes

Silence compiler Wmaybe-unitialized warnings in petscdt.h

Loosen tolerances for single precision tests

Fix C++ PetscD

Fix uninitialized, C++, and Fortran issues with PetscDTAltV, other small
fixes

Silence compiler Wmaybe-unitialized warnings in petscdt.h

Loosen tolerances for single precision tests

Fix C++ PetscDTAltv test ex7

Fix complex<->real loops in PetscDTJacobianInverse_Internal

Fix -Wmaybe-unitialized errors in ex7

Improve documentation of PetscDTAltV thanks to criticism from Barry Smith

Protect PetscDTFactorial against negative input

Unbracket SETERRQ

Make "the" -> "a"/"an" more input argument documentation in PetscDTAltV

Rename formIndex to formDegree for consistent notation

show more ...

Add PetscQuadraturePushForward to make it easier to use arbitrary node sets in mixed elements

Use getri in PetscDTJacobianInverse_Internal

Remove unneeded header from dt.c that was introduced in fe

Add PetscQuadraturePushForward to make it easier to use arbitrary node sets in mixed elements

Use getri in PetscDTJacobianInverse_Internal

Remove unneeded header from dt.c that was introduced in feature development

show more ...

Document static inlines in petscdt.h

Expose PetscDTFactorial and PetscDTFactorialInt, match with PetscDTBinomial and PetscDTBinomialInt

Document PetscDTAltV

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 ...

PetscQuadrature: Finally it gets a class id, and like Pinocchio, becomes a real boy

Merge branch 'maint'