rtp: add dipole calculation to linear scaling part

svn-origin-rev: 14529

src/rt_propagation_output.F

@@ -183,6 +183,7 @@ CONTAINS
     IF(rtp%converged)THEN
        IF(dft_control%qs_control%dftb)THEN
           IF (rtp%linear_scaling) THEN
+             CALL make_moment(qs_env,error)
              CALL get_rtp(rtp=rtp,rho_new=rho_new,error=error)
              IF (rtp%write_restart) CALL write_rt_p_to_restart(rho_new,particle_set,dft_section,&
                               atomic_kind_set,qs_kind_set,error)
@@ -196,6 +197,7 @@ CONTAINS
        ELSE
           dft_section =>  section_vals_get_subs_vals(input,"DFT",error=error)
           IF (rtp%linear_scaling) THEN
+             CALL make_moment(qs_env,error)
              CALL write_mo_free_results(qs_env,error)
              CALL get_rtp(rtp=rtp,rho_new=rho_new,error=error)
              IF (rtp%write_restart) CALL write_rt_p_to_restart(rho_new,particle_set,dft_section,&

tests/QS/regtest-ls-rtp/H2O-ls-emd.inp

@@ -43,6 +43,13 @@
       &XC_FUNCTIONAL PADE
       &END XC_FUNCTIONAL
     &END XC
+    &PRINT
+      &MOMENTS MEDIUM
+         ADD_LAST NUMERIC
+         FILENAME output
+         PERIODIC .FALSE.
+      &END MOMENTS
+    &END PRINT
   &END DFT
   &SUBSYS
     &CELL