Mercurial > octave

comparison libinterp/corefcn/dassl.cc @ 33395:a258493e726a

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Backed out changeset c714266d9f0d
author Arun Giridhar <arungiridhar@gmail.com>
date Sat, 13 Apr 2024 14:28:48 -0400
parents c714266d9f0d
children
comparison
equal deleted inserted replaced
33394:ee94862503e2 33395:a258493e726a
65 dassl_user_function (const ColumnVector& x, const ColumnVector& xdot, 65 dassl_user_function (const ColumnVector& x, const ColumnVector& xdot,
66 double t, octave_idx_type& ires) 66 double t, octave_idx_type& ires)
67 { 67 {
68 ColumnVector retval; 68 ColumnVector retval;
69 69
70 if (x.numel () != xdot.numel ()) 70 panic_unless (x.numel () == xdot.numel ());
71 error ("dassl_user_function: X and XDOT must have the same number of elements");
72 71
73 octave_value_list args; 72 octave_value_list args;
74 73
75 args(2) = t; 74 args(2) = t;
76 args(1) = xdot; 75 args(1) = xdot;
117 dassl_user_jacobian (const ColumnVector& x, const ColumnVector& xdot, 116 dassl_user_jacobian (const ColumnVector& x, const ColumnVector& xdot,
118 double t, double cj) 117 double t, double cj)
119 { 118 {
120 Matrix retval; 119 Matrix retval;
121 120
122 if (x.numel () != xdot.numel ()) 121 panic_unless (x.numel () == xdot.numel ());
123 error ("dassl_user_jacobian: X and XDOT must have the same number of elements");
124 122
125 octave_value_list args; 123 octave_value_list args;
126 124
127 args(3) = cj; 125 args(3) = cj;
128 args(2) = t; 126 args(2) = t;