LCPQ/qp2
9
1
Fork 0
mirror of https://github.com/QuantumPackage/qp2.git synced 2024-12-21 11:03:29 +01:00
Code Releases Activity

changed the test routines in dav_general_mat

Browse Source
This commit is contained in:
Emmanuel Giner 2021-07-02 16:48:20 +02:00
parent 9cdb127e33
commit ff87f67cb7
2 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
src/dav_general_mat/dav_ext_rout.irp.f
View File

@ -434,6 +434,9 @@ subroutine hcalc_template(v,u,N_st,sze)
v(i,istate) += H_matrix_all_dets(j,i) * u(j,istate)
enddo
enddo
do i = 1, sze
v(i,istate) += u(i,istate) * nuclear_repulsion
enddo
enddo
end

7
src/dav_general_mat/test_dav.irp.f
View File

@ -31,12 +31,13 @@ subroutine routine
do j = 1, sze
h_mat(j,i) = H_matrix_all_dets(j,i)
enddo
H_jj(i) = H_mat(i,i)
H_jj(i) = H_mat(i,i) + nuclear_repulsion
h_mat(i,i) = H_mat(i,i) + nuclear_repulsion
enddo
provide nthreads_davidson
call davidson_general(u_in,H_jj,energies,dim_in,sze,N_st,N_st_diag_in,converged,h_mat)
print*,'energies = ',energies + nuclear_repulsion
print*,'energies = ',energies
call davidson_general_ext_rout(u_in,H_jj,energies,dim_in,sze,N_st,N_st_diag_in,converged,hcalc_template)
print*,'energies = ',energies + nuclear_repulsion
print*,'energies = ',energies
end