Forgot to change the memalign calls

Makefile.am

@@ -549,6 +549,11 @@ elpa2_test_real_default_kernel_qr_decomposition@SUFFIX@.sh:
 	echo '$(wrapper)./elpa2_test_real_default_kernel_qr_decomposition@SUFFIX@' > $@
 	chmod +x $@
 
+if WANT_SINGLE_PRECISION_REAL
+elpa2_test_real_default_kernel_qr_decomposition_single_precision@SUFFIX@.sh:
+	echo '$(wrapper)./elpa2_test_real_default_kernel_qr_decomposition_single_precision@SUFFIX@' > $@
+	chmod +x $@
+endif
 # Preprocessed files (just used for manual inspection)
 mod_precision.i: $(top_srcdir)/src/mod_precision.F90
 	$(CPP) $(CPPFLAGS) -I$(top_builddir)/ -I$(top_srcdir)/ -c $(top_srcdir)/src/mod_precision.F90 -o $@

src/elpa2_compute_complex_template.X90

@@ -2999,6 +2999,7 @@
 #else
       complex(kind=COMPLEX_DATATYPE), pointer     :: a(:,:,:)
 #endif
+      complex(kind=COMPLEX_DATATYPE)              :: a_complex
       type(c_ptr)                                 :: a_ptr
       complex(kind=COMPLEX_DATATYPE), allocatable :: row(:)
       complex(kind=COMPLEX_DATATYPE), allocatable :: row_group(:,:)
@@ -3224,7 +3225,7 @@
       endif
 
       if (.not.(useGPU)) then
-         if (posix_memalign(a_ptr, 64_C_SIZE_T, stripe_width*a_dim2*stripe_count*max_threads*C_SIZEOF(a(1,1,1,1))) /= 0) then
+         if (posix_memalign(a_ptr, 64_C_SIZE_T, stripe_width*a_dim2*stripe_count*max_threads*C_SIZEOF(a_complex)) /= 0) then
           print *,"trans_ev_tridi_to_band_complex: error allocating a "//errorMessage
           stop
         endif
@@ -3238,7 +3239,7 @@
 #else /* OpenMP */
 
       if (.not.(useGPU)) then
-        if (posix_memalign(a_ptr, 64_C_SIZE_T, stripe_width*a_dim2*stripe_count*C_SIZEOF(a(1,1,1))) /= 0) then
+        if (posix_memalign(a_ptr, 64_C_SIZE_T, stripe_width*a_dim2*stripe_count*C_SIZEOF(a_complex)) /= 0) then
           print *,"trans_ev_tridi_to_band_complex: error allocating a "//errorMessage
           stop
         endif

src/elpa2_compute_real_template.X90

@@ -3414,6 +3414,7 @@
 #else
       real(kind=REAL_DATATYPE), pointer        :: a(:,:,:)
 #endif
+      real(kind=REAL_DATATYPE)                 :: a_real
       type(c_ptr)                              :: a_ptr
       real(kind=REAL_DATATYPE)   , allocatable :: row(:)
       real(kind=REAL_DATATYPE)   , allocatable :: row_group(:,:)
@@ -3608,7 +3609,7 @@
 #endif
 
 #ifdef WITH_OPENMP
-        if (posix_memalign(a_ptr, 64_C_SIZE_T, stripe_width*a_dim2*stripe_count*max_threads*C_SIZEOF(a(1,1,1,1))) /= 0) then
+        if (posix_memalign(a_ptr, 64_C_SIZE_T, stripe_width*a_dim2*stripe_count*max_threads*C_SIZEOF(a_real)) /= 0) then
           print *,"trans_ev_tridi_to_band_real: error when allocating a"//errorMessage
           stop
         endif
@@ -3618,7 +3619,7 @@
 
         ! a(:,:,:,:) should be set to 0 in a parallel region, not here!
 #else
-        if (posix_memalign(a_ptr, 64_C_SIZE_T, stripe_width*a_dim2*stripe_count*C_SIZEOF(a(1,1,1))) /= 0) then
+        if (posix_memalign(a_ptr, 64_C_SIZE_T, stripe_width*a_dim2*stripe_count*C_SIZEOF(a_real)) /= 0) then
           print *,"trans_ev_tridi_to_band_real: error when allocating a"//errorMessage
           stop
         endif