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

twoside_transform for complex and for non-symetric case

Browse Source
This commit is contained in:
Jiri Pittner
2025-11-26 16:57:21 +01:00
parent 3da3efa57c
commit b556b06442
1 changed files with 24 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

27
nonclass.h
View File

@@ -68,22 +68,43 @@ return positive_power(z,i);
template <class T> template <class T>
const NRSMat<T> twoside_transform(const NRSMat<T> &S, const NRMat<T> &C, bool transp=0) //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 for symmetric/hermitian S
{ {
char tchar='t';
if(LA_traits<T>::is_complex()) tchar='c';
if(transp) if(transp)
{ {
NRMat<T> tmp = C * S; NRMat<T> tmp = C * S;
NRMat<T> result(C.nrows(),C.nrows()); NRMat<T> result(C.nrows(),C.nrows());
result.gemm((T)0,tmp,'n',C,'t',(T)1); result.gemm((T)0,tmp,'n',C,tchar,(T)1);
return NRSMat<T>(result); return NRSMat<T>(result);
} }
NRMat<T> tmp = S * C; NRMat<T> tmp = S * C;
NRMat<T> result(C.ncols(),C.ncols()); NRMat<T> result(C.ncols(),C.ncols());
result.gemm((T)0,C,'t',tmp,'n',(T)1); result.gemm((T)0,C,tchar,tmp,'n',(T)1);
return NRSMat<T>(result); return NRSMat<T>(result);
} }
template <class T>
const NRMat<T> twoside_transform(const NRMat<T> &S, const NRMat<T> &C, bool transp=0) //calculate C^dagger S C for non-symmetric S
{
char tchar='t';
if(LA_traits<T>::is_complex()) tchar='c';
if(transp)
{
NRMat<T> tmp = C * S;
NRMat<T> result(C.nrows(),C.nrows());
result.gemm((T)0,tmp,'n',C,tchar,(T)1);
return result;
}
NRMat<T> tmp = S * C;
NRMat<T> result(C.ncols(),C.ncols());
result.gemm((T)0,C,tchar,tmp,'n',(T)1);
return result;
}
template <class T> template <class T>