portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-15 10:15:36 +00:00
Code Releases Activity

Added ImplicitlyRestartedLanczosCJ to Algorithms.h

Browse Source
This commit is contained in:
Chulwoo Jung 2017-07-28 15:33:59 -04:00
parent 93650f3a61
commit 191fbf85fc
2 changed files with 9 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
lib/algorithms/Algorithms.h
View File

@ -38,6 +38,10 @@ Author: Peter Boyle <paboyle@ph.ed.ac.uk>
#include <Grid/algorithms/approx/Remez.h> #include <Grid/algorithms/approx/Remez.h>
#include <Grid/algorithms/approx/MultiShiftFunction.h> #include <Grid/algorithms/approx/MultiShiftFunction.h>
#include <Grid/algorithms/densematrix/DenseMatrix.h>
#include <Grid/algorithms/densematrix/Francis.h>
#include <Grid/algorithms/densematrix/Householder.h>
#include <Grid/algorithms/iterative/ConjugateGradient.h> #include <Grid/algorithms/iterative/ConjugateGradient.h>
#include <Grid/algorithms/iterative/ConjugateResidual.h> #include <Grid/algorithms/iterative/ConjugateResidual.h>
#include <Grid/algorithms/iterative/NormalEquations.h> #include <Grid/algorithms/iterative/NormalEquations.h>
@ -48,6 +52,7 @@ Author: Peter Boyle <paboyle@ph.ed.ac.uk>
// Lanczos support // Lanczos support
//#include <Grid/algorithms/iterative/MatrixUtils.h> //#include <Grid/algorithms/iterative/MatrixUtils.h>
#include <Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h> #include <Grid/algorithms/iterative/ImplicitlyRestartedLanczos.h>
#include <Grid/algorithms/iterative/ImplicitlyRestartedLanczosCJ.h>
#include <Grid/algorithms/iterative/SimpleLanczos.h> #include <Grid/algorithms/iterative/SimpleLanczos.h>
#include <Grid/algorithms/CoarsenedMatrix.h> #include <Grid/algorithms/CoarsenedMatrix.h>
#include <Grid/algorithms/FFT.h> #include <Grid/algorithms/FFT.h>

8
lib/algorithms/iterative/ImplicitlyRestartedLanczosCJ.h
View File

@ -27,8 +27,8 @@ Author: Chulwoo Jung <chulwoo@bnl.gov>
See the full license in the file "LICENSE" in the top level distribution directory See the full license in the file "LICENSE" in the top level distribution directory
*************************************************************************************/ *************************************************************************************/
/* END LEGAL */ /* END LEGAL */
#ifndef GRID_IRL_H #ifndef GRID_IRL_CJ_H
#define GRID_IRL_H #define GRID_IRL_CJ_H
#include <string.h> //memset #include <string.h> //memset
@ -91,7 +91,7 @@ public:
void init(void){}; void init(void){};
void Abort(int ff, DenseVector<RealD> &evals, DenseVector<DenseVector<RealD> > &evecs); void Abort(int ff, DenseVector<RealD> &evals, DenseVector<DenseVector<RealD> > &evecs);
ImplicitlyRestartedLanczos( ImplicitlyRestartedLanczosCJ(
LinearOperatorBase<Field> &Linop, // op LinearOperatorBase<Field> &Linop, // op
OperatorFunction<Field> & poly, // polynmial OperatorFunction<Field> & poly, // polynmial
int _Nstop, // sought vecs int _Nstop, // sought vecs
@ -110,7 +110,7 @@ public:
Np = Nm-Nk; assert(Np>0); Np = Nm-Nk; assert(Np>0);
}; };
ImplicitlyRestartedLanczos( ImplicitlyRestartedLanczosCJ(
LinearOperatorBase<Field> &Linop, // op LinearOperatorBase<Field> &Linop, // op
OperatorFunction<Field> & poly, // polynmial OperatorFunction<Field> & poly, // polynmial
int _Nk, // sought vecs int _Nk, // sought vecs