octave-nkf: scripts/special-matrix/wilkinson.m annotate

annotate scripts/special-matrix/wilkinson.m @ 20818:9d2023d1a63c

binoinv.m: Implement binary search algorithm for 28X performance increase (bug #34363). * binoinv.m: Call new functions scalar_binoinv or vector_binoinv to calculate binoinv. If there are still uncalculated values then call bin_search_binoinv to perform binary search for remaining values. Add more BIST tests. * binoinv.m (scalar_binoinv): New subfunction to calculate binoinv for scalar x. Stops when x > 1000. * binoinv.m (vector_binoinv): New subfunction to calculate binoinv for scalar x. Stops when x > 1000.

author	Lachlan Andrew <lachlanbis@gmail.com>
date	Sun, 11 Oct 2015 19:49:40 -0700
parents	2645f9ef8c88
children

rev	line source
19898 4197fc428c7d maint: Update copyright notices for 2015. John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	1 ## Copyright (C) 1999-2015 Peter Ekberg
5827 1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	2 ##
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	3 ## This file is part of Octave.
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	4 ##
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	5 ## Octave is free software; you can redistribute it and/or modify it
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	6 ## under the terms of the GNU General Public License as published by
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5827 diff changeset	7 ## the Free Software Foundation; either version 3 of the License, or (at
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5827 diff changeset	8 ## your option) any later version.
5827 1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	9 ##
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	10 ## Octave is distributed in the hope that it will be useful, but
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	13 ## General Public License for more details.
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	14 ##
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	15 ## You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5827 diff changeset	16 ## along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5827 diff changeset	17 ## <http://www.gnu.org/licenses/>.
5827 1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	18
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	19 ## -- texinfo --
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	20 ## @deftypefn {Function File} {} wilkinson (@var{n})
20372 2645f9ef8c88 doc: Update more docstrings to have one sentence summary as first line. Rik <rik@octave.org> parents: 20038 diff changeset	21 ## Return the Wilkinson matrix of order @var{n}.
2645f9ef8c88 doc: Update more docstrings to have one sentence summary as first line. Rik <rik@octave.org> parents: 20038 diff changeset	22 ##
2645f9ef8c88 doc: Update more docstrings to have one sentence summary as first line. Rik <rik@octave.org> parents: 20038 diff changeset	23 ## Wilkinson matrices are symmetric and tridiagonal with pairs of nearly, but
2645f9ef8c88 doc: Update more docstrings to have one sentence summary as first line. Rik <rik@octave.org> parents: 20038 diff changeset	24 ## not exactly, equal eigenvalues. They are useful in testing the behavior and
2645f9ef8c88 doc: Update more docstrings to have one sentence summary as first line. Rik <rik@octave.org> parents: 20038 diff changeset	25 ## performance of eigenvalue solvers.
5827 1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	26 ##
12639 4d777e05d47c doc: Review and update documentation for "Matrix Manipulation" chapter. Rik <octave@nomad.inbox5.com> parents: 12551 diff changeset	27 ## @seealso{rosser, eig}
5827 1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	28 ## @end deftypefn
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	29
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	30 ## Author: Peter Ekberg
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	31 ## (peda)
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	32
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	33 function retval = wilkinson (n)
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	34
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	35 if (nargin != 1)
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	36 print_usage ();
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	37 endif
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	38
13890 3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	39 if (! (isscalar (n) && n >= 0 && (n == fix (n))))
12551 e4a1ede4e832 Update test for 'wilkinson(1)' to reflect changes to diag(). John W. Eaton <jwe@octave.org> parents: 12448 diff changeset	40 error ("wilkinson: N must be a non-negative integer");
5827 1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	41 endif
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	42
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	43 side = ones (n-1, 1);
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	44 center = abs (-(n-1)/2:(n-1)/2);
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	45 retval = diag (side, -1) + diag (center) + diag (side, 1);
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	46
1fe78adb91bc [project @ 2006-05-22 06:25:14 by jwe] jwe parents: diff changeset	47 endfunction
11098 dcde7c5a1d29 new tests for special-matrix functions John W. Eaton <jwe@octave.org> parents: 10791 diff changeset	48
13890 3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	49
3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	50 %!assert (wilkinson (0), [])
3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	51 %!assert (wilkinson (1), 0)
3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	52 %!assert (wilkinson (2), [0.5,1;1,0.5])
3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	53 %!assert (wilkinson (3), [1,1,0;1,0,1;0,1,1])
3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	54 %!assert (wilkinson (4), [1.5,1,0,0;1,0.5,1,0;0,1,0.5,1;0,0,1,1.5])
3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	55
20038 9fc020886ae9 maint: Clean up m-files to follow Octave coding conventions. Rik <rik@octave.org> parents: 19898 diff changeset	56 ## Test input validation
13890 3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	57 %!error wilkinson ()
3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	58 %!error wilkinson (1,2)
3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	59 %!error <N must be a non-negative integer> wilkinson (ones (2))
3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	60 %!error <N must be a non-negative integer> wilkinson (-1)
3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	61 %!error <N must be a non-negative integer> wilkinson (1.5)
3a2f28c08fbd wilkinson.m: Use Octave spacing conventions in code. Add more input validation tests. Rik <octave@nomad.inbox5.com> parents: 12639 diff changeset	62

Mercurial > hg > octave-nkf

annotate scripts/special-matrix/wilkinson.m @ 20818:9d2023d1a63c