3
0
mirror of https://github.com/triqs/dft_tools synced 2024-12-25 13:53:40 +01:00

gf: correct bug in matsubara_freq mesh point

- bool at_end() for a mesh point means AFTER the last point ,
as end in STL. It was not correct ( +1 missing).
- also the at_end was not computed in the mesh product.
- it slightly changes the test_fit_tail --> just changed the output.
--> did this bug affect other functions/codes ?
This commit is contained in:
Olivier Parcollet 2014-02-19 21:14:59 +01:00
parent dcc7ad2c8c
commit 73d81cc55a
4 changed files with 6 additions and 6 deletions

View File

@ -66,7 +66,7 @@
... Order 2 = ... Order 2 =
[[(1,0)]] [[(1,0)]]
... Order 3 = ... Order 3 =
[[(0.999288,0)]] [[(0.999247,0)]]
... Order 4 = ... Order 4 =
[[(0.998655,0)]] [[(0.998655,0)]]
@ -80,7 +80,7 @@
... Order 2 = ... Order 2 =
[[(1,0)]] [[(1,0)]]
... Order 3 = ... Order 3 =
[[(0.999209,0)]] [[(0.999186,0)]]
... Order 4 = ... Order 4 =
[[(0.998631,0)]] [[(0.998631,0)]]
@ -94,7 +94,7 @@
... Order 2 = ... Order 2 =
[[(1,0)]] [[(1,0)]]
... Order 3 = ... Order 3 =
[[(0.999251,0)]] [[(0.999214,0)]]
... Order 4 = ... Order 4 =
[[(0.998655,0)]] [[(0.998655,0)]]

View File

@ -92,7 +92,7 @@ namespace gfs {
void advance() { ++n; } void advance() { ++n; }
long linear_index() const { return n - first_index; } long linear_index() const { return n - first_index; }
long index() const { return n; } long index() const { return n; }
bool at_end() const { return (n == index_stop); } bool at_end() const { return (n == index_stop + 1); } // at_end means " one after the last one", as in STL
void reset() { n = first_index; } void reset() { n = first_index; }
private: private:

View File

@ -156,7 +156,7 @@ namespace gfs {
}; };
public: public:
void advance() { triqs::tuple::fold(_aux1(), _c, false); } void advance() { _atend = ! (triqs::tuple::fold(_aux1(), _c, false)) ; }
// index_t index() const { return _index;} // not implemented yet // index_t index() const { return _index;} // not implemented yet
bool at_end() const { return _atend; } bool at_end() const { return _atend; }

View File

@ -116,7 +116,7 @@ namespace gfs {
p.reset(); p.reset();
return false; return false;
}; };
triqs::tuple::fold(l, _c, false); _atend = !(triqs::tuple::fold(l, _c, false));
} }
// index_t index() const { return _index;} // not implemented yet // index_t index() const { return _index;} // not implemented yet