Skip to content
Snippets Groups Projects
Select Git revision
  • 5ee346e4fdb5605ab3f561fc8dc7b105774df274
  • master default protected
  • nlf-fixes
  • DSMH
  • OneStep2
  • SMC
  • online-filter-as-a-sampler
  • nonlinear-filter-fixes
  • SMCsamplers
  • smc-sampler
  • 4.5
  • dynamic-striated
  • occbin
  • exo_steady_state
  • filter_initial_state
  • declare_vars_in_model_block
  • exceptions
  • rmExtraExo
  • julia
  • error_msg_undeclared_model_vars
  • static_aux_vars
  • 4.5.6
  • 4.5.5
  • 4.5.4
  • 4.5.3
  • 4.5.2
  • 4.5.1
  • 4.5.0
  • 4.4.3
  • 4.4.2
  • 4.4.1
  • 4.4.0
  • 4.4-beta1
  • 4.3.3
  • 4.3.2
  • 4.3.1
  • 4.3.0
  • 4.2.5
  • 4.2.4
  • 4.2.3
  • 4.2.2
41 results

NumericalInitialization.cc

Blame
  • Forked from Dynare / dynare
    Source project has a limited visibility.
    NumericalInitialization.hh 6.76 KiB
    /*
     * Copyright (C) 2003-2017 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 _NUMERICALINITIALIZATION_HH
    #define _NUMERICALINITIALIZATION_HH
    
    using namespace std;
    
    #include <string>
    #include <vector>
    #include <map>
    
    #include "SymbolTable.hh"
    #include "ExprNode.hh"
    #include "Statement.hh"
    
    class InitParamStatement : public Statement
    {
    private:
      const int symb_id;
      const expr_t param_value;
      const SymbolTable &symbol_table;
    public:
      InitParamStatement(int symb_id_arg, const expr_t param_value_arg,
                         const SymbolTable &symbol_table_arg);
      virtual void checkPass(ModFileStructure &mod_file_struct, WarningConsolidation &warnings);
      virtual void writeOutput(ostream &output, const string &basename, bool minimal_workspace) const;
      virtual void writeJuliaOutput(ostream &output, const string &basename);
      virtual void writeCOutput(ostream &output, const string &basename);
      virtual void writeJsonOutput(ostream &output) const;
      //! Fill eval context with parameter value
      void fillEvalContext(eval_context_t &eval_context) const;
    };
    
    class InitOrEndValStatement : public Statement
    {
    public:
      /*!
        We use a vector instead of a map, since the order of declaration matters:
        an initialization can depend on a previously initialized variable inside the block
      */
      typedef vector<pair<int, expr_t> > init_values_t;
    protected:
      const init_values_t init_values;
      const SymbolTable &symbol_table;
      const bool all_values_required;
    public:
      InitOrEndValStatement(const init_values_t &init_values_arg,
                            const SymbolTable &symbol_table_arg,
                            const bool &all_values_required_arg);
      //! Return set of unused variables by type
      set<int> getUninitializedVariables(SymbolType type);
      //! Fill eval context with variables values
      void fillEvalContext(eval_context_t &eval_context) const;
    protected:
      void writeInitValues(ostream &output) const;
      void writeJsonInitValues(ostream &output) const;
    };
    
    class InitValStatement : public InitOrEndValStatement
    {
    public:
      InitValStatement(const init_values_t &init_values_arg,
                       const SymbolTable &symbol_table_arg,
                       const bool &all_values_required_arg);
      virtual void checkPass(ModFileStructure &mod_file_struct, WarningConsolidation &warnings);
      virtual void writeOutput(ostream &output, const string &basename, bool minimal_workspace) const;
      virtual void writeJsonOutput(ostream &output) const;
      //! Writes initializations for oo_.exo_simul and oo_.exo_det_simul
      void writeOutputPostInit(ostream &output) const;
    };
    
    class EndValStatement : public InitOrEndValStatement
    {
    public:
      EndValStatement(const init_values_t &init_values_arg,
                      const SymbolTable &symbol_table_arg,
                      const bool &all_values_required_arg);
      //! Workaround for trac ticket #35
      virtual void checkPass(ModFileStructure &mod_file_struct, WarningConsolidation &warnings);
      virtual void writeOutput(ostream &output, const string &basename, bool minimal_workspace) const;
      virtual void writeJsonOutput(ostream &output) const;
    };
    
    class HistValStatement : public Statement
    {
    public:
      /*!
        Contrary to Initval and Endval, we use a map, since it is impossible to reuse
        a given initialization value in a second initialization inside the block.
        Maps pairs (symbol_id, lag) to expr_t
      */
      typedef map<pair<int, int>, expr_t> hist_values_t;
    private:
      const hist_values_t hist_values;
      const SymbolTable &symbol_table;
      const bool all_values_required;
    public:
      HistValStatement(const hist_values_t &hist_values_arg,
                       const SymbolTable &symbol_table_arg,
                       const bool &all_values_required_arg);
      //! Workaround for trac ticket #157
      virtual void checkPass(ModFileStructure &mod_file_struct, WarningConsolidation &warnings);
      virtual void writeOutput(ostream &output, const string &basename, bool minimal_workspace) const;
      virtual void writeJsonOutput(ostream &output) const;
    };
    
    class InitvalFileStatement : public Statement
    {
    private:
      const string filename;
    public:
      InitvalFileStatement(const string &filename_arg);
      virtual void writeOutput(ostream &output, const string &basename, bool minimal_workspace) const;
      virtual void writeJsonOutput(ostream &output) const;
    };
    
    class HistvalFileStatement : public Statement
    {
    private:
      const string filename;
    public:
      HistvalFileStatement(const string &filename_arg);
      virtual void writeOutput(ostream &output, const string &basename, bool minimal_workspace) const;
      virtual void writeJsonOutput(ostream &output) const;
    };
    
    class HomotopyStatement : public Statement
    {
    public:
      //! Stores the declarations of homotopy_setup
      /*! Order matter so we use a vector. First expr_t can be NULL if no initial value given. */
      typedef vector<pair<int, pair<expr_t, expr_t> > > homotopy_values_t;
    private:
      const homotopy_values_t homotopy_values;
      const SymbolTable &symbol_table;
    public:
      HomotopyStatement(const homotopy_values_t &homotopy_values_arg,
                        const SymbolTable &symbol_table_arg);
      virtual void writeOutput(ostream &output, const string &basename, bool minimal_workspace) const;
      virtual void writeJsonOutput(ostream &output) const;
    };
    
    class SaveParamsAndSteadyStateStatement : public Statement
    {
    private:
      const string filename;
    public:
      SaveParamsAndSteadyStateStatement(const string &filename_arg);
      virtual void writeOutput(ostream &output, const string &basename, bool minimal_workspace) const;
      virtual void writeJsonOutput(ostream &output) const;
    };
    
    class LoadParamsAndSteadyStateStatement : public Statement
    {
    private:
      const SymbolTable &symbol_table;
      //! Content of the file
      /*! Maps symbol ID to numeric value (stored as string) */
      map<int, string> content;
    public:
      LoadParamsAndSteadyStateStatement(const string &filename,
                                        const SymbolTable &symbol_table_arg,
                                        WarningConsolidation &warnings);
      virtual void writeOutput(ostream &output, const string &basename, bool minimal_workspace) const;
      //! Fill eval context with parameters/variables values
      void fillEvalContext(eval_context_t &eval_context) const;
      virtual void writeJsonOutput(ostream &output) const;
    };
    
    #endif