| 8f1c130e | 26-Nov-2014 |
Stefano Zampini <stefano.zampini@gmail.com> |
PCBDDC: no need to reorder the constraints in local ordering after SVD
This also guaranties that the change of basis computed with QR will be coherent among processes
It also extends the reordering
PCBDDC: no need to reorder the constraints in local ordering after SVD
This also guaranties that the change of basis computed with QR will be coherent among processes
It also extends the reordering to the case of using QR on a single cc
show more ...
|
| 911cabfe | 26-Nov-2014 |
Stefano Zampini <stefano.zampini@gmail.com> |
PCBDDC: fixed bug when using more than 1 constraint per connected component via SVD
The bug originated from the (wrong!) assumption that the returned eigenmodes from the SVD would be consistent amon
PCBDDC: fixed bug when using more than 1 constraint per connected component via SVD
The bug originated from the (wrong!) assumption that the returned eigenmodes from the SVD would be consistent among sharing processes
This is not true since the computed SVD depends on the local ordering of dofs and each process computes its own SVD (separately for each connected component).
The fix consists in ordering the constrained indices consistently across processes by using global ordering
show more ...
|
| a6b551f4 | 24-Nov-2014 |
Stefano Zampini <stefano.zampini@gmail.com> |
PCBDDC: fix wrong computation of change of basis blocks
The last column of the block should be given by the orthonormalized coefficients of the quadrature rule
Note that the computed block is the sa
PCBDDC: fix wrong computation of change of basis blocks
The last column of the block should be given by the orthonormalized coefficients of the quadrature rule
Note that the computed block is the same as before if computed with respect the arithmetic average
show more ...
|
| fa434743 | 24-Nov-2014 |
Stefano Zampini <stefano.zampini@gmail.com> |
PCBDDC: QR decomposition for change of basis can now be also used for connected components with only 1 constraint
This partially solves the following issue with Raviart Thomas elements
Using standar
PCBDDC: QR decomposition for change of basis can now be also used for connected components with only 1 constraint
This partially solves the following issue with Raviart Thomas elements
Using standard change of basis, the coarse matrix produced is singular and CG doesn't converge to the exact solution (even if the eigenvalues of the preconditioned operator are the same of the unchanged basis case)
On the other hand, using QR, CG converges to the exact solution but the eigenvalues of the preconditioned operator are very different from not using the change of basis.
Still needs to be investigated...
show more ...
|