remap_bicubic: allocate and init srcGridMask (bug fix)

src/remap_bicubic.cc

@@ -61,12 +61,10 @@ renormalize_weights(const double (&srcLats)[4], double (&weights)[4][4])
 }
 
 static void
-bicubic_sort_weights(size_t (&indices)[4], double (&weights)[4][4])
+bicubic_sort_weights_by_index(size_t (&indices)[4], double (&weights)[4][4])
 {
   constexpr size_t numWeights = 4;
 
-  if (is_sorted_list(numWeights, indices)) return;
-
   struct IndexWeightX
   {
     size_t index;
@@ -90,6 +88,15 @@ bicubic_sort_weights(size_t (&indices)[4], double (&weights)[4][4])
     }
 }
 
+static void
+bicubic_sort_weights(size_t (&indices)[4], double (&weights)[4][4])
+{
+  constexpr size_t numWeights = 4;
+  if (is_sorted_list(numWeights, indices)) return;
+
+  bicubic_sort_weights_by_index(indices, weights);
+}
+
 static void
 bicubic_warning()
 {
@@ -238,7 +245,7 @@ remap_bicubic(RemapSearch &rsearch, const Varray<T> &srcArray, Varray<T> &tgtArr
   auto tgtGridSize = tgtGrid->size;
   auto srcGridSize = srcGrid->size;
 
-  Varray<short> srcGridMask;
+  Varray<short> srcGridMask(srcGridSize, 1);
   if (numMissVals) remap_set_mask(srcGridSize, srcArray, numMissVals, missval, srcGridMask);
 
   // Compute mappings from source to target grid

src/remap_bilinear.cc

@@ -135,7 +135,7 @@ bilinear_sort_weights_by_index_zip(size_t (&indices)[4], double (&weights)[4])
 }
 #endif
 
-void
+static void
 bilinear_sort_weights_by_index(size_t (&indices)[4], double (&weights)[4])
 {
   constexpr size_t numWeights = 4;