IjkGridParametricRepresentation.h

/*-----------------------------------------------------------------------
Licensed to the Apache Software Foundation (ASF) under one
or more contributor license agreements.  See the NOTICE file
distributed with this work for additional information
regarding copyright ownership.  The ASF licenses this file
to you under the Apache License, Version 2.0 (the
"License"; you may not use this file except in compliance
with the License.  You may obtain a copy of the License at

  http://www.apache.org/licenses/LICENSE-2.0

Unless required by applicable law or agreed to in writing,
software distributed under the License is distributed on an
"AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
KIND, either express or implied.  See the License for the
specific language governing permissions and limitations
under the License.
-----------------------------------------------------------------------*/
#pragma once

#include "resqml2_0_1/AbstractIjkGridRepresentation.h"

#include "tools/BSpline.h"

namespace RESQML2_0_1_NS
{
    class DLL_IMPORT_OR_EXPORT IjkGridParametricRepresentation : public AbstractIjkGridRepresentation
    {
    private:
        void getXyzPointsOfPatchFromParametricPoints(gsoap_resqml2_0_1::resqml2__Point3dParametricArray* parametricPoint3d, double * xyzPoints) const;
    
        class PillarInformation
        {
        public:

            unsigned int maxControlPointCount;
            unsigned int parametricLineCount;
            unsigned int splitLineCount;
            double * controlPoints;
            double * controlPointParameters;
            short * pillarKind;
            unsigned int* pillarOfSplitCoordLines;
            std::vector< std::vector< geometry::BSpline > > splines;

            PillarInformation():maxControlPointCount(0), parametricLineCount(0), splitLineCount(0),
                    controlPoints(nullptr), controlPointParameters(nullptr), pillarKind(nullptr), pillarOfSplitCoordLines(nullptr) {}

            ~PillarInformation() {
                cleanMemory();
            }

            void cleanMemory() {
                if (controlPoints != nullptr) delete[] controlPoints;
                if (controlPointParameters != nullptr) delete[] controlPointParameters;
                if (pillarKind != nullptr) delete[] pillarKind;
                if (pillarOfSplitCoordLines != nullptr) delete[] pillarOfSplitCoordLines;
            }
        };

        void loadPillarInformation(PillarInformation & pillarInfo) const;

        gsoap_resqml2_0_1::resqml2__KDirection computeKDirection(double * controlPoints, const unsigned int & controlPointCountPerPillar);
    protected:
        PillarInformation* pillarInformation;
    
    public:
        IjkGridParametricRepresentation(soap* soapContext, RESQML2_NS::AbstractLocal3dCrs * crs,
            const std::string & guid, const std::string & title,
            const unsigned int & iCount, const unsigned int & jCount, const unsigned int & kCount,
            bool withTruncatedPillars = false);

        IjkGridParametricRepresentation(RESQML2_NS::AbstractFeatureInterpretation* interp, RESQML2_NS::AbstractLocal3dCrs * crs,
            const std::string & guid, const std::string & title,
            const unsigned int & iCount, const unsigned int & jCount, const unsigned int & kCount,
            bool withTruncatedPillars = false);

        IjkGridParametricRepresentation(gsoap_resqml2_0_1::_resqml2__IjkGridRepresentation* fromGsoap) : AbstractIjkGridRepresentation(fromGsoap), pillarInformation(nullptr) {}
        IjkGridParametricRepresentation(gsoap_resqml2_0_1::_resqml2__TruncatedIjkGridRepresentation* fromGsoap) : AbstractIjkGridRepresentation(fromGsoap), pillarInformation(nullptr) {}

        virtual ~IjkGridParametricRepresentation() { 
            if (pillarInformation != nullptr) 
                delete pillarInformation; 
        }

        std::string getHdfProxyUuid() const;

        ULONG64 getXyzPointCountOfPatch(const unsigned int & patchIndex) const;

        void getXyzPointsOfKInterfaceSequenceOfPatch(const unsigned int & kInterfaceStart, const unsigned int & kInterfaceEnd, const unsigned int & patchIndex, double * xyzPoints);

        void getXyzPointsOfBlockOfPatch(const unsigned int & patchIndex, double * xyzPoints);

        void getXyzPointsOfPatch(const unsigned int & patchIndex, double * xyzPoints) const;

        unsigned int getControlPointMaxCountPerPillar() const;

        void getControlPoints(double * controlPoints, bool reverseIAxis = false, bool reverseJAxis= false, bool reverseKAxis= false) const;

        bool hasControlPointParameters() const;

        void getControlPointParameters(double * controlPointParameters, bool reverseIAxis = false, bool reverseJAxis= false, bool reverseKAxis= false) const;

        bool isParametricLineKindConstant() const;

        /*
        * Get the constant parametric line count in the grid.
        */
        short getConstantParametricLineKind() const;

        void getParametricLineKind(short * pillarKind, bool reverseIAxis = false, bool reverseJAxis= false) const;

        void getParametersOfNodes(double * parameters, bool reverseIAxis = false, bool reverseJAxis= false, bool reverseKAxis= false) const;

        void setGeometryAsParametricNonSplittedPillarNodes(
            const gsoap_resqml2_0_1::resqml2__PillarShape & mostComplexPillarGeometry, const bool & isRightHanded,
            double * parameters, double * controlPoints, double * controlPointParameters, const unsigned int & controlPointMaxCountPerPillar, short * pillarKind, COMMON_NS::AbstractHdfProxy* proxy);

        void setGeometryAsParametricNonSplittedPillarNodesUsingExistingDatasets(
            const gsoap_resqml2_0_1::resqml2__PillarShape & mostComplexPillarGeometry, const gsoap_resqml2_0_1::resqml2__KDirection & kDirectionKind, const bool & isRightHanded,
            const std::string & parameters, const std::string & controlPoints, const std::string & controlPointParameters, const unsigned int & controlPointMaxCountPerPillar, const std::string & pillarKind, const std::string & definedPillars, COMMON_NS::AbstractHdfProxy* proxy);

        void setGeometryAsParametricSplittedPillarNodes(
            const gsoap_resqml2_0_1::resqml2__PillarShape & mostComplexPillarGeometry, const bool & isRightHanded,
            double * parameters, double * controlPoints, double * controlPointParameters, const unsigned int & controlPointMaxCountPerPillar, short * pillarKind, COMMON_NS::AbstractHdfProxy* proxy,
            const unsigned long & splitCoordinateLineCount, unsigned int * pillarOfCoordinateLine,
            unsigned int * splitCoordinateLineColumnCumulativeCount, unsigned int * splitCoordinateLineColumns);

        void setGeometryAsParametricSplittedPillarNodesUsingExistingDatasets(
            const gsoap_resqml2_0_1::resqml2__PillarShape & mostComplexPillarGeometry, const gsoap_resqml2_0_1::resqml2__KDirection & kDirectionKind, const bool & isRightHanded,
            const std::string & parameters, const std::string & controlPoints, const std::string & controlPointParameters, const unsigned int & controlPointMaxCountPerPillar, const std::string & pillarKind, const std::string & definedPillars, COMMON_NS::AbstractHdfProxy* proxy,
            const unsigned long & splitCoordinateLineCount, const std::string & pillarOfCoordinateLine,
            const std::string & splitCoordinateLineColumnCumulativeCount, const std::string & splitCoordinateLineColumns);

        void setGeometryAsParametricSplittedPillarNodes(const bool & isRightHanded,
            double * parameters, double * controlPoints, double * controlPointParameters, const unsigned int & controlPointCountPerPillar, short pillarKind, COMMON_NS::AbstractHdfProxy* proxy,
            const unsigned long & splitCoordinateLineCount, unsigned int * pillarOfCoordinateLine,
            unsigned int * splitCoordinateLineColumnCumulativeCount, unsigned int * splitCoordinateLineColumns);

        void setGeometryAsParametricSplittedPillarNodesUsingExistingDatasets(
            const gsoap_resqml2_0_1::resqml2__KDirection & kDirectionKind, const bool & isRightHanded,
            const std::string & parameters, const std::string & controlPoints, const std::string & controlPointParameters, const unsigned int & controlPointCountPerPillar, short pillarKind, COMMON_NS::AbstractHdfProxy* proxy,
            const unsigned long & splitCoordinateLineCount, const std::string & pillarOfCoordinateLine,
            const std::string & splitCoordinateLineColumnCumulativeCount, const std::string & splitCoordinateLineColumns);

        geometryKind getGeometryKind() const;
    };
}