2008-01-11 14:42:14 +01:00
|
|
|
/*
|
2010-09-16 19:00:48 +02:00
|
|
|
* Copyright (C) 2003-2010 Dynare Team
|
2008-01-11 14:42:14 +01:00
|
|
|
*
|
|
|
|
* This file is part of Dynare.
|
|
|
|
*
|
|
|
|
* Dynare 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 3 of the License, or
|
|
|
|
* (at your option) any later version.
|
|
|
|
*
|
|
|
|
* Dynare 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.
|
|
|
|
*
|
|
|
|
* You should have received a copy of the GNU General Public License
|
|
|
|
* along with Dynare. If not, see <http://www.gnu.org/licenses/>.
|
|
|
|
*/
|
|
|
|
|
2006-11-05 00:31:17 +01:00
|
|
|
#ifndef _SIGMAEINITIALIZATION_HH
|
|
|
|
#define _SIGMAEINITIALIZATION_HH
|
2006-12-12 12:54:30 +01:00
|
|
|
|
2006-11-05 00:31:17 +01:00
|
|
|
#include <string>
|
|
|
|
#include <vector>
|
|
|
|
|
2007-03-09 18:27:46 +01:00
|
|
|
#include "ExprNode.hh"
|
2006-12-12 12:54:30 +01:00
|
|
|
#include "Statement.hh"
|
|
|
|
|
|
|
|
//! Stores a Sigma_e statement
|
|
|
|
class SigmaeStatement : public Statement
|
|
|
|
{
|
|
|
|
public:
|
|
|
|
//! Matrix form (lower or upper triangular) enum
|
2010-09-16 19:00:48 +02:00
|
|
|
enum matrix_form_t
|
2006-11-05 00:31:17 +01:00
|
|
|
{
|
2006-12-12 12:54:30 +01:00
|
|
|
eLower = 0, //!< Lower triangular matrix
|
|
|
|
eUpper = 1 //!< Upper triangular matrix
|
2006-11-05 00:31:17 +01:00
|
|
|
};
|
2006-12-12 12:54:30 +01:00
|
|
|
//! Type of a matrix row
|
2010-09-16 19:18:45 +02:00
|
|
|
typedef vector<expr_t> row_t;
|
2006-12-12 12:54:30 +01:00
|
|
|
//! Type of a complete matrix
|
2010-09-16 19:00:48 +02:00
|
|
|
typedef vector<row_t> matrix_t;
|
2006-12-12 12:54:30 +01:00
|
|
|
|
|
|
|
//! An exception indicating that a matrix is neither upper triangular nor lower triangular
|
|
|
|
class MatrixFormException
|
|
|
|
{
|
|
|
|
};
|
|
|
|
private:
|
|
|
|
//! The matrix
|
2010-09-16 19:00:48 +02:00
|
|
|
const matrix_t matrix;
|
2006-12-12 12:54:30 +01:00
|
|
|
//! Matrix form (lower or upper)
|
2010-09-16 19:00:48 +02:00
|
|
|
const matrix_form_t matrix_form;
|
2006-11-05 00:31:17 +01:00
|
|
|
|
2006-12-12 12:54:30 +01:00
|
|
|
//! Returns the type (upper or lower triangular) of a given matrix
|
|
|
|
/*! Throws an exception if it is neither upper triangular nor lower triangular */
|
2010-09-16 19:00:48 +02:00
|
|
|
static matrix_form_t determineMatrixForm(const matrix_t &matrix) throw (MatrixFormException);
|
2006-11-05 00:31:17 +01:00
|
|
|
|
2009-12-16 18:13:23 +01:00
|
|
|
public:
|
2010-09-16 19:00:48 +02:00
|
|
|
SigmaeStatement(const matrix_t &matrix_arg) throw (MatrixFormException);
|
2006-12-19 00:28:52 +01:00
|
|
|
virtual void writeOutput(ostream &output, const string &basename) const;
|
2006-11-05 00:31:17 +01:00
|
|
|
};
|
2006-12-12 12:54:30 +01:00
|
|
|
|
2006-11-05 00:31:17 +01:00
|
|
|
#endif
|