jiri/LA_library
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

*** empty log message ***

Browse Source
This commit is contained in:
jiri 2015-03-05 15:18:46 +00:00
parent 80fe44fab2
commit 8ed1014ddc
1 changed files with 23 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
fourindex.h
View File

@ -290,6 +290,7 @@ protected:
public: public:
void resize(I n) {nn=n;}
fourindex_ext(const int file, const fourindexsymtype s=undefined_symmetry, const I n=0, const unsigned int b=1024) :current(NULL),fd(file),nread(0),symmetry(s),nn(n) fourindex_ext(const int file, const fourindexsymtype s=undefined_symmetry, const I n=0, const unsigned int b=1024) :current(NULL),fd(file),nread(0),symmetry(s),nn(n)
{ {
struct statfs sfs; struct statfs sfs;
@ -648,6 +649,28 @@ std::istream& operator>>(std::istream &s, fourindex<I,T> &x)
return s; return s;
} }
template <class I, class T>
std::istream& operator>>(std::istream &s, fourindex_ext<I,T> &x)
{
int n;
s >> n;
x.resize(n);
typename LA_traits_io<I>::IOtype i,j,k,l;
typename LA_traits_io<T>::IOtype elem;
s >> i >> j >>k >>l;
while(i!= (typename LA_traits_io<I>::IOtype)-1 && j!= (typename LA_traits_io<I>::IOtype)-1 && k != (typename LA_traits_io<I>::IOtype)-1 && l!= (typename LA_traits_io<I>::IOtype)-1)
{
s>>elem;
x.put((I)i,(I)j,(I)k,(I)l,(T)elem);
s >> i >> j >>k >>l;
}
x.flush();
return s;
}
/////////////////////densely stored fourindex/////////////////////////////////// /////////////////////densely stored fourindex///////////////////////////////////
//not all symmetry cases implemented yet, but a general template declaration used //not all symmetry cases implemented yet, but a general template declaration used