Mercurial > forge
changeset 2374:83bd96c90380 octave-forge
Changed the directory structure of special-matrix to match the package system
author | hauberg |
---|---|
date | Sun, 20 Aug 2006 14:06:10 +0000 |
parents | ee27f189687e |
children | 3b5a47384193 |
files | main/special-matrix/COPYING main/special-matrix/DESCRIPTION main/special-matrix/inst/lauchli.m main/special-matrix/lauchli.m |
diffstat | 4 files changed, 36 insertions(+), 24 deletions(-) [+] |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/main/special-matrix/COPYING Sun Aug 20 14:06:10 2006 +0000 @@ -0,0 +1,1 @@ +This program is donated to the public domain
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/main/special-matrix/DESCRIPTION Sun Aug 20 14:06:10 2006 +0000 @@ -0,0 +1,11 @@ +Name: special-matrix +Version: 1.0.0 +Date: 2006-08-05 +Author: Paul Kienzle <pkienzle@users.sf.net> +Maintainer: Paul Kienzle <pkienzle@users.sf.net> +Title: Special Matrices. +Description: Add a description to this package! +Categories: Special Matrices +Depends: octave (>= 2.9.7) +License: GPL version 2 or later +Url: http://octave.sf.net
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/main/special-matrix/inst/lauchli.m Sun Aug 20 14:06:10 2006 +0000 @@ -0,0 +1,24 @@ +## A = lauchli (n [,mu]) +## Creates the matrix [ ones(1,n); mu*eye(n) ] +## The value mu defaults to sqrt(eps) +## This is an ill-conditioned system for testing the +## accuracy of the QR routine. +## E.g., +## A = lauchli(15); +## [Q, R] = qr(A); +## norm(Q*R - A) +## norm(Q'*Q - eye(rows(Q))) + +## This program is in the public domain +## Author: Paul Kienzle <pkienzle@users.sf.net> + +function A = lauchli(n,mu) + if (nargin < 1 || nargin > 2) + usage("A = lauchli(n [, mu])"); + endif + + if (nargin < 2), mu = sqrt(eps); endif + + A = [ ones(1,n); mu*eye(n) ]; + +endfunction
--- a/main/special-matrix/lauchli.m Sun Aug 20 14:04:05 2006 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,24 +0,0 @@ -## A = lauchli (n [,mu]) -## Creates the matrix [ ones(1,n); mu*eye(n) ] -## The value mu defaults to sqrt(eps) -## This is an ill-conditioned system for testing the -## accuracy of the QR routine. -## E.g., -## A = lauchli(15); -## [Q, R] = qr(A); -## norm(Q*R - A) -## norm(Q'*Q - eye(rows(Q))) - -## This program is in the public domain -## Author: Paul Kienzle <pkienzle@users.sf.net> - -function A = lauchli(n,mu) - if (nargin < 1 || nargin > 2) - usage("A = lauchli(n [, mu])"); - endif - - if (nargin < 2), mu = sqrt(eps); endif - - A = [ ones(1,n); mu*eye(n) ]; - -endfunction