portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-06-12 20:27:06 +01:00
Code Releases Activity

Added one flavour rational ratios (unprec)

Browse Source
This commit is contained in:
Peter Boyle
2015-08-19 04:58:40 +01:00
parent 570150f1d3
commit fc160eeccc
9 changed files with 315 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
lib/algorithms/approx/MultiShiftFunction.h
View File

@ -27,6 +27,8 @@ public:
if ( inverse ) remez.getIPFE (&residues[0],&poles[0],&norm);
else remez.getPFE (&residues[0],&poles[0],&norm);
}
// Allow deferred initialisation
MultiShiftFunction(void){};
MultiShiftFunction(AlgRemez & remez,double tol,bool inverse)
{
Init(remez,tol,inverse);

6
lib/algorithms/iterative/ConjugateGradientMultiShift.h
View File

@ -27,10 +27,14 @@ public:
void operator() (LinearOperatorBase<Field> &Linop, const Field &src, Field &psi)
{
GridBase *grid = src._grid;
int nshift = shifts.order;
std::vector<Field> results(nshift,grid);
(*this)(Linop,src,results,psi);
}
void operator() (LinearOperatorBase<Field> &Linop, const Field &src, std::vector<Field> &results, Field &psi)
{
int nshift = shifts.order;
(*this)(Linop,src,results);