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 2008-07-15 14:28:38 +00:00
parent 85913e8d69
commit f32be3195a
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
fourindex.h
View File

@ -680,12 +680,12 @@ public:
const T& operator() (unsigned int i, unsigned int j, unsigned int k, unsigned int l) const; const T& operator() (unsigned int i, unsigned int j, unsigned int k, unsigned int l) const;
void resize(const int n) {(*this).NRSMat<T>::resize(n*(n+1)/2);}; void resize(const int n) {(*this).NRSMat<T>::resize(n*(n+1)/2);};
void putext(int f, T thr=1e-15); void putext(int f, T thr=1e-15);
int nbas() const {return (int)sqrt(2*nn);}; int nbas() const {return (int)sqrt(2*(*this).nrows());};
}; };
template<class T, class I> template<class T, class I>
fourindex_dense<twoelectronrealmullikan,T,I>::putext(int f, T thr) void fourindex_dense<twoelectronrealmullikan,T,I>::putext(int f, T thr)
{ {
T y; T y;
for(int i=1; i<=nbas(); ++i) for(int j=1; j<=i; ++j) for(int i=1; i<=nbas(); ++i) for(int j=1; j<=i; ++j)
@ -698,7 +698,7 @@ for(int i=1; i<=nbas(); ++i) for(int j=1; j<=i; ++j)
x.index.indiv.j=j; x.index.indiv.j=j;
x.index.indiv.k=k; x.index.indiv.k=k;
x.index.indiv.l=l; x.index.indiv.l=l;
if(sizeof(matel4stored<I,T>) != write(f,&x,sizeof(matel4stored<I,T>)) if(sizeof(matel4stored<I,T>) != write(f,&x,sizeof(matel4stored<I,T>)) )
laerror("write error in putext"); laerror("write error in putext");
} }
} }