/***************************************************************************
                          generalmatrix.h  -  A matrix class to solve matrix equations
                             -------------------
    begin                : do feb 5 2004
    copyright            : (C) 2004 by CJP
    email                : cornware-cjp@users.sourceforge.net
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/

#ifndef GENERALMATRIX_H
#define GENERALMATRIX_H


/**
  *@author CJP
  */

#include "generalvector.h"

#define GM_GAUSS		1
#define GM_PIVOT		2
#define GM_SCALED		4
//#define GM_MODIFIED		8

#define GM_INEQUAL		16

#define GM_LU			32
#define GM_CHOLESKI		64
#define GM_CROUT		128

class CGeneralMatrix {
public:
	CGeneralMatrix(unsigned int rows, unsigned int cols);
	~CGeneralMatrix();

	int solve(CGeneralVector *b, unsigned int solvemethod=GM_GAUSS);
	void setElement(unsigned int i, unsigned int j, float value);
	float getElement(unsigned int i, unsigned int j) const;
	void removeRowCol(unsigned int n);

	CGeneralVector operator*(const CGeneralVector &val) const;
	CGeneralMatrix operator*(const CGeneralMatrix &val) const;
	CGeneralMatrix transpose();

	unsigned int getNumRows() const {return m_numRows;}
	unsigned int getNumCols() const {return m_numCols;}

	bool m_debug;
protected:
	vector<CGeneralVector> m_Matrix;

	//saved for matrix equation solving:
	CGeneralVector *m_Vector;

	unsigned int m_numRows, m_numCols;

	int solveGauss();

	//int solveModified();
	//int solveInequal();
	bool satisfiesInequalities(const CGeneralVector &m) const;
	void restrictPositive(CGeneralVector &v, float desiredSize=-1.0) const;

	void swaprow(unsigned int i, unsigned int j);
	void addrow(unsigned int i, float mul, unsigned int j);

	bool isSmall(float val);
};

#endif