GNU Octave  4.2.1
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
op-pm-pm.cc
Go to the documentation of this file.
1 /*
2 
3 Copyright (C) 2008-2017 Jaroslav Hajek
4 
5 This file is part of Octave.
6 
7 Octave is free software; you can redistribute it and/or modify it
8 under the terms of the GNU General Public License as published by the
9 Free Software Foundation; either version 3 of the License, or (at your
10 option) any later version.
11 
12 Octave is distributed in the hope that it will be useful, but WITHOUT
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
15 for more details.
16 
17 You should have received a copy of the GNU General Public License
18 along with Octave; see the file COPYING. If not, see
19 <http://www.gnu.org/licenses/>.
20 
21 */
22 
23 #if defined (HAVE_CONFIG_H)
24 # include "config.h"
25 #endif
26 
27 #include "errwarn.h"
28 #include "ovl.h"
29 #include "ov.h"
30 #include "ov-perm.h"
31 #include "ov-re-mat.h"
32 #include "ov-scalar.h"
33 #include "ov-typeinfo.h"
34 #include "ops.h"
35 #include "xpow.h"
36 
37 DEFUNOP (transpose, perm_matrix)
38 {
39  const octave_perm_matrix& v = dynamic_cast<const octave_perm_matrix&> (a);
40  return octave_value (v.perm_matrix_value ().transpose ());
41 }
42 
43 DEFBINOP_OP (mul, perm_matrix, perm_matrix, *)
44 
45 DEFBINOP (div, perm_matrix, perm_matrix)
46 {
47  const octave_perm_matrix& v1 = dynamic_cast<const octave_perm_matrix&> (a1);
48  const octave_perm_matrix& v2 = dynamic_cast<const octave_perm_matrix&> (a2);
49 
50  return (v1.perm_matrix_value () * v2.perm_matrix_value ().inverse ());
51 }
52 
53 DEFBINOP (ldiv, perm_matrix, perm_matrix)
54 {
55  const octave_perm_matrix& v1 = dynamic_cast<const octave_perm_matrix&> (a1);
56  const octave_perm_matrix& v2 = dynamic_cast<const octave_perm_matrix&> (a2);
57 
58  return (v1.perm_matrix_value ().inverse () * v2.perm_matrix_value ());
59 }
60 
61 DEFBINOP (pow, perm_matrix, scalar)
62 {
63  const octave_perm_matrix& v1 = dynamic_cast<const octave_perm_matrix&> (a1);
64  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
65 
66  return xpow (v1.perm_matrix_value (), v2.scalar_value ());
67 }
68 
69 CONVDECL (perm_matrix_to_matrix)
70 {
71  const octave_perm_matrix& v = dynamic_cast<const octave_perm_matrix&> (a);
72 
73  return new octave_matrix (v.matrix_value ());
74 }
75 
76 void
78 {
81 
86 
88  INSTALL_WIDENOP (octave_perm_matrix, octave_matrix, perm_matrix_to_matrix);
89 }
octave_value xpow(const SparseMatrix &a, double b)
Definition: sparse-xpow.cc:58
#define DEFBINOP(name, t1, t2)
Definition: ops.h:223
#define DEFUNOP(name, t)
Definition: ops.h:172
nd group nd example oindent but is performed more efficiently If only and it is a scalar
Definition: data.cc:5342
double scalar_value(bool=false) const
Definition: ov-scalar.h:150
static void transpose(octave_idx_type N, const octave_idx_type *ridx, const octave_idx_type *cidx, octave_idx_type *ridx2, octave_idx_type *cidx2)
Definition: symrcm.cc:382
PermMatrix transpose(void) const
Definition: PermMatrix.cc:110
#define INSTALL_WIDENOP(t1, t2, f)
Definition: ops.h:70
octave_value op_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1519
PermMatrix perm_matrix_value(void) const
Definition: ov-perm.h:138
Matrix matrix_value(bool=false) const
Definition: ov-perm.cc:227
#define DEFBINOP_OP(name, t1, t2, op)
Definition: ops.h:228
const octave_base_value & a2
calling an anonymous function involves an overhead quite comparable to the overhead of an m file function Passing a handle to a built in function is because the interpreter is not involved in the internal loop For a
Definition: cellfun.cc:398
octave_value op_div(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1517
void install_pm_pm_ops(void)
Definition: op-pm-pm.cc:77
#define INSTALL_BINOP(op, t1, t2, f)
Definition: ops.h:48
PermMatrix inverse(void) const
Definition: PermMatrix.cc:123
octave_value op_transpose(const octave_value &a)
Definition: ov.h:1488
octave_int< T > pow(const octave_int< T > &a, const octave_int< T > &b)
const octave_char_matrix & v2
#define INSTALL_UNOP(op, t, f)
Definition: ops.h:40
octave_value op_ldiv(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1520
octave_value op_hermitian(const octave_value &a)
Definition: ov.h:1489
#define INSTALL_ASSIGNCONV(t1, t2, tr)
Definition: ops.h:66
octave_value op_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1516
#define CONVDECL(name)
Definition: ops.h:161
return octave_value(v1.char_array_value().concat(v2.char_array_value(), ra_idx),((a1.is_sq_string()||a2.is_sq_string())? '\'': '"'))