Split vector products routines into new module

hiromatsui · hiromatsui · commit ca3be929dce8 · 2025-09-24T16:07:29.000+09:00
diff --git a/src/Fortran_libraries/MHD_src/sph_MHD/cal_ene_flux_by_sym_rtp.f90 b/src/Fortran_libraries/MHD_src/sph_MHD/cal_ene_flux_by_sym_rtp.f90
@@ -165,7 +165,7 @@ subroutine cal_work_of_lorentz_on_node                            &
      &          nnod, ntot_comp_fld, fld_rtp, ntot_comp_frc, frc_rtp,   &
      &          ntot_comp_flx, flx_rtp)
 !
-      use cal_products_smp
+      use cal_vector_products
 !
       type(base_field_address), intent(in) :: bs_trns_base
       type(base_force_address), intent(in) :: f_trns_frc
@@ -180,7 +180,7 @@ subroutine cal_work_of_lorentz_on_node                            &
       real(kind = kreal), intent(inout) :: flx_rtp(nnod,ntot_comp_flx)
 !
       if(fs_trns_eflux%i_ujb .gt. 0) then
-        call cal_dot_prod_no_coef_smp(nnod,                             &
+        call cal_dot_product_no_coef(nnod,                              &
      &      frc_rtp(1,f_trns_frc%i_lorentz),                            &
      &      fld_rtp(1,bs_trns_base%i_velo),                             &
      &      flx_rtp(1,fs_trns_eflux%i_ujb) )
@@ -195,7 +195,7 @@ subroutine cal_work_of_inertia_on_node                            &
      &          nnod, ntot_comp_fld, fld_rtp, ntot_comp_frc, frc_rtp,   &
      &          ntot_comp_flx, flx_rtp)
 !
-      use cal_products_smp
+      use cal_vector_products
 !
       type(base_field_address), intent(in) :: bs_trns_base
       type(base_force_address), intent(in) :: f_trns_frc
@@ -210,7 +210,7 @@ subroutine cal_work_of_inertia_on_node                            &
       real(kind = kreal), intent(inout) :: flx_rtp(nnod,ntot_comp_flx)
 !
       if(fs_trns_eflux%i_m_advect_work .gt. 0) then
-        call cal_dot_prod_no_coef_smp(nnod,                             &
+        call cal_dot_product_no_coef(nnod,                              &
      &      frc_rtp(1,f_trns_frc%i_m_advect),                           &
      &      fld_rtp(1,bs_trns_base%i_velo),                             &
      &      flx_rtp(1,fs_trns_eflux%i_m_advect_work) )
@@ -226,7 +226,7 @@ subroutine cal_ene_flux_by_induct_on_node                         &
      &          nnod, ntot_comp_fld, fld_rtp, ntot_comp_frc, frc_rtp,   &
      &          ntot_comp_flx, flx_rtp)
 !
-      use cal_products_smp
+      use cal_vector_products
 !
       type(base_field_address), intent(in) :: bs_trns_base
       type(base_force_address), intent(in) :: ipol_frc
@@ -241,7 +241,7 @@ subroutine cal_ene_flux_by_induct_on_node                         &
       real(kind = kreal), intent(inout) :: flx_rtp(nnod,ntot_comp_flx)
 !
       if(fs_trns_eflux%i_me_gen .gt. 0) then
-        call cal_dot_prod_no_coef_smp(nnod,                             &
+        call cal_dot_product_no_coef(nnod,                              &
      &      frc_rtp(1,ipol_frc%i_induction),                            &
      &      fld_rtp(1,bs_trns_base%i_magne),                            &
      &      flx_rtp(1,fs_trns_eflux%i_me_gen) )
diff --git a/src/Fortran_libraries/MHD_src/sph_MHD/cal_energy_flux_rtp.f90 b/src/Fortran_libraries/MHD_src/sph_MHD/cal_energy_flux_rtp.f90
@@ -128,7 +128,7 @@ subroutine cal_energy_fluxes_on_node                              &
      &          nnod, ntot_comp_fld, fld_rtp, ntot_comp_frc, frc_rtp,   &
      &          ntot_comp_uxb, fub_rtp, ntot_comp_flx, flx_rtp)
 !
-      use cal_products_smp
+      use cal_vector_products
 !
       type(base_field_address), intent(in) :: bs_trns_base
       type(base_force_address), intent(in) :: f_trns_frc
@@ -145,21 +145,21 @@ subroutine cal_energy_fluxes_on_node                              &
       real(kind = kreal), intent(inout) :: flx_rtp(nnod,ntot_comp_flx)
 !
       if(fs_trns_eflux%i_ujb .gt. 0) then
-        call cal_dot_prod_no_coef_smp(nnod,                             &
+        call cal_dot_product_no_coef(nnod,                              &
      &      frc_rtp(1,f_trns_frc%i_lorentz),                            &
      &      fld_rtp(1,bs_trns_base%i_velo),                             &
      &      flx_rtp(1,fs_trns_eflux%i_ujb) )
       end if
 !
       if(fs_trns_eflux%i_nega_ujb .gt. 0) then
-        call cal_dot_prod_w_coef_smp(nnod, dminus,                      &
+        call cal_dot_product_w_coef(nnod, dminus,                       &
      &      frc_rtp(1,f_trns_frc%i_lorentz),                            &
      &      fld_rtp(1,bs_trns_base%i_velo),                             &
      &      flx_rtp(1,fs_trns_eflux%i_nega_ujb))
       end if
 !
       if(fs_trns_eflux%i_me_gen .gt. 0) then
-        call cal_dot_prod_no_coef_smp(nnod,                             &
+        call cal_dot_product_no_coef(nnod,                              &
      &      fub_rtp(1,be_trns_frc%i_induction),                         &
      &      fld_rtp(1,bs_trns_base%i_magne),                            &
      &      flx_rtp(1,fs_trns_eflux%i_me_gen))
diff --git a/src/Fortran_libraries/MHD_src/sph_MHD/cal_helicities_rtp.f90 b/src/Fortran_libraries/MHD_src/sph_MHD/cal_helicities_rtp.f90
@@ -103,7 +103,7 @@ subroutine cal_helicity_on_node                                   &
      &         (bs_trns_base, fe_trns_prod, nnod,                       &
      &          ntot_comp_fld, fld_rtp, ntot_comp_hls, fhls_rtp)
 !
-      use cal_products_smp
+      use cal_vector_products
 !
       type(base_field_address), intent(in) :: bs_trns_base
       type(phys_products_address), intent(in) :: fe_trns_prod
@@ -116,19 +116,19 @@ subroutine cal_helicity_on_node                                   &
 !
 !
       if(fe_trns_prod%i_k_heli .gt. 0) then
-        call cal_dot_prod_no_coef_smp(nnod,                             &
+        call cal_dot_product_no_coef(nnod,                              &
      &      fld_rtp(1,bs_trns_base%i_velo),                             &
      &      fld_rtp(1,bs_trns_base%i_vort),                             &
      &      fhls_rtp(1,fe_trns_prod%i_k_heli))
       end if
       if(fe_trns_prod%i_c_heli .gt. 0) then
-        call cal_dot_prod_no_coef_smp(nnod,                             &
+        call cal_dot_product_no_coef(nnod,                              &
      &      fld_rtp(1,bs_trns_base%i_magne),                            &
      &      fld_rtp(1,bs_trns_base%i_current),                          &
      &      fhls_rtp(1,fe_trns_prod%i_c_heli))
       end if
       if(fe_trns_prod%i_x_heli .gt. 0) then
-        call cal_dot_prod_no_coef_smp(nnod,                             &
+        call cal_dot_product_no_coef(nnod,                              &
      &      fld_rtp(1,bs_trns_base%i_velo),                             &
      &      fld_rtp(1,bs_trns_base%i_magne),                            &
      &      fhls_rtp(1,fe_trns_prod%i_x_heli))
@@ -185,13 +185,13 @@ subroutine cal_square_vector_on_node                              &
      &      fmag_rtp(1,fe_trns_prod%i_square_j))
       end if
       if(fe_trns_prod%i_square_t .gt. 0) then
-        call cal_scalar_prod_no_coef_smp(nnod,                          &
+        call cal_scalar_product_no_coef(nnod,                           &
      &      fld_rtp(1,bs_trns_base%i_temp),                             &
      &      fld_rtp(1,bs_trns_base%i_temp),                             &
      &      fmag_rtp(1,fe_trns_prod%i_square_t))
       end if
       if(fe_trns_prod%i_square_c .gt. 0) then
-        call cal_scalar_prod_no_coef_smp(nnod,                          &
+        call cal_scalar_product_no_coef(nnod,                           &
      &      fld_rtp(1,bs_trns_base%i_light),                            &
      &      fld_rtp(1,bs_trns_base%i_light),                            &
      &      fmag_rtp(1,fe_trns_prod%i_square_c))
diff --git a/src/Fortran_libraries/MHD_src/sph_MHD/cal_lorentz_by_dipole.f90 b/src/Fortran_libraries/MHD_src/sph_MHD/cal_lorentz_by_dipole.f90
@@ -67,7 +67,7 @@ subroutine lorentz_frc_and_work_by_dipole                         &
      &         (fl_prop, bs_trns_base, bs_trns_prod, fs_trns_prod,      &
      &          nnod, ntot_comp_fld, fld_rtp, ntot_comp_flx, frc_rtp)
 !
-      use cal_products_smp
+      use cal_vector_products
 !
       type(fluid_property), intent(in) :: fl_prop
       type(base_field_address), intent(in) :: bs_trns_base
@@ -83,17 +83,19 @@ subroutine lorentz_frc_and_work_by_dipole                         &
 !
       if(fs_trns_prod%i_dipole_Lorentz .gt. 0) then
 ! 
-      if(iflag_debug.gt.0) write(*,*) 'cal_dot_prod_no_coef_smp for dipole lorentz'
-        call cal_cross_prod_w_coef_smp(nnod, fl_prop%coef_lor,          &
+      if(iflag_debug.gt.0) write(*,*)                                   &
+     &      'cal_cross_product_w_coef for dipole lorentz'
+        call cal_cross_product_w_coef(nnod, fl_prop%coef_lor,           &
      &      fld_rtp(1,bs_trns_prod%i_dipole_J),                         &
      &      fld_rtp(1,bs_trns_base%i_magne),                            &
      &      frc_rtp(1,fs_trns_prod%i_dipole_Lorentz) )
       end if
 !
       if(fs_trns_prod%i_dipole_ujb .gt. 0) then
 ! 
-        if(iflag_debug.gt.0) write(*,*) 'cal_dot_prod_no_coef_smp for dipole ujb'
-        call cal_dot_prod_no_coef_smp(nnod,                             &
+        if(iflag_debug.gt.0) write(*,*)                                 &
+     &      'cal_dot_product_no_coef for dipole ujb'
+        call cal_dot_product_no_coef(nnod,                              &
      &      frc_rtp(1,fs_trns_prod%i_dipole_Lorentz),                   &
      &      fld_rtp(1,bs_trns_base%i_velo),                             &
      &      frc_rtp(1,fs_trns_prod%i_dipole_ujb) )
diff --git a/src/Fortran_libraries/MHD_src/sph_MHD/cal_nonlinear_sph_MHD.f90 b/src/Fortran_libraries/MHD_src/sph_MHD/cal_nonlinear_sph_MHD.f90
@@ -158,6 +158,7 @@ subroutine nonlinear_terms_on_node                                &
      &          nnod, ntot_comp_fld, fld_rtp, ntot_comp_frc, frc_rtp)
 !
       use cal_products_smp
+      use cal_vector_products
 !
       type(MHD_evolution_param), intent(in) :: MHD_prop
       type(base_field_address), intent(in) :: b_trns_base
@@ -169,22 +170,22 @@ subroutine nonlinear_terms_on_node                                &
       real(kind = kreal), intent(inout) :: frc_rtp(nnod,ntot_comp_frc)
 !
       if(f_trns_frc%i_m_advect .gt. 0) then
-        call cal_cross_prod_w_coef_smp                                  &
+        call cal_cross_product_w_coef                                   &
      &     (nnod, MHD_prop%fl_prop%coef_velo,                           &
      &      fld_rtp(1,b_trns_base%i_vort),                              &
      &      fld_rtp(1,b_trns_base%i_velo),                              &
      &      frc_rtp(1,f_trns_frc%i_m_advect) )
       end if
 !
       if(f_trns_frc%i_lorentz .gt. 0) then
-        call cal_cross_prod_w_coef_smp                                  &
+        call cal_cross_product_w_coef                                   &
      &     (nnod, MHD_prop%fl_prop%coef_lor,                            &
      &      fld_rtp(1,b_trns_base%i_current),                           &
      &      fld_rtp(1,b_trns_base%i_magne),                             &
      &      frc_rtp(1,f_trns_frc%i_lorentz) )
 !
         if(b_trns_base%i_back_B .gt. 0) then
-          call add_cross_prod_w_coef_smp                                &
+          call add_cross_product_w_coef                                 &
      &       (nnod, MHD_prop%fl_prop%coef_lor,                          &
      &        fld_rtp(1,b_trns_base%i_current),                         &
      &        fld_rtp(1,b_trns_base%i_back_B),                          &
@@ -194,14 +195,14 @@ subroutine nonlinear_terms_on_node                                &
 !
 !
       if(f_trns_frc%i_vp_induct .gt. 0) then
-        call cal_cross_prod_w_coef_smp                                  &
+        call cal_cross_product_w_coef                                   &
      &     (nnod, MHD_prop%cd_prop%coef_induct,                         &
      &      fld_rtp(1,b_trns_base%i_velo),                              &
      &      fld_rtp(1,b_trns_base%i_magne),                             &
      &      frc_rtp(1,f_trns_frc%i_vp_induct) )
 !
         if(b_trns_base%i_back_B .gt. 0) then
-          call add_cross_prod_w_coef_smp                                &
+          call add_cross_product_w_coef                                 &
      &       (nnod, MHD_prop%cd_prop%coef_induct,                       &
      &        fld_rtp(1,b_trns_base%i_velo),                            &
      &        fld_rtp(1,b_trns_base%i_back_B),                          &
@@ -211,15 +212,15 @@ subroutine nonlinear_terms_on_node                                &
 !
 !
       if(f_trns_frc%i_h_flux .gt. 0) then
-        call cal_vec_scalar_prod_w_coef_smp                             &
+        call cal_vec_scalar_product_w_coef                              &
      &     (nnod, MHD_prop%ht_prop%coef_advect,                         &
      &      fld_rtp(1,b_trns_base%i_velo),                              &
      &      fld_rtp(1,b_trns_base%i_temp),                              &
      &      frc_rtp(1,f_trns_frc%i_h_flux) )
       end if
 !
       if(f_trns_frc%i_c_flux .gt. 0) then
-        call cal_vec_scalar_prod_w_coef_smp                             &
+        call cal_vec_scalar_product_w_coef                              &
      &     (nnod, MHD_prop%cp_prop%coef_advect,                         &
      &      fld_rtp(1,b_trns_base%i_velo),                              &
      &      fld_rtp(1,b_trns_base%i_light),                             &
@@ -236,6 +237,7 @@ subroutine nonlinear_terms_on_node_w_sym                          &
      &          nnod, ntot_comp_fld, fld_rtp, ntot_comp_frc, frc_rtp)
 !
       use cal_products_smp
+      use cal_vector_products
 !
       type(MHD_evolution_param), intent(in) :: MHD_prop
       type(base_field_address), intent(in) :: b_trns_base_1
@@ -249,22 +251,22 @@ subroutine nonlinear_terms_on_node_w_sym                          &
 !
 !
       if(f_trns_frc%i_m_advect .gt. 0) then
-        call cal_cross_prod_w_coef_smp                                  &
+        call cal_cross_product_w_coef                                   &
      &     (nnod, MHD_prop%fl_prop%coef_velo,                           &
      &      fld_rtp(1,b_trns_base_1%i_vort),                            &
      &      fld_rtp(1,b_trns_base_2%i_velo),                            &
      &      frc_rtp(1,f_trns_frc%i_m_advect) )
       end if
 !
       if(f_trns_frc%i_lorentz .gt. 0) then
-        call cal_cross_prod_w_coef_smp                                  &
+        call cal_cross_product_w_coef                                   &
      &     (nnod, MHD_prop%fl_prop%coef_lor,                            &
      &      fld_rtp(1,b_trns_base_1%i_current),                         &
      &      fld_rtp(1,b_trns_base_2%i_magne),                           &
      &      frc_rtp(1,f_trns_frc%i_lorentz) )
 !
         if(b_trns_base_2%i_back_B .gt. 0) then
-          call add_cross_prod_w_coef_smp                                &
+          call add_cross_product_w_coef                                 &
      &       (nnod, MHD_prop%fl_prop%coef_lor,                          &
      &        fld_rtp(1,b_trns_base_1%i_current),                       &
      &        fld_rtp(1,b_trns_base_2%i_back_B),                        &
@@ -274,14 +276,14 @@ subroutine nonlinear_terms_on_node_w_sym                          &
 !
 !
       if(f_trns_frc%i_vp_induct .gt. 0) then
-        call cal_cross_prod_w_coef_smp                                  &
+        call cal_cross_product_w_coef                                   &
      &     (nnod, MHD_prop%cd_prop%coef_induct,                         &
      &      fld_rtp(1,b_trns_base_1%i_velo),                            &
      &      fld_rtp(1,b_trns_base_2%i_magne),                           &
      &      frc_rtp(1,f_trns_frc%i_vp_induct) )
 !
         if(b_trns_base_2%i_back_B .gt. 0) then
-          call add_cross_prod_w_coef_smp                                &
+          call add_cross_product_w_coef                                 &
      &       (nnod, MHD_prop%cd_prop%coef_induct,                       &
      &        fld_rtp(1,b_trns_base_1%i_velo),                          &
      &        fld_rtp(1,b_trns_base_2%i_back_B),                        &
@@ -291,15 +293,15 @@ subroutine nonlinear_terms_on_node_w_sym                          &
 !
 !
       if(f_trns_frc%i_h_flux .gt. 0) then
-        call cal_vec_scalar_prod_w_coef_smp                             &
+        call cal_vec_scalar_product_w_coef                              &
      &     (nnod, MHD_prop%ht_prop%coef_advect,                         &
      &      fld_rtp(1,b_trns_base_1%i_velo),                            &
      &      fld_rtp(1,b_trns_base_2%i_temp),                            &
      &      frc_rtp(1,f_trns_frc%i_h_flux) )
       end if
 !
       if(f_trns_frc%i_c_flux .gt. 0) then
-        call cal_vec_scalar_prod_w_coef_smp                             &
+        call cal_vec_scalar_product_w_coef                              &
      &     (nnod, MHD_prop%cp_prop%coef_advect,                         &
      &      fld_rtp(1,b_trns_base_1%i_velo),                            &
      &      fld_rtp(1,b_trns_base_2%i_light),                           &
diff --git a/src/Fortran_libraries/SERIAL_src/BASE/cal_vector_products.f90 b/src/Fortran_libraries/SERIAL_src/BASE/cal_vector_products.f90

Original file line number	Diff line number	Diff line change
`@@ -165,7 +165,7 @@ subroutine cal_work_of_lorentz_on_node &`
`165`	`165`	`& nnod, ntot_comp_fld, fld_rtp, ntot_comp_frc, frc_rtp, &`
`166`	`166`	`& ntot_comp_flx, flx_rtp)`
`167`	`167`	`!`
`168`		`- use cal_products_smp`
	`168`	`+ use cal_vector_products`
`169`	`169`	`!`
`170`	`170`	`type(base_field_address), intent(in) :: bs_trns_base`
`171`	`171`	`type(base_force_address), intent(in) :: f_trns_frc`
`@@ -180,7 +180,7 @@ subroutine cal_work_of_lorentz_on_node &`
`180`	`180`	`real(kind = kreal), intent(inout) :: flx_rtp(nnod,ntot_comp_flx)`
`181`	`181`	`!`
`182`	`182`	`if(fs_trns_eflux%i_ujb .gt. 0) then`
`183`		`- call cal_dot_prod_no_coef_smp(nnod, &`
	`183`	`+ call cal_dot_product_no_coef(nnod, &`
`184`	`184`	`& frc_rtp(1,f_trns_frc%i_lorentz), &`
`185`	`185`	`& fld_rtp(1,bs_trns_base%i_velo), &`
`186`	`186`	`& flx_rtp(1,fs_trns_eflux%i_ujb) )`
`@@ -195,7 +195,7 @@ subroutine cal_work_of_inertia_on_node &`
`195`	`195`	`& nnod, ntot_comp_fld, fld_rtp, ntot_comp_frc, frc_rtp, &`
`196`	`196`	`& ntot_comp_flx, flx_rtp)`
`197`	`197`	`!`
`198`		`- use cal_products_smp`
	`198`	`+ use cal_vector_products`
`199`	`199`	`!`
`200`	`200`	`type(base_field_address), intent(in) :: bs_trns_base`
`201`	`201`	`type(base_force_address), intent(in) :: f_trns_frc`
`@@ -210,7 +210,7 @@ subroutine cal_work_of_inertia_on_node &`
`210`	`210`	`real(kind = kreal), intent(inout) :: flx_rtp(nnod,ntot_comp_flx)`
`211`	`211`	`!`
`212`	`212`	`if(fs_trns_eflux%i_m_advect_work .gt. 0) then`
`213`		`- call cal_dot_prod_no_coef_smp(nnod, &`
	`213`	`+ call cal_dot_product_no_coef(nnod, &`
`214`	`214`	`& frc_rtp(1,f_trns_frc%i_m_advect), &`
`215`	`215`	`& fld_rtp(1,bs_trns_base%i_velo), &`
`216`	`216`	`& flx_rtp(1,fs_trns_eflux%i_m_advect_work) )`
`@@ -226,7 +226,7 @@ subroutine cal_ene_flux_by_induct_on_node &`
`226`	`226`	`& nnod, ntot_comp_fld, fld_rtp, ntot_comp_frc, frc_rtp, &`
`227`	`227`	`& ntot_comp_flx, flx_rtp)`
`228`	`228`	`!`
`229`		`- use cal_products_smp`
	`229`	`+ use cal_vector_products`
`230`	`230`	`!`
`231`	`231`	`type(base_field_address), intent(in) :: bs_trns_base`
`232`	`232`	`type(base_force_address), intent(in) :: ipol_frc`
`@@ -241,7 +241,7 @@ subroutine cal_ene_flux_by_induct_on_node &`
`241`	`241`	`real(kind = kreal), intent(inout) :: flx_rtp(nnod,ntot_comp_flx)`
`242`	`242`	`!`
`243`	`243`	`if(fs_trns_eflux%i_me_gen .gt. 0) then`
`244`		`- call cal_dot_prod_no_coef_smp(nnod, &`
	`244`	`+ call cal_dot_product_no_coef(nnod, &`
`245`	`245`	`& frc_rtp(1,ipol_frc%i_induction), &`
`246`	`246`	`& fld_rtp(1,bs_trns_base%i_magne), &`
`247`	`247`	`& flx_rtp(1,fs_trns_eflux%i_me_gen) )`
Original file line number	Diff line number	Diff line change
`@@ -67,7 +67,7 @@ subroutine lorentz_frc_and_work_by_dipole &`
`67`	`67`	`& (fl_prop, bs_trns_base, bs_trns_prod, fs_trns_prod, &`
`68`	`68`	`& nnod, ntot_comp_fld, fld_rtp, ntot_comp_flx, frc_rtp)`
`69`	`69`	`!`
`70`		`- use cal_products_smp`
	`70`	`+ use cal_vector_products`
`71`	`71`	`!`
`72`	`72`	`type(fluid_property), intent(in) :: fl_prop`
`73`	`73`	`type(base_field_address), intent(in) :: bs_trns_base`
`@@ -83,17 +83,19 @@ subroutine lorentz_frc_and_work_by_dipole &`
`83`	`83`	`!`
`84`	`84`	`if(fs_trns_prod%i_dipole_Lorentz .gt. 0) then`
`85`	`85`	`!`
`86`		`- if(iflag_debug.gt.0) write(,) 'cal_dot_prod_no_coef_smp for dipole lorentz'`
`87`		`- call cal_cross_prod_w_coef_smp(nnod, fl_prop%coef_lor, &`
	`86`	`+ if(iflag_debug.gt.0) write(,) &`
	`87`	`+ & 'cal_cross_product_w_coef for dipole lorentz'`
	`88`	`+ call cal_cross_product_w_coef(nnod, fl_prop%coef_lor, &`
`88`	`89`	`& fld_rtp(1,bs_trns_prod%i_dipole_J), &`
`89`	`90`	`& fld_rtp(1,bs_trns_base%i_magne), &`
`90`	`91`	`& frc_rtp(1,fs_trns_prod%i_dipole_Lorentz) )`
`91`	`92`	`end if`
`92`	`93`	`!`
`93`	`94`	`if(fs_trns_prod%i_dipole_ujb .gt. 0) then`
`94`	`95`	`!`
`95`		`- if(iflag_debug.gt.0) write(,) 'cal_dot_prod_no_coef_smp for dipole ujb'`
`96`		`- call cal_dot_prod_no_coef_smp(nnod, &`
	`96`	`+ if(iflag_debug.gt.0) write(,) &`
	`97`	`+ & 'cal_dot_product_no_coef for dipole ujb'`
	`98`	`+ call cal_dot_product_no_coef(nnod, &`
`97`	`99`	`& frc_rtp(1,fs_trns_prod%i_dipole_Lorentz), &`
`98`	`100`	`& fld_rtp(1,bs_trns_base%i_velo), &`
`99`	`101`	`& frc_rtp(1,fs_trns_prod%i_dipole_ujb) )`