grid_search_point_unstruct: added m_srcCoords for KnnMethod::gaussWeighted

src/pointsearch_reg2d.cc

@@ -81,8 +81,9 @@ PointsearchReg2d::store_distance_reg2d(double plon, double plat, KnnData &knnDat
 
   if (knnData.m_method == KnnMethod::gaussWeighted)
     {
-      if (numWeights > knnData.m_srcCoords.size()) knnData.m_srcCoords.resize(numWeights);
-      for (size_t i = 0; i < numWeights; ++i) { compute_point(indices[i], knnData.m_srcCoords[i]); }
+      auto numNeighbors = knnData.m_numNeighbors;
+      if (numNeighbors > knnData.m_srcCoords.size()) knnData.m_srcCoords.resize(numNeighbors);
+      for (size_t i = 0; i < numNeighbors; ++i) { compute_point(knnData.m_indices[i], knnData.m_srcCoords[i]); }
     }
 }
 

src/remap_point_search.cc

@@ -130,6 +130,16 @@ grid_search_point_unstruct(GridPointsearch &gps, const PointLonLat &pointLL, Knn
   for (size_t i = 0; i < numDist; ++i) knnData.store_distance(indices[i], dist[i], numNeighbors);
 
   knnData.check_distance();
+
+  if (knnData.m_method == KnnMethod::gaussWeighted)
+    {
+      numNeighbors = knnData.m_numNeighbors;
+      if (numNeighbors > knnData.m_srcCoords.size()) knnData.m_srcCoords.resize(numNeighbors);
+      for (size_t i = 0; i < numNeighbors; ++i)
+        {
+          gcLLtoXYZ(gps.plons[knnData.m_indices[i]], gps.plats[knnData.m_indices[i]], knnData.m_srcCoords[i]);
+        }
+    }
 }
 
 void