Changeset [6afc8368186cc737c9d4a9498ac4c9858c6db5bd] 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/...
Committed by Jed Brown
- 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
Create new ticket
Create your profile
Help contribute to this project by taking a few moments to create your personal profile. Create your profile ยป
An implementation of the ``dual order hp'' version of the finite element method. This project targets parallel domain-decomposition methods for strongly coupled nonlinear problems with PDE constraints.