old_matrix.h

//  Percy++ Copyright 2014 Ian Goldberg <iang@cs.uwaterloo.ca>,
//  Wouter Lueks <lueks@cs.ru.nl>
//
//  This program is free software; you can redistribute it and/or modify
//  it under the terms of version 2 of the GNU General Public License as
//  published by the Free Software Foundation.
//
//  This program is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.
//
//  There is a copy of the GNU General Public License in the COPYING file
//  packaged with this plugin; if you cannot find it, write to the Free
//  Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
//  02110-1301 USA

#ifndef __MATRIX_H__
#define __MATRIX_H__

#include "gf2e.h"
#include <iostream>
#include <stdint.h>
#include <string.h>
#include "xor.h"

template <typename Element> class SubMatrix;
template <typename Element> class Col;
template <typename Element> class Row;
template <typename Element> class Elem;

template <typename Element>
class Matrix {
public:
    Element *data;
    dbsize_t num_rows;
    dbsize_t num_cols;

    // Constructor and destructor
    Matrix ();
    Matrix (Element *data, dbsize_t num_rows, dbsize_t num_cols);

    void print();

    void strassen_submatrices(SubMatrix<Element> &m11,
        SubMatrix<Element> &m12,
        SubMatrix<Element> &m21,
        SubMatrix<Element> &m22);

    void clear();

    void is_sum_of(SubMatrix<const Element> &m1,
        SubMatrix<const Element> &m2);
    void copy_from(SubMatrix<const Element> &src);
};

template <typename Element>
class SubMatrix {
public:
    Matrix<Element> matrix;
    dbsize_t row;
    dbsize_t col;
    dbsize_t num_rows;
    dbsize_t num_cols;

    // Constructors
    SubMatrix ();
    SubMatrix (Matrix<Element> matrix,
        dbsize_t row, dbsize_t col,
        dbsize_t num_rows, dbsize_t num_cols);

    Element* get(dbsize_t row, dbsize_t col);
    Element* first();

    void add(Matrix<Element> &m);

    void print();

    inline void is_mult_of(
    Col<const Element> col_a,
    Row<const Element> row_b);
};

template <typename Element>
class Row : public virtual SubMatrix<Element> {
public:
    Row (Matrix<Element> matrix,
        dbsize_t row, dbsize_t col, dbsize_t num_cols);

    Element* get(dbsize_t row);


    inline void is_mult_of(Row<const Element> row_a,
    SubMatrix<const Element> mat_b);
};

template <typename Element>
class Col : public virtual SubMatrix<Element> {
public:
    Col (Matrix<Element> matrix,
        dbsize_t row, dbsize_t col, dbsize_t num_rows);

    Element* get(dbsize_t col);

    inline void is_mult_of(
    SubMatrix<const Element> mat_a,
    Col<const Element> col_b);

    inline void is_mult_of(
    Col<const Element> col_a,
    Elem<const Element> elem_b);
};

template <typename Element>
class Elem : public Row<Element>, public Col<Element> {
public:
    Elem (Matrix<Element> matrix,
        dbsize_t row, dbsize_t col);

    Element* get(dbsize_t nothing);

    inline void is_mult_of(
    Row<const Element> row_a,
    Col<const Element> col_b);

    inline void is_mult_of(
    Elem<const Element> elem_a,
    Elem<const Element> elem_b);
};

#include "matrix_impl.h"

#endif