implement version of STOREVALS for MatFDColoringApply_SeqBAIJ, but it does not improve the version of !STOREVALS

bugfix of MatFDColoringApply_BAIJ

mv MatFDColoringApply_SeqAIJ() from aij.c to fdaij.c; rm MatFDColoringApply_SeqAIJ_old; MatFDColoringCreate_SeqAIJ_den2sp() supports SeqBAIJ matrix.

mv redundant valications in MatFDColoringApply_xxx()

Merge branch 'barry/improve-memory-logging'

Conflicts:
src/ksp/pc/impls/redundant/redundant.c
src/mat/interface/matnull.c

removed Mat option for skipping checking of compressed row format since it is now free using the nonzerorowcnt previously introduced

changed seqbaij to count nonzero rows in assembly and use those values to compute flops in MatMult

Merge branch 'barry/fix-PetscObjectPrintClassNamePrefixType'

* barry/fix-PetscObjectPrintClassNamePrefixType:
update tests of src/mat/examples/tests/ex136.c
Use classname in object for object na

Merge branch 'barry/fix-PetscObjectPrintClassNamePrefixType'

* barry/fix-PetscObjectPrintClassNamePrefixType:
update tests of src/mat/examples/tests/ex136.c
Use classname in object for object name when calling PetscObjectPrintClassNamePrefixType() instead of some string with less precise name passed in

Conflicts:
src/snes/interface/snes.c

show more ...

baijfact2.c: splitup the big file into multiple files - so asto
reduce load on some compilers [which break in c++/complex/optimzied mode]

baijsolv.c baijsolvnat.c baijsolvtran.c baijsolvtrannat.c

logging memory now credits to all ancestors

Use classname in object for object name when calling PetscObjectPrintClassNamePrefixType() instead of some string with less precise name passed in

Merge branch 'jed/stack-call-blas'

Fix performance bug that calling BLAS takes an amount of time dependent
on the number of memory allocations that the program has made.

* jed/stack-call-blas:
Pe

Merge branch 'jed/stack-call-blas'

Fix performance bug that calling BLAS takes an amount of time dependent
on the number of memory allocations that the program has made.

* jed/stack-call-blas:
PetscStackCallBLAS: new macro instead of PetscStackCall (skip CHKMEMQ)
Stack: consolidate handling of PetscStack

show more ...

Merge branch 'toby/fix-sor-baij'

Shell portability and silence bad compiler warning.

* toby/fix-sor-baij:
MatSOR_SeqBAIJ: initialize work for pedantic compilers.
KSP test: avoid bash4 case synt

Merge branch 'toby/fix-sor-baij'

Shell portability and silence bad compiler warning.

* toby/fix-sor-baij:
MatSOR_SeqBAIJ: initialize work for pedantic compilers.
KSP test: avoid bash4 case syntax ';&' in favor of simpler alternative

show more ...

MatSOR_SeqBAIJ: initialize work for pedantic compilers.

The work variable is initialized in a macro within a do{}while(0) construct:
some compilers don't realize that this initialization always happ

MatSOR_SeqBAIJ: initialize work for pedantic compilers.

The work variable is initialized in a macro within a do{}while(0) construct:
some compilers don't realize that this initialization always happens, so they
complain about uninitialized values. This protects against those warnings.

show more ...

Merge branch 'toby/fix-sor-baij'

Consolidates code and adds support for BAIJ SOR with nonzero initial
guess (required for use as a MG smoother).

* toby/fix-sor-baij:
Mat: add MatSOR/BAIJ test to

Merge branch 'toby/fix-sor-baij'

Consolidates code and adds support for BAIJ SOR with nonzero initial
guess (required for use as a MG smoother).

* toby/fix-sor-baij:
Mat: add MatSOR/BAIJ test to nightlies
Mat: move new MatSOR macros to a common header blockmatmult.h
MatSOR_SeqBAIJ: refactor, enable nonzero initial guess.
MatSOR, test ex43: removed BlockMat from test
MatSOR: added test ex43, currently fails
MatSOR_SeqBAIJ: protect against zero vector length

Conflicts:
src/ksp/ksp/examples/tests/makefile

show more ...

PetscStackCallBLAS: new macro instead of PetscStackCall (skip CHKMEMQ)

When CHKMEMQ occurs on every call to BLAS or Lapack, we sometimes see
unacceptable performance degradation (orders of magnitude

PetscStackCallBLAS: new macro instead of PetscStackCall (skip CHKMEMQ)

When CHKMEMQ occurs on every call to BLAS or Lapack, we sometimes see
unacceptable performance degradation (orders of magnitude).
PetscStackCallBLAS is currently identical to PetscStackCall except for
skipping CHKMEMQ, but such checking could be made optional in the
future.

show more ...

Mat: move new MatSOR macros to a common header blockmatmult.h

These macros would be useful for condensing the BAIJ triangular solve
code, for example.

MatSOR_SeqBAIJ: refactor, enable nonzero initial guess.

At the cost of a couple of conditionals per iteration, the various
MatSOR_SeqAIJ_# functions are combined.

Passes make testexamples, includin

MatSOR_SeqBAIJ: refactor, enable nonzero initial guess.

At the cost of a couple of conditionals per iteration, the various
MatSOR_SeqAIJ_# functions are combined.

Passes make testexamples, including recently added ksp test ex43

show more ...

MatSOR_SeqBAIJ: protect against zero vector length

removed string version of function name for XXXRegister(), PetscFunctionListAdd() and PetscObjectComposeFunction()

removed excessive use of CHKMEMQ;