Mercurial > octave-nkf
diff src/OPERATORS/op-sm-m.cc @ 5164:57077d0ddc8e
[project @ 2005-02-25 19:55:24 by jwe]
author | jwe |
---|---|
date | Fri, 25 Feb 2005 19:55:28 +0000 |
parents | |
children | 4c8a2e4e0717 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/src/OPERATORS/op-sm-m.cc Fri Feb 25 19:55:28 2005 +0000 @@ -0,0 +1,151 @@ +/* + +Copyright (C) 2004 David Bateman +Copyright (C) 1998-2004 Andy Adler + +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 2, 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 this program; see the file COPYING. If not, write to the Free +Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. + +*/ + +#ifdef HAVE_CONFIG_H +#include <config.h> +#endif + +#include "gripes.h" +#include "oct-obj.h" +#include "ov.h" +#include "ov-typeinfo.h" +#include "ov-re-mat.h" +#include "ops.h" +#include "xdiv.h" + +#include "sparse-xpow.h" +#include "sparse-xdiv.h" +#include "smx-sm-m.h" +#include "smx-m-sm.h" +#include "ov-re-sparse.h" + +// sparse matrix by matrix ops. + +DEFBINOP_OP (add, sparse_matrix, matrix, +) +DEFBINOP_OP (sub, sparse_matrix, matrix, -) + +DEFBINOP (mul, sparse_matrix, matrix) +{ + CAST_BINOP_ARGS (const octave_sparse_matrix&, const octave_matrix&); + + Matrix tmp (v1.matrix_value ()); + + return octave_value ( tmp * v2.matrix_value()); +} + +DEFBINOP (div, sparse_matrix, matrix) +{ + CAST_BINOP_ARGS (const octave_sparse_matrix&, const octave_matrix&); + + return xdiv (v1.matrix_value (), v2.matrix_value ()); +} + +DEFBINOPX (pow, sparse_matrix, matrix) +{ + error ("can't do A ^ B for A and B both matrices"); + return octave_value (); +} + +DEFBINOP (ldiv, sparse_matrix, matrix) +{ + CAST_BINOP_ARGS (const octave_sparse_matrix&, const octave_matrix&); + + return xleftdiv (v1.sparse_matrix_value (), v2.matrix_value ()); +} + + +DEFBINOP_FN (lt, sparse_matrix, matrix, mx_el_lt) +DEFBINOP_FN (le, sparse_matrix, matrix, mx_el_le) +DEFBINOP_FN (eq, sparse_matrix, matrix, mx_el_eq) +DEFBINOP_FN (ge, sparse_matrix, matrix, mx_el_ge) +DEFBINOP_FN (gt, sparse_matrix, matrix, mx_el_gt) +DEFBINOP_FN (ne, sparse_matrix, matrix, mx_el_ne) + +DEFBINOP_FN (el_mul, sparse_matrix, matrix, product) +DEFBINOP_FN (el_div, sparse_matrix, matrix, quotient) + +DEFBINOP (el_pow, sparse_matrix, matrix) +{ + CAST_BINOP_ARGS (const octave_sparse_matrix&, const octave_matrix&); + + return octave_value (elem_xpow (v1.sparse_matrix_value (), + SparseMatrix (v2.matrix_value ()))); +} + +DEFBINOP (el_ldiv, sparse_matrix, matrix) +{ + CAST_BINOP_ARGS (const octave_sparse_matrix&, const octave_matrix&); + + return octave_value + (quotient (v2.matrix_value (), v1.sparse_matrix_value ())); +} + +DEFBINOP_FN (el_and, sparse_matrix, matrix, mx_el_and) +DEFBINOP_FN (el_or, sparse_matrix, matrix, mx_el_or) + +DEFCATOP (sm_m, sparse_matrix, matrix) +{ + CAST_BINOP_ARGS (octave_sparse_matrix&, const octave_matrix&); + SparseMatrix tmp (v2.matrix_value ()); + return octave_value (v1.sparse_matrix_value (). concat (tmp, ra_idx)); +} + +DEFASSIGNOP (assign, sparse_matrix, matrix) +{ + CAST_BINOP_ARGS (octave_sparse_matrix&, const octave_matrix&); + + SparseMatrix tmp (v2.matrix_value ()); + v1.assign (idx, tmp); + return octave_value (); +} + +void +install_sm_m_ops (void) +{ + INSTALL_BINOP (op_add, octave_sparse_matrix, octave_matrix, add); + INSTALL_BINOP (op_sub, octave_sparse_matrix, octave_matrix, sub); + INSTALL_BINOP (op_mul, octave_sparse_matrix, octave_matrix, mul); + INSTALL_BINOP (op_div, octave_sparse_matrix, octave_matrix, div); + INSTALL_BINOP (op_pow, octave_sparse_matrix, octave_matrix, pow); + INSTALL_BINOP (op_ldiv, octave_sparse_matrix, octave_matrix, ldiv); + INSTALL_BINOP (op_lt, octave_sparse_matrix, octave_matrix, lt); + INSTALL_BINOP (op_le, octave_sparse_matrix, octave_matrix, le); + INSTALL_BINOP (op_eq, octave_sparse_matrix, octave_matrix, eq); + INSTALL_BINOP (op_ge, octave_sparse_matrix, octave_matrix, ge); + INSTALL_BINOP (op_gt, octave_sparse_matrix, octave_matrix, gt); + INSTALL_BINOP (op_ne, octave_sparse_matrix, octave_matrix, ne); + INSTALL_BINOP (op_el_mul, octave_sparse_matrix, octave_matrix, el_mul); + INSTALL_BINOP (op_el_div, octave_sparse_matrix, octave_matrix, el_div); + INSTALL_BINOP (op_el_pow, octave_sparse_matrix, octave_matrix, el_pow); + INSTALL_BINOP (op_el_ldiv, octave_sparse_matrix, octave_matrix, el_ldiv); + INSTALL_BINOP (op_el_and, octave_sparse_matrix, octave_matrix, el_and); + INSTALL_BINOP (op_el_or, octave_sparse_matrix, octave_matrix, el_or); + + INSTALL_CATOP (octave_sparse_matrix, octave_matrix, sm_m); + + INSTALL_ASSIGNOP (op_asn_eq, octave_sparse_matrix, octave_matrix, assign); +} + +/* +;;; Local Variables: *** +;;; mode: C++ *** +;;; End: *** +*/