Skip to content
Snippets Groups Projects
Select Git revision
  • 85dbb649b78bd14f6cedcede0e06fcb5475aebab
  • 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.hh

Blame
  • NumericalInitialization.cc 19.78 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/>.
     */
    
    #include <iostream>
    #include <fstream>
    #include <sstream>
    #include <cstdlib>
    
    #include "NumericalInitialization.hh"
    
    InitParamStatement::InitParamStatement(int symb_id_arg,
                                           const expr_t param_value_arg,
                                           const SymbolTable &symbol_table_arg) :
      symb_id(symb_id_arg),
      param_value(param_value_arg),
      symbol_table(symbol_table_arg)
    {
    }
    
    void
    InitParamStatement::checkPass(ModFileStructure &mod_file_struct, WarningConsolidation &warnings)
    {
      if (symbol_table.getName(symb_id) == "dsge_prior_weight")
        mod_file_struct.dsge_prior_weight_initialized = true;
    }
    
    void
    InitParamStatement::writeOutput(ostream &output, const string &basename, bool minimal_workspace) const
    {
      int id = symbol_table.getTypeSpecificID(symb_id) + 1;
      output << "M_.params( " << id << " ) = ";
      param_value->writeOutput(output);
      output << ";" << endl;
      if (!minimal_workspace)
        output << symbol_table.getName(symb_id) << " = M_.params( " << id << " );" << endl;
    }
    
    void
    InitParamStatement::writeJuliaOutput(ostream &output, const string &basename)
    {
      int id = symbol_table.getTypeSpecificID(symb_id) + 1;
      output << "model_.params[ " << id << " ] = ";
      param_value->writeOutput(output);
      output << endl;
      // Do we really need this?
      // if (!minimal_workspace)
      //   output << symbol_table.getName(symb_id) << " = model_.params[ " << id << " ]" << endl;
    }
    
    void
    InitParamStatement::writeJsonOutput(ostream &output) const
    {
      deriv_node_temp_terms_t tef_terms;
      output << "{\"statementName\": \"param_init\", \"name\": \"" << symbol_table.getName(symb_id) << "\", " << "\"value\": \"";