Method InsertKnot

src/LNLib/Geometry/Volume/NurbsVolume.cpp:339–484

}

void LNLib::NurbsVolume::InsertKnot(const LN_NurbsVolume& volume, double insertKnot, int times, VolumeDirection direction, LN_NurbsVolume& result)
{
	int degreeU = volume.DegreeU;
	int degreeV = volume.DegreeV;
	int degreeW = volume.DegreeW;

	std::vector<double> kvU = volume.KnotVectorU;
	std::vector<double> kvV = volume.KnotVectorV;
	std::vector<double> kvW = volume.KnotVectorW;

	std::vector<std::vector<std::vector<LNLib::XYZW>>> controlPoints = volume.ControlPoints;

	int n = controlPoints.size() - 1;
	int m = controlPoints[0].size() - 1;
	int l = controlPoints[0][0].size() - 1;

	result.DegreeU = degreeU;
	result.DegreeV = degreeV;
	result.DegreeW = degreeW;

	std::vector<std::vector<std::vector<LNLib::XYZW>>> newControlPoints;
	if (direction == VolumeDirection::U)
	{
		newControlPoints.resize(n + times + 1, std::vector<std::vector<LNLib::XYZW>>(m + 1, std::vector<LNLib::XYZW>(l + 1)));
	}
	else if (direction == VolumeDirection::V)
	{
		newControlPoints.resize(n + 1, std::vector<std::vector<LNLib::XYZW>>(m + times + 1, std::vector<LNLib::XYZW>(l + 1)));
	}
	else if (direction == VolumeDirection::W)
	{
		newControlPoints.resize(n + 1, std::vector<std::vector<LNLib::XYZW>>(m + 1, std::vector<LNLib::XYZW>(l + times + 1)));
	}

	if (direction == VolumeDirection::U)
	{
		result.KnotVectorV = kvV;
		result.KnotVectorW = kvW;

		std::vector<double> newKvU;

		for (int sep = 0; sep <= l; sep++)
		{
			for (int col = 0; col <= m; col++)
			{
				std::vector<LNLib::XYZW> uControlPoints(n + 1);
				for (int i = 0; i <= n; ++i) {
					uControlPoints[i] = controlPoints[i][col][sep];
				}

				LNLib::LN_NurbsCurve curve;
				curve.Degree = degreeU;
				curve.KnotVector = kvU;
				curve.ControlPoints = uControlPoints;

				LNLib::LN_NurbsCurve resultCurve;
				LNLib::NurbsCurve::InsertKnot(curve, insertKnot, times, resultCurve);