gallo/atrip
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Fix some complex pointer errors

Browse Source
This commit is contained in:
Alejandro Gallo 2022-07-18 19:27:27 +02:00
parent e8b733cb1a
commit 1392e8dc36
1 changed files with 8 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
atrip.org
View File

@ -2996,7 +2996,10 @@ V^{{\color{blue}ab}}_{{\color{red}e}i} T^{{\color{blue}c}{\color{red}e}}_{ij} \
} \ } \
) )
#if defined(HAVE_CUDA) #if defined(HAVE_CUDA)
#define __TO_DEVICEPTR(_v) (DataFieldType<F>*)(CUdeviceptr)thrust::raw_pointer_cast((_v)) #define __TO_DEVICEPTR(_v) \
((DataFieldType<F>*) \
(CUdeviceptr) \
thrust::raw_pointer_cast((_v)))
#define DGEMM_PARTICLES(__A, __B) \ #define DGEMM_PARTICLES(__A, __B) \
atrip::xgemm<F>("T", \ atrip::xgemm<F>("T", \
"N", \ "N", \
@ -3028,8 +3031,8 @@ V^{{\color{blue}ab}}_{{\color{red}e}i} T^{{\color{blue}c}{\color{red}e}}_{ij} \
); );
#define MAYBE_CONJ(_conj, _buffer) \ #define MAYBE_CONJ(_conj, _buffer) \
for (size_t __i = 0; __i < NoNoNo; ++__i) \ for (size_t __i = 0; __i < NoNoNo; ++__i) \
_conj[__i] = \ __TO_DEVICEPTR(_conj.data())[__i] = \
maybeConjugate<DataFieldType<F>>(((DataFieldType<F>*)_buffer)[__i]); maybeConjugate<DataFieldType<F>>(((DataFieldType<F>*)(_buffer))[__i]);
#else #else
#define __TO_DEVICEPTR(_v) (_v) #define __TO_DEVICEPTR(_v) (_v)
#define DGEMM_PARTICLES(__A, __B) \ #define DGEMM_PARTICLES(__A, __B) \
@ -3087,7 +3090,8 @@ V^{{\color{blue}ab}}_{{\color{red}e}i} T^{{\color{blue}c}{\color{red}e}}_{ij} \
{ // Holes part %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% { // Holes part %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#ifdef HAVE_CUDA #ifdef HAVE_CUDA
thrust::device_vector< DataFieldType<F> > _vhhh(NoNoNo); thrust::device_vector< DataFieldType<F> > _vhhh;
_vhhh.reserve(NoNoNo);
#else #else
std::vector<F> _vhhh(NoNoNo); std::vector<F> _vhhh(NoNoNo);
#endif #endif