diff -r 756c7a550542 scripts/optimization/fsolve.m --- a/scripts/optimization/fsolve.m Tue Feb 28 16:54:08 2017 -0800 +++ b/scripts/optimization/fsolve.m Mon Mar 06 08:59:54 2017 +0100 @@ -235,6 +235,10 @@ ## Calculate function value and Jacobian (possibly via FD). if (has_jac) [fvec, fjac] = fcn (reshape (x, xsiz)); + if ~all(size(fjac) == [m, n]) + error ("fsolve: jacobian should have size (%d,%d), has (%d,%d)", ... + m, n, size(fjac)); + endif ## If the Jacobian is sparse, disable Broyden updating. if (issparse (fjac)) updating = false;