GRIDEDIT-247 Minor update

Deltares · Jul 17, 2023 · da067b2 · da067b2
1 parent b5e8c82
commit da067b2
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/libs/MeshKernel/src/Splines.cpp b/libs/MeshKernel/src/Splines.cpp
@@ -711,11 +711,14 @@ void Splines::snapSpline(const LandBoundaries& landBoundary,
         {
             Point point(xf(i), yf(i));
             landBoundary.nearestPointOnLandBoundary(point, 2, nearestPoint, smallestDistance, segmentIndex, scaledDistance);
-            // toLand(xf[i], 1, mxlan, 2, nearestPoint, smallestDistance, segmentIndex, scaledDistance);
             xbVec(i) = nearestPoint.x;
             ybVec(i) = nearestPoint.y;
         }
 
+        // TODO how to get a pointwise multiply here?
+        // xbVec = weights.colwise() * xbVec.colwise();
+        // Can then replace the two nested loops below with a gemv.
+
         // Could be replaced with gemv.
         for (size_t i = 0; i < splinePoints.size(); ++i)
         {
@@ -729,6 +732,7 @@ void Splines::snapSpline(const LandBoundaries& landBoundary,
             }
         }
 
+        // Compute constraints.
         lambda = eMatrix * (bMatrix * atwaInverse * atwxb + cMatrix * atwaInverse * atwyb - dVector);
 
         rhsx = atwxb - bMatrix.transpose() * lambda;