view liboctave/numeric/lo-slatec-proto.h @ 31256:fa4bb329a51a

Allow empty sparse matrix inputs to qr() (bug #63069) * qr.cc (Fqr): Add BIST tests for bug #63069. * sparse-qr.cc (sparse_qr<SparseMatrix>::sparse_qr_rep::sparse_qr_rep, sparse_qr<SparseMatrix>::sparse_qr_rep::C, sparse_qr<SparseMatrix>::sparse_qr_rep::tall_solve<MArray<double>, Matrix>, sparse_qr<SparseComplexMatrix>::sparse_qr_rep::sparse_qr_rep, sparse_qr<SparseComplexMatrix>::sparse_qr_rep::C, sparse_qr<SPARSE_T>::solve): Change input validation from "<= 0" to "< 0". Change error message to "matrix dimension with negative size" from "matrix dimension with negative or zero size".
author Rik <rik@octave.org>
date Tue, 04 Oct 2022 17:05:21 -0700
parents 796f54d4ddbf
children 597f3ee61a48
line wrap: on
line source

////////////////////////////////////////////////////////////////////////
//
// Copyright (C) 2016-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_lo_slatec_proto_h)
#define octave_lo_slatec_proto_h 1

#include "octave-config.h"

#include "f77-fcn.h"

extern "C"
{
  // PCHIM

  F77_RET_T
  F77_FUNC (dpchim, DPCHIM) (const F77_INT& n, const F77_DBLE *x,
                             const F77_DBLE *f, F77_DBLE *d,
                             const F77_INT& incfd,
                             F77_INT& ierr);

  F77_RET_T
  F77_FUNC (pchim, PCHIM) (const F77_INT& n, const F77_REAL *x,
                           const F77_REAL *f, F77_REAL *d,
                           const F77_INT& incfd,
                           F77_INT& ierr);

  // PSIFN

  F77_RET_T
  F77_FUNC (psifn, PSIFN) (const F77_REAL&, const F77_INT&,
                           const F77_INT&, const F77_INT&,
                           F77_REAL&, F77_INT&, F77_INT&);

  F77_RET_T
  F77_FUNC (dpsifn, DPSIFN) (const F77_DBLE&, const F77_INT&,
                             const F77_INT&, const F77_INT&,
                             F77_DBLE&, F77_INT&, F77_INT&);
}

#endif