Mercurial > octave-nkf
view libinterp/corefcn/xnorm.cc @ 18518:0bdecd41b2dd stable
correctly size fread result (bug #41648)
* oct-stream.cc (octave_base_stream::read): When reading to EOF, don't
add extra column to the result matrix if the number of elements found
is an exact multiple of the number of rows requested.
Avoid mixed signed/unsigned comparisons.
* io.tst: New tests.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Sat, 22 Feb 2014 13:06:18 -0500 |
parents | d63878346099 |
children | 7bf70dc7f06c |
line wrap: on
line source
/* Copyright (C) 2008-2013 VZLU Prague, a.s. This file is part of Octave. Octave is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 3 of the License, or (at your option) any later version. Octave is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with Octave; see the file COPYING. If not, see <http://www.gnu.org/licenses/>. */ // author: Jaroslav Hajek <highegg@gmail.com> #ifdef HAVE_CONFIG_H #include <config.h> #endif #include <cassert> #include <cfloat> #include <cmath> #include "oct-norm.h" #include "error.h" #include "xnorm.h" #include "ov.h" #include "gripes.h" octave_value xnorm (const octave_value& x, const octave_value& p) { octave_value retval; bool isvector = (x.columns () == 1 || x.rows () == 1); bool iscomplex = x.is_complex_type (); bool issparse = x.is_sparse_type (); bool isfloat = x.is_single_type (); if (isfloat || x.is_double_type ()) { if (isvector) { if (isfloat & iscomplex) retval = xnorm (x.float_complex_column_vector_value (), p.float_value ()); else if (isfloat) retval = xnorm (x.float_column_vector_value (), p.float_value ()); else if (iscomplex) retval = xnorm (x.complex_column_vector_value (), p.double_value ()); else retval = xnorm (x.column_vector_value (), p.double_value ()); } else if (issparse) { if (iscomplex) retval = xnorm (x.sparse_complex_matrix_value (), p.double_value ()); else retval = xnorm (x.sparse_matrix_value (), p.double_value ()); } else { if (isfloat & iscomplex) retval = xnorm (x.float_complex_matrix_value (), p.float_value ()); else if (isfloat) retval = xnorm (x.float_matrix_value (), p.float_value ()); else if (iscomplex) retval = xnorm (x.complex_matrix_value (), p.double_value ()); else retval = xnorm (x.matrix_value (), p.double_value ()); } } else gripe_wrong_type_arg ("xnorm", x, true); return retval; } octave_value xcolnorms (const octave_value& x, const octave_value& p) { octave_value retval; bool iscomplex = x.is_complex_type (); bool issparse = x.is_sparse_type (); bool isfloat = x.is_single_type (); if (isfloat || x.is_double_type ()) { if (issparse) { if (iscomplex) retval = xcolnorms (x.sparse_complex_matrix_value (), p.double_value ()); else retval = xcolnorms (x.sparse_matrix_value (), p.double_value ()); } else { if (isfloat & iscomplex) retval = xcolnorms (x.float_complex_matrix_value (), p.float_value ()); else if (isfloat) retval = xcolnorms (x.float_matrix_value (), p.float_value ()); else if (iscomplex) retval = xcolnorms (x.complex_matrix_value (), p.double_value ()); else retval = xcolnorms (x.matrix_value (), p.double_value ()); } } else gripe_wrong_type_arg ("xcolnorms", x, true); return retval; } octave_value xrownorms (const octave_value& x, const octave_value& p) { octave_value retval; bool iscomplex = x.is_complex_type (); bool issparse = x.is_sparse_type (); bool isfloat = x.is_single_type (); if (isfloat || x.is_double_type ()) { if (issparse) { if (iscomplex) retval = xrownorms (x.sparse_complex_matrix_value (), p.double_value ()); else retval = xrownorms (x.sparse_matrix_value (), p.double_value ()); } else { if (isfloat & iscomplex) retval = xrownorms (x.float_complex_matrix_value (), p.float_value ()); else if (isfloat) retval = xrownorms (x.float_matrix_value (), p.float_value ()); else if (iscomplex) retval = xrownorms (x.complex_matrix_value (), p.double_value ()); else retval = xrownorms (x.matrix_value (), p.double_value ()); } } else gripe_wrong_type_arg ("xrownorms", x, true); return retval; } octave_value xfrobnorm (const octave_value& x) { octave_value retval; bool iscomplex = x.is_complex_type (); bool issparse = x.is_sparse_type (); bool isfloat = x.is_single_type (); if (isfloat || x.is_double_type ()) { if (issparse) { if (iscomplex) retval = xfrobnorm (x.sparse_complex_matrix_value ()); else retval = xfrobnorm (x.sparse_matrix_value ()); } else { if (isfloat & iscomplex) retval = xfrobnorm (x.float_complex_matrix_value ()); else if (isfloat) retval = xfrobnorm (x.float_matrix_value ()); else if (iscomplex) retval = xfrobnorm (x.complex_matrix_value ()); else retval = xfrobnorm (x.matrix_value ()); } } else gripe_wrong_type_arg ("xfrobnorm", x, true); return retval; }