Skip to content
Snippets Groups Projects
Select Git revision
  • d9f78ae33649572bed698ed476ef181a40b0705d
  • master default protected
  • julia protected
  • 6.x protected
  • python-codegen
  • llvm-15
  • 5.x protected
  • 4.6 protected
  • uop
  • rework_pac
  • aux_vars_fix
  • julia-7.0.0
  • julia-6.4.0
  • julia-6.3.0
  • julia-6.2.0
15 results

ComputingTasks.cc

Blame
  • DynamicModel.cc 275.81 KiB
    /*
     * Copyright © 2003-2021 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 <https://www.gnu.org/licenses/>.
     */
    
    #include <iostream>
    #include <cmath>
    #include <cstdlib>
    #include <cassert>
    #include <algorithm>
    #include <numeric>
    #include <regex>
    #include <sstream>
    
    #include "DynamicModel.hh"
    
    void
    DynamicModel::copyHelper(const DynamicModel &m)
    {
      auto f = [this](const ExprNode *e) { return e->clone(*this); };
    
      for (const auto &it : m.static_only_equations)
        static_only_equations.push_back(dynamic_cast<BinaryOpNode *>(f(it)));
    
        auto convert_block_derivative = [f](const map<tuple<int, int, int>, expr_t> &dt)
                                        {
                                          map<tuple<int, int, int>, expr_t> dt2;
                                          for (const auto &it : dt)
                                            dt2[it.first] = f(it.second);
                                          return dt2;
                                        };
      for (const auto &it : m.blocks_derivatives_other_endo)
        blocks_derivatives_other_endo.emplace_back(convert_block_derivative(it));
      for (const auto &it : m.blocks_derivatives_exo)
        blocks_derivatives_exo.emplace_back(convert_block_derivative(it));
      for (const auto &it : m.blocks_derivatives_exo_det)
        blocks_derivatives_exo_det.emplace_back(convert_block_derivative(it));
    
      for (const auto &[key, expr] : m.pac_expectation_substitution)
        pac_expectation_substitution.emplace(key, f(expr));
    }
    
    DynamicModel::DynamicModel(SymbolTable &symbol_table_arg,
                               NumericalConstants &num_constants_arg,
                               ExternalFunctionsTable &external_functions_table_arg,
                               TrendComponentModelTable &trend_component_model_table_arg,
                               VarModelTable &var_model_table_arg) :
      ModelTree{symbol_table_arg, num_constants_arg, external_functions_table_arg, true},
      trend_component_model_table{trend_component_model_table_arg},
      var_model_table{var_model_table_arg}
    {
    }
    
    DynamicModel::DynamicModel(const DynamicModel &m) :
      ModelTree{m},
      trend_component_model_table{m.trend_component_model_table},