Mercurial > octave
view libinterp/corefcn/__dsearchn__.cc @ 28217:87554d9ac6f8 stable
Warn if -v7.3 format is used for load/save (bug #45706).
* load-save.cc (Fload, Fsave): Document that -v7.3 format is not implemented.
* load-save.cc (load_save_system::parse_save_options): Detect '-v7.3' option
and issue an error() that format is not implemented.
* load-save.cc (load_save_system::load): Detect '-v7.3' option and issue an
error() that format is not implemented. Re-order if/elseif tree for clarity.
author | Rik <rik@octave.org> |
---|---|
date | Tue, 14 Apr 2020 11:09:42 -0700 |
parents | bd51beb6205e |
children | 7bc983bc2cbd 0a5b15007766 |
line wrap: on
line source
//////////////////////////////////////////////////////////////////////// // // Copyright (C) 2007-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 <cmath> #include <string> #include "defun.h" #include "error.h" #include "ovl.h" DEFUN (__dsearchn__, args, , doc: /* -*- texinfo -*- @deftypefn {} {[@var{idx}, @var{d}] =} dsearch (@var{x}, @var{xi}) Undocumented internal function. @end deftypefn */) { if (args.length () != 2) print_usage (); Matrix x = args(0).matrix_value ().transpose (); Matrix xi = args(1).matrix_value ().transpose (); if (x.rows () != xi.rows () || x.columns () < 1) error ("__dsearch__: number of rows of X and XI must match"); octave_idx_type n = x.rows (); octave_idx_type nx = x.columns (); octave_idx_type nxi = xi.columns (); ColumnVector idx (nxi); double *pidx = idx.fortran_vec (); ColumnVector dist (nxi); double *pdist = dist.fortran_vec (); #define DIST(dd, y, yi, m) \ dd = 0.0; \ for (octave_idx_type k = 0; k < m; k++) \ { \ double yd = y[k] - yi[k]; \ dd += yd * yd; \ } \ dd = sqrt (dd) const double *pxi = xi.fortran_vec (); for (octave_idx_type i = 0; i < nxi; i++) { double d0; const double *px = x.fortran_vec (); DIST(d0, px, pxi, n); *pidx = 1.; for (octave_idx_type j = 1; j < nx; j++) { px += n; double d; DIST (d, px, pxi, n); if (d < d0) { d0 = d; *pidx = static_cast<double> (j + 1); } octave_quit (); } *pdist++ = d0; pidx++; pxi += n; } return ovl (idx, dist); } /* ## No test needed for internal helper function. %!assert (1) */