10
0
mirror of https://github.com/QuantumPackage/qp2.git synced 2024-06-26 15:12:19 +02:00

fixed floating point exception in AOs

This commit is contained in:
Emmanuel Giner 2020-06-08 11:50:41 +02:00
parent 125294a05a
commit 47528188b5

View File

@ -97,6 +97,7 @@ subroutine give_all_aos_at_r(r,aos_array)
dz2 = dz**power_ao(3)
do l = 1,ao_prim_num(k)
beta = ao_expo_ordered_transp_per_nucl(l,j,i)
if(dabs(beta*r2).gt.40.d0)cycle
aos_array(k)+= ao_coef_normalized_ordered_transp_per_nucl(l,j,i) * dexp(-beta*r2)
enddo
aos_array(k) = aos_array(k) * dx2 * dy2 * dz2