octave: libinterp/corefcn/__eigs_

annotate libinterp/corefcn/eigs.cc @ 29958:32c3a5805893

move DEFUN and DEFMETHOD functions inside octave namespace * mk-opts.pl: Surround emitted function definitions with OCTAVE_NAMESPACE_BEGIN and OCTAVE_NAMESPACE_END tags. * mk-builtins.pl: Surround emitted function declarations with OCTAVE_NAMESPACE_BEGIN and OCTAVE_NAMESPACE_END tags. Emit deprecated global inline functions. * __betainc__.cc, __contourc__.cc, __dsearchn__.cc, __eigs__.cc, __expint__.cc, __ftp__.cc, __gammainc__.cc, __ichol__.cc, __ilu__.cc, __lin_interpn__.cc, __magick_read__.cc, __pchip_deriv__.cc, __qp__.cc, amd.cc, balance.cc, besselj.cc, bitfcns.cc, bsxfun.cc, call-stack.cc, ccolamd.cc, cellfun.cc, chol.cc, colamd.cc, colloc.cc, conv2.cc, daspk.cc, dasrt.cc, dassl.cc, data.cc, debug.cc, defaults.cc, defun.cc, det.cc, dirfns.cc, display.cc, dlmread.cc, dmperm.cc, dot.cc, eig.cc, ellipj.cc, environment.cc, error.cc, event-manager.cc, fcn-info.cc, fft.cc, fft2.cc, fftn.cc, file-io.cc, filter.cc, find.cc, gcd.cc, getgrent.cc, getpwent.cc, getrusage.cc, givens.cc, graphics.cc, gsvd.cc, hash.cc, help.cc, hess.cc, hex2num.cc, input.cc, interpreter.cc, inv.cc, jsondecode.cc, jsonencode.cc, kron.cc, load-path.cc, load-save.cc, lookup.cc, ls-oct-text.cc, lsode.cc, lu.cc, mappers.cc, matrix_type.cc, max.cc, mgorth.cc, nproc.cc, oct-hist.cc, ordqz.cc, ordschur.cc, pager.cc, pinv.cc, pr-flt-fmt.cc, pr-output.cc, psi.cc, qr.cc, quad.cc, quadcc.cc, qz.cc, rand.cc, rcond.cc, regexp.cc, schur.cc, settings.cc, sighandlers.cc, sparse.cc, spparms.cc, sqrtm.cc, stream-euler.cc, strfind.cc, strfns.cc, sub2ind.cc, svd.cc, sylvester.cc, symbfact.cc, symrcm.cc, symtab.cc, syscalls.cc, sysdep.cc, time.cc, toplev.cc, tril.cc, tsearch.cc, typecast.cc, urlwrite.cc, utils.cc, variables.cc, __delaunayn__.cc, __fltk_uigetfile__.cc, __glpk__.cc, __init_gnuplot__.cc, __ode15__.cc, __voronoi__.cc, audiodevinfo.cc, audioread.cc, convhulln.cc, fftw.cc, gzip.cc, ov-base.cc, ov-bool-mat.cc, ov-cell.cc, ov-class.cc, ov-classdef.cc, ov-fcn-handle.cc, ov-java.cc, ov-null-mat.cc, ov-oncleanup.cc, ov-struct.cc, ov-typeinfo.cc, ov-usr-fcn.cc, ov.cc, octave.cc, lex.ll, oct-parse.yy, profiler.cc, pt-eval.cc: Surround DEFUN and DEFMETHOD function defnitions with OCTAVE_NAMESPACE_BEGIN and OCTAVE_NAMESPACE_END tags.

author	John W. Eaton <jwe@octave.org>
date	Fri, 13 Aug 2021 21:53:51 -0400
parents	7b160bf6b897
children	7d6709900da7

rev	line source
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	1 ////////////////////////////////////////////////////////////////////////
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	2 //
29358 0a5b15007766 update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 28024 diff changeset	3 // Copyright (C) 2005-2021 The Octave Project Developers
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	4 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	5 // See the file COPYRIGHT.md in the top-level directory of this
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	6 // distribution or <https://octave.org/copyright/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	7 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	8 // This file is part of Octave.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	9 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	10 // Octave is free software: you can redistribute it and/or modify it
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	11 // under the terms of the GNU General Public License as published by
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	12 // the Free Software Foundation, either version 3 of the License, or
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	13 // (at your option) any later version.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	14 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	15 // Octave is distributed in the hope that it will be useful, but
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	16 // WITHOUT ANY WARRANTY; without even the implied warranty of
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	18 // GNU General Public License for more details.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	19 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	20 // You should have received a copy of the GNU General Public License
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	21 // along with Octave; see the file COPYING. If not, see
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	22 // <https://www.gnu.org/licenses/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	23 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	24 ////////////////////////////////////////////////////////////////////////
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	25
21724 aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef" John W. Eaton <jwe@octave.org> parents: 21691 diff changeset	26 #if defined (HAVE_CONFIG_H)
21301 40de9f8f23a6 Use '#include "config.h"' rather than <config.h>. Rik <rik@octave.org> parents: 21200 diff changeset	27 # include "config.h"
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	28 #endif
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	29
23024 a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	30 #include <limits>
a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	31 #include <string>
a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	32
a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	33 #include "Matrix.h"
a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	34 #include "eigs-base.h"
a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	35 #include "unwind-prot.h"
a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	36
28024 c28b8ba841fb move sparse functions from dldfcn to corefcn directory (bug #57459) John W. Eaton <jwe@octave.org> parents: 27923 diff changeset	37 #include "defun.h"
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	38 #include "error.h"
21100 e39e05d90788 Switch gripe_XXX to either err_XXX or warn_XXX naming scheme. Rik <rik@octave.org> parents: 20956 diff changeset	39 #include "errwarn.h"
26864 041caa61ed34 use get_function_handle instead of extract_function John W. Eaton <jwe@octave.org> parents: 26495 diff changeset	40 #include "interpreter-private.h"
21190 342764537e5a don't install eigs-base.cc John W. Eaton <jwe@octave.org> parents: 21109 diff changeset	41 #include "oct-map.h"
342764537e5a don't install eigs-base.cc John W. Eaton <jwe@octave.org> parents: 21109 diff changeset	42 #include "ov.h"
23024 a6a7b054e4ba Rationalize #includes in libinterp/dldfcn directory. Rik <rik@octave.org> parents: 22755 diff changeset	43 #include "ovl.h"
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	44 #include "pager.h"
23473 55916f99b8b6 prefer feval over do_multi_index_op John W. Eaton <jwe@octave.org> parents: 23457 diff changeset	45 #include "parse.h"
21190 342764537e5a don't install eigs-base.cc John W. Eaton <jwe@octave.org> parents: 21109 diff changeset	46 #include "variables.h"
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	47
21691 263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21575 diff changeset	48 #if defined (HAVE_ARPACK)
263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21575 diff changeset	49
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	50 struct eigs_callback {
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	51 // Pointer for user defined function.
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	52 octave_value eigs_fcn;
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	53
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	54 // Have we warned about imaginary values returned from user function?
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	55 bool warned_imaginary = false;
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	56
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	57 ColumnVector
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	58 eigs_func (const ColumnVector& x, int& eigs_error);
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	59
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	60 ComplexColumnVector
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	61 eigs_complex_func (const ComplexColumnVector& x, int& eigs_error);
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	62 };
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	63
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	64 // Is this a recursive call?
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	65 static int call_depth = 0;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	66
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	67 ColumnVector
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	68 eigs_callback::eigs_func (const ColumnVector& x, int& eigs_error)
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	69 {
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	70 ColumnVector retval;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	71 octave_value_list args;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	72 args(0) = x;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	73
26864 041caa61ed34 use get_function_handle instead of extract_function John W. Eaton <jwe@octave.org> parents: 26495 diff changeset	74 if (eigs_fcn.is_defined ())
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	75 {
20743 b6408331bfa2 eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20703 diff changeset	76 octave_value_list tmp;
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	77
20743 b6408331bfa2 eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20703 diff changeset	78 try
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	79 {
23478 bce43f1a3cce * __eigs__.cc: Avoid deprecated function warnings. John W. Eaton <jwe@octave.org> parents: 23473 diff changeset	80 tmp = octave::feval (eigs_fcn, args, 1);
20743 b6408331bfa2 eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20703 diff changeset	81 }
29163 8f67ad8b3103 maint: Updating naming conventions for exceptions and use const where possible. Rik <rik@octave.org> parents: 28817 diff changeset	82 catch (octave::execution_exception& ee)
20743 b6408331bfa2 eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20703 diff changeset	83 {
29163 8f67ad8b3103 maint: Updating naming conventions for exceptions and use const where possible. Rik <rik@octave.org> parents: 28817 diff changeset	84 err_user_supplied_eval (ee, "eigs");
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	85 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	86
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	87 if (tmp.length () && tmp(0).is_defined ())
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	88 {
23581 c3075ae020e1 maint: Deprecate is_complex_type and replace with iscomplex. Rik <rik@octave.org> parents: 23478 diff changeset	89 if (! warned_imaginary && tmp(0).iscomplex ())
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	90 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	91 warning ("eigs: ignoring imaginary part returned from user-supplied function");
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	92 warned_imaginary = true;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	93 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	94
20743 b6408331bfa2 eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20703 diff changeset	95 retval = tmp(0).xvector_value ("eigs: evaluation of user-supplied function failed");
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	96 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	97 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	98 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	99 eigs_error = 1;
21100 e39e05d90788 Switch gripe_XXX to either err_XXX or warn_XXX naming scheme. Rik <rik@octave.org> parents: 20956 diff changeset	100 err_user_supplied_eval ("eigs");
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	101 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	102 }
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	103
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	104 return retval;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	105 }
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	106
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	107 ComplexColumnVector
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	108 eigs_callback::eigs_complex_func (const ComplexColumnVector& x,
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	109 int& eigs_error)
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	110 {
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	111 ComplexColumnVector retval;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	112 octave_value_list args;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	113 args(0) = x;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	114
26864 041caa61ed34 use get_function_handle instead of extract_function John W. Eaton <jwe@octave.org> parents: 26495 diff changeset	115 if (eigs_fcn.is_defined ())
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	116 {
20743 b6408331bfa2 eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20703 diff changeset	117 octave_value_list tmp;
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	118
20743 b6408331bfa2 eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20703 diff changeset	119 try
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	120 {
23478 bce43f1a3cce * __eigs__.cc: Avoid deprecated function warnings. John W. Eaton <jwe@octave.org> parents: 23473 diff changeset	121 tmp = octave::feval (eigs_fcn, args, 1);
20743 b6408331bfa2 eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20703 diff changeset	122 }
29163 8f67ad8b3103 maint: Updating naming conventions for exceptions and use const where possible. Rik <rik@octave.org> parents: 28817 diff changeset	123 catch (octave::execution_exception& ee)
20743 b6408331bfa2 eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20703 diff changeset	124 {
29163 8f67ad8b3103 maint: Updating naming conventions for exceptions and use const where possible. Rik <rik@octave.org> parents: 28817 diff changeset	125 err_user_supplied_eval (ee, "eigs");
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	126 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	127
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	128 if (tmp.length () && tmp(0).is_defined ())
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	129 {
26344 acfd683ea00b __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26164 diff changeset	130 retval = tmp(0).xcomplex_vector_value ("eigs: evaluation of user-supplied function failed");
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	131 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	132 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	133 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	134 eigs_error = 1;
21100 e39e05d90788 Switch gripe_XXX to either err_XXX or warn_XXX naming scheme. Rik <rik@octave.org> parents: 20956 diff changeset	135 err_user_supplied_eval ("eigs");
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	136 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	137 }
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	138
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	139 return retval;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	140 }
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	141
21691 263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21575 diff changeset	142 #endif
263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21575 diff changeset	143
29958 32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace John W. Eaton <jwe@octave.org> parents: 29481 diff changeset	144 OCTAVE_NAMESPACE_BEGIN
32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace John W. Eaton <jwe@octave.org> parents: 29481 diff changeset	145
28024 c28b8ba841fb move sparse functions from dldfcn to corefcn directory (bug #57459) John W. Eaton <jwe@octave.org> parents: 27923 diff changeset	146 DEFMETHOD (__eigs__, interp, args, nargout,
c28b8ba841fb move sparse functions from dldfcn to corefcn directory (bug #57459) John W. Eaton <jwe@octave.org> parents: 27923 diff changeset	147 doc: /* -- texinfo --
21966 112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	148 @deftypefn {} {@var{d} =} __eigs__ (@var{A})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	149 @deftypefnx {} {@var{d} =} __eigs__ (@var{A}, @var{k})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	150 @deftypefnx {} {@var{d} =} __eigs__ (@var{A}, @var{k}, @var{sigma})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	151 @deftypefnx {} {@var{d} =} __eigs__ (@var{A}, @var{k}, @var{sigma}, @var{opts})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	152 @deftypefnx {} {@var{d} =} __eigs__ (@var{A}, @var{B})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	153 @deftypefnx {} {@var{d} =} __eigs__ (@var{A}, @var{B}, @var{k})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	154 @deftypefnx {} {@var{d} =} __eigs__ (@var{A}, @var{B}, @var{k}, @var{sigma})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	155 @deftypefnx {} {@var{d} =} __eigs__ (@var{A}, @var{B}, @var{k}, @var{sigma}, @var{opts})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	156 @deftypefnx {} {@var{d} =} __eigs__ (@var{af}, @var{n})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	157 @deftypefnx {} {@var{d} =} __eigs__ (@var{af}, @var{n}, @var{B})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	158 @deftypefnx {} {@var{d} =} __eigs__ (@var{af}, @var{n}, @var{k})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	159 @deftypefnx {} {@var{d} =} __eigs__ (@var{af}, @var{n}, @var{B}, @var{k})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	160 @deftypefnx {} {@var{d} =} __eigs__ (@var{af}, @var{n}, @var{k}, @var{sigma})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	161 @deftypefnx {} {@var{d} =} __eigs__ (@var{af}, @var{n}, @var{B}, @var{k}, @var{sigma})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	162 @deftypefnx {} {@var{d} =} __eigs__ (@var{af}, @var{n}, @var{k}, @var{sigma}, @var{opts})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	163 @deftypefnx {} {@var{d} =} __eigs__ (@var{af}, @var{n}, @var{B}, @var{k}, @var{sigma}, @var{opts})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	164 @deftypefnx {} {[@var{V}, @var{d}] =} __eigs__ (@var{A}, @dots{})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	165 @deftypefnx {} {[@var{V}, @var{d}] =} __eigs__ (@var{af}, @var{n}, @dots{})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	166 @deftypefnx {} {[@var{V}, @var{d}, @var{flag}] =} __eigs__ (@var{A}, @dots{})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	167 @deftypefnx {} {[@var{V}, @var{d}, @var{flag}] =} __eigs__ (@var{af}, @var{n}, @dots{})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	168 Undocumented internal function.
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21743 diff changeset	169 @end deftypefn */)
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	170 {
21724 aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef" John W. Eaton <jwe@octave.org> parents: 21691 diff changeset	171 #if defined (HAVE_ARPACK)
21691 263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21575 diff changeset	172
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	173 int nargin = args.length ();
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	174
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	175 if (nargin == 0)
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	176 print_usage ();
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	177
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	178 octave_value_list retval;
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	179
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	180 std::string fcn_name;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	181 octave_idx_type n = 0;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	182 octave_idx_type k = 6;
26344 acfd683ea00b __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26164 diff changeset	183 Complex sigma = 0.0;
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	184 double sigmar, sigmai;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	185 bool have_sigma = false;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	186 std::string typ = "LM";
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	187 Matrix amm, bmm, bmt;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	188 ComplexMatrix acm, bcm, bct;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	189 SparseMatrix asmm, bsmm, bsmt;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	190 SparseComplexMatrix ascm, bscm, bsct;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	191 int b_arg = 0;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	192 bool have_b = false;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	193 bool have_a_fun = false;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	194 bool a_is_complex = false;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	195 bool b_is_complex = false;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	196 bool symmetric = false;
12649 8579e7a6e6e0 Improve performance (2X) of svds by avoiding issymmetric calculation (bug #33073) Rik <octave@nomad.inbox5.com> parents: 12329 diff changeset	197 bool sym_tested = false;
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	198 bool cholB = false;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	199 bool a_is_sparse = false;
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	200 bool b_is_sparse = false;
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	201 ColumnVector permB;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	202 int arg_offset = 0;
15220 61822c866ba1 use std::numeric_limits<T>::epsilon in C++ code John W. Eaton <jwe@octave.org> parents: 15195 diff changeset	203 double tol = std::numeric_limits<double>::epsilon ();
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	204 int maxit = 300;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	205 int disp = 0;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	206 octave_idx_type p = -1;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	207 ColumnVector resid;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	208 ComplexColumnVector cresid;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	209 octave_idx_type info = 1;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	210
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	211 eigs_callback callback;
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	212
28817 aef839f40da2 Replace unwind_protect with more efficient constructs (bug #59192). Rik <rik@octave.org> parents: 28024 diff changeset	213 octave::unwind_protect_var<int> restore_var (call_depth);
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	214 call_depth++;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	215
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	216 if (call_depth > 1)
20956 850e3d2533d4 maint: Eliminate more useless statements after error(). Rik <rik@octave.org> parents: 20955 diff changeset	217 error ("eigs: invalid recursive call");
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	218
20790 c2d9556d51d0 eliminate return statements after calls to print_usage John W. Eaton <jwe@octave.org> parents: 20785 diff changeset	219 if (args(0).is_function_handle () \|\| args(0).is_inline_function ()
c2d9556d51d0 eliminate return statements after calls to print_usage John W. Eaton <jwe@octave.org> parents: 20785 diff changeset	220 \|\| args(0).is_string ())
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	221 {
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	222 callback.eigs_fcn = octave::get_function_handle (interp, args(0), "x");
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	223
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	224 if (callback.eigs_fcn.is_undefined ())
20825 66cd994d1b79 eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20790 diff changeset	225 error ("eigs: unknown function");
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	226
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	227 if (nargin < 2)
20825 66cd994d1b79 eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20790 diff changeset	228 error ("eigs: incorrect number of arguments");
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	229
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	230 n = args(1).nint_value ();
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	231 arg_offset = 1;
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	232 have_a_fun = true;
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	233 }
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	234 else
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	235 {
23581 c3075ae020e1 maint: Deprecate is_complex_type and replace with iscomplex. Rik <rik@octave.org> parents: 23478 diff changeset	236 if (args(0).iscomplex ())
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	237 {
23583 b7747a2c88b2 maint: Deprecate is_sparse_type and replace with issparse. Rik <rik@octave.org> parents: 23581 diff changeset	238 if (args(0).issparse ())
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	239 {
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14844 diff changeset	240 ascm = (args(0).sparse_complex_matrix_value ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	241 a_is_sparse = true;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	242 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	243 else
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14844 diff changeset	244 acm = (args(0).complex_matrix_value ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	245 a_is_complex = true;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	246 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	247 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	248 {
23583 b7747a2c88b2 maint: Deprecate is_sparse_type and replace with issparse. Rik <rik@octave.org> parents: 23581 diff changeset	249 if (args(0).issparse ())
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	250 {
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14844 diff changeset	251 asmm = (args(0).sparse_matrix_value ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	252 a_is_sparse = true;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	253 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	254 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	255 {
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14844 diff changeset	256 amm = (args(0).matrix_value ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	257 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	258 }
9593 cf8403208c43 Fix nesting error in options parsing of eigs David Bateman <dbateman@free.fr> parents: 9377 diff changeset	259 }
cf8403208c43 Fix nesting error in options parsing of eigs David Bateman <dbateman@free.fr> parents: 9377 diff changeset	260
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	261 // Note hold off reading B until later to avoid issues of double
12649 8579e7a6e6e0 Improve performance (2X) of svds by avoiding issymmetric calculation (bug #33073) Rik <octave@nomad.inbox5.com> parents: 12329 diff changeset	262 // copies of the matrix if B is full/real while A is complex.
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	263 if (nargin > 1 + arg_offset
19864 17d647821d61 maint: More cleanup of C++ code to follow Octave coding conventions. John W. Eaton <jwe@octave.org> parents: 19697 diff changeset	264 && ! (args(1 + arg_offset).is_real_scalar ()))
9593 cf8403208c43 Fix nesting error in options parsing of eigs David Bateman <dbateman@free.fr> parents: 9377 diff changeset	265 {
23581 c3075ae020e1 maint: Deprecate is_complex_type and replace with iscomplex. Rik <rik@octave.org> parents: 23478 diff changeset	266 if (args(1+arg_offset).iscomplex ())
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	267 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	268 b_arg = 1+arg_offset;
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	269 if (args(b_arg).issparse ())
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	270 {
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	271 bscm = (args(b_arg).sparse_complex_matrix_value ());
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	272 b_is_sparse = true;
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	273 }
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	274 else
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	275 bcm = (args(b_arg).complex_matrix_value ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	276 have_b = true;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	277 b_is_complex = true;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	278 arg_offset++;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	279 }
9593 cf8403208c43 Fix nesting error in options parsing of eigs David Bateman <dbateman@free.fr> parents: 9377 diff changeset	280 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	281 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	282 b_arg = 1+arg_offset;
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	283 if (args(b_arg).issparse ())
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	284 {
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	285 bsmm = (args(b_arg).sparse_matrix_value ());
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	286 b_is_sparse = true;
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	287 }
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	288 else
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	289 bmm = (args(b_arg).matrix_value ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	290 have_b = true;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	291 arg_offset++;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	292 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	293 }
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	294
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	295 if (nargin > (1+arg_offset))
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	296 k = args(1+arg_offset).nint_value ();
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	297
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	298 if (nargin > (2+arg_offset))
8683 e4f89f4a7cf8 eigs.cc (Feigs): avoid warnings; style fixes; return after errors John W. Eaton <jwe@octave.org> parents: 8585 diff changeset	299 {
8845 5a6db6bd1a02 eigs.cc (Feigs): fix handling of sigma arg John W. Eaton <jwe@octave.org> parents: 8702 diff changeset	300 if (args(2+arg_offset).is_string ())
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	301 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	302 typ = args(2+arg_offset).string_value ();
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	303
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	304 // Use STL function to convert to upper case
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	305 transform (typ.begin (), typ.end (), typ.begin (), toupper);
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	306
26344 acfd683ea00b __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26164 diff changeset	307 sigma = 0.0;
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	308 }
8683 e4f89f4a7cf8 eigs.cc (Feigs): avoid warnings; style fixes; return after errors John W. Eaton <jwe@octave.org> parents: 8585 diff changeset	309 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	310 {
20703 85e5efae848a eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20687 diff changeset	311 sigma = args(2+arg_offset).xcomplex_value ("eigs: SIGMA must be a scalar or a string");
8845 5a6db6bd1a02 eigs.cc (Feigs): fix handling of sigma arg John W. Eaton <jwe@octave.org> parents: 8702 diff changeset	312
20703 85e5efae848a eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20687 diff changeset	313 have_sigma = true;
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	314 }
8683 e4f89f4a7cf8 eigs.cc (Feigs): avoid warnings; style fixes; return after errors John W. Eaton <jwe@octave.org> parents: 8585 diff changeset	315 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	316
22274 edcad35c364f don't rely on std:: overloads of arg, conj, imag, and real for real args John W. Eaton <jwe@octave.org> parents: 22022 diff changeset	317 sigmar = sigma.real ();
edcad35c364f don't rely on std:: overloads of arg, conj, imag, and real for real args John W. Eaton <jwe@octave.org> parents: 22022 diff changeset	318 sigmai = sigma.imag ();
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	319
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	320 if (nargin > (3+arg_offset))
8683 e4f89f4a7cf8 eigs.cc (Feigs): avoid warnings; style fixes; return after errors John W. Eaton <jwe@octave.org> parents: 8585 diff changeset	321 {
23584 7ed6b258db91 maint: Deprecate is_map and replace with isstruct. Rik <rik@octave.org> parents: 23583 diff changeset	322 if (! args(3+arg_offset).isstruct ())
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	323 error ("eigs: OPTS argument must be a structure");
20703 85e5efae848a eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20687 diff changeset	324
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	325 octave_scalar_map map = args(3+arg_offset).xscalar_map_value ("eigs: OPTS argument must be a scalar structure");
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	326
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	327 octave_value tmp;
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	328
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	329 // issym is ignored for complex matrix inputs
27692 e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	330 tmp = map.getfield ("issym");
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	331 if (tmp.is_defined ())
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	332 {
27692 e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	333 if (tmp.numel () != 1)
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	334 error ("eigs: OPTS.issym must be a scalar value");
26495 f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	335
27692 e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	336 symmetric = tmp.xbool_value ("eigs: OPTS.issym must be a logical value");
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	337 sym_tested = true;
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	338 }
11044 757249bb4888 eigs.cc: use octave_scalar_map instead of Octave_map John W. Eaton <jwe@octave.org> parents: 10846 diff changeset	339
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	340 // isreal is ignored if A is not a function
26495 f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	341 if (have_a_fun)
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	342 {
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	343 tmp = map.getfield ("isreal");
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	344 if (tmp.is_defined ())
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	345 {
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	346 if (tmp.numel () != 1)
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	347 error ("eigs: OPTS.isreal must be a scalar value");
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	348
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	349 a_is_complex = ! tmp.xbool_value ("eigs: OPTS.isreal must be a logical value");
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	350 }
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	351 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	352
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	353 tmp = map.getfield ("tol");
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	354 if (tmp.is_defined ())
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	355 tol = tmp.double_value ();
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	356
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	357 tmp = map.getfield ("maxit");
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	358 if (tmp.is_defined ())
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	359 maxit = tmp.nint_value ();
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	360
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	361 tmp = map.getfield ("p");
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	362 if (tmp.is_defined ())
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	363 p = tmp.nint_value ();
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	364
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	365 tmp = map.getfield ("v0");
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	366 if (tmp.is_defined ())
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	367 {
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	368 if (a_is_complex \|\| b_is_complex)
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	369 cresid = ComplexColumnVector (tmp.complex_vector_value ());
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	370 else
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	371 resid = ColumnVector (tmp.vector_value ());
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	372 }
20703 85e5efae848a eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20687 diff changeset	373
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	374 tmp = map.getfield ("disp");
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	375 if (tmp.is_defined ())
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	376 disp = tmp.nint_value ();
20703 85e5efae848a eliminate more uses of error_state John W. Eaton <jwe@octave.org> parents: 20687 diff changeset	377
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	378 tmp = map.getfield ("cholB");
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	379 if (tmp.is_defined ())
26495 f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	380 {
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	381 if (tmp.numel () != 1)
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	382 error ("eigs: OPTS.cholB must be a scalar value");
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	383
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	384 cholB = tmp.xbool_value ("eigs: OPTS.cholB must be a logical value");
f0065630ad9d __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	385 }
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	386
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	387 tmp = map.getfield ("permB");
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	388 if (tmp.is_defined ())
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	389 permB = ColumnVector (tmp.vector_value ()) - 1.0;
8683 e4f89f4a7cf8 eigs.cc (Feigs): avoid warnings; style fixes; return after errors John W. Eaton <jwe@octave.org> parents: 8585 diff changeset	390 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	391
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	392 if (nargin > (4+arg_offset))
20825 66cd994d1b79 eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20790 diff changeset	393 error ("eigs: incorrect number of arguments");
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	394
12649 8579e7a6e6e0 Improve performance (2X) of svds by avoiding issymmetric calculation (bug #33073) Rik <octave@nomad.inbox5.com> parents: 12329 diff changeset	395 // Test undeclared (no issym) matrix inputs for symmetry
20955 77f5591878bf maint: Use '! expr' rather than '!expr' to conform to coding guidelines. Rik <rik@octave.org> parents: 20939 diff changeset	396 if (! sym_tested && ! have_a_fun)
12649 8579e7a6e6e0 Improve performance (2X) of svds by avoiding issymmetric calculation (bug #33073) Rik <octave@nomad.inbox5.com> parents: 12329 diff changeset	397 {
27692 e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	398 if (a_is_complex)
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	399 {
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	400 if (a_is_sparse)
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	401 symmetric = ascm.ishermitian ();
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	402 else
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	403 symmetric = acm.ishermitian ();
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	404 }
12649 8579e7a6e6e0 Improve performance (2X) of svds by avoiding issymmetric calculation (bug #33073) Rik <octave@nomad.inbox5.com> parents: 12329 diff changeset	405 else
27692 e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	406 {
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	407 if (a_is_sparse)
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	408 symmetric = asmm.issymmetric ();
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	409 else
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	410 symmetric = amm.issymmetric ();
e557774a2604 Return real values for Hermitian inputs to eigs (bug #57196). Rik <rik@octave.org> parents: 27683 diff changeset	411 }
12649 8579e7a6e6e0 Improve performance (2X) of svds by avoiding issymmetric calculation (bug #33073) Rik <octave@nomad.inbox5.com> parents: 12329 diff changeset	412 }
8579e7a6e6e0 Improve performance (2X) of svds by avoiding issymmetric calculation (bug #33073) Rik <octave@nomad.inbox5.com> parents: 12329 diff changeset	413
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	414 if (have_b)
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	415 {
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	416 if (a_is_complex \|\| b_is_complex)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	417 {
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	418 if (b_is_sparse)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	419 bscm = args(b_arg).sparse_complex_matrix_value ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	420 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	421 bcm = args(b_arg).complex_matrix_value ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	422 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	423 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	424 {
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	425 if (b_is_sparse)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	426 bsmm = args(b_arg).sparse_matrix_value ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	427 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	428 bmm = args(b_arg).matrix_value ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	429 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	430 }
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	431
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11553 diff changeset	432 // Mode 1 for SM mode seems unstable for some reason.
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	433 // Use Mode 3 instead, with sigma = 0.
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	434 if (! have_sigma && typ == "SM")
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	435 have_sigma = true;
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	436
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	437 octave_idx_type nconv;
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	438 if (a_is_complex \|\| b_is_complex)
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	439 {
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	440 EigsComplexFunc
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	441 eigs_complex_fcn = [&callback] (const ComplexColumnVector& x,
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	442 int& eigs_error)
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	443 {
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	444 return callback.eigs_complex_func (x, eigs_error);
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	445 };
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	446
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	447 ComplexMatrix eig_vec;
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	448 ComplexColumnVector eig_val;
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	449
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	450 if (have_a_fun)
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	451 {
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	452 if (b_is_sparse)
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	453 nconv = EigsComplexNonSymmetricFunc
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	454 (eigs_complex_fcn, n, typ, sigma, k, p, info, eig_vec,
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	455 eig_val, bscm, permB, cresid, octave_stdout, tol,
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	456 (nargout > 1), cholB, disp, maxit);
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	457 else
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	458 nconv = EigsComplexNonSymmetricFunc
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	459 (eigs_complex_fcn, n, typ, sigma, k, p, info, eig_vec,
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	460 eig_val, bcm, permB, cresid, octave_stdout, tol,
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	461 (nargout > 1), cholB, disp, maxit);
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	462 }
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	463 else if (have_sigma)
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	464 {
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	465 if (a_is_sparse)
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	466 nconv = EigsComplexNonSymmetricMatrixShift
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	467 (ascm, sigma, k, p, info, eig_vec, eig_val, bscm, permB,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	468 cresid, octave_stdout, tol, (nargout > 1), cholB, disp,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	469 maxit);
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	470 else
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	471 nconv = EigsComplexNonSymmetricMatrixShift
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	472 (acm, sigma, k, p, info, eig_vec, eig_val, bcm, permB,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	473 cresid, octave_stdout, tol, (nargout > 1), cholB, disp,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	474 maxit);
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	475 }
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	476 else
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	477 {
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	478 if (a_is_sparse)
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	479 nconv = EigsComplexNonSymmetricMatrix
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	480 (ascm, typ, k, p, info, eig_vec, eig_val, bscm, permB,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	481 cresid, octave_stdout, tol, (nargout > 1), cholB, disp,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	482 maxit);
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	483 else
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	484 nconv = EigsComplexNonSymmetricMatrix
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	485 (acm, typ, k, p, info, eig_vec, eig_val, bcm, permB,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	486 cresid, octave_stdout, tol, (nargout > 1), cholB, disp,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	487 maxit);
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	488 }
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	489
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	490 if (nargout < 2)
27683 6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	491 {
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	492 if (symmetric)
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	493 retval(0) = real (eig_val);
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	494 else
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	495 retval(0) = eig_val;
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	496 }
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	497 else
27683 6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	498 {
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	499 if (symmetric)
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	500 retval = ovl (eig_vec, DiagMatrix (real (eig_val)), double (info));
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	501 else
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	502 retval = ovl (eig_vec, ComplexDiagMatrix (eig_val), double (info));
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	503 }
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	504 }
26344 acfd683ea00b __eigs__.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26164 diff changeset	505 else if (sigmai != 0.0)
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	506 {
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	507 EigsComplexFunc
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	508 eigs_complex_fcn = [&callback] (const ComplexColumnVector& x,
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	509 int& eigs_error)
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	510 {
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	511 return callback.eigs_complex_func (x, eigs_error);
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	512 };
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	513
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	514 // Promote real problem to a complex one.
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	515 ComplexMatrix eig_vec;
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	516 ComplexColumnVector eig_val;
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	517
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	518 if (have_a_fun)
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	519 {
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	520 if (b_is_sparse)
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	521 nconv = EigsComplexNonSymmetricFunc
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	522 (eigs_complex_fcn, n, typ, sigma, k, p, info, eig_vec,
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	523 eig_val, bscm, permB, cresid, octave_stdout, tol,
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	524 (nargout > 1), cholB, disp, maxit);
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	525 else
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	526 nconv = EigsComplexNonSymmetricFunc
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	527 (eigs_complex_fcn, n, typ, sigma, k, p, info, eig_vec,
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	528 eig_val, bcm, permB, cresid, octave_stdout, tol,
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	529 (nargout > 1), cholB, disp, maxit);
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	530 }
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	531 else
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	532 {
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	533 if (a_is_sparse)
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	534 nconv = EigsComplexNonSymmetricMatrixShift
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	535 (SparseComplexMatrix (asmm), sigma, k, p, info, eig_vec,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	536 eig_val, SparseComplexMatrix (bsmm), permB, cresid,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	537 octave_stdout, tol, (nargout > 1), cholB, disp, maxit);
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	538 else
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	539 nconv = EigsComplexNonSymmetricMatrixShift
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	540 (ComplexMatrix (amm), sigma, k, p, info, eig_vec,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	541 eig_val, ComplexMatrix (bmm), permB, cresid,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	542 octave_stdout, tol, (nargout > 1), cholB, disp, maxit);
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	543 }
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	544
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	545 if (nargout < 2)
27683 6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	546 {
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	547 if (symmetric)
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	548 retval(0) = real (eig_val);
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	549 else
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	550 retval(0) = eig_val;
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	551 }
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	552 else
27683 6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	553 {
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	554 if (symmetric)
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	555 retval = ovl (eig_vec, DiagMatrix (real (eig_val)), double (info));
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	556 else
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	557 retval = ovl (eig_vec, ComplexDiagMatrix (eig_val), double (info));
6604b41ca63f Return only real eigenvalues from eigs() for symmetric matrices (bug #57196). Rik <rik@octave.org> parents: 26864 diff changeset	558 }
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	559 }
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	560 else
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	561 {
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	562 EigsFunc eigs_fcn = [&callback] (const ColumnVector& x, int& eigs_error)
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	563 {
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	564 return callback.eigs_func (x, eigs_error);
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	565 };
7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	566
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	567 if (symmetric)
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	568 {
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	569 Matrix eig_vec;
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	570 ColumnVector eig_val;
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	571
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	572 if (have_a_fun)
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	573 {
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	574 if (b_is_sparse)
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	575 nconv = EigsRealSymmetricFunc
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	576 (eigs_fcn, n, typ, sigmar, k, p, info, eig_vec,
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	577 eig_val, bsmm, permB, resid, octave_stdout, tol,
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	578 (nargout > 1), cholB, disp, maxit);
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	579 else
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	580 nconv = EigsRealSymmetricFunc
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	581 (eigs_fcn, n, typ, sigmar, k, p, info, eig_vec,
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	582 eig_val, bmm, permB, resid, octave_stdout, tol,
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	583 (nargout > 1), cholB, disp, maxit);
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	584 }
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	585 else if (have_sigma)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	586 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	587 if (a_is_sparse)
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	588 nconv = EigsRealSymmetricMatrixShift
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	589 (asmm, sigmar, k, p, info, eig_vec, eig_val, bsmm,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	590 permB, resid, octave_stdout, tol, (nargout > 1),
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	591 cholB, disp, maxit);
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	592 else
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	593 nconv = EigsRealSymmetricMatrixShift
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	594 (amm, sigmar, k, p, info, eig_vec, eig_val, bmm,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	595 permB, resid, octave_stdout, tol, (nargout > 1),
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	596 cholB, disp, maxit);
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	597 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	598 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	599 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	600 if (a_is_sparse)
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	601 nconv = EigsRealSymmetricMatrix
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	602 (asmm, typ, k, p, info, eig_vec, eig_val, bsmm,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	603 permB, resid, octave_stdout, tol, (nargout > 1),
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	604 cholB, disp, maxit);
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	605 else
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	606 nconv = EigsRealSymmetricMatrix
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	607 (amm, typ, k, p, info, eig_vec, eig_val, bmm, permB,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	608 resid, octave_stdout, tol, (nargout > 1), cholB,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	609 disp, maxit);
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	610 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	611
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	612 if (nargout < 2)
14844 5bc9b9cb4362 maint: Use Octave coding conventions for cuddled parenthesis in retval assignments. Rik <octave@nomad.inbox5.com> parents: 14501 diff changeset	613 retval(0) = eig_val;
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	614 else
20871 1e4edab2d22f 2015 Code Sprint: __eigs__.cc: use ovl (). Rik <rik@octave.org> parents: 20853 diff changeset	615 retval = ovl (eig_vec, DiagMatrix (eig_val), double (info));
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	616 }
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	617 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	618 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	619 ComplexMatrix eig_vec;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	620 ComplexColumnVector eig_val;
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	621
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	622 if (have_a_fun)
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	623 {
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	624 if (b_is_sparse)
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	625 nconv = EigsRealNonSymmetricFunc
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	626 (eigs_fcn, n, typ, sigmar, k, p, info, eig_vec,
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	627 eig_val, bsmm, permB, resid, octave_stdout, tol,
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	628 (nargout > 1), cholB, disp, maxit);
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	629 else
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	630 nconv = EigsRealNonSymmetricFunc
29481 7b160bf6b897 eigs: Remove static variables (bug #60311). Markus Mützel <markus.muetzel@gmx.de> parents: 29480 diff changeset	631 (eigs_fcn, n, typ, sigmar, k, p, info, eig_vec,
25696 3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	632 eig_val, bmm, permB, resid, octave_stdout, tol,
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	633 (nargout > 1), cholB, disp, maxit);
3b6691ff0f59 Fix eigs for generalized problems with user defined function (bug #54167). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	634 }
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	635 else if (have_sigma)
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	636 {
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	637 if (a_is_sparse)
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	638 nconv = EigsRealNonSymmetricMatrixShift
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	639 (asmm, sigmar, k, p, info, eig_vec, eig_val, bsmm,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	640 permB, resid, octave_stdout, tol, (nargout > 1),
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	641 cholB, disp, maxit);
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	642 else
4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	643 nconv = EigsRealNonSymmetricMatrixShift
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	644 (amm, sigmar, k, p, info, eig_vec, eig_val, bmm,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	645 permB, resid, octave_stdout, tol, (nargout > 1),
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	646 cholB, disp, maxit);
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	647 }
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	648 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	649 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	650 if (a_is_sparse)
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	651 nconv = EigsRealNonSymmetricMatrix
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	652 (asmm, typ, k, p, info, eig_vec, eig_val, bsmm,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	653 permB, resid, octave_stdout, tol, (nargout > 1),
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	654 cholB, disp, maxit);
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	655 else
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	656 nconv = EigsRealNonSymmetricMatrix
22407 34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	657 (amm, typ, k, p, info, eig_vec, eig_val, bmm, permB,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	658 resid, octave_stdout, tol, (nargout > 1), cholB,
34ce5be04942 maint: Style check C++ code in libinterp/. Rik <rik@octave.org> parents: 22327 diff changeset	659 disp, maxit);
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	660 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	661
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	662 if (nargout < 2)
14844 5bc9b9cb4362 maint: Use Octave coding conventions for cuddled parenthesis in retval assignments. Rik <octave@nomad.inbox5.com> parents: 14501 diff changeset	663 retval(0) = eig_val;
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	664 else
20871 1e4edab2d22f 2015 Code Sprint: __eigs__.cc: use ovl (). Rik <rik@octave.org> parents: 20853 diff changeset	665 retval = ovl (eig_vec, ComplexDiagMatrix (eig_val), double (info));
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 10066 diff changeset	666 }
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	667 }
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	668
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	669 if (nconv <= 0)
24334 d3d38a443df8 Add new warning ID "Octave:eigs:UnconvergedEigenvalues". Rik <rik@octave.org> parents: 24119 diff changeset	670 warning_with_id ("Octave:eigs:UnconvergedEigenvalues",
26164 7f6a50f73625 Silence compiler warnings about format identifier for octave_idx_type (bug #55046). Markus Mützel <markus.muetzel@gmx.de> parents: 25696 diff changeset	671 "eigs: None of the %" OCTAVE_IDX_TYPE_FORMAT
7f6a50f73625 Silence compiler warnings about format identifier for octave_idx_type (bug #55046). Markus Mützel <markus.muetzel@gmx.de> parents: 25696 diff changeset	672 " requested eigenvalues converged", k);
20556 4bed806ee3d4 eliminate more simple uses of error_state John W. Eaton <jwe@octave.org> parents: 19864 diff changeset	673 else if (nconv < k)
24334 d3d38a443df8 Add new warning ID "Octave:eigs:UnconvergedEigenvalues". Rik <rik@octave.org> parents: 24119 diff changeset	674 warning_with_id ("Octave:eigs:UnconvergedEigenvalues",
26164 7f6a50f73625 Silence compiler warnings about format identifier for octave_idx_type (bug #55046). Markus Mützel <markus.muetzel@gmx.de> parents: 25696 diff changeset	675 "eigs: Only %" OCTAVE_IDX_TYPE_FORMAT
7f6a50f73625 Silence compiler warnings about format identifier for octave_idx_type (bug #55046). Markus Mützel <markus.muetzel@gmx.de> parents: 25696 diff changeset	676 " of the %" OCTAVE_IDX_TYPE_FORMAT
7f6a50f73625 Silence compiler warnings about format identifier for octave_idx_type (bug #55046). Markus Mützel <markus.muetzel@gmx.de> parents: 25696 diff changeset	677 " requested eigenvalues converged",
24334 d3d38a443df8 Add new warning ID "Octave:eigs:UnconvergedEigenvalues". Rik <rik@octave.org> parents: 24119 diff changeset	678 nconv, k);
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	679
654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	680 if (! fcn_name.empty ())
23699 b29904962d2d deprecate some global functions that access the symbol table John W. Eaton <jwe@octave.org> parents: 23596 diff changeset	681 {
b29904962d2d deprecate some global functions that access the symbol table John W. Eaton <jwe@octave.org> parents: 23596 diff changeset	682 octave::symbol_table& symtab = interp.get_symbol_table ();
b29904962d2d deprecate some global functions that access the symbol table John W. Eaton <jwe@octave.org> parents: 23596 diff changeset	683
b29904962d2d deprecate some global functions that access the symbol table John W. Eaton <jwe@octave.org> parents: 23596 diff changeset	684 symtab.clear_function (fcn_name);
b29904962d2d deprecate some global functions that access the symbol table John W. Eaton <jwe@octave.org> parents: 23596 diff changeset	685 }
20825 66cd994d1b79 eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20790 diff changeset	686
20939 b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	687 return retval;
b17fda023ca6 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20871 diff changeset	688
14144 834df9f10963 remove ARPACK files from sources and restore configure checks for external ARPACK library John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	689 #else
21691 263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21575 diff changeset	690
24119 3cdd0705745b maint: silence compiler warning when building without ARPACK Mike Miller <mtmiller@octave.org> parents: 23699 diff changeset	691 octave_unused_parameter (interp);
21691 263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21575 diff changeset	692 octave_unused_parameter (args);
263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21575 diff changeset	693 octave_unused_parameter (nargout);
263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21575 diff changeset	694
21109 bd1752782e56 Use err_disabled_feature, warn_disabled_feature throughout code base. Rik <rik@octave.org> parents: 21100 diff changeset	695 err_disabled_feature ("eigs", "ARPACK");
21691 263d18409fdf Eliminate unused variable warnings for conditionally compiled code. John W. Eaton <jwe@octave.org> parents: 21575 diff changeset	696
14144 834df9f10963 remove ARPACK files from sources and restore configure checks for external ARPACK library John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	697 #endif
8417 654bcfb937bf Add the eigs and svds functions David Bateman <dbateman@free.fr> parents: diff changeset	698 }
21575 bc9aa534bc29 maint: Clean up BIST tests so they don't produce warnings. Rik <rik@octave.org> parents: 21301 diff changeset	699
bc9aa534bc29 maint: Clean up BIST tests so they don't produce warnings. Rik <rik@octave.org> parents: 21301 diff changeset	700 /*
bc9aa534bc29 maint: Clean up BIST tests so they don't produce warnings. Rik <rik@octave.org> parents: 21301 diff changeset	701 ## No test needed for internal helper function.
bc9aa534bc29 maint: Clean up BIST tests so they don't produce warnings. Rik <rik@octave.org> parents: 21301 diff changeset	702 %!assert (1)
bc9aa534bc29 maint: Clean up BIST tests so they don't produce warnings. Rik <rik@octave.org> parents: 21301 diff changeset	703 */
29958 32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace John W. Eaton <jwe@octave.org> parents: 29481 diff changeset	704
32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace John W. Eaton <jwe@octave.org> parents: 29481 diff changeset	705 OCTAVE_NAMESPACE_END

Mercurial > octave

annotate libinterp/corefcn/__eigs__.cc @ 29958:32c3a5805893

annotate libinterp/corefcn/eigs.cc @ 29958:32c3a5805893