DeloneType.h

#ifndef DELONETYPE_H
#define DELONETYPE_H

#include <cstdlib>
#include <cstdio>
#include <iostream>
#include <string>
#include <tuple>
#include <vector>

#include "MatS6.h"
#include "LabeledSellaMatrices.h"
#include "LabeledDeloneTypeMatrices.h"


class MatricesForOneDeloneType {
   typedef std::vector<MatS6> vmat;
public:
   MatricesForOneDeloneType() {}
   MatricesForOneDeloneType(const vmat& prj, const vmat& perp, const vmat& toCanon)
      : m_prjs(prj)
      , m_perps(perp)
      , m_toCanons(toCanon)
   {}

   MatS6 GetPrj( const size_t n ) const { return m_prjs[n]; }
   MatS6 GetPerp( const size_t n ) const { return m_perps[n]; }
   MatS6 GetToCanon( const size_t n ) const { return m_toCanons[n]; }
   std::vector<MatS6> GetPrjs    ( void ) const { return m_prjs; }
   std::vector<MatS6> GetPerps   ( void ) const { return m_perps; }
   std::vector<MatS6> GetToCanons( void ) const { return m_toCanons; }

   size_t size(void) const { return m_prjs.size(); }

protected:
   vmat m_prjs;
   vmat m_perps;
   vmat m_toCanons;

};

class DeloneType {
public:
   friend std::ostream& operator<< (std::ostream&, const DeloneType&);

   DeloneType() {}
   DeloneType(const std::string& deloneType, const std::string& bravaisType, const std::string& character, const std::pair<std::string, MatS6 >&);
   DeloneType(
      const std::string& deloneType,
      const std::string& bravaisType,
      const std::string& character,
      const std::string& toCentered_E3,
      const MatS6& ToCentered,
      const MatricesForOneDeloneType& vm);

   std::tuple<double, S6, MatS6> GetFit(const S6& s6) const;

   std::string GetName(void) const { return m_deloneName; }
   std::string GetBravaisType(void) const { return m_bravaisType; }
   std::string GetCharacter(void) const { return m_character; }
   std::string GetToCenteredString(void) const { return m_toCentered_E3; }
   MatS6 GetToCentered(void) const { return m_toCentered_S6; }
   MatricesForOneDeloneType GetFamilyOfMattrices() const  { return m_matrices; }
   std::vector<MatS6> GetPrjs( void ) const { return m_matrices.GetPrjs( ); }
   std::vector<MatS6> GetPerps( void ) const { return m_matrices.GetPerps( ); }

private:
   std::string m_deloneName;
   std::string m_bravaisType;
   std::string m_character;
   std::string m_toCentered_E3;
   MatS6       m_toCentered_S6;

   MatricesForOneDeloneType m_matrices;


};


#endif  // DELONETYPE_H