Commit d6321012 authored by Gregor Michalicek's avatar Gregor Michalicek

Fixed segmentation fault in force calculation

parent bbabe6a3
......@@ -75,7 +75,7 @@ CONTAINS
IF (nlo(n)>0) THEN
mlo=sum(nlo(:n-1))+1
mlolo=dot_product(nlo(:n-1),nlo(:n-1)+2)/2+1
mlolo=dot_product(nlo(:n-1),nlo(:n-1)+1)/2+1
tuulo(:,:,mlo:(mlo+nlo(n)-1))=&
td_stored%tuulo(:size(tuulo,1),:size(tuulo,2),mlo:mlo+nlo(n)-1,jspin)
tdulo(:,:,mlo:(mlo+nlo(n)-1))=&
......
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