stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

- Correction of a bug in the Jacobian Matrix computation for a non block-decomposed model using bytecode

time-shift
Ferhat Mihoubi 2010-11-20 16:06:07 +01:00 committed by Sébastien Villemot
parent 8d79ed5db9
commit 88cdf0193c
1 changed files with 15 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
preprocessor/StaticModel.cc
View File

@ -501,23 +501,26 @@ StaticModel::writeModelEquationsCode(const string file_name, const string bin_ba
{ {
FLDR_ fldr(i); FLDR_ fldr(i);
fldr.write(code_file, instruction_number); fldr.write(code_file, instruction_number);
for(vector<pair<int, int> >::const_iterator it = derivatives[i].begin(); if (derivatives[i].size())
it != derivatives[i].end(); it++)
{ {
FLDSU_ fldsu(it->second); for(vector<pair<int, int> >::const_iterator it = derivatives[i].begin();
fldsu.write(code_file, instruction_number); it != derivatives[i].end(); it++)
FLDSV_ fldsv(eEndogenous, it->first);
fldsv.write(code_file, instruction_number);
FBINARY_ fbinary(oTimes);
fbinary.write(code_file, instruction_number);
if (it != derivatives[i].begin())
{ {
FBINARY_ fbinary(oPlus); FLDSU_ fldsu(it->second);
fldsu.write(code_file, instruction_number);
FLDSV_ fldsv(eEndogenous, it->first);
fldsv.write(code_file, instruction_number);
FBINARY_ fbinary(oTimes);
fbinary.write(code_file, instruction_number); fbinary.write(code_file, instruction_number);
if (it != derivatives[i].begin())
{
FBINARY_ fbinary(oPlus);
fbinary.write(code_file, instruction_number);
}
} }
FBINARY_ fbinary(oMinus);
fbinary.write(code_file, instruction_number);
} }
FBINARY_ fbinary(oMinus);
fbinary.write(code_file, instruction_number);
FSTPSU_ fstpsu(i); FSTPSU_ fstpsu(i);
fstpsu.write(code_file, instruction_number); fstpsu.write(code_file, instruction_number);
} }