Mercurial > octave-nkf
view src/OPERATORS/op-fcs-fcm.cc @ 8968:91d53dc37f79
Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations.
Nothing terribly fancy in any of this. There probably is some
mechanism for using the permutation vectors and some assign or index
method in the sparse classes, but I've never understood all the
intricacies. I'm opting for a simple implementation at the cost of
possibly duplicating some functionality.
author | Jason Riedy <jason@acm.org> |
---|---|
date | Tue, 10 Mar 2009 21:54:44 -0400 |
parents | eb63fbe60fab |
children | 319e2ab9b8ae |
line wrap: on
line source
/* Copyright (C) 1996, 1997, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008 John W. Eaton 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/>. */ #ifdef HAVE_CONFIG_H #include <config.h> #endif #include "gripes.h" #include "oct-obj.h" #include "ov.h" #include "ov-complex.h" #include "ov-flt-complex.h" #include "ov-cx-mat.h" #include "ov-flt-cx-mat.h" #include "ov-typeinfo.h" #include "ops.h" #include "xdiv.h" #include "xpow.h" // complex scalar by complex matrix ops. DEFNDBINOP_OP (add, float_complex, float_complex_matrix, float_complex, float_complex_array, +) DEFNDBINOP_OP (sub, float_complex, float_complex_matrix, float_complex, float_complex_array, -) DEFNDBINOP_OP (mul, float_complex, float_complex_matrix, float_complex, float_complex_array, *) DEFBINOP (div, float_complex, float_complex_matrix) { CAST_BINOP_ARGS (const octave_float_complex&, const octave_float_complex_matrix&); FloatComplexMatrix m1 = v1.float_complex_matrix_value (); FloatComplexMatrix m2 = v2.float_complex_matrix_value (); MatrixType typ = v2.matrix_type (); FloatComplexMatrix ret = xdiv (m1, m2, typ); v2.matrix_type (typ); return ret; } DEFBINOP_FN (pow, float_complex, float_complex_matrix, xpow) DEFBINOP (ldiv, float_complex, float_complex_matrix) { CAST_BINOP_ARGS (const octave_float_complex&, const octave_float_complex_matrix&); FloatComplex d = v1.float_complex_value (); if (d == static_cast<float>(0.0)) gripe_divide_by_zero (); return octave_value (v2.float_complex_array_value () / d); } DEFNDBINOP_FN (lt, float_complex, float_complex_matrix, float_complex, float_complex_array, mx_el_lt) DEFNDBINOP_FN (le, float_complex, float_complex_matrix, float_complex, float_complex_array, mx_el_le) DEFNDBINOP_FN (eq, float_complex, float_complex_matrix, float_complex, float_complex_array, mx_el_eq) DEFNDBINOP_FN (ge, float_complex, float_complex_matrix, float_complex, float_complex_array, mx_el_ge) DEFNDBINOP_FN (gt, float_complex, float_complex_matrix, float_complex, float_complex_array, mx_el_gt) DEFNDBINOP_FN (ne, float_complex, float_complex_matrix, float_complex, float_complex_array, mx_el_ne) DEFNDBINOP_OP (el_mul, float_complex, float_complex_matrix, float_complex, float_complex_array, *) DEFNDBINOP_FN (el_div, float_complex, float_complex_matrix, float_complex, float_complex_array, x_el_div) DEFNDBINOP_FN (el_pow, float_complex, float_complex_matrix, float_complex, float_complex_array, elem_xpow) DEFBINOP (el_ldiv, float_complex, float_complex_matrix) { CAST_BINOP_ARGS (const octave_float_complex&, const octave_float_complex_matrix&); FloatComplex d = v1.float_complex_value (); if (d == static_cast<float>(0.0)) gripe_divide_by_zero (); return octave_value (v2.float_complex_array_value () / d); } DEFNDBINOP_FN (el_and, float_complex, float_complex_matrix, float_complex, float_complex_array, mx_el_and) DEFNDBINOP_FN (el_or, float_complex, float_complex_matrix, float_complex, float_complex_array, mx_el_or) DEFNDCATOP_FN (fcs_fcm, float_complex, float_complex_matrix, float_complex_array, float_complex_array, concat) DEFNDCATOP_FN (cs_fcm, complex, float_complex_matrix, float_complex_array, float_complex_array, concat) DEFNDCATOP_FN (fcs_cm, float_complex, complex_matrix, float_complex_array, float_complex_array, concat) DEFCONV (float_complex_matrix_conv, float_complex, float_complex_matrix) { CAST_CONV_ARG (const octave_float_complex&); return new octave_float_complex_matrix (v.float_complex_matrix_value ()); } void install_fcs_fcm_ops (void) { INSTALL_BINOP (op_add, octave_float_complex, octave_float_complex_matrix, add); INSTALL_BINOP (op_sub, octave_float_complex, octave_float_complex_matrix, sub); INSTALL_BINOP (op_mul, octave_float_complex, octave_float_complex_matrix, mul); INSTALL_BINOP (op_div, octave_float_complex, octave_float_complex_matrix, div); INSTALL_BINOP (op_pow, octave_float_complex, octave_float_complex_matrix, pow); INSTALL_BINOP (op_ldiv, octave_float_complex, octave_float_complex_matrix, ldiv); INSTALL_BINOP (op_lt, octave_float_complex, octave_float_complex_matrix, lt); INSTALL_BINOP (op_le, octave_float_complex, octave_float_complex_matrix, le); INSTALL_BINOP (op_eq, octave_float_complex, octave_float_complex_matrix, eq); INSTALL_BINOP (op_ge, octave_float_complex, octave_float_complex_matrix, ge); INSTALL_BINOP (op_gt, octave_float_complex, octave_float_complex_matrix, gt); INSTALL_BINOP (op_ne, octave_float_complex, octave_float_complex_matrix, ne); INSTALL_BINOP (op_el_mul, octave_float_complex, octave_float_complex_matrix, el_mul); INSTALL_BINOP (op_el_div, octave_float_complex, octave_float_complex_matrix, el_div); INSTALL_BINOP (op_el_pow, octave_float_complex, octave_float_complex_matrix, el_pow); INSTALL_BINOP (op_el_ldiv, octave_float_complex, octave_float_complex_matrix, el_ldiv); INSTALL_BINOP (op_el_and, octave_float_complex, octave_float_complex_matrix, el_and); INSTALL_BINOP (op_el_or, octave_float_complex, octave_float_complex_matrix, el_or); INSTALL_CATOP (octave_float_complex, octave_float_complex_matrix, fcs_fcm); INSTALL_CATOP (octave_complex, octave_float_complex_matrix, cs_fcm); INSTALL_CATOP (octave_float_complex, octave_complex_matrix, fcs_cm); INSTALL_ASSIGNCONV (octave_float_complex, octave_float_complex_matrix, octave_float_complex_matrix); INSTALL_ASSIGNCONV (octave_complex, octave_float_complex_matrix, octave_complex_matrix); INSTALL_WIDENOP (octave_float_complex, octave_float_complex_matrix, float_complex_matrix_conv); } /* ;;; Local Variables: *** ;;; mode: C++ *** ;;; End: *** */