Mercurial > octave
view libinterp/corefcn/xnorm.cc @ 25054:6652d3823428 stable
maint: Update copyright dates in all source files.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Fri, 30 Mar 2018 09:19:05 -0400 |
parents | 194eb4bd202b |
children | 00f796120a6d |
line wrap: on
line source
/* Copyright (C) 2008-2018 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 <https://www.gnu.org/licenses/>. */ // author: Jaroslav Hajek <highegg@gmail.com> #if defined (HAVE_CONFIG_H) # include "config.h" #endif #include "oct-norm.h" #include "errwarn.h" #include "ov.h" #include "xnorm.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.iscomplex (); bool issparse = x.issparse (); bool isfloat = x.is_single_type (); if (! isfloat && ! x.is_double_type ()) err_wrong_type_arg ("xnorm", x); if (x.isempty ()) retval = octave_value (0); else 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 ()); } return retval; } octave_value xcolnorms (const octave_value& x, const octave_value& p) { octave_value retval; bool iscomplex = x.iscomplex (); bool issparse = x.issparse (); bool isfloat = x.is_single_type (); if (! isfloat && ! x.is_double_type ()) err_wrong_type_arg ("xcolnorms", x); 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 ()); } return retval; } octave_value xrownorms (const octave_value& x, const octave_value& p) { octave_value retval; bool iscomplex = x.iscomplex (); bool issparse = x.issparse (); bool isfloat = x.is_single_type (); if (! isfloat && ! x.is_double_type ()) err_wrong_type_arg ("xrownorms", x); 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 ()); } return retval; } octave_value xfrobnorm (const octave_value& x) { octave_value retval; bool iscomplex = x.iscomplex (); bool issparse = x.issparse (); bool isfloat = x.is_single_type (); if (! isfloat && ! x.is_double_type ()) err_wrong_type_arg ("xfrobnorm", x); 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 ()); } return retval; }