Commit 04a5197c authored by Florian Goth's avatar Florian Goth
Browse files

remove some more instances of old MMULTS

parent 2ee180f4
Pipeline #9583 passed with stage
in 57 seconds
......@@ -151,9 +151,11 @@
!Local
Integer :: nc, I, n
Complex (Kind=Kind(0.D0)) :: alpha, beta
Out = In
alpha = 1.D0
beta = 0.D0
do nc = 1,Ncheck
If ( dble( Op_T(nc,nf)%g*conjg(Op_T(nc,nf)%g) ) > Zero ) then
do I = 1,Ndim
......@@ -161,7 +163,8 @@
V_Hlp(n,I) = Out(Op_T(nc,nf)%P(n),I)
enddo
enddo
Call mmult(V_HLP1,Exp_T_m1(:,:,nc,nf, tstage), V_Hlp)
CALL ZGEMM('N', 'N', Ndim_hop, Ndim, Ndim_hop, alpha, Exp_T_m1(:,:,nc,nf, tstage), Ndim_hop,V_hlp,Ndim_hop,beta,V_HLP1,Ndim_hop)
! Call mmult(V_HLP1,Exp_T_m1(:,:,nc,nf, tstage), V_Hlp)
DO I = 1,Ndim
do n = 1,Ndim_hop
OUT(OP_T(nc,nf)%P(n),I) = V_hlp1(n,I)
......@@ -187,14 +190,18 @@
!Local
Integer :: nc, n
Complex (Kind=Kind(0.D0)) :: alpha, beta
Out = In
alpha = 1.D0
beta = 0.D0
do nc = 1, Ncheck
If ( dble( Op_T(nc,nf)%g*conjg(Op_T(nc,nf)%g) ) > Zero ) then
do n = 1,Ndim_hop
call zcopy(Ndim, Out(1, Op_T(nc,nf)%P(n)), 1, U_Hlp(1, n), 1)
enddo
Call mmult(U_Hlp1,U_Hlp,Exp_T(:,:,nc, nf, tstage))
call ZGEMM('N', 'N', Ndim, Ndim_hop, Ndim_hop, alpha, U_hlp, Ndim, Exp_T(:,:,nc, nf, tstage), Ndim_hop, beta, U_hlp1, Ndim)
! Call mmult(U_Hlp1,U_Hlp,Exp_T(:,:,nc, nf, tstage))
do n = 1,Ndim_hop
call zcopy(Ndim, U_hlp1(1, n), 1, OUT(1,OP_T(nc,nf)%P(n)), 1)
Enddo
......@@ -218,14 +225,18 @@
!Local
Integer :: nc, n
Complex (Kind=Kind(0.D0)) :: alpha, beta
Out = In
alpha = 1.D0
beta = 0.D0
do nc = Ncheck,1,-1
If ( dble( Op_T(nc,nf)%g*conjg(Op_T(nc,nf)%g) ) > Zero ) then
do n = 1,Ndim_hop
call zcopy(Ndim, Out(1, Op_T(nc,nf)%P(n)), 1, U_Hlp(1, n), 1)
enddo
Call mmult(U_Hlp1,U_Hlp,Exp_T_M1(:,:, nc, nf, tstage))
call ZGEMM('N', 'N', Ndim, Ndim_hop, Ndim_hop, alpha, U_Hlp, Ndim, Exp_T_M1(:,:, nc, nf, tstage), Ndim_hop, beta, U_hlp1, Ndim)
! Call mmult(U_Hlp1,U_Hlp,Exp_T_M1(:,:, nc, nf, tstage))
do n = 1, Ndim_hop
call zcopy(Ndim, U_Hlp1(1, n), 1, Out(1, Op_T(nc,nf)%P(n)), 1)
Enddo
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment