# HG changeset patch # User jwe # Date 943067572 0 # Node ID 7ed630f7b7be0c282363c192e3638bb0c6cb0c1f # Parent fac05a83b4c5b99abaca3f2f6ab70de2844ffe37 [project @ 1999-11-20 03:11:17 by jwe] diff -r fac05a83b4c5 -r 7ed630f7b7be scripts/ChangeLog --- a/scripts/ChangeLog Sat Nov 20 03:05:45 1999 +0000 +++ b/scripts/ChangeLog Sat Nov 20 03:12:52 1999 +0000 @@ -1,5 +1,7 @@ 1999-11-19 John W. Eaton + * control/pinv.m: Delete. + * Makefile.in (DOCSTRINGS): Don't echo $(ALL_M_FILES). * strings/blanks.m: Texinfoize doc string. diff -r fac05a83b4c5 -r 7ed630f7b7be scripts/control/pinv.m --- a/scripts/control/pinv.m Sat Nov 20 03:05:45 1999 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,53 +0,0 @@ -# Copyright (C) 1994 John W. Eaton -# -# This file is part of Octave. -# -# Octave is free software; you can redistribute it and/or modify it -# under the terms of the GNU General Public License as published by the -# Free Software Foundation; either version 2, or (at your option) any -# later version. -# -# Octave is distributed in the hope that it will be useful, but WITHOUT -# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or -# FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License -# for more details. -# -# You should have received a copy of the GNU General Public License -# along with Octave; see the file COPYING. If not, write to the Free -# Software Foundation, 59 Temple Place, Suite 330, Boston, MA 02111 USA. - -## -*- texinfo -*- -## @deftypefn {Function File } { } pinv ( @var{X}@{,@var{tol}@} ) -## Returns the pseudoinverse of X; singular values less than tol are ignored. -## -## If the second arguement is ommited , it is assummed that -## @example -## tol = max (size (X)) * sigma_max (X) * eps, -## @end example -## where sigma_max(X) is the maximal singular value of X. -## @end deftypefn - -function retval = pinv (X, tol) -# Written by Kurt Hornik (hornik@neuro.tuwien.ac.at) March 1993. -# Dept of Probability Theory and Statistics TU Wien, Austria. - - if (nargin < 1 || nargin > 2) - error ("usage: pinv (X [, tol])"); - endif - - [U, S, V] = svd(X); - s = diag(S); - - if (nargin == 1) - tol = max (size (X)) * s (1) * eps; - endif - - r = sum (s > tol); - if (r == 0) - retval = zeros (X'); - else - D = diag (ones (r, 1) ./ s (1:r)); - retval = V (:, 1:r) * D * U (:, 1:r)'; - endif - -endfunction diff -r fac05a83b4c5 -r 7ed630f7b7be src/ChangeLog --- a/src/ChangeLog Sat Nov 20 03:05:45 1999 +0000 +++ b/src/ChangeLog Sat Nov 20 03:12:52 1999 +0000 @@ -1,5 +1,7 @@ 1999-11-19 John W. Eaton + * DLD-FUNCTIONS/pinv.cc (Fpinv): Texinfoize doc string. + * defun-int.h (DEFUN_DLD_INTERNAL): New macro. * defun-dld.h [MAKE_BUILTINS] (DEFUN_DLD): Simply expand to DEFUN_DLD_INTERNAL. diff -r fac05a83b4c5 -r 7ed630f7b7be src/DLD-FUNCTIONS/pinv.cc --- a/src/DLD-FUNCTIONS/pinv.cc Sat Nov 20 03:05:45 1999 +0000 +++ b/src/DLD-FUNCTIONS/pinv.cc Sat Nov 20 03:12:52 1999 +0000 @@ -31,8 +31,16 @@ #include "utils.h" DEFUN_DLD (pinv, args, , - "pinv ( [, tol])\n\ -Returns the pseudoinverse of X; singular values less than tol are ignored.") + "-*- texinfo -*-\n\ +@deftypefn {Function File } { } pinv (@var{X}, @var{tol})\n\ +Returns the pseudoinverse of X; singular values less than tol are ignored.\n\ +\n\ +If the second arguement is ommited , it is assummed that\n\ +@example\n\ +tol = max (size (X)) * sigma_max (X) * eps,\n\ +@end example\n\ +where sigma_max(X) is the maximal singular value of X.\n\ +@end deftypefn") { octave_value_list retval;