Mercurial > octave
view libinterp/corefcn/__dsearchn__.cc @ 23738:8acd390d16c9
don't use singleton for stream_list object
* oct-stream.h, oct-stream.cc (stream_list): Don't use singleton
pattern. Change all uses.
(stream_list::stream_list): Initialize stdin, stdout, and stderr
streams here.
(stream_list::~stream_list): Close files here.
(stream_list::m_stdin_file, stream_list::m_stdout_file,
stream_list::m_stderr_file): New data members.
(stream_list::stdin_file, stream_list::stdout_file,
stream_list::stderr_file): New functions.
* interpreter.h, interpreter.cc (interpreter::m_stream_list):
New data member.
(interpreter::interpreter): Initialize it. Don't call
initialize_file_io.
(interpreter::~interpreter): Don't call close_files.
(interpreter::get_stream_list): New function.
* file-io.h, file-io.cc (stdin_file, stdout_file, stderr_file,
stdin_stream, stdout_stream, stderr_stream): Delete static variables.
(initialize_file_io, close_files): Delete. These actions are now
handled by the stream_list constructor and destructor.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Thu, 06 Jul 2017 13:34:31 -0400 |
parents | bd77ab816e43 |
children | 194eb4bd202b |
line wrap: on
line source
/* Copyright (C) 2007-2017 David Bateman 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 <http://www.gnu.org/licenses/>. */ #if defined (HAVE_CONFIG_H) # include "config.h" #endif #include <cmath> #include <iostream> #include <fstream> #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) */