discretization_interface.h

/* Copyright (C) 2004-2009 MBSim Development Team
 *
 * 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
 *
 * Contact: martin.o.foerg@googlemail.com
 */
 
#ifndef _DISCRETIZATION_INTERFACE_H_
#define _DISCRETIZATION_INTERFACE_H_
 
#include <fmatvec/fmatvec.h>
 
namespace MBSimFlexibleBody {
 
  class DiscretizationInterface {
    public:
      DiscretizationInterface() = default;    
 
      virtual ~DiscretizationInterface() = default;    
 
      virtual const fmatvec::SymMat& getM() const = 0;    
 
      virtual const fmatvec::Vec& geth() const = 0;
 
      virtual const fmatvec::SqrMat& getdhdq() const = 0;    
 
      virtual const fmatvec::SqrMat& getdhdu() const = 0;
 
      virtual int getqSize() const = 0;
 
      virtual int getuSize() const = 0;
 
      virtual void computeM(const fmatvec::Vec& q) = 0;
 
      virtual void computeh(const fmatvec::Vec& q,const fmatvec::Vec& u) = 0;
 
      virtual void computedhdz(const fmatvec::Vec& q,const fmatvec::Vec& u) = 0;
 
      virtual double computeKineticEnergy(const fmatvec::Vec& q,const fmatvec::Vec& u) = 0;
 
      virtual double computeGravitationalEnergy(const fmatvec::Vec& q) = 0;    
 
      virtual double computeElasticEnergy(const fmatvec::Vec& q) = 0;
 
//      /*!
//       * \brief compute position of contour in physical representation
//       * \param generalised positions
//       * \param contour location
//       */
//      virtual fmatvec::Vec computePosition(const fmatvec::Vec& q, const fmatvec::Vec2 &zeta) = 0;
//
//      /*!
//       * \brief compute orientation of contour in physical representation
//       * \param generalised coordiantes
//       * \param contour location
//       */
//      virtual fmatvec::SqrMat computeOrientation(const fmatvec::Vec& q, const fmatvec::Vec2 &zeta) = 0;
//
//      /*!
//       * \brief compute translational velocity of contour in physical representation
//       * \param generalised positions
//       * \param generalised velocities
//       * \param contour location
//       */
//      virtual fmatvec::Vec computeVelocity(const fmatvec::Vec& q, const fmatvec::Vec& u, const fmatvec::Vec2 &zeta) = 0;
//
//      /*!
//       * \brief compute angular velocity of contour in physical representation
//       * \param generalised positions
//       * \param generalised velocities
//       * \param contour location
//       */
//      virtual fmatvec::Vec computeAngularVelocity(const fmatvec::Vec& q, const fmatvec::Vec& u, const fmatvec::Vec2 &zeta) = 0;
//
//      /*!
//       * \brief compute Jacobian of minimal representation regarding physical representation
//       * \param generalised positions
//       * \param contour location
//       */
//      virtual fmatvec::Mat computeJacobianOfMotion(const fmatvec::Vec& q, const fmatvec::Vec2 &zeta) = 0;
  };
 
}
 
#endif /* _DISCRETIZATION_INTERFACE_H_ */