- Adding missing files: StaticDllModel.cc and StaticDllModel.hh
git-svn-id: https://www.dynare.org/svn/dynare/trunk@2867 ac1d8469-bf42-47a9-8791-bf33cf982152issue#70
parent
4897b602c6
commit
081cf9accf
File diff suppressed because it is too large
Load Diff
|
@ -0,0 +1,192 @@
|
|||
/*
|
||||
* Copyright (C) 2003-2009 Dynare Team
|
||||
*
|
||||
* 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/>.
|
||||
*/
|
||||
|
||||
#ifndef _STATICDLLMODEL_HH
|
||||
#define _STATICDLLMODEL_HH
|
||||
|
||||
using namespace std;
|
||||
|
||||
#include <fstream>
|
||||
|
||||
#include "StaticModel.hh"
|
||||
#include "BlockTriangular.hh"
|
||||
|
||||
//! Stores a static model
|
||||
class StaticDllModel : public ModelTree
|
||||
{
|
||||
public:
|
||||
//! The modes in which StaticDllModel can work
|
||||
enum mode_t
|
||||
{
|
||||
eStandardMode, //!< Standard mode (static file in Matlab)
|
||||
eSparseMode, //!< Sparse mode (static file in Matlab with block decomposition)
|
||||
eDLLMode, //!< DLL mode (static file in C)
|
||||
eSparseDLLMode //!< Sparse DLL mode (static file in C with block decomposition plus a binary file)
|
||||
};
|
||||
private:
|
||||
typedef map<pair<int, int>, int> deriv_id_table_t;
|
||||
//! Maps a pair (symbol_id, lag) to a deriv ID
|
||||
deriv_id_table_t deriv_id_table;
|
||||
//! Maps a deriv ID to a pair (symbol_id, lag)
|
||||
vector<pair<int, int> > inv_deriv_id_table;
|
||||
|
||||
//! Maps a deriv_id to the column index of the static Jacobian
|
||||
/*! Contains only endogenous, exogenous and exogenous deterministic */
|
||||
map<int, int> dyn_jacobian_cols_table;
|
||||
|
||||
//! Maximum lag and lead over all types of variables (positive values)
|
||||
/*! Set by computeDerivID() */
|
||||
int max_lag, max_lead;
|
||||
//! Maximum lag and lead over endogenous variables (positive values)
|
||||
/*! Set by computeDerivID() */
|
||||
int max_endo_lag, max_endo_lead;
|
||||
//! Maximum lag and lead over exogenous variables (positive values)
|
||||
/*! Set by computeDerivID() */
|
||||
int max_exo_lag, max_exo_lead;
|
||||
//! Maximum lag and lead over deterministic exogenous variables (positive values)
|
||||
/*! Set by computeDerivID() */
|
||||
int max_exo_det_lag, max_exo_det_lead;
|
||||
|
||||
//! Number of columns of static jacobian
|
||||
/*! Set by computeDerivID() and computeDynJacobianCols() */
|
||||
int dynJacobianColsNbr;
|
||||
|
||||
//! Derivatives of the jacobian w.r. to parameters
|
||||
/*! First index is equation number, second is endo/exo/exo_det variable, and third is parameter.
|
||||
Only non-null derivatives are stored in the map.
|
||||
Variable and parameter indices are those of the getDerivID() method.
|
||||
*/
|
||||
second_derivatives_type params_derivatives;
|
||||
|
||||
//! Temporary terms for the file containing parameters dervicatives
|
||||
temporary_terms_type params_derivs_temporary_terms;
|
||||
|
||||
typedef map< pair< int, pair< int, int> >, NodeID> first_chain_rule_derivatives_type;
|
||||
first_chain_rule_derivatives_type first_chain_rule_derivatives;
|
||||
|
||||
|
||||
//! Writes static model file (Matlab version)
|
||||
void writeStaticMFile(const string &static_basename) const;
|
||||
//! Writes static model file (C version)
|
||||
/*! \todo add third derivatives handling */
|
||||
void writeStaticCFile(const string &static_basename) const;
|
||||
//! Writes static model file when SparseDLL option is on
|
||||
void writeSparseStaticMFile(const string &static_basename, const string &basename, const int mode) const;
|
||||
//! Writes the Block reordred structure of the model in M output
|
||||
void writeModelEquationsOrdered_M(Model_Block *ModelBlock, const string &static_basename) const;
|
||||
//! Writes the code of the Block reordred structure of the model in virtual machine bytecode
|
||||
void writeModelEquationsCodeOrdered(const string file_name, const Model_Block *ModelBlock, const string bin_basename, map_idx_type map_idx) const;
|
||||
//! Computes jacobian and prepares for equation normalization
|
||||
/*! Using values from initval/endval blocks and parameter initializations:
|
||||
- computes the jacobian for the model w.r. to contemporaneous variables
|
||||
- removes edges of the incidence matrix when derivative w.r. to the corresponding variable is too close to zero (below the cutoff)
|
||||
*/
|
||||
void evaluateJacobian(const eval_context_type &eval_context, jacob_map *j_m, bool dynamic);
|
||||
void BlockLinear(Model_Block *ModelBlock);
|
||||
map_idx_type map_idx;
|
||||
//! Build The incidence matrix form the modeltree
|
||||
void BuildIncidenceMatrix();
|
||||
|
||||
void computeTemporaryTermsOrdered(Model_Block *ModelBlock);
|
||||
//! Write derivative code of an equation w.r. to a variable
|
||||
void compileDerivative(ofstream &code_file, int eq, int symb_id, int lag, map_idx_type &map_idx) const;
|
||||
//! Write chain rule derivative code of an equation w.r. to a variable
|
||||
void compileChainRuleDerivative(ofstream &code_file, int eq, int var, int lag, map_idx_type &map_idx) const;
|
||||
|
||||
virtual int computeDerivID(int symb_id, int lag);
|
||||
//! Get the type corresponding to a derivation ID
|
||||
SymbolType getTypeByDerivID(int deriv_id) const throw (UnknownDerivIDException);
|
||||
//! Get the lag corresponding to a derivation ID
|
||||
int getLagByDerivID(int deriv_id) const throw (UnknownDerivIDException);
|
||||
//! Get the symbol ID corresponding to a derivation ID
|
||||
int getSymbIDByDerivID(int deriv_id) const throw (UnknownDerivIDException);
|
||||
//! Compute the column indices of the static Jacobian
|
||||
void computeStatJacobianCols();
|
||||
//! Computes chain rule derivatives of the Jacobian w.r. to endogenous variables
|
||||
void computeChainRuleJacobian(Model_Block *ModelBlock);
|
||||
//! Computes derivatives of the Jacobian w.r. to parameters
|
||||
void computeParamsDerivatives();
|
||||
//! Computes temporary terms for the file containing parameters derivatives
|
||||
void computeParamsDerivativesTemporaryTerms();
|
||||
//! Collect only the first derivatives
|
||||
map<pair<int, pair<int, int> >, NodeID> collect_first_order_derivatives_endogenous();
|
||||
|
||||
//! Helper for writing the Jacobian elements in MATLAB and C
|
||||
/*! Writes either (i+1,j+1) or [i+j*no_eq] */
|
||||
void jacobianHelper(ostream &output, int eq_nb, int col_nb, ExprNodeOutputType output_type) const;
|
||||
|
||||
//! Helper for writing the sparse Hessian elements in MATLAB and C
|
||||
/*! Writes either (i+1,j+1) or [i+j*NNZDerivatives[1]] */
|
||||
void hessianHelper(ostream &output, int row_nb, int col_nb, ExprNodeOutputType output_type) const;
|
||||
|
||||
//! Write chain rule derivative of a recursive equation w.r. to a variable
|
||||
void writeChainRuleDerivative(ostream &output, int eq, int var, int lag, ExprNodeOutputType output_type, const temporary_terms_type &temporary_terms) const;
|
||||
|
||||
|
||||
public:
|
||||
StaticDllModel(SymbolTable &symbol_table_arg, NumericalConstants &num_constants);
|
||||
//! Mode in which the ModelTree is supposed to work (Matlab, DLL or SparseDLL)
|
||||
mode_t mode;
|
||||
//! Adds a variable node
|
||||
/*! This implementation allows for non-zero lag */
|
||||
virtual NodeID AddVariable(const string &name, int lag = 0);
|
||||
//! Absolute value under which a number is considered to be zero
|
||||
double cutoff;
|
||||
//! The weight of the Markowitz criteria to determine the pivot in the linear solver (simul_NG1 and simul_NG from simulate.cc)
|
||||
double markowitz;
|
||||
//! Compute the minimum feedback set in the static model:
|
||||
/*! 0 : all endogenous variables are considered as feedback variables
|
||||
1 : the variables belonging to a non linear equation are considered as feedback variables
|
||||
2 : the variables belonging to a non normalizable non linear equation are considered as feedback variables
|
||||
default value = 0 */
|
||||
int mfs;
|
||||
//! the file containing the model and the derivatives code
|
||||
ofstream code_file;
|
||||
//! Execute computations (variable sorting + derivation)
|
||||
/*!
|
||||
\param jacobianExo whether derivatives w.r. to exo and exo_det should be in the Jacobian (derivatives w.r. to endo are always computed)
|
||||
\param hessian whether 2nd derivatives w.r. to exo, exo_det and endo should be computed (implies jacobianExo = true)
|
||||
\param thirdDerivatives whether 3rd derivatives w.r. to endo/exo/exo_det should be computed (implies jacobianExo = true)
|
||||
\param paramsDerivatives whether 2nd derivatives w.r. to a pair (endo/exo/exo_det, parameter) should be computed (implies jacobianExo = true)
|
||||
\param eval_context evaluation context for normalization
|
||||
\param no_tmp_terms if true, no temporary terms will be computed in the static files
|
||||
*/
|
||||
void computingPass(const eval_context_type &eval_context, bool no_tmp_terms);
|
||||
//! Writes model initialization and lead/lag incidence matrix to output
|
||||
void writeOutput(ostream &output, const string &basename) const;
|
||||
//! Write statements to be added to the main M-file, after computational tasks
|
||||
void writeOutputPostComputing(ostream &output, const string &basename) const;
|
||||
//! Complete set to block decompose the model
|
||||
BlockTriangular block_triangular;
|
||||
//! Adds informations for simulation in a binary file
|
||||
void Write_Inf_To_Bin_File(const string &static_basename, const string &bin_basename,
|
||||
const int &num, int &u_count_int, bool &file_open) const;
|
||||
//! Writes static model file
|
||||
void writeStaticFile(const string &basename) const;
|
||||
//! Writes file containing parameters derivatives
|
||||
void writeParamsDerivativesFile(const string &basename) const;
|
||||
|
||||
//! Writes LaTeX file with the equations of the static model
|
||||
void writeLatexFile(const string &basename) const;
|
||||
|
||||
virtual int getDerivID(int symb_id, int lag) const throw (UnknownDerivIDException);
|
||||
virtual int getDynJacobianCol(int deriv_id) const throw (UnknownDerivIDException);
|
||||
};
|
||||
|
||||
#endif
|
Loading…
Reference in New Issue