Mercurial > hg > octave-nkf
diff liboctave/LSODE.cc @ 5552:bcc328604953
[project @ 2005-11-30 05:48:59 by jwe]
| author | jwe |
|---|---|
| date | Wed, 30 Nov 2005 05:48:59 +0000 |
| parents | 4c8a2e4e0717 |
| children | 87c9641d938f |
line wrap: on
line diff
--- a/liboctave/LSODE.cc +++ b/liboctave/LSODE.cc @@ -134,8 +134,8 @@ else method_flag = 22; - liw = 20 + n; - lrw = 22 + n * (9 + n); + liw = 20 + n + 10000; + lrw = 22 + n * (9 + n) + 10000; } else { @@ -149,6 +149,16 @@ maxord = maximum_order (); + iwork.resize (liw); + + for (octave_idx_type i = 4; i < 9; i++) + iwork(i) = 0; + + rwork.resize (lrw); + + for (octave_idx_type i = 4; i < 9; i++) + rwork(i) = 0; + if (maxord >= 0) { if (maxord > 0 && maxord <= max_maxord) @@ -165,16 +175,6 @@ } } - iwork.resize (liw); - - for (octave_idx_type i = 4; i < 9; i++) - iwork(i) = 0; - - rwork.resize (lrw); - - for (octave_idx_type i = 4; i < 9; i++) - rwork(i) = 0; - if (stop_time_set) { itask = 4;