trunk mex/octave/rcond.m: simplification + header update (merge needed)

git-svn-id: https://www.dynare.org/svn/dynare/trunk@2618 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
sebastien 2009-04-23 13:26:34 +00:00
parent 225a23d8a7
commit e58ef63b96
1 changed files with 27 additions and 5 deletions

View File

@ -1,5 +1,27 @@
function r = rcond(A)
% Computes reciprocal condition number
% part of DYNARE, copyright Dynare Team (2008)
% Gnu Public License.
r = 1/(norm(A,1) * norm(inv(A), 1));
## Copyright (C) 2008-2009 Dynare Team
##
## This file is part of Dynare.
##
## Dynare is free software: you can redistribute it and/or modify
## it under the terms of the GNU General Public License as published by
## the Free Software Foundation, either version 3 of the License, or
## (at your option) any later version.
##
## Dynare 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 Dynare. If not, see <http://www.gnu.org/licenses/>.
## -*- texinfo -*-
## @deftypefn {Function File} {@var{rcn} =} rcond (@var{A})
##
## Computes an estimate for the reciprocal of the condition number
## of matrix @var{A}, using the 1-norm.
## @end deftypefn
function rcn = rcond(A)
[dummy, rcn] = inv(A);
endfunction