Moved remap_set_mask() to remaplib.c

src/remap.h

@@ -216,40 +216,8 @@ LonLatPoint remapgrid_get_lonlat(const RemapGrid *grid, size_t index);
 void remap_check_area(size_t grid_size, const Varray<double> &cell_area, const char *name);
 
 template <typename T>
-void
-remap_set_mask(const Varray<T> &v, size_t gridsize, size_t numMissVals, double mv, Varray<int8_t> &mask)
-{
-  T missval = mv;
-  mask.resize(gridsize);
-
-  if (numMissVals)
-    {
-      if (std::isnan(missval))
-        {
-#ifdef _OPENMP
-#pragma omp parallel for default(shared) schedule(static)
-#endif
-          for (size_t i = 0; i < gridsize; ++i) mask[i] = fp_is_not_equal(v[i], missval);
-        }
-      else
-        {
-#ifdef _OPENMP
-#pragma omp parallel for default(shared) schedule(static)
-#endif
-          for (size_t i = 0; i < gridsize; ++i) mask[i] = !is_equal(v[i], missval);
-        }
-    }
-  else
-    {
-      for (size_t i = 0; i < gridsize; ++i) mask[i] = 1;
-    }
-}
-
-inline void
-remap_set_mask(const Field &field1, size_t gridsize, size_t numMissVals, double missval, Varray<int8_t> &imask)
-{
-  auto func = [&](auto &v) { remap_set_mask(v, gridsize, numMissVals, missval, imask); };
-  field_operation(func, field1);
-}
+void remap_set_mask(const Varray<T> &v, size_t gridsize, size_t numMissVals, double mv, Varray<int8_t> &mask);
+
+void remap_set_mask(const Field &field1, size_t gridsize, size_t numMissVals, double missval, Varray<int8_t> &imask);
 
 #endif /* REMAP_H */

src/remaplib.cc

@@ -1145,3 +1145,44 @@ remap_check_area(size_t gridSize, const Varray<double> &cell_area, const char *n
       if (cell_area[i] < -0.01) cdo_print("%s grid area error: %zu %g", name, i, cell_area[i]);
     }
 }
+
+template <typename T>
+void
+remap_set_mask(const Varray<T> &v, size_t gridsize, size_t numMissVals, double mv, Varray<int8_t> &mask)
+{
+  T missval = mv;
+  mask.resize(gridsize);
+
+  if (numMissVals)
+    {
+      if (std::isnan(missval))
+        {
+#ifdef _OPENMP
+#pragma omp parallel for default(shared) schedule(static)
+#endif
+          for (size_t i = 0; i < gridsize; ++i) mask[i] = fp_is_not_equal(v[i], missval);
+        }
+      else
+        {
+#ifdef _OPENMP
+#pragma omp parallel for default(shared) schedule(static)
+#endif
+          for (size_t i = 0; i < gridsize; ++i) mask[i] = !is_equal(v[i], missval);
+        }
+    }
+  else
+    {
+      for (size_t i = 0; i < gridsize; ++i) mask[i] = 1;
+    }
+}
+
+// Explicit instantiation
+template void remap_set_mask(const Varray<float> &v, size_t gridsize, size_t numMissVals, double mv, Varray<int8_t> &mask);
+template void remap_set_mask(const Varray<double> &v, size_t gridsize, size_t numMissVals, double mv, Varray<int8_t> &mask);
+
+void
+remap_set_mask(const Field &field1, size_t gridsize, size_t numMissVals, double missval, Varray<int8_t> &imask)
+{
+  auto func = [&](auto &v) { remap_set_mask(v, gridsize, numMissVals, missval, imask); };
+  field_operation(func, field1);
+}