Commit 3fa14031 authored by Sébastien Villemot's avatar Sébastien Villemot
Browse files

Fix bug in *_{dynamic,static}_resid_g1_g2_g3.m functions

parent 732a3170
......@@ -2185,7 +2185,7 @@ DynamicModel::writeWrapperFunctions(const string &basename, const string &ending
<< " g2 = " << basename + "_g2(T, y, x, params, steady_state, it_, false);" << endl;
else if (ending == "g3")
output << " T = " << basename + "_" + ending + "_tt(T, y, x, params, steady_state, it_);" << endl
<< " [residual, g1, g2] = " << basename + "_resid_g1(T, y, x, params, steady_state, it_, false);" << endl
<< " [residual, g1, g2] = " << basename + "_resid_g1_g2(T, y, x, params, steady_state, it_, false);" << endl
<< " g3 = " << basename + "_g3(T, y, x, params, steady_state, it_, false);" << endl;
output << endl << "end" << endl;
......
......@@ -1209,7 +1209,7 @@ StaticModel::writeWrapperFunctions(const string &basename, const string &ending)
<< " g2 = " << basename + "_g2(T, y, x, params, false);" << endl;
else if (ending == "g3")
output << " T = " << basename + "_" + ending + "_tt(T, y, x, params);" << endl
<< " [residual, g1, g2] = " << basename + "_resid_g1(T, y, x, params, false);" << endl
<< " [residual, g1, g2] = " << basename + "_resid_g1_g2(T, y, x, params, false);" << endl
<< " g3 = " << basename + "_g3(T, y, x, params, false);" << endl;
output << endl << "end" << endl;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment