FFLASMatrix.h

Go to the documentation of this file.

#pragma once


#include <iostream> 
#include <string> 

#include <vector>
#include <assert.h>
#include <cassert>

#include "modular-balanced.h"
#include "ffpack.h"


using namespace std;

template <class Field >
class FFpackMatrix
{

        public:

                typedef         typename Field::Element         ElementType;
        private:

                int   characteristic_m;
                // zwei Moeglichkeiten: ffpack-Feld intern und extern zu definieren     
                const Field* F; 
                
                int     rows_m; 
                int     cols_m; 
                
                int     blocksPerRow;
        
        
                typename Field::Element * data_m; 
                
                string          name;


        protected:

                        //FFpackMatrix(typename Field::Element *data_m,int _rows, int _cols, int characteristic, string name);

                        FFpackMatrix(typename Field::Element *data_m,int _rows, int _cols, int characteristic, string name);
        public:


                        FFpackMatrix(int _rows,int _cols,  int characteristic, string name="unnamedFFPACKMatrix");
                        
                        
                        //template <class TMatrix>
                        //FFpackMatrix(const TMatrix &mat,  int characteristic);
        
                        //FFpackMatrix(int _rows,int _cols, int characteristic, string name="unnamedFFPACKMatrix");
        
                        template <class TMatrix>
                        FFpackMatrix(const TMatrix &mat, int characteristic, bool transpose=false);
        
                        ~FFpackMatrix();
                
                         void           checkBounds(int row, int col) const;
                         inline void    checkDimensions() const;
                         typename Field::Element        getVal(int row, int col) const;
                         void                           setVal(int row, int col, typename Field::Element elem);
                        string                  getName() const { return name; };

                        inline int              getColNum() const { return cols_m; };
                        inline int              getRowNum() const { return rows_m; };

                         const Field*           getField() const  { return F;   };

                        bool isZero()   const;
                        bool operator==(const FFpackMatrix& mat) const;

                        FFpackMatrix*           getRightKernel(int &Rkdim);

                        unsigned int            getRank();
                        FFpackMatrix*           getLeftInverse();

                        typename Field::Element*                solveLGS(const typename Field::Element* rightHandSide,int & rankRef);
                        //FFpackMatrix*                 invertInPlace();
                        void  printValue(std::ostream & os) const;
};


        #include "FFLASMatrix.cpp"