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

sparse_hessian_times_B_kronecker_C MEX: bugfix relating to number of rows of B and C 

In practice however, it seems that this bug never hits, because B and C always
have the same number of rows when the MEX is called.
time-shift
Sébastien Villemot 2021-06-08 17:58:40 +02:00
parent 3bd3c78e0e
commit a40a9b079c
No known key found for this signature in database
GPG Key ID: 2CECE9350ECEBE4A
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
mex/sources/kronecker/sparse_hessian_times_B_kronecker_C.cc
View File

@ -1,5 +1,5 @@
/*
* Copyright © 2007-2020 Dynare Team
* Copyright © 2007-2021 Dynare Team
*
* This file is part of Dynare.
*
@ -114,8 +114,8 @@ sparse_hessian_times_B_kronecker_C(const mwIndex *isparseA, const mwIndex *jspar
if (k1 < k2) // otherwise column ii of A does not have non zero elements (and there is nothing to compute).
{
++nz_in_column_ii_of_A;
mwIndex iC = ii % mB;
mwIndex iB = ii / mB;
mwIndex iC = ii % mC;
mwIndex iB = ii / mC;
double cb = C[jC*mC+iC]*B[jB*mB+iB];
/*
** Loop over the non zero entries of A(:,ii).