diff --git a/ELPA_2014.06/src/elpa2.F90 b/ELPA_2014.06/src/elpa2.F90 index 5791a625ea682f36214055ee423d0fa0b2b10521..6e056fe521f21d5720ac0766d4beedc40fa19b84 100644 --- a/ELPA_2014.06/src/elpa2.F90 +++ b/ELPA_2014.06/src/elpa2.F90 @@ -3884,7 +3884,8 @@ subroutine tridiag_band_complex(na, nb, nblk, a, lda, d, e, mpi_comm_rows, mpi_c h(2:nb) = h(2:nb) - x*hv(2:nb) ! Unfortunately there is no BLAS routine like DSYR2 for a nonsymmetric rank 2 update ("DGER2") do i=2,nb - ab(2+nb-i:1+nb+nr-i,i+ns-1) = ab(2+nb-i:1+nb+nr-i,i+ns-1) - hv_t(1:nr,my_thread)*conjg(h(i)) - hs(1:nr)*conjg(hv(i)) + ab(2+nb-i:1+nb+nr-i,i+ns-1) = ab(2+nb-i:1+nb+nr-i,i+ns-1) & + - hv_t(1:nr,my_thread)*conjg(h(i)) - hs(1:nr)*conjg(hv(i)) enddo else