Commit 129cb555 authored by Gregor Michalicek's avatar Gregor Michalicek

Allocate vCoul%pw_w such that its number of spin dimensions is the number of spin dimensions

Due to the former reduced dimension of this array we had a nasty bug in writing it out.
parent d0c0e47e
......@@ -66,9 +66,10 @@ CONTAINS
#ifndef CPP_OLDINTEL
ALLOCATE(vTot%pw_w,mold=vTot%pw)
#else
ALLOCATE( vTot%pw_w(size(vTot%pw,1),size(vTot%pw,2)) )
ALLOCATE( vTot%pw_w(size(vTot%pw,1),size(vTot%pw,2)))
#endif
ALLOCATE(vCoul%pw_w(SIZE(den%pw,1),1))
ALLOCATE(vCoul%pw_w(SIZE(vCoul%pw,1),size(vCoul%pw,2)))
vCoul%pw_w = CMPLX(0.0,0.0)
CALL workDen%init(stars,atoms,sphhar,vacuum,noco,input%jspins,0)
......
Markdown is supported
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