view libinterp/operators/ops.h @ 30564:796f54d4ddbf stable

update Octave Project Developers copyright for the new year In files that have the "Octave Project Developers" copyright notice, update for 2021. In all .txi and .texi files except gpl.txi and gpl.texi in the doc/liboctave and doc/interpreter directories, change the copyright to "Octave Project Developers", the same as used for other source files. Update copyright notices for 2022 (not done since 2019). For gpl.txi and gpl.texi, change the copyright notice to be "Free Software Foundation, Inc." and leave the date at 2007 only because this file only contains the text of the GPL, not anything created by the Octave Project Developers. Add Paul Thomas to contributors.in.
author John W. Eaton <jwe@octave.org>
date Tue, 28 Dec 2021 18:22:40 -0500
parents b260322f6730
children 08b08b7f05b2
line wrap: on
line source

////////////////////////////////////////////////////////////////////////
//
// Copyright (C) 1996-2022 The Octave Project Developers
//
// See the file COPYRIGHT.md in the top-level directory of this
// distribution or <https://octave.org/copyright/>.
//
// 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/>.
//
////////////////////////////////////////////////////////////////////////

#if ! defined (octave_ops_h)
#define octave_ops_h 1

#include "octave-config.h"

#include "Array-util.h"

namespace octave
{
  class type_info;
}

// Concatenation macros that enforce argument prescan
#define CONCAT2X(x, y) x ## y
#define CONCAT2(x, y) CONCAT2X (x, y)

#define CONCAT3X(x, y, z) x ## y ## z
#define CONCAT3(x, y, z) CONCAT3X (x, y, z)

#define INSTALL_UNOP_TI(ti, op, t, f)                                   \
  ti.install_unary_op                                                   \
  (octave_value::op, t::static_type_id (), CONCAT2 (oct_unop_, f));

#define INSTALL_NCUNOP_TI(ti, op, t, f)                                 \
  ti.install_non_const_unary_op                                         \
  (octave_value::op, t::static_type_id (), CONCAT2 (oct_unop_, f));

#define INSTALL_BINOP_TI(ti, op, t1, t2, f)                             \
  ti.install_binary_op                                                  \
  (octave_value::op, t1::static_type_id (), t2::static_type_id (),      \
   CONCAT2 (oct_binop_, f));

#define INSTALL_CATOP_TI(ti, t1, t2, f)                                 \
  ti.install_cat_op                                                     \
  (t1::static_type_id (), t2::static_type_id (), CONCAT2 (oct_catop_, f));

#define INSTALL_ASSIGNOP_TI(ti, op, t1, t2, f)                          \
  ti.install_assign_op                                                  \
  (octave_value::op, t1::static_type_id (), t2::static_type_id (),      \
   CONCAT2 (oct_assignop_, f));

#define INSTALL_ASSIGNANYOP_TI(ti, op, t1, f)                           \
  ti.install_assignany_op                                               \
  (octave_value::op, t1::static_type_id (), CONCAT2 (oct_assignop_, f));

#define INSTALL_ASSIGNCONV_TI(ti, t1, t2, tr)                           \
  ti.install_pref_assign_conv                                           \
  (t1::static_type_id (), t2::static_type_id (), tr::static_type_id ());

#define INSTALL_WIDENOP_TI(ti, t1, t2, f)                               \
  ti.install_widening_op                                                \
  (t1::static_type_id (), t2::static_type_id (), CONCAT2 (oct_conv_, f));

#define DEFASSIGNOP(name, t1, t2)                               \
  static octave_value                                           \
  CONCAT2 (oct_assignop_, name) (octave_base_value& a1,         \
                                 const octave_value_list& idx,  \
                                 const octave_base_value& a2)

#define DEFASSIGNOP_FN(name, t1, t2, f)                                 \
  static octave_value                                                   \
  CONCAT2 (oct_assignop_, name) (octave_base_value& a1,                 \
                                 const octave_value_list& idx,          \
                                 const octave_base_value& a2)           \
  {                                                                     \
    CONCAT2 (octave_, t1)& v1 = dynamic_cast<CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    v1.f (idx, v2.CONCAT2 (t1, _value) ());                             \
    return octave_value ();                                             \
  }

#define DEFNULLASSIGNOP_FN(name, t, f)                                  \
  static octave_value                                                   \
  CONCAT2 (oct_assignop_, name) (octave_base_value& a,                  \
                                 const octave_value_list& idx,          \
                                 const octave_base_value&)              \
  {                                                                     \
    CONCAT2 (octave_, t)& v = dynamic_cast<CONCAT2 (octave_, t)&> (a);  \
                                                                        \
    v.f (idx);                                                          \
    return octave_value ();                                             \
  }

#define DEFNDASSIGNOP_FN(name, t1, t2, e, f)                            \
  static octave_value                                                   \
  CONCAT2 (oct_assignop_, name) (octave_base_value& a1,                 \
                                 const octave_value_list& idx,          \
                                 const octave_base_value& a2)           \
  {                                                                     \
    CONCAT2 (octave_, t1)& v1 = dynamic_cast<CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    v1.f (idx, v2.CONCAT2 (e, _value) ());                              \
    return octave_value ();                                             \
  }

// FIXME: the following currently don't handle index.
#define DEFNDASSIGNOP_OP(name, t1, t2, f, op)                           \
  static octave_value                                                   \
  CONCAT2 (oct_assignop_, name) (octave_base_value& a1,                 \
                                 const octave_value_list& idx,          \
                                 const octave_base_value& a2)           \
  {                                                                     \
    CONCAT2 (octave_, t1)& v1 = dynamic_cast<CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    assert (idx.empty ());                                              \
    v1.matrix_ref () op v2.CONCAT2 (f, _value) ();                      \
                                                                        \
    return octave_value ();                                             \
  }

#define DEFNDASSIGNOP_FNOP(name, t1, t2, f, fnop)                       \
  static octave_value                                                   \
  CONCAT2 (oct_assignop_, name) (octave_base_value& a1,                 \
                                 const octave_value_list& idx,          \
                                 const octave_base_value& a2)           \
  {                                                                     \
    CONCAT2 (octave_, t1)& v1 = dynamic_cast<CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    assert (idx.empty ());                                              \
    fnop (v1.matrix_ref (), v2.CONCAT2 (f, _value) ());                 \
                                                                        \
    return octave_value ();                                             \
  }

#define DEFASSIGNANYOP_FN(name, t1, f)                                  \
  static octave_value                                                   \
  CONCAT2 (oct_assignop_, name) (octave_base_value& a1,                 \
                                 const octave_value_list& idx,          \
                                 const octave_value& a2)                \
  {                                                                     \
    CONCAT2 (octave_, t1)& v1 = dynamic_cast<CONCAT2 (octave_, t1)&> (a1); \
                                                                        \
    v1.f (idx, a2);                                                     \
    return octave_value ();                                             \
  }

#define CONVDECL(name)                                          \
  static octave_base_value *                                    \
  CONCAT2 (oct_conv_, name) (const octave_base_value& a)

#define DEFCONV(name, a_dummy, b_dummy)         \
  CONVDECL (name)

#define DEFUNOPX(name, t)                               \
  static octave_value                                   \
  CONCAT2 (oct_unop_, name) (const octave_base_value&)

#define DEFUNOP(name, t)                                        \
  static octave_value                                           \
  CONCAT2 (oct_unop_, name) (const octave_base_value& a)

#define DEFUNOP_OP(name, t, op)                                         \
  static octave_value                                                   \
  CONCAT2 (oct_unop_, name) (const octave_base_value& a)                \
  {                                                                     \
    const CONCAT2 (octave_, t)& v = dynamic_cast<const CONCAT2 (octave_, t)&> (a); \
    return octave_value (op v.CONCAT2 (t, _value) ());                  \
  }

#define DEFNDUNOP_OP(name, t, e, op)                                    \
  static octave_value                                                   \
  CONCAT2 (oct_unop_, name) (const octave_base_value& a)                \
  {                                                                     \
    const CONCAT2 (octave_, t)& v = dynamic_cast<const CONCAT2 (octave_, t)&> (a); \
    return octave_value (op v.CONCAT2 (e, _value) ());                  \
  }

// FIXME: in some cases, the constructor isn't necessary.

#define DEFUNOP_FN(name, t, f)                                          \
  static octave_value                                                   \
  CONCAT2 (oct_unop_, name) (const octave_base_value& a)                \
  {                                                                     \
    const CONCAT2 (octave_, t)& v = dynamic_cast<const CONCAT2 (octave_, t)&> (a); \
    return octave_value (f (v.CONCAT2 (t, _value) ()));                 \
  }

#define DEFNDUNOP_FN(name, t, e, f)                                     \
  static octave_value                                                   \
  CONCAT2 (oct_unop_, name) (const octave_base_value& a)                \
  {                                                                     \
    const CONCAT2 (octave_, t)& v = dynamic_cast<const CONCAT2 (octave_, t)&> (a); \
    return octave_value (f (v.CONCAT2 (e, _value) ()));                 \
  }

#define DEFNCUNOP_METHOD(name, t, method)                               \
  static void                                                           \
  CONCAT2 (oct_unop_, name) (octave_base_value& a)                      \
  {                                                                     \
    CONCAT2 (octave_, t)& v = dynamic_cast<CONCAT2 (octave_, t)&> (a);  \
    v.method ();                                                        \
  }

#define DEFBINOPX(name, t1, t2)                         \
  static octave_value                                   \
  CONCAT2 (oct_binop_, name) (const octave_base_value&, \
                              const octave_base_value&)

#define DEFBINOP(name, t1, t2)                                  \
  static octave_value                                           \
  CONCAT2 (oct_binop_, name) (const octave_base_value& a1,      \
                              const octave_base_value& a2)

#define DEFBINOP_OP(name, t1, t2, op)                                   \
  static octave_value                                                   \
  CONCAT2 (oct_binop_, name) (const octave_base_value& a1,              \
                              const octave_base_value& a2)              \
  {                                                                     \
    const CONCAT2 (octave_, t1)& v1 = dynamic_cast<const CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    return octave_value                                                 \
      (v1.CONCAT2 (t1, _value) () op v2.CONCAT2 (t2, _value) ());       \
  }

#define DEFCMPLXCMPOP_OP(name, t1, t2, op)                              \
  static octave_value                                                   \
  CONCAT2 (oct_binop_, name) (const octave_base_value& a1,              \
                              const octave_base_value& a2)              \
  {                                                                     \
    const CONCAT2 (octave_, t1)& v1 = dynamic_cast<const CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    warn_complex_cmp ();                                                \
                                                                        \
    return octave_value                                                 \
      (v1.CONCAT2 (t1, _value) () op v2.CONCAT2 (t2, _value) ());       \
  }

#define DEFSCALARBOOLOP_OP(name, t1, t2, op)                            \
  static octave_value                                                   \
  CONCAT2 (oct_binop_, name) (const octave_base_value& a1,              \
                              const octave_base_value& a2)              \
  {                                                                     \
    const CONCAT2 (octave_, t1)& v1 = dynamic_cast<const CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    if (octave::math::isnan (v1.CONCAT2 (t1, _value) ()) || octave::math::isnan (v2.CONCAT2 (t2, _value) ())) \
      octave::err_nan_to_logical_conversion ();                                 \
                                                                        \
    return octave_value                                                 \
      (v1.CONCAT2 (t1, _value) () op v2.CONCAT2 (t2, _value) ());       \
  }

#define DEFNDBINOP_OP(name, t1, t2, e1, e2, op)                         \
  static octave_value                                                   \
  CONCAT2 (oct_binop_, name) (const octave_base_value& a1,              \
                              const octave_base_value& a2)              \
  {                                                                     \
    const CONCAT2 (octave_, t1)& v1 = dynamic_cast<const CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    return octave_value                                                 \
      (v1.CONCAT2 (e1, _value) () op v2.CONCAT2 (e2, _value) ());       \
  }

// FIXME: in some cases, the constructor isn't necessary.

#define DEFBINOP_FN(name, t1, t2, f)                                    \
  static octave_value                                                   \
  CONCAT2 (oct_binop_, name) (const octave_base_value& a1,              \
                              const octave_base_value& a2)              \
  {                                                                     \
    const CONCAT2 (octave_, t1)& v1 = dynamic_cast<const CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    return octave_value (f (v1.CONCAT2 (t1, _value) (), v2.CONCAT2 (t2, _value) ())); \
  }

#define DEFNDBINOP_FN(name, t1, t2, e1, e2, f)                          \
  static octave_value                                                   \
  CONCAT2 (oct_binop_, name) (const octave_base_value& a1,              \
                              const octave_base_value& a2)              \
  {                                                                     \
    const CONCAT2 (octave_, t1)& v1 = dynamic_cast<const CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    return octave_value (f (v1.CONCAT2 (e1, _value) (), v2.CONCAT2 (e2, _value) ())); \
  }

#define DEFNDCMPLXCMPOP_FN(name, t1, t2, e1, e2, f)                     \
  static octave_value                                                   \
  CONCAT2 (oct_binop_, name) (const octave_base_value& a1,              \
                              const octave_base_value& a2)              \
  {                                                                     \
    const CONCAT2 (octave_, t1)& v1 = dynamic_cast<const CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    return octave_value (f (v1.CONCAT2 (e1, _value) (), v2.CONCAT2 (e2, _value) ())); \
  }

#define DEFCATOPX(name, t1, t2)                                         \
  static octave_value                                                   \
  CONCAT2 (oct_catop_, name) (const octave_base_value&,                 \
                              const octave_base_value&,                 \
                              const Array<octave_idx_type>& ra_idx)

#define DEFCATOP(name, t1, t2)                                          \
  static octave_value                                                   \
  CONCAT2 (oct_catop_, name) (const octave_base_value& a1,              \
                              const octave_base_value& a2,              \
                              const Array<octave_idx_type>& ra_idx)

// FIXME: in some cases, the constructor isn't necessary.

#define DEFCATOP_FN(name, t1, t2, f)                                    \
  static octave_value                                                   \
  CONCAT2 (oct_catop_, name) (const octave_base_value& a1,              \
                              const octave_base_value& a2,              \
                              const Array<octave_idx_type>& ra_idx)     \
  {                                                                     \
    const CONCAT2 (octave_, t1)& v1 = dynamic_cast<const CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    return octave_value (v1.CONCAT2 (t1, _value) () . f (v2.CONCAT2 (t2, _value) (), ra_idx)); \
  }

#define DEFNDCATOP_FN(name, t1, t2, e1, e2, f)                          \
  static octave_value                                                   \
  CONCAT2 (oct_catop_, name) (const octave_base_value& a1,              \
                              const octave_base_value& a2,              \
                              const Array<octave_idx_type>& ra_idx)     \
  {                                                                     \
    const CONCAT2 (octave_, t1)& v1 = dynamic_cast<const CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    return octave_value (v1.CONCAT2 (e1, _value) () . f (v2.CONCAT2 (e2, _value) (), ra_idx)); \
  }

#define DEFNDCHARCATOP_FN(name, t1, t2, f)                              \
  static octave_value                                                   \
  CONCAT2 (oct_catop_, name) (const octave_base_value& a1,              \
                              const octave_base_value& a2,              \
                              const Array<octave_idx_type>& ra_idx)     \
  {                                                                     \
    const CONCAT2 (octave_, t1)& v1 = dynamic_cast<const CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    return octave_value (v1.char_array_value () . f (v2.char_array_value (), ra_idx), \
                         ((a1.is_sq_string () || a2.is_sq_string ())    \
                          ? '\'' : '"'));                               \
  }

// For compatibility, the second arg is always converted to the type
// of the first.  Hmm.

#define DEFNDCATOP_FN2(name, t1, t2, tc1, tc2, e1, e2, f)               \
  static octave_value                                                   \
  CONCAT2 (oct_catop_, name) (const octave_base_value& a1,              \
                              const octave_base_value& a2,              \
                              const Array<octave_idx_type>& ra_idx)     \
  {                                                                     \
    const CONCAT2 (octave_, t1)& v1 = dynamic_cast<const CONCAT2 (octave_, t1)&> (a1); \
    const CONCAT2 (octave_, t2)& v2 = dynamic_cast<const CONCAT2 (octave_, t2)&> (a2); \
                                                                        \
    return octave_value (tc1 (v1.CONCAT2 (e1, _value) ()) . f (tc2 (v2.CONCAT2 (e2, _value) ()), ra_idx)); \
  }

#endif