psc_bnd_fields_impl: fix vector t

Author: JamesMcClung

src/libpsc/psc_bnd_fields/psc_bnd_fields_impl.hxx

@@ -557,10 +557,10 @@ struct BndFields_ : BndFieldsBase
       real_t s = 0.0;
       real_t p = 0.0;
       if (radiation) {
-        Real3 x3_s =
-          (Real3(edge_idx) + Real3::unit(d1) * real_t(0.5)) * grid.domain.dx;
-        Real3 x3_p =
-          (Real3(edge_idx) + Real3::unit(d2) * real_t(0.5)) * grid.domain.dx;
+        Real3 x3_s = (Real3(edge_idx) + Real3::unit(d1) * real_t(0.5)) *
+                     Real3(grid.domain.dx);
+        Real3 x3_p = (Real3(edge_idx) + Real3::unit(d2) * real_t(0.5)) *
+                     Real3(grid.domain.dx);
         s = radiation->pulse_s_lower(grid.time(), d0, p, x3_s);
         p = radiation->pulse_p_lower(grid.time(), d0, p, x3_p);
       }
@@ -608,10 +608,10 @@ struct BndFields_ : BndFieldsBase
       real_t s = 0.0;
       real_t p = 0.0;
       if (radiation) {
-        Real3 x3_s =
-          (Real3(edge_idx) + Real3::unit(d1) * real_t(0.5)) * grid.domain.dx;
-        Real3 x3_p =
-          (Real3(edge_idx) + Real3::unit(d2) * real_t(0.5)) * grid.domain.dx;
+        Real3 x3_s = (Real3(edge_idx) + Real3::unit(d1) * real_t(0.5)) *
+                     Real3(grid.domain.dx);
+        Real3 x3_p = (Real3(edge_idx) + Real3::unit(d2) * real_t(0.5)) *
+                     Real3(grid.domain.dx);
         s = radiation->pulse_s_upper(grid.time(), d0, p, x3_s);
         p = radiation->pulse_p_upper(grid.time(), d0, p, x3_p);
       }