mirror of
https://github.com/paboyle/Grid.git
synced 2024-11-10 15:55:37 +00:00
Use 2 passes of GS in coarse operator construction
This commit is contained in:
parent
194e4b94bb
commit
faf4278019
@ -334,10 +334,9 @@ namespace Grid {
|
||||
|
||||
CoarseVector iProj(Grid());
|
||||
CoarseVector oProj(Grid());
|
||||
CoarseScalar InnerProd(Grid());
|
||||
|
||||
// Orthogonalise the subblocks over the basis
|
||||
blockOrthogonalise(InnerProd,Subspace.subspace);
|
||||
Subspace.Orthogonalise();
|
||||
|
||||
// Compute the matrix elements of linop between this orthonormal
|
||||
// set of vectors.
|
||||
|
Loading…
Reference in New Issue
Block a user