view liboctave/util/oct-sparse.cc @ 29274:fed765133d53

sparse-qr.cc: Use correct integer types (bug #57033). * liboctave/numeric/sparse-qr.cc (sparse_qr_rep): Use correct integer types for members. Use cholmod types with SuiteSparse_long for SPQR. Check when converting to narrower integer types. Do not static_cast pointer types. (ros2rcs, cos2ccs, rod2ccd, cod2ccd, rcs2ros, ccs2cos, ros2ccs): Move functions from oct-sparse.cc and use cholmod types with SuiteSparse_long. Check when converting to narrower integer types. Do not static_cast pointer types. Copy arrays to correct integer pointer type when necessary. * liboctave/util/oct-sparse.h (from_suitesparse_long, from_size_t): New functions. * liboctave/util/oct-sparse.cc (ros2rcs, cos2ccs, rod2ccd, cod2ccd, rcs2ros, ccs2cos, ros2ccs): Remove functions that were moved to sparse-qr.cc.
author Markus Mützel <markus.muetzel@gmx.de>
date Thu, 07 Jan 2021 00:19:21 +0100
parents 66f162b6fa03
children 7854d5752dd2
line wrap: on
line source

////////////////////////////////////////////////////////////////////////
//
// Copyright (C) 2017-2020 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 (HAVE_CONFIG_H)
#  include "config.h"
#endif

#include "lo-error.h"
#include "oct-sparse.h"

#if (defined (HAVE_AMD) || defined (HAVE_CCOLAMD)               \
     || defined (HAVE_CHOLMOD) || defined (HAVE_COLAMD)         \
     || defined (HAVE_CXSPARSE) || defined (HAVE_UMFPACK))

namespace octave
{
  static inline void
  check_suitesparse_integer_size (void)
  {
    // FIXME: maybe it would be better to make this a configure check and
    // disable suitesparse if it fails?

    if (sizeof (suitesparse_integer) != sizeof (octave_idx_type))
      (*current_liboctave_error_handler)
        ("size of suitesparse integer does not match octave_idx_type!");
  }

  suitesparse_integer *
  to_suitesparse_intptr (octave_idx_type *i)
  {
    check_suitesparse_integer_size ();

    return reinterpret_cast<suitesparse_integer *> (i);
  }

  const suitesparse_integer *
  to_suitesparse_intptr (const octave_idx_type *i)
  {
    check_suitesparse_integer_size ();

    return reinterpret_cast<const suitesparse_integer *> (i);
  }

  octave_idx_type*
  to_octave_idx_type_ptr (suitesparse_integer *i)
  {
    check_suitesparse_integer_size ();

    return reinterpret_cast<octave_idx_type *> (i);
  }

  const octave_idx_type*
  to_octave_idx_type_ptr (const suitesparse_integer *i)
  {
    check_suitesparse_integer_size ();

    return reinterpret_cast<const octave_idx_type *> (i);
  }
}

#endif