nurbsdisk.h

/*
   OpenMBV - Open Multi Body Viewer.
   Copyright (C) 2009 Markus Friedrich
 
  This library is free software; you can redistribute it and/or 
  modify it under the terms of the GNU Lesser General Public 
  License as published by the Free Software Foundation; either 
  version 2.1 of the License, or (at your option) any later version. 
   
  This library is distributed in the hope that it will be useful, 
  but WITHOUT ANY WARRANTY; without even the implied warranty of 
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU 
  Lesser General Public License for more details. 
   
  You should have received a copy of the GNU Lesser General Public 
  License along with this library; if not, write to the Free Software 
  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301 USA
   */
 
#ifndef _OPENMBV_NURBSDISK_H_
#define _OPENMBV_NURBSDISK_H_
 
#include <openmbvcppinterface/dynamiccoloredbody.h>
#include <hdf5serie/vectorserie.h>
#include <vector>
#include <cassert>
#include <stdexcept>
 
namespace OpenMBV {
 
  class NurbsDisk : public DynamicColoredBody {
    friend class ObjectFactory;
    public:
      void setLocalFrame(bool f) { localFrameStr=(f)?"true":"false"; }
 
      bool getLocalFrame() { return localFrameStr=="true"?true:false; }
 
      void setScaleFactor(const double scale) {
        scaleFactor=scale;
      }
 
      double getScaleFactor() { return scaleFactor; }
 
      void setDrawDegree(int drawDegree_) {
        drawDegree=drawDegree_;
      }
 
      int getDrawDegree() { return drawDegree; }
 
      void setRadii(double Ri_, double Ro_) {
        Ri=Ri_;
        Ro=Ro_;
      }
 
      void setRi(double Ri_) {
        Ri=Ri_;
      }
 
      void setRo(double Ro_) {
        Ro=Ro_;
      }
 
      double getRi() { return Ri; }
      double getRo() { return Ro; }
 
      void setKnotVecAzimuthal(const std::vector<double> &KnotVecAzimuthal_) {
        KnotVecAzimuthal=KnotVecAzimuthal_;
      }
 
      std::vector<double> getKnotVecAzimuthal() { return KnotVecAzimuthal; }
 
      void setKnotVecRadial(const std::vector<double> &KnotVecRadial_) {
        KnotVecRadial=KnotVecRadial_;
      }
 
      std::vector<double> getKnotVecRadial() { return KnotVecRadial; }
 
      void setElementNumberAzimuthal(int ElementNumberAzimuthal_) {
        ElementNumberAzimuthal=ElementNumberAzimuthal_;
      }
 
      int getElementNumberAzimuthal() { return ElementNumberAzimuthal; }
 
      void setElementNumberRadial(int ElementNumberRadial_) {
        ElementNumberRadial=ElementNumberRadial_;
      }
 
      int getElementNumberRadial() { return ElementNumberRadial; }
 
      void setInterpolationDegreeRadial(int InterpolationDegreeRadial_) {
        InterpolationDegreeRadial=InterpolationDegreeRadial_;
      }
 
      int getInterpolationDegreeRadial() { return InterpolationDegreeRadial; }
 
      void setInterpolationDegreeAzimuthal(int InterpolationDegreeAzimuthal_) {
        InterpolationDegreeAzimuthal=InterpolationDegreeAzimuthal_;
      }
 
      int getInterpolationDegreeAzimuthal() { return InterpolationDegreeAzimuthal; }
      void setDiskNormal(float *DiskNormal_) {
        DiskNormal=DiskNormal_;
      }
 
      float* getDiskNormal() { return DiskNormal; }
 
      void setDiskPoint(float *DiskPoint_) {
        DiskPoint=DiskPoint_;
      }
 
      float* getDiskPoint() { return DiskPoint; }
 
      template<typename T>
      void append(const T& row) { 
        if(data==nullptr) throw std::runtime_error("can not append data to an environment object");
        data->append(row);
      }
 
      int getRows() override { return data?data->getRows():0; }
      std::vector<double> getRow(int i) override {
        int NodeDofs = (getElementNumberRadial() + 1) * (getElementNumberAzimuthal() + getInterpolationDegreeAzimuthal());
        return data?data->getRow(i):std::vector<double>(7+3*NodeDofs+3*getElementNumberAzimuthal()*drawDegree*2);
      }
 
      void initializeUsingXML(xercesc::DOMElement *element) override;
 
      xercesc::DOMElement* writeXMLFile(xercesc::DOMNode *parent) override;
    protected:
      NurbsDisk(); 
      ~NurbsDisk() override;
 
      H5::VectorSerie<double>* data{nullptr};
 
      std::string localFrameStr;
 
      double scaleFactor{1};
 
      int drawDegree{1};
 
      double Ri{0.};
      double Ro{0.};
 
      int ElementNumberAzimuthal{0};
      int ElementNumberRadial{0};
 
      int InterpolationDegreeAzimuthal{8};
      int InterpolationDegreeRadial{3};
 
      std::vector<double> KnotVecAzimuthal, KnotVecRadial;
 
      float *DiskNormal{nullptr};
 
      float *DiskPoint{nullptr};
 
      void createHDF5File() override;
      void openHDF5File() override;
  };
 
}
 
#endif /* _OPENMBV_NURBSDISK_H_ */