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 2006-08-15 20:29:01 +00:00
parent 2c99357806
commit 47a1a24807
1 changed files with 8 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
nonclass.h
View File

@ -9,23 +9,20 @@
//MISC //MISC
export template <class T> export template <class T>
const NRSMat<T> twoside_transform(const NRSMat<T> &S, const NRMat<T> &C) //calculate C^dagger S C const NRSMat<T> twoside_transform(const NRSMat<T> &S, const NRMat<T> &C, bool transp=0) //calculate C^dagger S C
{ {
NRMat<T> tmp = S * C; if(transp)
NRMat<T> result(C.ncols(),C.ncols());
result.gemm((T)0,C,'t',tmp,'n',(T)1);
return NRSMat<T>(result);
}
export template <class T>
const NRSMat<T> twoside_transform_t(const NRSMat<T> &S, const NRMat<T> &C) //calculate C S C^dagger
{ {
NRMat<T> tmp = C * S; NRMat<T> tmp = C * S;
NRMat<T> result(C.rows(),C.nrows()); NRMat<T> result(C.rows(),C.nrows());
result.gemm((T)0,tmp,'n',C,'t',(T)1); result.gemm((T)0,tmp,'n',C,'t',(T)1);
return NRSMat<T>(result); return NRSMat<T>(result);
} }
NRMat<T> tmp = S * C;
NRMat<T> result(C.ncols(),C.ncols());
result.gemm((T)0,C,'t',tmp,'n',(T)1);
return NRSMat<T>(result);
}