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(FMatrix<Data, M, N> const & other);

      ~FMatrix();

      FMatrix<Data, M, N>& operator = (FMatrix<Data, M, N> const & 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(FMatrix<Data, M, N> const & 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 = (FMatrix<Data, M, N> const & other)
   {
      // Check for self assignment
      if (this == &other) return *this;
 
      // Copy all 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