tdhf_gradient Subroutine

public subroutine tdhf_gradient(infos)

Uses

  • proc~~tdhf_gradient~~UsesGraph proc~tdhf_gradient tdhf_gradient module~basis_tools basis_tools proc~tdhf_gradient->module~basis_tools module~dft dft proc~tdhf_gradient->module~dft module~grd1 grd1 proc~tdhf_gradient->module~grd1 module~io_constants io_constants proc~tdhf_gradient->module~io_constants module~mathlib mathlib proc~tdhf_gradient->module~mathlib module~messages messages proc~tdhf_gradient->module~messages module~mod_dft_gridint_tdxc_grad mod_dft_gridint_tdxc_grad proc~tdhf_gradient->module~mod_dft_gridint_tdxc_grad module~mod_dft_molgrid mod_dft_molgrid proc~tdhf_gradient->module~mod_dft_molgrid module~oqp_tagarray_driver oqp_tagarray_driver proc~tdhf_gradient->module~oqp_tagarray_driver module~printing printing proc~tdhf_gradient->module~printing module~strings strings proc~tdhf_gradient->module~strings module~tdhf_lib tdhf_lib proc~tdhf_gradient->module~tdhf_lib module~types types proc~tdhf_gradient->module~types module~util util proc~tdhf_gradient->module~util module~basis_tools->module~io_constants iso_fortran_env iso_fortran_env module~basis_tools->iso_fortran_env module~atomic_structure_m atomic_structure_m module~basis_tools->module~atomic_structure_m module~constants constants module~basis_tools->module~constants module~parallel parallel module~basis_tools->module~parallel module~precision precision module~basis_tools->module~precision module~dft->module~basis_tools module~dft->module~io_constants module~dft->module~messages module~dft->module~mod_dft_molgrid module~dft->module~precision module~grd1->module~basis_tools module~grd1->module~io_constants module~grd1->module~mathlib module~grd1->module~types module~grd1->module~atomic_structure_m module~ecp_tool ecp_tool module~grd1->module~ecp_tool module~mod_1e_primitives mod_1e_primitives module~grd1->module~mod_1e_primitives module~mod_shell_tools mod_shell_tools module~grd1->module~mod_shell_tools module~grd1->module~precision module~oqp_linalg oqp_linalg module~mathlib->module~oqp_linalg module~mathlib->module~precision module~messages->module~io_constants comm_IOFILE comm_IOFILE module~messages->comm_IOFILE comm_PAR comm_PAR module~messages->comm_PAR module~messages->module~precision module~mod_dft_gridint mod_dft_gridint module~mod_dft_gridint_tdxc_grad->module~mod_dft_gridint module~mod_dft_gridint_tdxc_grad->module~precision module~bragg_slater_radii bragg_slater_radii module~mod_dft_molgrid->module~bragg_slater_radii module~lebedev lebedev module~mod_dft_molgrid->module~lebedev module~mod_grid_storage mod_grid_storage module~mod_dft_molgrid->module~mod_grid_storage module~mod_dft_molgrid->module~precision iso_c_binding iso_c_binding module~oqp_tagarray_driver->iso_c_binding tagarray tagarray module~oqp_tagarray_driver->tagarray module~printing->module~precision module~strings->iso_c_binding module~tdhf_lib->module~basis_tools module~int2_compute int2_compute module~tdhf_lib->module~int2_compute module~tdhf_lib->module~oqp_linalg module~tdhf_lib->module~precision module~types->module~basis_tools module~types->iso_c_binding module~types->module~atomic_structure_m module~functionals functionals module~types->module~functionals module~types->module~parallel module~types->module~precision module~types->tagarray module~util->module~precision module~atomic_structure_m->iso_c_binding module~bragg_slater_radii->module~precision module~constants->module~precision module~ecp_tool->module~basis_tools module~ecp_tool->iso_c_binding module~ecp_tool->iso_fortran_env module~ecp_tool->module~precision module~libecp_result libecp_result module~ecp_tool->module~libecp_result module~libecpint_wrapper libecpint_wrapper module~ecp_tool->module~libecpint_wrapper module~functionals->iso_c_binding module~functionals->module~precision xc_f03_lib_m xc_f03_lib_m module~functionals->xc_f03_lib_m module~int2_compute->module~basis_tools module~int2_compute->module~messages module~int2_compute->module~atomic_structure_m module~int2_compute->module~parallel module~int2_compute->module~precision module~int2_pairs int2_pairs module~int2_compute->module~int2_pairs module~int2e_libint int2e_libint module~int2_compute->module~int2e_libint module~int2e_rys int2e_rys module~int2_compute->module~int2e_rys module~lebedev->module~precision module~mod_1e_primitives->iso_fortran_env module~mod_1e_primitives->module~constants module~mod_1e_primitives->module~mod_shell_tools module~mod_gauss_hermite mod_gauss_hermite module~mod_1e_primitives->module~mod_gauss_hermite module~rys rys module~mod_1e_primitives->module~rys module~xyz_order xyz_order module~mod_1e_primitives->module~xyz_order module~mod_dft_gridint->module~basis_tools module~mod_dft_gridint->module~io_constants module~mod_dft_gridint->module~mod_dft_molgrid module~mod_dft_gridint->module~functionals module~mod_dft_gridint->module~oqp_linalg module~mod_dft_gridint->module~parallel module~mod_dft_gridint->module~precision module~mod_dft_xc_libxc mod_dft_xc_libxc module~mod_dft_gridint->module~mod_dft_xc_libxc module~mod_grid_storage->module~precision module~mod_shell_tools->module~basis_tools module~mod_shell_tools->module~precision module~blas_wrap blas_wrap module~oqp_linalg->module~blas_wrap module~lapack_wrap lapack_wrap module~oqp_linalg->module~lapack_wrap module~parallel->iso_c_binding module~parallel->iso_fortran_env module~parallel->module~precision mpi mpi module~parallel->mpi module~precision->iso_fortran_env module~blas_wrap->module~messages module~blas_wrap->module~precision module~mathlib_types mathlib_types module~blas_wrap->module~mathlib_types module~int2_pairs->module~precision module~int2e_libint->iso_c_binding module~int2e_libint->module~constants module~int2e_libint->module~precision module~int2e_libint->module~int2_pairs module~libint_f libint_f module~int2e_libint->module~libint_f module~int2e_rys->module~basis_tools module~int2e_rys->module~constants module~int2e_rys->module~precision module~lapack_wrap->module~messages module~lapack_wrap->module~precision module~lapack_wrap->module~mathlib_types module~libecp_result->iso_c_binding module~libecpint_wrapper->iso_c_binding module~mod_dft_xc_libxc->module~precision module~mod_dft_xclib mod_dft_xclib module~mod_dft_xc_libxc->module~mod_dft_xclib module~mod_gauss_hermite->module~precision module~rys->module~constants module~rys->module~precision module~rys_lut rys_lut module~rys->module~rys_lut module~libint_f->iso_c_binding module~mod_dft_xclib->module~functionals module~mod_dft_xclib->module~precision

Arguments

Type IntentOptional Attributes Name
type(information), intent(inout), target :: infos

Calls

proc~~tdhf_gradient~~CallsGraph proc~tdhf_gradient tdhf_gradient interface~data_has_tags data_has_tags proc~tdhf_gradient->interface~data_has_tags interface~show_message show_message proc~tdhf_gradient->interface~show_message interface~tagarray_get_data tagarray_get_data proc~tdhf_gradient->interface~tagarray_get_data interface~unpack_matrix unpack_matrix proc~tdhf_gradient->interface~unpack_matrix proc~dft_initialize dft_initialize proc~tdhf_gradient->proc~dft_initialize proc~dftclean dftclean proc~tdhf_gradient->proc~dftclean proc~iatogen iatogen proc~tdhf_gradient->proc~iatogen proc~measure_time measure_time proc~tdhf_gradient->proc~measure_time proc~orthogonal_transform orthogonal_transform proc~tdhf_gradient->proc~orthogonal_transform proc~print_gradient print_gradient proc~tdhf_gradient->proc~print_gradient proc~print_module_info print_module_info proc~tdhf_gradient->proc~print_module_info proc~symmetrize_matrix symmetrize_matrix proc~tdhf_gradient->proc~symmetrize_matrix proc~tddft_xc_gradient tddft_xc_gradient proc~tdhf_gradient->proc~tddft_xc_gradient proc~tdhf_1e_grad tdhf_1e_grad proc~tdhf_gradient->proc~tdhf_1e_grad proc~tdhf_2e_grad tdhf_2e_grad proc~tdhf_gradient->proc~tdhf_2e_grad proc~unpack_f90 UNPACK_F90 interface~unpack_matrix->proc~unpack_f90 proc~dft_initialize->interface~show_message none~add_atomic_grid dft_grid_t%add_atomic_grid proc~dft_initialize->none~add_atomic_grid none~add_grid sorted_grid_t%add_grid proc~dft_initialize->none~add_grid none~compress dft_grid_t%compress proc~dft_initialize->none~compress none~find_neighbours dft_grid_t%find_neighbours proc~dft_initialize->none~find_neighbours none~reset dft_grid_t%reset proc~dft_initialize->none~reset none~set_screening basis_set%set_screening proc~dft_initialize->none~set_screening proc~c_f_char c_f_char proc~dft_initialize->proc~c_f_char proc~dft_fc_blk dft_fc_blk proc~dft_initialize->proc~dft_fc_blk proc~get_radial_grid get_radial_grid proc~dft_initialize->proc~get_radial_grid proc~libxc_input libxc_input proc~dft_initialize->proc~libxc_input proc~set_bragg_slater set_bragg_slater proc~dft_initialize->proc~set_bragg_slater proc~libxc_destroy libxc_destroy proc~dftclean->proc~libxc_destroy proc~orthogonal_transform->interface~show_message proc~oqp_dgemm_i64 oqp_dgemm_i64 proc~orthogonal_transform->proc~oqp_dgemm_i64 none~clean~23 xc_consumer_tdg_t%clean proc~tddft_xc_gradient->none~clean~23 none~init~14 par_env_t%init proc~tddft_xc_gradient->none~init~14 proc~run_xc run_xc proc~tddft_xc_gradient->proc~run_xc proc~tdhf_1e_grad->interface~data_has_tags proc~tdhf_1e_grad->interface~show_message proc~tdhf_1e_grad->interface~tagarray_get_data proc~eijden eijden proc~tdhf_1e_grad->proc~eijden proc~grad_1e_ecp grad_1e_ecp proc~tdhf_1e_grad->proc~grad_1e_ecp proc~grad_ee_kinetic grad_ee_kinetic proc~tdhf_1e_grad->proc~grad_ee_kinetic proc~grad_ee_overlap grad_ee_overlap proc~tdhf_1e_grad->proc~grad_ee_overlap proc~grad_en_hellman_feynman grad_en_hellman_feynman proc~tdhf_1e_grad->proc~grad_en_hellman_feynman proc~grad_en_pulay grad_en_pulay proc~tdhf_1e_grad->proc~grad_en_pulay proc~grad_nn grad_nn proc~tdhf_1e_grad->proc~grad_nn proc~tdhf_2e_grad->interface~show_message clean clean proc~tdhf_2e_grad->clean init init proc~tdhf_2e_grad->init proc~grd2_driver grd2_driver proc~tdhf_2e_grad->proc~grd2_driver none~add_slices dft_grid_t%add_slices none~add_atomic_grid->none~add_slices none~get~5 list_grid_3d_t%get none~add_grid->none~get~5 none~set~4 list_grid_3d_t%set none~add_grid->none~set~4 proc~lebedev_get_grid lebedev_get_grid none~add_grid->proc~lebedev_get_grid none~set~2 partition_function%set none~find_neighbours->none~set~2 mpi_comm_rank mpi_comm_rank none~init~14->mpi_comm_rank mpi_comm_size mpi_comm_size none~init~14->mpi_comm_size none~init~21 sorted_grid_t%init none~reset->none~init~21 none~getbyid list_grid_3d_t%getByID proc~dft_fc_blk->none~getbyid proc~dft_fc_blk->none~set~2 proc~eijden->interface~data_has_tags proc~eijden->interface~tagarray_get_data proc~eijden->interface~unpack_matrix proc~orthogonal_transform_sym orthogonal_transform_sym proc~eijden->proc~orthogonal_transform_sym set set proc~get_radial_grid->set transform transform proc~get_radial_grid->transform proc~grad_1e_ecp->none~init~14 none~bcast par_env_t%bcast proc~grad_1e_ecp->none~bcast proc~add_ecpder add_ecpder proc~grad_1e_ecp->proc~add_ecpder proc~grad_ee_kinetic->interface~unpack_matrix interface~bas_norm_matrix bas_norm_matrix proc~grad_ee_kinetic->interface~bas_norm_matrix none~alloc~2 shpair_t%alloc proc~grad_ee_kinetic->none~alloc~2 none~fetch_by_id shell_t%fetch_by_id proc~grad_ee_kinetic->none~fetch_by_id none~shell_pair shpair_t%shell_pair proc~grad_ee_kinetic->none~shell_pair proc~comp_kinetic_der1 comp_kinetic_der1 proc~grad_ee_kinetic->proc~comp_kinetic_der1 proc~grad_ee_overlap->interface~unpack_matrix proc~grad_ee_overlap->interface~bas_norm_matrix proc~grad_ee_overlap->none~alloc~2 proc~grad_ee_overlap->none~fetch_by_id proc~grad_ee_overlap->none~shell_pair proc~comp_overlap_der1 comp_overlap_der1 proc~grad_ee_overlap->proc~comp_overlap_der1 proc~grad_en_hellman_feynman->interface~unpack_matrix proc~grad_en_hellman_feynman->interface~bas_norm_matrix proc~grad_en_hellman_feynman->none~alloc~2 proc~grad_en_hellman_feynman->none~fetch_by_id proc~grad_en_hellman_feynman->none~shell_pair proc~comp_coulomb_helfeyder1 comp_coulomb_helfeyder1 proc~grad_en_hellman_feynman->proc~comp_coulomb_helfeyder1 proc~grad_en_pulay->interface~unpack_matrix proc~grad_en_pulay->interface~bas_norm_matrix proc~grad_en_pulay->none~alloc~2 proc~grad_en_pulay->none~fetch_by_id proc~grad_en_pulay->none~shell_pair proc~comp_coulomb_der1 comp_coulomb_der1 proc~grad_en_pulay->proc~comp_coulomb_der1 proc~grd2_driver->interface~show_message proc~grd2_driver->none~init~14 get_density get_density proc~grd2_driver->get_density none~alloc int2_pair_storage%alloc proc~grd2_driver->none~alloc none~allreduce par_env_t%allreduce proc~grd2_driver->none~allreduce none~clean~24 grd2_int_data_t%clean proc~grd2_driver->none~clean~24 none~compute int2_pair_storage%compute proc~grd2_driver->none~compute none~set int2_cutoffs_t%set proc~grd2_driver->none~set none~set_ids~2 grd2_int_data_t%set_ids proc~grd2_driver->none~set_ids~2 proc~grd2_rys_compute grd2_rys_compute proc~grd2_driver->proc~grd2_rys_compute proc~ints_exchange ints_exchange proc~grd2_driver->proc~ints_exchange none~destroy~2 functional_t%destroy proc~libxc_destroy->none~destroy~2 proc~libxc_input->interface~show_message none~add_functional functional_t%add_functional proc~libxc_input->none~add_functional none~can_calculate functional_t%can_calculate proc~libxc_input->none~can_calculate xc_f03_reference xc_f03_reference proc~libxc_input->xc_f03_reference xc_f03_reference_doi xc_f03_reference_doi proc~libxc_input->xc_f03_reference_doi xc_f03_version_string xc_f03_version_string proc~libxc_input->xc_f03_version_string proc~oqp_dgemm_i64->interface~show_message dgemm dgemm proc~oqp_dgemm_i64->dgemm proc~run_xc->none~allreduce none~compaos xc_engine_t%compAOs proc~run_xc->none~compaos none~compxc xc_engine_t%compXC proc~run_xc->none~compxc none~getslicenonzero dft_grid_t%getSliceNonZero proc~run_xc->none~getslicenonzero none~getstats xc_engine_t%getStats proc~run_xc->none~getstats none~init~18 xc_engine_t%init proc~run_xc->none~init~18 none~pruneaos xc_engine_t%pruneAOs proc~run_xc->none~pruneaos none~resetpointers xc_engine_t%resetPointers proc~run_xc->none~resetpointers parallel_start parallel_start proc~run_xc->parallel_start parallel_stop parallel_stop proc~run_xc->parallel_stop postUpdate postUpdate proc~run_xc->postUpdate update update proc~run_xc->update proc~unpack_f90->interface~show_message proc~oqp_dtpttr_i64 oqp_dtpttr_i64 proc~unpack_f90->proc~oqp_dtpttr_i64 none~add_functional->interface~show_message xc_f03_func_get_info xc_f03_func_get_info none~add_functional->xc_f03_func_get_info xc_f03_func_info_get_family xc_f03_func_info_get_family none~add_functional->xc_f03_func_info_get_family xc_f03_func_info_get_flags xc_f03_func_info_get_flags none~add_functional->xc_f03_func_info_get_flags xc_f03_func_info_get_kind xc_f03_func_info_get_kind none~add_functional->xc_f03_func_info_get_kind xc_f03_func_info_get_name xc_f03_func_info_get_name none~add_functional->xc_f03_func_info_get_name xc_f03_func_info_get_references xc_f03_func_info_get_references none~add_functional->xc_f03_func_info_get_references xc_f03_func_init xc_f03_func_init none~add_functional->xc_f03_func_init xc_f03_func_reference_get_doi xc_f03_func_reference_get_doi none~add_functional->xc_f03_func_reference_get_doi xc_f03_func_reference_get_ref xc_f03_func_reference_get_ref none~add_functional->xc_f03_func_reference_get_ref xc_f03_func_set_ext_params xc_f03_func_set_ext_params none~add_functional->xc_f03_func_set_ext_params xc_f03_hyb_cam_coef xc_f03_hyb_cam_coef none~add_functional->xc_f03_hyb_cam_coef xc_f03_hyb_exx_coef xc_f03_hyb_exx_coef none~add_functional->xc_f03_hyb_exx_coef none~add_slices->none~get~5 none~extend dft_grid_t%extend none~add_slices->none~extend none~setslice dft_grid_t%setSlice none~add_slices->none~setslice none~par_env_t_allreduce_byte par_env_t%par_env_t_allreduce_byte none~allreduce->none~par_env_t_allreduce_byte none~par_env_t_allreduce_c_bool par_env_t%par_env_t_allreduce_c_bool none~allreduce->none~par_env_t_allreduce_c_bool none~par_env_t_allreduce_dp_1d par_env_t%par_env_t_allreduce_dp_1d none~allreduce->none~par_env_t_allreduce_dp_1d none~par_env_t_allreduce_dp_2d par_env_t%par_env_t_allreduce_dp_2d none~allreduce->none~par_env_t_allreduce_dp_2d none~par_env_t_allreduce_dp_3d par_env_t%par_env_t_allreduce_dp_3d none~allreduce->none~par_env_t_allreduce_dp_3d none~par_env_t_allreduce_dp_4d par_env_t%par_env_t_allreduce_dp_4d none~allreduce->none~par_env_t_allreduce_dp_4d none~par_env_t_allreduce_dp_scalar par_env_t%par_env_t_allreduce_dp_scalar none~allreduce->none~par_env_t_allreduce_dp_scalar none~par_env_t_allreduce_int32_1d par_env_t%par_env_t_allreduce_int32_1d none~allreduce->none~par_env_t_allreduce_int32_1d none~par_env_t_allreduce_int32_scalar par_env_t%par_env_t_allreduce_int32_scalar none~allreduce->none~par_env_t_allreduce_int32_scalar none~par_env_t_allreduce_int64_1d par_env_t%par_env_t_allreduce_int64_1d none~allreduce->none~par_env_t_allreduce_int64_1d none~par_env_t_allreduce_int64_scalar par_env_t%par_env_t_allreduce_int64_scalar none~allreduce->none~par_env_t_allreduce_int64_scalar none~par_env_t_bcast_byte par_env_t%par_env_t_bcast_byte none~bcast->none~par_env_t_bcast_byte none~par_env_t_bcast_c_bool par_env_t%par_env_t_bcast_c_bool none~bcast->none~par_env_t_bcast_c_bool none~par_env_t_bcast_dp_1d par_env_t%par_env_t_bcast_dp_1d none~bcast->none~par_env_t_bcast_dp_1d none~par_env_t_bcast_dp_2d par_env_t%par_env_t_bcast_dp_2d none~bcast->none~par_env_t_bcast_dp_2d none~par_env_t_bcast_dp_3d par_env_t%par_env_t_bcast_dp_3d none~bcast->none~par_env_t_bcast_dp_3d none~par_env_t_bcast_dp_4d par_env_t%par_env_t_bcast_dp_4d none~bcast->none~par_env_t_bcast_dp_4d none~par_env_t_bcast_dp_scalar par_env_t%par_env_t_bcast_dp_scalar none~bcast->none~par_env_t_bcast_dp_scalar none~par_env_t_bcast_int32_1d par_env_t%par_env_t_bcast_int32_1d none~bcast->none~par_env_t_bcast_int32_1d none~par_env_t_bcast_int32_scalar par_env_t%par_env_t_bcast_int32_scalar none~bcast->none~par_env_t_bcast_int32_scalar none~par_env_t_bcast_int64_1d par_env_t%par_env_t_bcast_int64_1d none~bcast->none~par_env_t_bcast_int64_1d none~par_env_t_bcast_int64_scalar par_env_t%par_env_t_bcast_int64_scalar none~bcast->none~par_env_t_bcast_int64_scalar none~aoval basis_set%aoval none~compaos->none~aoval none~compmos xc_engine_t%compMOs none~compxc->none~compmos none~comprhoall xc_engine_t%compRhoAll none~compxc->none~comprhoall proc~compute~2 xc_libxc_t%compute none~compxc->proc~compute~2 xc_f03_func_end xc_f03_func_end none~destroy~2->xc_f03_func_end none~getslicenonzero->none~getbyid none~getenergy xc_lib_t%getEnergy none~getstats->none~getenergy none~check grid_3d_t%check none~get~5->none~check proc~init xc_libxc_t%init none~init~18->proc~init none~init~22 list_grid_3d_t%init none~init~21->none~init~22 none~resetprunedpointers xc_engine_t%resetPrunedPointers none~pruneaos->none~resetprunedpointers none~resetorbpointers~2 xc_engine_t%resetOrbPointers none~resetpointers->none~resetorbpointers~2 none~resetxcpointers xc_engine_t%resetXCPointers none~resetpointers->none~resetxcpointers proc~setpts xc_libxc_t%setPts none~resetpointers->proc~setpts none~set~4->none~get~5 none~push list_grid_3d_t%push none~set~4->none~push interface~compute_first_derivs compute_first_derivs proc~add_ecpder->interface~compute_first_derivs interface~free_integrator free_integrator proc~add_ecpder->interface~free_integrator interface~free_result free_result proc~add_ecpder->interface~free_result interface~init_integrator init_integrator proc~add_ecpder->interface~init_integrator interface~init_integrator_instance init_integrator_instance proc~add_ecpder->interface~init_integrator_instance interface~set_ecp_basis set_ecp_basis proc~add_ecpder->interface~set_ecp_basis none~bf_to_shell basis_set%bf_to_shell proc~add_ecpder->none~bf_to_shell none~evaluate rys_root_t%evaluate proc~comp_coulomb_der1->none~evaluate proc~comp_coulomb_helfeyder1->none~evaluate proc~doquadgausshermite doQuadGaussHermite proc~comp_kinetic_der1->proc~doquadgausshermite proc~comp_overlap_der1->proc~doquadgausshermite proc~grd2_rys_compute->none~evaluate proc~ints_exchange->interface~show_message proc~ints_exchange->none~alloc proc~ints_exchange->none~compute proc~ints_exchange->none~set libint2_cleanup_eri libint2_cleanup_eri proc~ints_exchange->libint2_cleanup_eri libint2_init_eri libint2_init_eri proc~ints_exchange->libint2_init_eri none~clean~2 int2_rys_data_t%clean proc~ints_exchange->none~clean~2 none~init~2 int2_rys_data_t%init proc~ints_exchange->none~init~2 none~set_ids int2_rys_data_t%set_ids proc~ints_exchange->none~set_ids proc~genr22 genr22 proc~ints_exchange->proc~genr22 proc~int2_rys_compute int2_rys_compute proc~ints_exchange->proc~int2_rys_compute proc~libint_compute_eri libint_compute_eri proc~ints_exchange->proc~libint_compute_eri proc~lebedev_get_grid->interface~show_message proc~oqp_dtpttr_i64->interface~show_message dtpttr dtpttr proc~oqp_dtpttr_i64->dtpttr proc~orthogonal_transform_sym->interface~show_message proc~orthogonal_transform_sym->proc~oqp_dgemm_i64 proc~orthogonal_transform_sym->proc~oqp_dtpttr_i64 proc~oqp_dsymm_i64 oqp_dsymm_i64 proc~orthogonal_transform_sym->proc~oqp_dsymm_i64 proc~oqp_dtrttp_i64 oqp_dtrttp_i64 proc~orthogonal_transform_sym->proc~oqp_dtrttp_i64 none~compaov basis_set%compAOv none~aoval->none~compaov none~compaovg basis_set%compAOvg none~aoval->none~compaovg none~compaovgg basis_set%compAOvgg none~aoval->none~compaovgg proc~mo_tran_gemm_ mo_tran_gemm_ none~compmos->proc~mo_tran_gemm_ proc~mo_tran_symm_ mo_tran_symm_ none~compmos->proc~mo_tran_symm_ none~clear~2 list_grid_3d_t%clear none~init~22->none~clear~2 mpi_allreduce mpi_allreduce none~par_env_t_allreduce_byte->mpi_allreduce none~par_env_t_allreduce_c_bool->mpi_allreduce none~par_env_t_allreduce_dp_1d->mpi_allreduce none~par_env_t_allreduce_dp_2d->mpi_allreduce none~par_env_t_allreduce_dp_3d->mpi_allreduce none~par_env_t_allreduce_dp_4d->mpi_allreduce none~par_env_t_allreduce_dp_scalar->mpi_allreduce none~par_env_t_allreduce_int32_1d->mpi_allreduce none~par_env_t_allreduce_int32_scalar->mpi_allreduce none~par_env_t_allreduce_int64_1d->mpi_allreduce none~par_env_t_allreduce_int64_scalar->mpi_allreduce mpi_bcast mpi_bcast none~par_env_t_bcast_byte->mpi_bcast none~par_env_t_bcast_c_bool->mpi_bcast none~par_env_t_bcast_dp_1d->mpi_bcast none~par_env_t_bcast_dp_2d->mpi_bcast none~par_env_t_bcast_dp_3d->mpi_bcast none~par_env_t_bcast_dp_4d->mpi_bcast none~par_env_t_bcast_dp_scalar->mpi_bcast none~par_env_t_bcast_int32_1d->mpi_bcast none~par_env_t_bcast_int32_scalar->mpi_bcast none~par_env_t_bcast_int64_1d->mpi_bcast none~par_env_t_bcast_int64_scalar->mpi_bcast none~calc_evfxc functional_t%calc_evfxc proc~compute~2->none~calc_evfxc none~calc_evxc functional_t%calc_evxc proc~compute~2->none~calc_evxc none~calc_xc functional_t%calc_xc proc~compute~2->none~calc_xc none~scalexc xc_lib_t%scalexc proc~compute~2->none~scalexc abrt abrt proc~doquadgausshermite->abrt fgrid fgrid proc~genr22->fgrid rfinc rfinc proc~genr22->rfinc rmr rmr proc~genr22->rmr xgrid xgrid proc~genr22->xgrid none~clean~13 xc_lib_t%clean proc~init->none~clean~13 none~resetenergy xc_lib_t%resetEnergy proc~init->none~resetenergy proc~int2_rys_compute->none~evaluate ab_x ab_x proc~libint_compute_eri->ab_x ab_y ab_y proc~libint_compute_eri->ab_y ab_z ab_z proc~libint_compute_eri->ab_z alpha1over_zetapluseta alpha1over_zetapluseta proc~libint_compute_eri->alpha1over_zetapluseta alpha1rho_over_zeta2 alpha1rho_over_zeta2 proc~libint_compute_eri->alpha1rho_over_zeta2 alpha2over_zetapluseta alpha2over_zetapluseta proc~libint_compute_eri->alpha2over_zetapluseta alpha2rho_over_zeta2 alpha2rho_over_zeta2 proc~libint_compute_eri->alpha2rho_over_zeta2 alpha3over_zetapluseta alpha3over_zetapluseta proc~libint_compute_eri->alpha3over_zetapluseta alpha3rho_over_eta2 alpha3rho_over_eta2 proc~libint_compute_eri->alpha3rho_over_eta2 alpha4over_zetapluseta alpha4over_zetapluseta proc~libint_compute_eri->alpha4over_zetapluseta alpha4rho_over_eta2 alpha4rho_over_eta2 proc~libint_compute_eri->alpha4rho_over_eta2 cd_x cd_x proc~libint_compute_eri->cd_x cd_y cd_y proc~libint_compute_eri->cd_y cd_z cd_z proc~libint_compute_eri->cd_z f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_0 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_0 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_0 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_1 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_1 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_1 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_10 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_10 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_10 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_11 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_11 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_11 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_12 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_12 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_12 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_13 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_13 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_13 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_14 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_14 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_14 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_15 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_15 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_15 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_16 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_16 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_16 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_17 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_17 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_17 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_18 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_18 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_18 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_19 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_19 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_19 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_2 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_2 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_2 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_20 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_20 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_20 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_3 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_3 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_3 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_4 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_4 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_4 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_5 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_5 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_5 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_6 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_6 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_6 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_7 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_7 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_7 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_8 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_8 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_8 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_9 f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_9 proc~libint_compute_eri->f_ab_s___0__s___1___twoprep_s___0__s___1___ab__up_9 proc~libint_compute_eri->fgrid igrid igrid proc~libint_compute_eri->igrid irgrd irgrd proc~libint_compute_eri->irgrd oo2e oo2e proc~libint_compute_eri->oo2e oo2z oo2z proc~libint_compute_eri->oo2z oo2ze oo2ze proc~libint_compute_eri->oo2ze pa_x pa_x proc~libint_compute_eri->pa_x pa_y pa_y proc~libint_compute_eri->pa_y pa_z pa_z proc~libint_compute_eri->pa_z qc_x qc_x proc~libint_compute_eri->qc_x qc_y qc_y proc~libint_compute_eri->qc_y qc_z qc_z proc~libint_compute_eri->qc_z proc~libint_compute_eri->rfinc rho12_over_alpha1 rho12_over_alpha1 proc~libint_compute_eri->rho12_over_alpha1 rho12_over_alpha2 rho12_over_alpha2 proc~libint_compute_eri->rho12_over_alpha2 rho34_over_alpha3 rho34_over_alpha3 proc~libint_compute_eri->rho34_over_alpha3 rho34_over_alpha4 rho34_over_alpha4 proc~libint_compute_eri->rho34_over_alpha4 proc~libint_compute_eri->rmr roe roe proc~libint_compute_eri->roe roz roz proc~libint_compute_eri->roz two_alpha0_bra two_alpha0_bra proc~libint_compute_eri->two_alpha0_bra two_alpha0_ket two_alpha0_ket proc~libint_compute_eri->two_alpha0_ket two_alpha1bra two_alpha1bra proc~libint_compute_eri->two_alpha1bra two_alpha1ket two_alpha1ket proc~libint_compute_eri->two_alpha1ket wp_x wp_x proc~libint_compute_eri->wp_x wp_y wp_y proc~libint_compute_eri->wp_y wp_z wp_z proc~libint_compute_eri->wp_z wq_x wq_x proc~libint_compute_eri->wq_x wq_y wq_y proc~libint_compute_eri->wq_y wq_z wq_z proc~libint_compute_eri->wq_z proc~oqp_dsymm_i64->interface~show_message dsymm dsymm proc~oqp_dsymm_i64->dsymm proc~oqp_dtrttp_i64->interface~show_message dtrttp dtrttp proc~oqp_dtrttp_i64->dtrttp none~calc_evfxc->interface~show_message none~calc_evfxc->xc_f03_func_info_get_family none~calc_evfxc->xc_f03_func_info_get_kind none~calc_evfxc->xc_f03_func_info_get_name xc_f03_gga_exc_vxc_fxc xc_f03_gga_exc_vxc_fxc none~calc_evfxc->xc_f03_gga_exc_vxc_fxc xc_f03_lda_exc_vxc_fxc xc_f03_lda_exc_vxc_fxc none~calc_evfxc->xc_f03_lda_exc_vxc_fxc xc_f03_mgga_exc_vxc_fxc xc_f03_mgga_exc_vxc_fxc none~calc_evfxc->xc_f03_mgga_exc_vxc_fxc none~calc_evxc->interface~show_message none~calc_evxc->xc_f03_func_info_get_family none~calc_evxc->xc_f03_func_info_get_kind none~calc_evxc->xc_f03_func_info_get_name xc_f03_gga_exc_vxc xc_f03_gga_exc_vxc none~calc_evxc->xc_f03_gga_exc_vxc xc_f03_lda_exc_vxc xc_f03_lda_exc_vxc none~calc_evxc->xc_f03_lda_exc_vxc xc_f03_mgga_exc_vxc xc_f03_mgga_exc_vxc none~calc_evxc->xc_f03_mgga_exc_vxc none~calc_xc->interface~show_message none~calc_xc->xc_f03_func_info_get_family none~calc_xc->xc_f03_func_info_get_kind none~calc_xc->xc_f03_func_info_get_name xc_f03_gga_exc_vxc_fxc_kxc xc_f03_gga_exc_vxc_fxc_kxc none~calc_xc->xc_f03_gga_exc_vxc_fxc_kxc xc_f03_lda_exc_vxc_fxc_kxc xc_f03_lda_exc_vxc_fxc_kxc none~calc_xc->xc_f03_lda_exc_vxc_fxc_kxc xc_f03_mgga_exc_vxc_fxc_kxc xc_f03_mgga_exc_vxc_fxc_kxc none~calc_xc->xc_f03_mgga_exc_vxc_fxc_kxc proc~mo_tran_gemm_->proc~oqp_dgemm_i64 proc~mo_tran_symm_->proc~oqp_dsymm_i64

Called by

proc~~tdhf_gradient~~CalledByGraph proc~tdhf_gradient tdhf_gradient proc~tdhf_gradient_c tdhf_gradient_C proc~tdhf_gradient_c->proc~tdhf_gradient

Source Code

  subroutine tdhf_gradient(infos)
    use io_constants, only: iw
    use oqp_tagarray_driver

    use types, only: information
    use strings, only: Cstring, fstring
    use basis_tools, only: basis_set
    use messages, only: show_message, with_abort

    use grd1, only: eijden, print_gradient
    use util, only: measure_time
    use tdhf_lib, only: &
      iatogen, mntoia
    use dft, only: dft_initialize, dftclean
    use mathlib, only: symmetrize_matrix, orthogonal_transform
    use mod_dft_molgrid, only: dft_grid_t
    use mod_dft_gridint_tdxc_grad, only: tddft_xc_gradient
    use mathlib, only: unpack_matrix
    use printing, only: print_module_info

    implicit none

    character(len=*), parameter :: subroutine_name = "tdhf_gradient"

    type(basis_set), pointer :: basis
    type(information), target, intent(inout) :: infos

    integer :: nbf, nocc

    type(dft_grid_t) :: molGrid

  ! General data
    logical :: dft
    integer :: scf_type, mol_mult

    real(kind=dp), allocatable :: p(:,:,:), d(:,:,:), xpy2(:,:,:), xmy2(:,:,:), wrk1(:,:), wrk2(:,:)

    ! tagarray
    real(kind=dp), contiguous, pointer :: &
      dmat_a(:), mo_a(:,:), td_p(:,:), xpy(:), xmy(:)
    character(len=*), parameter :: tags_general(*) = [ character(len=80) :: &
      OQP_TD_XPY, OQP_TD_XMY, OQP_DM_A, OQP_VEC_MO_A, OQP_TD_P ]

    mol_mult = infos%mol_prop%mult
    if (mol_mult/=1) call show_message(&
            'RPA-TDDFT are only available for RHF reference', with_abort)

    scf_type = infos%control%scftype
    if (scf_type/=1) error stop

    dft = infos%control%hamilton == 20

  ! Files open
    open (unit=IW, file=infos%log_filename, position="append")
  !
    call print_module_info('TDHF_Grad','Computing Grdient of TDDFT')
!
    write(iw,'(/5X,"Gradient options"/&
                &5X,18("-")/&
                &5X,"Target State: ",I8/)') infos%tddft%target_state
   ! Load basis set
    basis => infos%basis
    basis%atoms => infos%atoms

    call data_has_tags(infos%dat, tags_general, module_name, subroutine_name, WITH_ABORT)
    call tagarray_get_data(infos%dat, OQP_DM_A, dmat_a)
    call tagarray_get_data(infos%dat, OQP_VEC_MO_A, mo_a)
    call tagarray_get_data(infos%dat, OQP_td_p, td_p)
    call tagarray_get_data(infos%dat, OQP_td_xpy, xpy)
    call tagarray_get_data(infos%dat, OQP_td_xmy, xmy)

  ! Allocate H, S ,T and D matrices
    nbf = basis%nbf
    nocc = infos%mol_prop%nocc

!   Compute 1e gradient
    call flush(iw)

    call tdhf_1e_grad(infos, basis)
    !call print_gradient(infos)

    write(iw,"(' ..... End Of 1-Eelectron Gradient ......')")
    call measure_time(print_total=1, log_unit=iw)
    call flush(iw)

    allocate(d(nbf,nbf,2), &
             p(nbf,nbf,2), &
             xpy2(nbf,nbf,2), &
             xmy2(nbf,nbf,2), &
             wrk1(nbf,nbf), &
             wrk2(nbf,nbf), &
             source=0.0d0)

    call iatogen(xpy, wrk1, nocc, nocc)
    call symmetrize_matrix(wrk1, nbf)
    wrk1 = wrk1*0.5
    call orthogonal_transform('t', nbf, mo_a, wrk1, xpy2(:,:,1), wrk2)

    call iatogen(xmy, wrk1, nocc, nocc)
    call orthogonal_transform('t', nbf, mo_a, wrk1, xmy2(:,:,1), wrk2)


    call unpack_matrix(td_p(:,1), p(:,:,1))
    call unpack_matrix(dmat_a, d(:,:,1))

!   Compute xc gradient
    if (dft) then
      call dft_initialize(infos, basis, molGrid)

      d(:,:,2) = d(:,:,1)
      p(:,:,2) = p(:,:,1)
      xpy2(:,:,2) = xpy2(:,:,1)
      call tddft_xc_gradient(basis=basis, &
             molGrid=molGrid, &
             dedft=infos%atoms%grad, &
             da=d(:,:,1), &
             pa=p(:,:,1:1), &
             xa=xpy2(:,:,1:1), &
             nmtx=1, &
             threshold=1.0d-14, &
             infos=infos)
      call dftclean(infos)
      call measure_time(print_total=1, log_unit=iw)
      call flush(iw)
    end if

!   Compute 2e gradient
    call tdhf_2e_grad(basis, infos, d(:,:,1:1), p(:,:,1:1), xpy2(:,:,1:1), xmy2(:,:,1:1))

    call print_gradient(infos)

!   Print timings
    call measure_time(print_total=1, log_unit=iw)

    close(iw)

  end subroutine tdhf_gradient