Mercurial > hg > octave-nkf
view examples/code/@polynomial/subsref.m @ 19631:aee5fea8a03e
compute condition number when solving triangular systems (bug #43875)
* dMatrix.cc (Matrix::solve): Pass true for calc_rcond in calls to
utsolve and ltsolve.
* CMatrix.cc (ComplexMatrix::solve): Likewise.
* fCMatrix.cc (FloatComplexMatrix::solve): Likewise.
* fMatrix.cc (FloatMatrix::solve): Likeiwse.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Fri, 26 Dec 2014 14:35:46 -0500 |
parents | c8240a60dd01 |
children |
line wrap: on
line source
function b = subsref (a, s) if (isempty (s)) error ("polynomial: missing index"); endif switch (s(1).type) case "()" ind = s(1).subs; if (numel (ind) != 1) error ("polynomial: need exactly one index"); else b = polyval (fliplr (a.poly), ind{1}); endif case "{}" ind = s(1).subs; if (numel (ind) != 1) error ("polynomial: need exactly one index"); else if (isnumeric (ind{1})) b = a.poly(ind{1}+1); else b = a.poly(ind{1}); endif endif case "." fld = s.subs; if (strcmp (fld, "poly")) b = a.poly; else error ("@polynomial/subsref: invalid property \"%s\"", fld); endif otherwise error ("invalid subscript type"); endswitch if (numel (s) > 1) b = subsref (b, s(2:end)); endif endfunction