Mercurial > octave-libgccjit
view src/OPERATORS/op-fpm-fpm.cc @ 8958:6ccc12cc65ef
implement raising a permutation matrix to integer power
author | Jaroslav Hajek <highegg@gmail.com> |
---|---|
date | Wed, 11 Mar 2009 14:07:24 +0100 |
parents | eb63fbe60fab |
children |
line wrap: on
line source
/* Copyright (C) 2008, 2009 Jaroslav Hajek 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-perm.h" #include "ov-flt-perm.h" #include "ov-flt-re-mat.h" #include "ov-float.h" #include "ov-typeinfo.h" #include "ops.h" #include "xpow.h" DEFUNOP (transpose, float_perm_matrix) { CAST_UNOP_ARG (const octave_float_perm_matrix&); return octave_value (v.perm_matrix_value().transpose (), true); } DEFBINOP (mul, float_perm_matrix, float_perm_matrix) { CAST_BINOP_ARGS (const octave_float_perm_matrix&, const octave_float_perm_matrix&); return octave_value (v1.perm_matrix_value () * v2.perm_matrix_value (), true); } DEFBINOP (div, float_perm_matrix, float_perm_matrix) { CAST_BINOP_ARGS (const octave_float_perm_matrix&, const octave_float_perm_matrix&); return octave_value (v1.perm_matrix_value () * v2.perm_matrix_value ().inverse (), false); } DEFBINOP (ldiv, float_perm_matrix, float_perm_matrix) { CAST_BINOP_ARGS (const octave_float_perm_matrix&, const octave_float_perm_matrix&); return octave_value (v1.perm_matrix_value ().inverse () * v2.perm_matrix_value (), false); } DEFBINOP (pow, float_perm_matrix, float_scalar) { CAST_BINOP_ARGS (const octave_float_perm_matrix&, const octave_float_scalar&); return xpow (v1.perm_matrix_value (), v2.float_scalar_value ()); } CONVDECL (float_perm_matrix_to_float_matrix) { CAST_CONV_ARG (const octave_float_perm_matrix&); return new octave_float_matrix (v.float_matrix_value ()); } CONVDECL (float_perm_matrix_to_perm_matrix) { CAST_CONV_ARG (const octave_float_perm_matrix&); return new octave_perm_matrix (v.perm_matrix_value ()); } void install_fpm_fpm_ops (void) { INSTALL_UNOP (op_transpose, octave_float_perm_matrix, transpose); INSTALL_UNOP (op_hermitian, octave_float_perm_matrix, transpose); INSTALL_BINOP (op_mul, octave_float_perm_matrix, octave_float_perm_matrix, mul); INSTALL_BINOP (op_div, octave_float_perm_matrix, octave_float_perm_matrix, div); INSTALL_BINOP (op_ldiv, octave_float_perm_matrix, octave_float_perm_matrix, ldiv); INSTALL_BINOP (op_pow, octave_float_perm_matrix, octave_float_scalar, pow); INSTALL_CONVOP (octave_float_perm_matrix, octave_float_matrix, float_perm_matrix_to_float_matrix); INSTALL_CONVOP (octave_float_perm_matrix, octave_perm_matrix, float_perm_matrix_to_perm_matrix); INSTALL_ASSIGNCONV (octave_float_perm_matrix, octave_float_matrix, octave_float_matrix); INSTALL_WIDENOP (octave_float_perm_matrix, octave_float_matrix, float_perm_matrix_to_float_matrix); }