Static cast in dot product.

JCSDA-internal · Jan 6, 2024 · c52e50f · c52e50f
1 parent f865b49
commit c52e50f
Showing 1 changed file with 4 additions and 6 deletions.
diff --git a/src/tests/interpolation/test_interpolation_spherical_vector.cc b/src/tests/interpolation/test_interpolation_spherical_vector.cc
@@ -141,9 +141,10 @@ template <int Rank>
 double dotProduct(const array::ArrayView<double, Rank>& a,
                   const array::ArrayView<double, Rank>& b) {
   auto dotProd = 0.;
-  arrayForEachDim(
-      std::make_integer_sequence<int, Rank>{}, std::tie(a, b),
-      [&](auto&& aElem, auto&& bElem) { dotProd += aElem * bElem; });
+  arrayForEachDim(std::make_integer_sequence<int, Rank>{}, std::tie(a, b),
+                  [&](auto&& aElem, auto&& bElem) {
+    dotProd += static_cast<double>(aElem) * static_cast<double>(bElem);
+  });
   return dotProd;
 }
 
@@ -175,9 +176,6 @@ void testInterpolation(const Config& config) {
   auto sourceView = array::make_view<double, Rank>(sourceField);
   auto targetView = array::make_view<double, Rank>(targetField);
 
-  // MacOS test failing if I don't pre-assign field.
-  targetView.assign(0.);
-
   ArrayForEach<0>::apply(std::tie(sourceLonLat, sourceView),
                          [](auto&& lonLat, auto&& sourceColumn) {