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

- bug fix in Simulate_NG 

git-svn-id: https://www.dynare.org/svn/dynare/trunk@2918 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
ferhat 2009-09-11 10:19:39 +00:00
parent 992370edff
commit e627f62441
1 changed files with 2 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
mex/sources/bytecode/SparseMatrix.cc
View File

@ -1335,11 +1335,11 @@ SparseMatrix::simulate_NG(int blck, int y_size, int it_, int y_kmin, int y_kmax,
int nb_var_piva = At_Row(pivj, &first_piva); int nb_var_piva = At_Row(pivj, &first_piva);
int nb_eq_todo = 0; int nb_eq_todo = 0;
for (j = 0; j < nb_eq && first; j++) for (j = 0; j < nb_eq && first; j++)
{ {
if (!line_done[first->r_index]) if (!line_done[first->r_index])
bc[nb_eq_todo++] = first; bc[nb_eq_todo++] = first;
first = first->NZE_C_N; first = first->NZE_C_N;
} }
//#pragma omp parallel for num_threads(atoi(getenv("DYNARE_NUM_THREADS"))) //#pragma omp parallel for num_threads(atoi(getenv("DYNARE_NUM_THREADS")))
for (j = 0; j < nb_eq_todo; j++) for (j = 0; j < nb_eq_todo; j++)
{ {
@ -1418,7 +1418,6 @@ SparseMatrix::simulate_NG(int blck, int y_size, int it_, int y_kmin, int y_kmax,
} }
} }
u[b[row]] -= u[b[pivj]]*first_elem; u[b[row]] -= u[b[pivj]]*first_elem;
first = first->NZE_C_N;
} }
} }
double slowc_lbx = slowc, res1bx; double slowc_lbx = slowc, res1bx;