| 08122e43 | 27-Jan-2015 |
Stefano Zampini <stefano.zampini@gmail.com> |
PCBDDC: added initial implementation for adaptive selection of constraints
positive definite case covered only
untested for complex
generalized eigenvalue problem is inv(sum(inv(S_F))) = lambda inv
PCBDDC: added initial implementation for adaptive selection of constraints
positive definite case covered only
untested for complex
generalized eigenvalue problem is inv(sum(inv(S_F))) = lambda inv(sum(inv(St_F)))
this is the right one for 2D problems and for faces in 3D; it is an heuristic for edges in 3D
tested with H(div) and H(curl) (3D) and works like a charm
cholesky support for matrix inversion should still be added
show more ...
|
| 1996caf8 | 22-Jan-2015 |
Toby Isaac <tisaac@ices.utexas.edu> |
Move MG and GAMG implementation headers to petsc-private
I want to be able to register GAMG plugins. I considered a few options:
a) Move the implementation headers to petsc-private so they can be
Move MG and GAMG implementation headers to petsc-private
I want to be able to register GAMG plugins. I considered a few options:
a) Move the implementation headers to petsc-private so they can be
installed
b) Make Get/Set methods for all of the GAMGOps: graph, coarsen,
optprol, etc.
c) Add a PCGAMGGet/SetOperation like MatShellGet/SetOperation for each
of these operations.
The easiest was (a), and the GAMG methods are so context-specific I'm
not sure who would benefit from their being part of the public petsc
interface. If GAMG get's refactored to where the operations are easier
to understand on their own, then maybe (b) or (c) would be a better
option.
show more ...
|