FroudeKrylov.h

#ifndef CFroudeKrylov_h__
#define CFroudeKrylov_h__
 
#include <complex>
#include <map>
#include <string>
#include <vector>
 
using std::string;
using std::vector;
using std::map;
using std::complex;
 
using Eigen::Matrix;
 
namespace Ship{
    class FroudeKrylovForce
    {
    public:
        virtual ~FroudeKrylovForce() = default;
 
        virtual void Read( std::string fname )=0;
        
        virtual std::complex<double> GetRAOValue(unsigned int DOF,  double ShipDirection, double ShipSpeed, double WaveFrequency ) = 0;
        //virtual void GetRAOValues( unsigned int DOF, double ShipDirection, double ShipSpeed, Eigen::Matrix< std::complex<double>,Eigen::Dynamic,1 > &outRAO)=0;
    
        //virtual void GetRAOValue2(unsigned int DOF,  double ShipDirection, double ShipSpeed, double WaveFrequency, double &Magnitude, double &Phase ) = 0;
        //virtual void GetRAOValues2( unsigned int DOF, double ShipDirection, double ShipSpeed, Eigen::Matrix< double,Eigen::Dynamic,1 > &outRAOmag,Eigen::Matrix< double,Eigen::Dynamic,1 > &outRAOphase)=0;
    
    
    protected:
        int m_NumDirections;    
        int m_NumSpeed;         
        int m_NumFreq;          
        int m_NumDOF;           
    
        std::vector<double> Direction;  
        std::vector<double> Speed;      
        std::vector<double> Freq;       
        std::vector<std::vector<std::vector<Eigen::Matrix< std::complex<double>, Eigen::Dynamic,1 > > > > RAO;  
        std::vector<std::vector<std::vector<Eigen::Matrix< double, Eigen::Dynamic,1 > > > > RAOmagnitude;
        std::vector<std::vector<std::vector<Eigen::Matrix< double, Eigen::Dynamic,1 > > > > RAOphase;
    };
}
#endif // CFroudeKrylov_h__