コード例 #1
0
ファイル: Flat.cs プロジェクト: olesar/Altaxo
        /// <inheritdoc />
        public void AddGeometry(
            Action <PointD3D, VectorD3D> AddPositionAndNormal,
            Action <int, int, int, bool> AddIndices,
            ref int vertexIndexOffset,
            bool isStartCap,
            PointD3D basePoint,
            VectorD3D eastVector,
            VectorD3D northVector,
            VectorD3D forwardVectorNormalized,
            ICrossSectionOfLine lineCrossSection,
            PointD3D[] baseCrossSectionPositions,
            VectorD3D[] baseCrossSectionNormals,
            ref object temporaryStorageSpace)
        {
            var crossSectionVertexCount = lineCrossSection.NumberOfVertices;
            var crossSectionNormalCount = lineCrossSection.NumberOfNormals;

            var capCrossSectionPositions = baseCrossSectionPositions ?? (PointD3D[])temporaryStorageSpace ?? (PointD3D[])(temporaryStorageSpace = new PointD3D[crossSectionVertexCount]);

            if (null == baseCrossSectionPositions) // if null the positions were not provided
            {
                var matrix = Math3D.Get2DProjectionToPlane(eastVector, northVector, basePoint);
                for (int i = 0, j = 0; i < crossSectionVertexCount; ++i, ++j)
                {
                    capCrossSectionPositions[i] = matrix.Transform(lineCrossSection.Vertices(i));
                }
            }

            AddGeometry(AddPositionAndNormal, AddIndices, ref vertexIndexOffset, isStartCap, basePoint, forwardVectorNormalized, capCrossSectionPositions);
        }
コード例 #2
0
ファイル: CrossSectionBase.cs プロジェクト: olesar/Altaxo
        /// <summary>
        /// Gets the vertices of a cross section from a certain vertex index to another vertex index (including the last vertex index).
        /// </summary>
        /// <param name="crossSection">The cross section.</param>
        /// <param name="firstVertexIndex">The first vertex index to include. Must be greater than or equal to zero.</param>
        /// <param name="lastIncludedVertexIndex">The last vertex index to include. Must be greater than <paramref name="firstVertexIndex"/>, but can be greater than the number of vertices of the cross section in order to describe a polygon that 'wraps around'.</param>
        /// <returns>The vertices from the first index up to and including the last index.</returns>
        /// <exception cref="ArgumentOutOfRangeException"></exception>
        /// <exception cref="ArgumentException"></exception>
        public static IEnumerable <PointD2D> GetVerticesFromToIncluding(ICrossSectionOfLine crossSection, int firstVertexIndex, int lastIncludedVertexIndex)
        {
            if (!(firstVertexIndex <= lastIncludedVertexIndex))
            {
                throw new ArgumentOutOfRangeException(nameof(lastIncludedVertexIndex) + " must be greater then or equal to " + nameof(firstVertexIndex));
            }
            if (!(firstVertexIndex >= 0))
            {
                throw new ArgumentException(nameof(firstVertexIndex) + "must be >=0");
            }

            int vertexCount = crossSection.NumberOfVertices;

            for (int i = firstVertexIndex; i <= lastIncludedVertexIndex; ++i)
            {
                yield return(crossSection.Vertices(i % vertexCount));
            }
        }
コード例 #3
0
ファイル: Flat.cs プロジェクト: olesar/Altaxo
        /// <summary>
        /// Adds the triangle geometry for this cap.
        /// </summary>
        /// <param name="AddPositionAndNormal">The procedure to add a vertex position and normal.</param>
        /// <param name="AddIndices">The procedure to add vertex indices for one triangle.</param>
        /// <param name="vertexIndexOffset">The vertex index offset. Must be actualized during this call.</param>
        /// <param name="isStartCap">If set to <c>true</c>, a start cap is drawn; otherwise, an end cap is drawn.</param>
        /// <param name="basePoint">The location of the middle point of the line at the cap's location.</param>
        /// <param name="westVector">The west vector for orientation of the cross section.</param>
        /// <param name="northVector">The north vector for orientation of the cross section.</param>
        /// <param name="forwardVector">The forward vector for orientation of the cross section.</param>
        /// <param name="crossSection">The cross section of the line.</param>
        public static void AddGeometry(
            Action <PointD3D, VectorD3D> AddPositionAndNormal,
            Action <int, int, int, bool> AddIndices,
            ref int vertexIndexOffset,
            bool isStartCap,
            PointD3D basePoint,
            VectorD3D westVector,
            VectorD3D northVector,
            VectorD3D forwardVector,
            ICrossSectionOfLine crossSection
            )
        {
            if (isStartCap)
            {
                forwardVector = -forwardVector;
            }

            var matrix = Matrix4x3.NewFromBasisVectorsAndLocation(westVector, northVector, forwardVector, basePoint);

            int currIndex = vertexIndexOffset;
            int crossSectionPositionCount = crossSection.NumberOfVertices;

            // Add the midpoint
            // add the middle point of the end cap and the normal of the end cap
            AddPositionAndNormal(basePoint, forwardVector);
            ++currIndex;

            for (int i = 0; i < crossSectionPositionCount; ++i)
            {
                var sp = matrix.Transform(crossSection.Vertices(i));
                AddPositionAndNormal(sp, forwardVector);

                AddIndices(
                    currIndex - 1, // mid point of the end cap
                    currIndex + i,
                    currIndex + (1 + i) % crossSectionPositionCount,
                    isStartCap);
            }

            currIndex        += crossSectionPositionCount;
            vertexIndexOffset = currIndex;
        }
コード例 #4
0
        public static void Add(
            Action <PointD3D, VectorD3D> AddPositionAndNormal,
            Action <int, int, int, bool> AddIndices,
            ref int vertexIndexOffset,
            bool isStartCap,
            PointD3D basePoint,
            VectorD3D westVector,
            VectorD3D northVector,
            VectorD3D forwardVectorNormalized,
            ICrossSectionOfLine lineCrossSection,
            PointD3D[] crossSectionPositions,
            VectorD3D[] crossSectionNormals,
            ref object temporaryStorageSpace,
            ILineCapContour capContour
            )
        {
            var crossSectionVertexCount = lineCrossSection.NumberOfVertices;
            var crossSectionNormalCount = lineCrossSection.NumberOfNormals;
            var contourZScale           = 0.5 * Math.Max(lineCrossSection.Size1, lineCrossSection.Size2);

            // do we need a flat end at the beginning of the cap?
            if (null == crossSectionPositions &&         // if lineCrossSectionPositions are null, it means that our cap is not connected to the line and needs a flat end
                capContour.Vertices(0) == _pointD2D_0_1) // furthermore the cap assumes to be started at the cross section
            {
                // the parameter isStartCap must be negated, because this flat cap is the "counterpart" of our cap to draw
                Flat.AddGeometry(
                    AddPositionAndNormal,
                    AddIndices,
                    ref vertexIndexOffset,
                    !isStartCap,
                    basePoint,
                    westVector,
                    northVector,
                    forwardVectorNormalized,
                    lineCrossSection);
            }

            if (isStartCap)
            {
                forwardVectorNormalized = -forwardVectorNormalized;
            }

            var contourVertexCount = capContour.NumberOfVertices;
            var contourNormalCount = capContour.NumberOfNormals;

            // now the calculation can start

            CrossSectionCases previousCrossSectionType = CrossSectionCases.MiddlePointSmooth;
            int  previousGeneratedPoints            = 0;
            int  previousContourVertexIndex         = 0;
            bool isOnSecondSideOfContourVertexSharp = true;

            for (int contourVertexIndex = 0, contourNormalIndex = 0; contourVertexIndex < contourVertexCount; ++contourVertexIndex, ++contourNormalIndex)
            {
                // we have 4 different situations here:
                // 1st) the crossSection.Y is zero, thus this is the middle point (the normal should then go in z-direction) -> we need only one single vertex and normal for that
                // 2nd) the countour normal is in x-direction  -> we need only a point for each crossSection vertex, but not for each crossSectionNormal
                // 3rd) the regular case -> we need a point for each crossSection normal

                var capContourVertex = capContour.Vertices(contourVertexIndex);
                var capContourNormal = capContour.Normals(contourNormalIndex);

                CrossSectionCases currentCrossSectionType;
                if (capContourVertex.Y == 0 && capContourNormal.Y == 0)
                {
                    currentCrossSectionType = CrossSectionCases.MiddlePointSmooth;
                }
                else if (capContourVertex.Y == 0)
                {
                    currentCrossSectionType = CrossSectionCases.MiddlePointSharp;
                }
                else if (0 == capContourNormal.Y)
                {
                    currentCrossSectionType = CrossSectionCases.VerticesOnly;
                }
                else
                {
                    currentCrossSectionType = CrossSectionCases.Regular;
                }

                var currentLocation = basePoint + forwardVectorNormalized * capContourVertex.X * contourZScale;
                var matrix          = Matrix4x3.NewFromBasisVectorsAndLocation(westVector, northVector, forwardVectorNormalized, currentLocation);

                int currentGeneratedPoints = 0;
                switch (currentCrossSectionType)
                {
                case CrossSectionCases.MiddlePointSmooth:
                {
                    var position = matrix.Transform(PointD2D.Empty);
                    var normal   = matrix.Transform(new VectorD3D(0, 0, capContourNormal.X));
                    AddPositionAndNormal(position, normal);
                    currentGeneratedPoints = 1;
                }
                break;

                case CrossSectionCases.MiddlePointSharp:
                {
                    for (int i = 0, j = 0; i < crossSectionVertexCount; ++i, ++j)
                    {
                        var normal1  = (i == 0) ? lineCrossSection.Normals(crossSectionNormalCount - 1) : lineCrossSection.Normals(j - 1);
                        var normal2  = lineCrossSection.Normals(j);
                        var sn       = (normal1 + normal2).Normalized;
                        var utNormal = GetNormalVector(lineCrossSection.Vertices(i), sn, capContourNormal, contourZScale);
                        AddPositionAndNormal(currentLocation, matrix.Transform(utNormal)); // store the tip point with the averaged normal
                        if (lineCrossSection.IsVertexSharp(i))
                        {
                            ++j;
                        }
                    }
                    currentGeneratedPoints = crossSectionVertexCount;
                }
                break;

                case CrossSectionCases.VerticesOnly:
                {
                    var commonNormal = matrix.Transform(new VectorD3D(0, 0, capContourNormal.X));
                    for (int i = 0; i < crossSectionVertexCount; ++i)
                    {
                        var position = matrix.Transform(lineCrossSection.Vertices(i) * capContourVertex.Y);
                        AddPositionAndNormal(position, commonNormal);
                    }
                    currentGeneratedPoints = crossSectionVertexCount;
                }
                break;

                case CrossSectionCases.Regular:
                {
                    for (int i = 0, j = 0; i < crossSectionVertexCount; ++i, ++j)
                    {
                        var sp       = lineCrossSection.Vertices(i);
                        var sn       = lineCrossSection.Normals(j);
                        var utNormal = GetNormalVector(sp, sn, capContourNormal, contourZScale);
                        var position = matrix.Transform(sp * capContourVertex.Y);
                        var normal   = matrix.Transform(utNormal);

                        AddPositionAndNormal(position, normal);

                        if (lineCrossSection.IsVertexSharp(i))
                        {
                            ++j;
                            sn       = lineCrossSection.Normals(j);
                            utNormal = GetNormalVector(sp, sn, capContourNormal, contourZScale);
                            normal   = matrix.Transform(utNormal);
                            AddPositionAndNormal(position, normal);
                        }
                    }
                    currentGeneratedPoints = crossSectionNormalCount;
                }
                break;

                default:
                    throw new NotImplementedException();
                }
                vertexIndexOffset += currentGeneratedPoints;

                // now we start generating triangles

                if (contourVertexIndex > previousContourVertexIndex)
                {
                    int voffset1 = vertexIndexOffset - currentGeneratedPoints;
                    int voffset0 = voffset1 - previousGeneratedPoints;
                    switch (previousCrossSectionType)
                    {
                    case CrossSectionCases.MiddlePointSmooth:
                    {
                        switch (currentCrossSectionType)
                        {
                        case CrossSectionCases.MiddlePointSmooth: // Middle point to middle point
                        {
                            // no triangles, since from middle point to middle point we have an infinity thin line
                        }
                        break;

                        case CrossSectionCases.MiddlePointSharp: // Middle point to middle point
                        {
                            // no triangles, since from middle point to middle point we have an infinity thin line
                        }
                        break;

                        case CrossSectionCases.VerticesOnly: // Middle point to vertices only
                        {
                            for (int i = 0; i < crossSectionVertexCount; ++i)
                            {
                                AddIndices(voffset0, voffset1 + i, voffset1 + (i + 1) % crossSectionVertexCount, isStartCap);
                            }
                        }
                        break;

                        case CrossSectionCases.Regular: // Middle point to regular
                        {
                            for (int i = 0; i < crossSectionNormalCount; ++i)
                            {
                                AddIndices(voffset0, voffset1 + i, voffset1 + (i + 1) % crossSectionNormalCount, isStartCap);
                            }
                        }
                        break;

                        default:
                            throw new NotImplementedException();
                        }
                    }
                    break;

                    case CrossSectionCases.MiddlePointSharp:
                    {
                        switch (currentCrossSectionType)
                        {
                        case CrossSectionCases.MiddlePointSmooth: // Middle point to middle point
                        {
                            // no triangles, since from middle point to middle point we have an infinity thin line
                        }
                        break;

                        case CrossSectionCases.MiddlePointSharp: // Middle point to middle point
                        {
                            // no triangles, since from middle point to middle point we have an infinity thin line
                        }
                        break;

                        case CrossSectionCases.VerticesOnly: // MiddlePointSharp to VerticesOnly
                        {
                            for (int i = 0; i < crossSectionVertexCount; ++i)
                            {
                                AddIndices(voffset0, voffset1 + i, voffset1 + (i + 1) % crossSectionVertexCount, isStartCap);
                            }
                        }
                        break;

                        case CrossSectionCases.Regular: // MiddlePointSharp to Regular
                        {
                            for (int i = 0, j = 0; i < crossSectionNormalCount; ++i, ++j)
                            {
                                AddIndices(voffset0 + i, voffset1 + i, voffset1 + (i + 1) % crossSectionNormalCount, isStartCap);
                            }
                        }
                        break;

                        default:
                            throw new NotImplementedException();
                        }
                    }
                    break;

                    case CrossSectionCases.VerticesOnly:
                    {
                        switch (currentCrossSectionType)
                        {
                        case CrossSectionCases.MiddlePointSmooth: // VerticesOnly to MiddlePoint
                        {
                            for (int i = 0; i < crossSectionVertexCount; ++i)
                            {
                                AddIndices(voffset1, voffset0 + i, voffset0 + (i + 1) % crossSectionVertexCount, isStartCap);
                            }
                        }
                        break;

                        case CrossSectionCases.VerticesOnly: // VerticesOnly to VerticesOnly
                        {
                            for (int i = 0; i < crossSectionVertexCount; ++i)
                            {
                                AddIndices(voffset0 + ((i == 0) ? crossSectionVertexCount - 1 : i - 1), voffset0 + i, voffset1 + i, isStartCap);
                                AddIndices(voffset0 + ((i == 0) ? crossSectionVertexCount - 1 : i - 1), voffset1 + i, voffset1 + ((i == 0) ? crossSectionVertexCount - 1 : i - 1), isStartCap);
                            }
                        }
                        break;

                        case CrossSectionCases.Regular: // VerticesOnly to regular
                        {
                            throw new NotImplementedException();
                        }
                        //break;

                        default:
                            throw new NotImplementedException();
                        }
                    }
                    break;

                    case CrossSectionCases.Regular:
                    {
                        switch (currentCrossSectionType)
                        {
                        case CrossSectionCases.MiddlePointSmooth: // Regular to MiddlePointOnly
                        {
                            for (int i = 0, j = 0; i < crossSectionVertexCount; ++i, ++j)
                            {
                                AddIndices(voffset0 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), voffset0 + j, voffset1, isStartCap);

                                if (lineCrossSection.IsVertexSharp(i))
                                {
                                    ++j;
                                }
                            }
                        }
                        break;

                        case CrossSectionCases.MiddlePointSharp: // Regular to MiddlePointSharp
                        {
                            for (int i = 0, j = 0; i < crossSectionVertexCount; ++i, ++j)
                            {
                                AddIndices(voffset0 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), voffset0 + j, voffset1 + i, isStartCap);

                                if (lineCrossSection.IsVertexSharp(i))
                                {
                                    ++j;
                                }
                            }
                        }
                        break;

                        case CrossSectionCases.VerticesOnly: // Regular to VerticesOnly
                        {
                            for (int i = 0, j = 0; i < crossSectionVertexCount; ++i, ++j)
                            {
                                AddIndices(voffset0 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), voffset0 + j, voffset1 + i, isStartCap);
                                AddIndices(voffset0 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), voffset1 + i, voffset1 + ((i == 0) ? crossSectionVertexCount - 1 : i - 1), isStartCap);

                                if (lineCrossSection.IsVertexSharp(i))
                                {
                                    ++j;
                                }
                            }
                        }
                        break;

                        case CrossSectionCases.Regular: // Regular to Regular
                        {
                            for (int i = 0, j = 0; i < crossSectionVertexCount; ++i, ++j)
                            {
                                AddIndices(voffset0 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), voffset0 + j, voffset1 + j, isStartCap);
                                AddIndices(voffset0 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), voffset1 + j, voffset1 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), isStartCap);

                                if (lineCrossSection.IsVertexSharp(i))
                                {
                                    ++j;
                                }
                            }
                        }
                        break;

                        default:
                            throw new NotImplementedException();
                        }
                    }
                    break;

                    default:
                        throw new NotImplementedException();
                    }
                }

                if (!isOnSecondSideOfContourVertexSharp && capContour.IsVertexSharp(contourVertexIndex) && contourVertexIndex < (contourVertexCount - 1))
                {
                    previousContourVertexIndex = contourVertexIndex;
                    --contourVertexIndex; // trick: decrement the vertex index, it is incremented then again in the following for loop, so that contourVertexIndex stays constant
                    isOnSecondSideOfContourVertexSharp = true;
                    continue;
                }
                isOnSecondSideOfContourVertexSharp = false;

                // now we switch the current calculated positions and normals with the old ones
                previousCrossSectionType   = currentCrossSectionType;
                previousGeneratedPoints    = currentGeneratedPoints;
                previousContourVertexIndex = contourVertexIndex;
            }
        }
コード例 #5
0
        /// <summary>
        /// Adds the triangle geometry. Here, the position of the startcap base and of the endcap base is already calculated and provided in the arguments.
        /// </summary>
        /// <param name="AddPositionAndNormal">The procedure to add a vertex position and normal.</param>
        /// <param name="AddIndices">The procedure to add vertex indices for one triangle.</param>
        /// <param name="vertexIndexOffset">The vertex index offset.</param>
        /// <param name="lineStart">The line start. This is the precalculated base of the start line cap.</param>
        /// <param name="lineEnd">The line end. Here, this is the precalculated base of the end line cap.</param>
        /// <param name="drawLine">If this parameter is true, the line segment between lineStart and lineEnd is drawn. If false, the line segment itself is not drawn, but the start end end caps are drawn.</param>
        /// <param name="overrideStartCap">If not null, this parameter override the start cap that is stored in this class.</param>
        /// <param name="overrideEndCap">If not null, this parameter overrides the end cap that is stored in this class.</param>
        /// <exception cref="System.InvalidProgramException">The structure is not initialized yet. Call Initialize before using it!</exception>
        public void AddGeometry(
            Action <PointD3D, VectorD3D> AddPositionAndNormal,
            Action <int, int, int, bool> AddIndices,
            ref int vertexIndexOffset,
            PointD3D lineStart,
            PointD3D lineEnd,
            bool drawLine,
            ILineCap overrideStartCap,
            ILineCap overrideEndCap
            )
        {
            if (null == _lastNormalsTransformed)
            {
                throw new InvalidProgramException("The structure is not initialized yet. Call Initialize before using it!");
            }

            var resultingStartCap = overrideStartCap ?? _dashStartCap;
            var resultingEndCap   = overrideEndCap ?? _dashEndCap;

            // draw the straight line if the remaining line length is >0
            if (drawLine)
            {
                // Get the matrix for the start plane
                var matrix = Math3D.Get2DProjectionToPlane(_westVector, _northVector, lineStart);
                for (int i = 0; i < _crossSectionVertexCount; ++i)
                {
                    _lastPositionsTransformedStart[i] = matrix.Transform(_crossSection.Vertices(i));
                }

                matrix = Math3D.Get2DProjectionToPlane(_westVector, _northVector, lineEnd);
                for (int i = 0; i < _crossSectionVertexCount; ++i)
                {
                    _lastPositionsTransformedEnd[i] = matrix.Transform(_crossSection.Vertices(i));
                }

                // draw the line segment now
                var currIndex = vertexIndexOffset;
                for (int i = 0, j = 0; i < _crossSectionVertexCount; ++i, ++j)
                {
                    if (j == 0)
                    {
                        AddIndices(currIndex, currIndex + 1, currIndex + 2 * _crossSectionNormalCount - 2, false);
                        AddIndices(currIndex + 2 * _crossSectionNormalCount - 2, currIndex + 1, currIndex + 2 * _crossSectionNormalCount - 1, false);
                    }
                    else
                    {
                        AddIndices(currIndex, currIndex + 1, currIndex - 2, false);
                        AddIndices(currIndex - 2, currIndex + 1, currIndex - 1, false);
                    }

                    AddPositionAndNormal(_lastPositionsTransformedStart[i], _lastNormalsTransformed[j]);
                    AddPositionAndNormal(_lastPositionsTransformedEnd[i], _lastNormalsTransformed[j]);
                    currIndex += 2;

                    if (_crossSection.IsVertexSharp(i))
                    {
                        ++j;
                        AddPositionAndNormal(_lastPositionsTransformedStart[i], _lastNormalsTransformed[j]);
                        AddPositionAndNormal(_lastPositionsTransformedEnd[i], _lastNormalsTransformed[j]);
                        currIndex += 2;
                    }
                }
                vertexIndexOffset = currIndex;
            }

            // now the start cap
            if (null != resultingStartCap)
            {
                resultingStartCap.AddGeometry(
                    AddPositionAndNormal,
                    AddIndices,
                    ref vertexIndexOffset,
                    true,
                    lineStart,
                    _westVector,
                    _northVector,
                    _forwardVector,
                    _crossSection,
                    drawLine ? _lastPositionsTransformedStart : null,
                    _lastNormalsTransformed,
                    ref _startCapTemporaryStorageSpace);
            }
            else if (drawLine)
            {
                LineCaps.Flat.AddGeometry(
                    AddPositionAndNormal,
                    AddIndices,
                    ref vertexIndexOffset,
                    true,
                    lineStart,
                    _forwardVector,
                    _lastPositionsTransformedStart
                    );
            }

            if (null != resultingEndCap)
            {
                resultingEndCap.AddGeometry(
                    AddPositionAndNormal,
                    AddIndices,
                    ref vertexIndexOffset,
                    false,
                    lineEnd,
                    _westVector,
                    _northVector,
                    _forwardVector,
                    _crossSection,
                    drawLine ? _lastPositionsTransformedEnd : null,
                    _lastNormalsTransformed,
                    ref _endCapTemporaryStorageSpace);
            }
            else if (drawLine)
            {
                LineCaps.Flat.AddGeometry(
                    AddPositionAndNormal,
                    AddIndices,
                    ref vertexIndexOffset,
                    false,
                    lineEnd,
                    _forwardVector,
                    _lastPositionsTransformedEnd
                    );
            }
        }
コード例 #6
0
		public static void Add(
			Action<PointD3D, VectorD3D> AddPositionAndNormal,
			Action<int, int, int, bool> AddIndices,
			ref int vertexIndexOffset,
			bool isStartCap,
			PointD3D basePoint,
			VectorD3D westVector,
			VectorD3D northVector,
			VectorD3D forwardVectorNormalized,
			ICrossSectionOfLine lineCrossSection,
			PointD3D[] crossSectionPositions,
			VectorD3D[] crossSectionNormals,
			ref object temporaryStorageSpace,
			ILineCapContour capContour
			)
		{
			var crossSectionVertexCount = lineCrossSection.NumberOfVertices;
			var crossSectionNormalCount = lineCrossSection.NumberOfNormals;
			var contourZScale = 0.5 * Math.Max(lineCrossSection.Size1, lineCrossSection.Size2);

			// do we need a flat end at the beginning of the cap?
			if (null == crossSectionPositions && // if lineCrossSectionPositions are null, it means that our cap is not connected to the line and needs a flat end
						capContour.Vertices(0) == _pointD2D_0_1) // furthermore the cap assumes to be started at the cross section
			{
				// the parameter isStartCap must be negated, because this flat cap is the "counterpart" of our cap to draw
				Flat.AddGeometry(
					AddPositionAndNormal,
					AddIndices,
					ref vertexIndexOffset,
					!isStartCap,
					basePoint,
					westVector,
					northVector,
					forwardVectorNormalized,
					lineCrossSection);
			}

			if (isStartCap)
				forwardVectorNormalized = -forwardVectorNormalized;

			var contourVertexCount = capContour.NumberOfVertices;
			var contourNormalCount = capContour.NumberOfNormals;

			// now the calculation can start

			CrossSectionCases previousCrossSectionType = CrossSectionCases.MiddlePointSmooth;
			int previousGeneratedPoints = 0;
			int previousContourVertexIndex = 0;
			bool isOnSecondSideOfContourVertexSharp = true;
			for (int contourVertexIndex = 0, contourNormalIndex = 0; contourVertexIndex < contourVertexCount; ++contourVertexIndex, ++contourNormalIndex)
			{
				// we have 4 different situations here:
				// 1st) the crossSection.Y is zero, thus this is the middle point (the normal should then go in z-direction) -> we need only one single vertex and normal for that
				// 2nd) the countour normal is in x-direction  -> we need only a point for each crossSection vertex, but not for each crossSectionNormal
				// 3rd) the regular case -> we need a point for each crossSection normal

				var capContourVertex = capContour.Vertices(contourVertexIndex);
				var capContourNormal = capContour.Normals(contourNormalIndex);

				CrossSectionCases currentCrossSectionType;
				if (capContourVertex.Y == 0 && capContourNormal.Y == 0)
					currentCrossSectionType = CrossSectionCases.MiddlePointSmooth;
				else if (capContourVertex.Y == 0)
					currentCrossSectionType = CrossSectionCases.MiddlePointSharp;
				else if (0 == capContourNormal.Y)
					currentCrossSectionType = CrossSectionCases.VerticesOnly;
				else
					currentCrossSectionType = CrossSectionCases.Regular;

				var currentLocation = basePoint + forwardVectorNormalized * capContourVertex.X * contourZScale;
				var matrix = Matrix4x3.NewFromBasisVectorsAndLocation(westVector, northVector, forwardVectorNormalized, currentLocation);

				int currentGeneratedPoints = 0;
				switch (currentCrossSectionType)
				{
					case CrossSectionCases.MiddlePointSmooth:
						{
							var position = matrix.Transform(PointD2D.Empty);
							var normal = matrix.Transform(new VectorD3D(0, 0, capContourNormal.X));
							AddPositionAndNormal(position, normal);
							currentGeneratedPoints = 1;
						}
						break;

					case CrossSectionCases.MiddlePointSharp:
						{
							for (int i = 0, j = 0; i < crossSectionVertexCount; ++i, ++j)
							{
								var normal1 = (i == 0) ? lineCrossSection.Normals(crossSectionNormalCount - 1) : lineCrossSection.Normals(j - 1);
								var normal2 = lineCrossSection.Normals(j);
								var sn = (normal1 + normal2).Normalized;
								var utNormal = GetNormalVector(lineCrossSection.Vertices(i), sn, capContourNormal, contourZScale);
								AddPositionAndNormal(currentLocation, matrix.Transform(utNormal)); // store the tip point with the averaged normal
								if (lineCrossSection.IsVertexSharp(i))
								{
									++j;
								}
							}
							currentGeneratedPoints = crossSectionVertexCount;
						}
						break;

					case CrossSectionCases.VerticesOnly:
						{
							var commonNormal = matrix.Transform(new VectorD3D(0, 0, capContourNormal.X));
							for (int i = 0; i < crossSectionVertexCount; ++i)
							{
								var position = matrix.Transform(lineCrossSection.Vertices(i) * capContourVertex.Y);
								AddPositionAndNormal(position, commonNormal);
							}
							currentGeneratedPoints = crossSectionVertexCount;
						}
						break;

					case CrossSectionCases.Regular:
						{
							for (int i = 0, j = 0; i < crossSectionVertexCount; ++i, ++j)
							{
								var sp = lineCrossSection.Vertices(i);
								var sn = lineCrossSection.Normals(j);
								var utNormal = GetNormalVector(sp, sn, capContourNormal, contourZScale);
								var position = matrix.Transform(sp * capContourVertex.Y);
								var normal = matrix.Transform(utNormal);

								AddPositionAndNormal(position, normal);

								if (lineCrossSection.IsVertexSharp(i))
								{
									++j;
									sn = lineCrossSection.Normals(j);
									utNormal = GetNormalVector(sp, sn, capContourNormal, contourZScale);
									normal = matrix.Transform(utNormal);
									AddPositionAndNormal(position, normal);
								}
							}
							currentGeneratedPoints = crossSectionNormalCount;
						}
						break;

					default:
						throw new NotImplementedException();
				}
				vertexIndexOffset += currentGeneratedPoints;

				// now we start generating triangles

				if (contourVertexIndex > previousContourVertexIndex)
				{
					int voffset1 = vertexIndexOffset - currentGeneratedPoints;
					int voffset0 = voffset1 - previousGeneratedPoints;
					switch (previousCrossSectionType)
					{
						case CrossSectionCases.MiddlePointSmooth:
							{
								switch (currentCrossSectionType)
								{
									case CrossSectionCases.MiddlePointSmooth: // Middle point to middle point
										{
											// no triangles, since from middle point to middle point we have an infinity thin line
										}
										break;

									case CrossSectionCases.MiddlePointSharp: // Middle point to middle point
										{
											// no triangles, since from middle point to middle point we have an infinity thin line
										}
										break;

									case CrossSectionCases.VerticesOnly: // Middle point to vertices only
										{
											for (int i = 0; i < crossSectionVertexCount; ++i)
											{
												AddIndices(voffset0, voffset1 + i, voffset1 + (i + 1) % crossSectionVertexCount, isStartCap);
											}
										}
										break;

									case CrossSectionCases.Regular: // Middle point to regular
										{
											for (int i = 0; i < crossSectionNormalCount; ++i)
											{
												AddIndices(voffset0, voffset1 + i, voffset1 + (i + 1) % crossSectionNormalCount, isStartCap);
											}
										}
										break;

									default:
										throw new NotImplementedException();
								}
							}
							break;

						case CrossSectionCases.MiddlePointSharp:
							{
								switch (currentCrossSectionType)
								{
									case CrossSectionCases.MiddlePointSmooth: // Middle point to middle point
										{
											// no triangles, since from middle point to middle point we have an infinity thin line
										}
										break;

									case CrossSectionCases.MiddlePointSharp: // Middle point to middle point
										{
											// no triangles, since from middle point to middle point we have an infinity thin line
										}
										break;

									case CrossSectionCases.VerticesOnly: // MiddlePointSharp to VerticesOnly
										{
											for (int i = 0; i < crossSectionVertexCount; ++i)
											{
												AddIndices(voffset0, voffset1 + i, voffset1 + (i + 1) % crossSectionVertexCount, isStartCap);
											}
										}
										break;

									case CrossSectionCases.Regular: // MiddlePointSharp to Regular
										{
											for (int i = 0, j = 0; i < crossSectionNormalCount; ++i, ++j)
											{
												AddIndices(voffset0 + i, voffset1 + i, voffset1 + (i + 1) % crossSectionNormalCount, isStartCap);
											}
										}
										break;

									default:
										throw new NotImplementedException();
								}
							}
							break;

						case CrossSectionCases.VerticesOnly:
							{
								switch (currentCrossSectionType)
								{
									case CrossSectionCases.MiddlePointSmooth: // VerticesOnly to MiddlePoint
										{
											for (int i = 0; i < crossSectionVertexCount; ++i)
											{
												AddIndices(voffset1, voffset0 + i, voffset0 + (i + 1) % crossSectionVertexCount, isStartCap);
											}
										}
										break;

									case CrossSectionCases.VerticesOnly: // VerticesOnly to VerticesOnly
										{
											for (int i = 0; i < crossSectionVertexCount; ++i)
											{
												AddIndices(voffset0 + ((i == 0) ? crossSectionVertexCount - 1 : i - 1), voffset0 + i, voffset1 + i, isStartCap);
												AddIndices(voffset0 + ((i == 0) ? crossSectionVertexCount - 1 : i - 1), voffset1 + i, voffset1 + ((i == 0) ? crossSectionVertexCount - 1 : i - 1), isStartCap);
											}
										}
										break;

									case CrossSectionCases.Regular: // VerticesOnly to regular
										{
											throw new NotImplementedException();
										}
									//break;

									default:
										throw new NotImplementedException();
								}
							}
							break;

						case CrossSectionCases.Regular:
							{
								switch (currentCrossSectionType)
								{
									case CrossSectionCases.MiddlePointSmooth: // Regular to MiddlePointOnly
										{
											for (int i = 0, j = 0; i < crossSectionVertexCount; ++i, ++j)
											{
												AddIndices(voffset0 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), voffset0 + j, voffset1, isStartCap);

												if (lineCrossSection.IsVertexSharp(i))
													++j;
											}
										}
										break;

									case CrossSectionCases.MiddlePointSharp: // Regular to MiddlePointSharp
										{
											for (int i = 0, j = 0; i < crossSectionVertexCount; ++i, ++j)
											{
												AddIndices(voffset0 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), voffset0 + j, voffset1 + i, isStartCap);

												if (lineCrossSection.IsVertexSharp(i))
													++j;
											}
										}
										break;

									case CrossSectionCases.VerticesOnly: // Regular to VerticesOnly
										{
											for (int i = 0, j = 0; i < crossSectionVertexCount; ++i, ++j)
											{
												AddIndices(voffset0 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), voffset0 + j, voffset1 + i, isStartCap);
												AddIndices(voffset0 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), voffset1 + i, voffset1 + ((i == 0) ? crossSectionVertexCount - 1 : i - 1), isStartCap);

												if (lineCrossSection.IsVertexSharp(i))
													++j;
											}
										}
										break;

									case CrossSectionCases.Regular: // Regular to Regular
										{
											for (int i = 0, j = 0; i < crossSectionVertexCount; ++i, ++j)
											{
												AddIndices(voffset0 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), voffset0 + j, voffset1 + j, isStartCap);
												AddIndices(voffset0 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), voffset1 + j, voffset1 + ((j == 0) ? crossSectionNormalCount - 1 : j - 1), isStartCap);

												if (lineCrossSection.IsVertexSharp(i))
													++j;
											}
										}
										break;

									default:
										throw new NotImplementedException();
								}
							}
							break;

						default:
							throw new NotImplementedException();
					}
				}

				if (!isOnSecondSideOfContourVertexSharp && capContour.IsVertexSharp(contourVertexIndex) && contourVertexIndex < (contourVertexCount - 1))
				{
					previousContourVertexIndex = contourVertexIndex;
					--contourVertexIndex; // trick: decrement the vertex index, it is incremented then again in the following for loop, so that contourVertexIndex stays constant
					isOnSecondSideOfContourVertexSharp = true;
					continue;
				}
				isOnSecondSideOfContourVertexSharp = false;

				// now we switch the current calculated positions and normals with the old ones
				previousCrossSectionType = currentCrossSectionType;
				previousGeneratedPoints = currentGeneratedPoints;
				previousContourVertexIndex = contourVertexIndex;
			}
		}