FMatrix.h

#ifndef UTIL_F_MATRIX_H
#define UTIL_F_MATRIX_H

/*
* Util Package - C++ Utilities for Scientific Computation
*
* Copyright 2010 - 2017, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/

#include <util/containers/Matrix.h>

namespace Util
{

   template <typename Data, int M, int N>
   class FMatrix : public Matrix<Data>
   { 
   
      using Matrix<Data>::data_;
      using Matrix<Data>::capacity1_;
      using Matrix<Data>::capacity2_;
   
   public: 
   
      FMatrix();

      FMatrix(const FMatrix<Data, M, N>& other);

      ~FMatrix();

      FMatrix<Data, M, N>& operator = (const FMatrix<Data, M, N>& other);

      template <class Archive>
      void serialize(Archive& ar, const unsigned int version);

   private:

      Data fixedArray_[M*N];
   
   };

   /*
   * Default constructor.
   */
   template <typename Data, int M, int N>
   FMatrix<Data,M,N>::FMatrix() 
     : Matrix<Data>()
   {
      data_ = fixedArray_;
      capacity1_ = M;
      capacity2_ = N;
   }
   
   /*
   * Copy constructor.
   */
   template <typename Data, int M, int N>
   FMatrix<Data,M,N>::FMatrix(const FMatrix<Data, M, N>& other) 
     : Matrix<Data>()
   {
      data_ = fixedArray_;
      capacity1_ = M;
      capacity2_ = N;
      for (int i = 0; i < M*N; ++i) {
         fixedArray_[i] = other.fixedArray_[i];
      }
   }
   
   /*
   * Destructor.
   */
   template <typename Data, int M, int N>
   FMatrix<Data,M,N>::~FMatrix()
   {}

   /*
   * Assignment.
   */
   template <typename Data, int M, int N>
   FMatrix<Data,M,N>& 
   FMatrix<Data,M,N>::operator = (const FMatrix<Data, M, N>& other) 
   {
      //Check for self assignment
      if (this == &other) return *this;

      // Copy elements
      for (int i = 0; i < M*N; ++i) {
         fixedArray_[i] = other.fixedArray_[i];
      }

      return *this;
   }
   
   /*
   * Serialize a FMatrix to/from an Archive.
   */
   template <class Data, int M, int N>
   template <class Archive>
   void FMatrix<Data,M,N>::serialize(Archive& ar, const unsigned int version)
   {
      for (int i = 0; i < M*N; ++i) {
         ar & fixedArray_[i];
      }
   }

} 
#endif