Changeset [6afc8368186cc737c9d4a9498ac4c9858c6db5bd] by Jed Brown
Revision 6afc8368186cc737c9d4a9498ac4c9858c6db5bd,
committed by Jed Brown
· May 10th, 2009 @ 05:35 PM
Add ability to apply left- and right-scaling to the sparse Jacobian. The scaling was optimized in the 1D case to minimize the condition number of the preconditioned spectral operator. Unfortunately, it does not seem to be a clear win on real problems (it is usually worse). It can be activated by -djac_tensor{mass,laplace}scale TRUE if supported by the assembly function (see ellip.c for example). http://github.com/jedbrown/dohp/...
11 changed files
-
M
include/dohpfs.h -
M
include/dohpjacobi.h -
M
include/private/jacimpl.h -
M
src/fs/mesh/interface/mesh.c -
M
src/fs/tests/ellip.c -
M
src/jacobi/impls/tensor/efstopo.c -
M
src/jacobi/impls/tensor/optimalscale.h -
M
src/jacobi/impls/tensor/tensor.c -
M
src/jacobi/impls/tensor/tensor.h -
M
src/jacobi/interface/jacobi.c -
M
src/jacobi/tests/ex1.c