Mercurial > octave

comparison liboctave/numeric/LSODE.cc @ 31549:ed7b17c7ddf3 stable

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
maint: Strip trailing spaces and add missing EOL to all files. * README, NEWS.6.md, RELEASE_CHECKLIST.md, README.md, besselj.cc, Array.h, LSODE.cc, set.m, audiorecorder.m, play.m, set.m, patch.m: Strip trailing spaces. * command-widget.cc, led-indicator.cc, led-indicator.h: Add missing EOL to files.
author Rik <rik@octave.org>
date Fri, 25 Nov 2022 21:38:22 -0800
parents 18a6c1408626
children 597f3ee61a48
comparison
equal deleted inserted replaced
31548:c8ad083a5802 31549:ed7b17c7ddf3
79 79
80 return 0; 80 return 0;
81 } 81 }
82 82
83 static F77_INT 83 static F77_INT
84 lsode_j (const F77_INT& neq, const double& time, double *, 84 lsode_j (const F77_INT& neq, const double& time, double *,
85 const F77_INT& ml, const F77_INT& mu, 85 const F77_INT& ml, const F77_INT& mu,
86 double *pd, const F77_INT& nrowpd) 86 double *pd, const F77_INT& nrowpd)
87 { 87 {
88 Matrix tmp_jac (neq, neq); 88 Matrix tmp_jac (neq, neq);
89 89
100 else 100 else
101 // upper left ends of subdiagonals in tmp_jac 101 // upper left ends of subdiagonals in tmp_jac
102 for (F77_INT i = 0, j = mu; i <= ml; j == 0 ? i++ : j--) 102 for (F77_INT i = 0, j = mu; i <= ml; j == 0 ? i++ : j--)
103 // walk down the subdiagonal in tmp_jac 103 // walk down the subdiagonal in tmp_jac
104 for (F77_INT k = i, l = j; k < neq && l < neq; k++, l++) 104 for (F77_INT k = i, l = j; k < neq && l < neq; k++, l++)
105 pd[nrowpd * l + k + mu - l] = tmp_jac (k, l); 105 pd[nrowpd * l + k + mu - l] = tmp_jac (k, l);
106 106
107 return 0; 107 return 0;
108 } 108 }
109 109
110 ColumnVector 110 ColumnVector
111 LSODE::do_integrate (double tout) 111 LSODE::do_integrate (double tout)