diff --git a/cuda_la.cc b/cuda_la.cc index b47f7ac..ea960fa 100644 --- a/cuda_la.cc +++ b/cuda_la.cc @@ -48,15 +48,15 @@ double *gpuputdouble(const double &x){ return (double *)ptr; } -complex *gpuputcomplex(const complex &x){ +std::complex *gpuputcomplex(const std::complex &x){ void *ptr = NULL; - cublasAlloc(1, sizeof(complex), &ptr); + cublasAlloc(1, sizeof(std::complex), &ptr); TEST_CUBLAS("cublasAlloc"); - cublasSetVector(1, sizeof(complex), &x, 1, ptr, 1); + cublasSetVector(1, sizeof(std::complex), &x, 1, ptr, 1); TEST_CUBLAS("cublasSetVector"); - return (complex *)ptr; + return (std::complex *)ptr; } diff --git a/noncblas.cc b/noncblas.cc index 320efe6..32e4489 100644 --- a/noncblas.cc +++ b/noncblas.cc @@ -255,8 +255,8 @@ void cblas_zgerc(const CBLAS_ORDER Order, const int M, const int N, { if(Order!=CblasRowMajor) laerror("CblasRowMajor order asserted"); //conjugate y -complex *p; -p= (complex *) Y; for(int i=0; iimag() = -p->imag(); p+= incY;} +std::complex *p; +p= (std::complex *) Y; for(int i=0; iimag(-p->imag()); p+= incY;} #ifdef FORINT const FINT mtmp=M; const FINT ntmp=N; @@ -268,7 +268,7 @@ p= (complex *) Y; for(int i=0; iimag() = -p->imag(); p+= in FORNAME(zgerc)(&N, &M, alpha, Y, &incY, X, &incX, A, &lda); #endif //conjugate y back -p= (complex *) Y; for(int i=0; iimag() = -p->imag(); p+= incY;} +p= (std::complex *) Y; for(int i=0; iimag(--p->imag()); p+= incY;} }