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

realpart and imagpart for Vec and SMat

Browse Source
This commit is contained in:
Jiri Pittner 2025-08-04 17:09:29 +02:00
parent 4fdb1b3ab7
commit 92f7815eef
1 changed files with 72 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

72
nonclass.cc
View File

@ -1407,6 +1407,78 @@ const NRMat<double> imagpart<NRMat< std::complex<double> > >(const NRMat< std::c
return result; return result;
} }
template<>
const NRSMat<double> realpart<NRSMat< std::complex<double> > >(const NRSMat< std::complex<double> > &a)
{
NRSMat<double> result(a.nrows());
#ifdef CUDALA
if(a.location == cpu){
#endif
cblas_dcopy(a.nrows()*(a.ncols()+1)/2, (const double *)(&a(0,0)), 2, result, 1);
#ifdef CUDALA
}else{
laerror("not implemented for cuda yet");
}
#endif
return result;
}
template<>
const NRSMat<double> imagpart<NRSMat< std::complex<double> > >(const NRSMat< std::complex<double> > &a)
{
NRSMat<double> result(a.nrows());
#ifdef CUDALA
if(a.location == cpu){
#endif
cblas_dcopy(a.nrows()*(a.ncols()+1)/2, (const double *)(&a(0,0))+1, 2, result, 1);
#ifdef CUDALA
}else{
laerror("not implemented for cuda yet");
}
#endif
return result;
}
template<>
const NRVec<double> realpart<NRVec< std::complex<double> > >(const NRVec< std::complex<double> > &a)
{
NRVec<double> result(a.size());
#ifdef CUDALA
if(a.location == cpu){
#endif
cblas_dcopy(a.size(), (const double *)&a[0], 2, result, 1);
#ifdef CUDALA
}else{
laerror("not implemented for cuda yet");
}
#endif
return result;
}
template<>
const NRVec<double> imagpart<NRVec< std::complex<double> > >(const NRVec< std::complex<double> > &a)
{
NRVec<double> result(a.size());
#ifdef CUDALA
if(a.location == cpu){
#endif
cblas_dcopy(a.size(), (const double *)&a[0]+1, 2, result, 1);
#ifdef CUDALA
}else{
laerror("not implemented for cuda yet");
}
#endif
return result;
}
template<> template<>
const NRMat< std::complex<double> > realmatrix<NRMat<double> > (const NRMat<double> &a) const NRMat< std::complex<double> > realmatrix<NRMat<double> > (const NRMat<double> &a)
{ {