diff --git a/libs/MeshKernel/include/MeshKernel/CurvilinearGrid/CurvilinearGridSplineToGrid.hpp b/libs/MeshKernel/include/MeshKernel/CurvilinearGrid/CurvilinearGridSplineToGrid.hpp index 6b4757277..a89bdf16b 100644 --- a/libs/MeshKernel/include/MeshKernel/CurvilinearGrid/CurvilinearGridSplineToGrid.hpp +++ b/libs/MeshKernel/include/MeshKernel/CurvilinearGrid/CurvilinearGridSplineToGrid.hpp @@ -48,7 +48,7 @@ namespace meshkernel private: /// @brief Maximum number of spline points allowed when doubling of the spline points. - static const UInt MaximumNumberOfSplinePoints = 10000; + static const UInt MaximumNumberOfSplinePoints = 2000; /// @brief The maximum number of checks for unlabeled splines static const UInt MaximumCumulativeUnlabeledSplineCount = 1000; diff --git a/libs/MeshKernel/src/CurvilinearGrid/CurvilinearGridSplineToGrid.cpp b/libs/MeshKernel/src/CurvilinearGrid/CurvilinearGridSplineToGrid.cpp index de17f9941..7ce6f5aa7 100644 --- a/libs/MeshKernel/src/CurvilinearGrid/CurvilinearGridSplineToGrid.cpp +++ b/libs/MeshKernel/src/CurvilinearGrid/CurvilinearGridSplineToGrid.cpp @@ -78,7 +78,9 @@ void meshkernel::CurvilinearGridSplineToGrid::DetermineIntersection(Splines& spl void meshkernel::CurvilinearGridSplineToGrid::DoubleSplinePoints(Splines& splines) const { std::vector doubledSplinePoints; - doubledSplinePoints.reserve(2 * splines.MaxSize() - 1); + doubledSplinePoints.reserve(2 * splines.MaxSizeIndex() - 1); + + std::cout << "maximum size: " << splines.Size(splines.MaxSizeIndex()) << std::endl; for (UInt splineIndex = 0; splineIndex < splines.GetNumSplines(); ++splineIndex) { @@ -108,7 +110,7 @@ bool meshkernel::CurvilinearGridSplineToGrid::ComputeAndCheckIntersection(Spline std::vector& splineType, VectorOfDoubleVectors& splineIntersections) const { - UInt maxSplineSize = splines.Size(splines.MaxSize()); + UInt maxSplineSize = splines.Size(splines.MaxSizeIndex()); double crossProduct; double normalisedIntersectionSplineI;