Mercurial > hg > octave-nkf
view examples/code/paramdemo.cc @ 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 | 2f8500ca91d3 |
line wrap: on
line source
#include <octave/oct.h> DEFUN_DLD (paramdemo, args, nargout, "Parameter Check Demo") { octave_value retval; int nargin = args.length (); if (nargin != 1) print_usage (); else if (nargout != 0) error ("paramdemo: OUTPUT argument required"); else { NDArray m = args(0).array_value (); double min_val = -10.0; double max_val = 10.0; octave_stdout << "Properties of input array:\n"; if (m.any_element_is_negative ()) octave_stdout << " includes negative values\n"; if (m.any_element_is_inf_or_nan ()) octave_stdout << " includes Inf or NaN values\n"; if (m.any_element_not_one_or_zero ()) octave_stdout << " includes other values than 1 and 0\n"; if (m.all_elements_are_int_or_inf_or_nan ()) octave_stdout << " includes only int, Inf or NaN values\n"; if (m.all_integers (min_val, max_val)) octave_stdout << " includes only integers in [-10,10]\n"; } return retval; }