Mercurial > hg > octave-nkf

comparison liboctave/LSODE.cc @ 3954:8194e0b10a9a

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2002-05-23 18:40:13 by jwe]
author jwe
date Thu, 23 May 2002 18:40:13 +0000
parents 3a8de88328af
children dc970f5941ee
comparison
equal deleted inserted replaced
3953:3a8de88328af 3954:8194e0b10a9a
64 integration_error = 0; 64 integration_error = 0;
65 restart = 1; 65 restart = 1;
66 66
67 istate = 1; 67 istate = 1;
68 itask = 1; 68 itask = 1;
69 iopt = 0;
70 69
71 liw = 20 + n; 70 liw = 20 + n;
72 lrw = 22 + n * (9 + n); 71 lrw = 22 + n * (9 + n);
73 72
74 sanity_checked = 0; 73 sanity_checked = 0;
85 integration_error = 0; 84 integration_error = 0;
86 restart = 1; 85 restart = 1;
87 86
88 istate = 1; 87 istate = 1;
89 itask = 1; 88 itask = 1;
90 iopt = 0;
91 89
92 liw = 20 + n; 90 liw = 20 + n;
93 lrw = 22 + n * (9 + n); 91 lrw = 22 + n * (9 + n);
94 92
95 sanity_checked = 0; 93 sanity_checked = 0;
181 179
182 for (int i = 4; i < 9; i++) 180 for (int i = 4; i < 9; i++)
183 rwork.elem (i) = 0; 181 rwork.elem (i) = 0;
184 } 182 }
185 183
186 if (jac) 184 int method_flag = jac ? 21 : 22;
187 method_flag = 21; 185 int iopt = 0;
188 else
189 method_flag = 22;
190 186
191 integration_error = 0; 187 integration_error = 0;
192 188
193 double *xp = x.fortran_vec (); 189 double *xp = x.fortran_vec ();
194 190