Mercurial > hg > octave-nkf
annotate scripts/general/prepad.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 | 7503499a252b |
children |
rev | line source |
---|---|
19898
4197fc428c7d
maint: Update copyright notices for 2015.
John W. Eaton <jwe@octave.org>
parents:
19855
diff
changeset
|
1 ## Copyright (C) 1994-2015 John W. Eaton |
2313 | 2 ## |
3 ## This file is part of Octave. | |
4 ## | |
5 ## Octave is free software; you can redistribute it and/or modify it | |
6 ## under the terms of the GNU General Public License as published by | |
7016 | 7 ## the Free Software Foundation; either version 3 of the License, or (at |
8 ## your option) any later version. | |
2313 | 9 ## |
10 ## Octave is distributed in the hope that it will be useful, but | |
11 ## WITHOUT ANY WARRANTY; without even the implied warranty of | |
12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
13 ## General Public License for more details. | |
14 ## | |
15 ## You should have received a copy of the GNU General Public License | |
7016 | 16 ## along with Octave; see the file COPYING. If not, see |
17 ## <http://www.gnu.org/licenses/>. | |
2303 | 18 |
3428 | 19 ## -*- texinfo -*- |
11539
6bac61388876
Add undocumented postpad function to documentation.
Rik <octave@nomad.inbox5.com>
parents:
11523
diff
changeset
|
20 ## @deftypefn {Function File} {} prepad (@var{x}, @var{l}) |
6bac61388876
Add undocumented postpad function to documentation.
Rik <octave@nomad.inbox5.com>
parents:
11523
diff
changeset
|
21 ## @deftypefnx {Function File} {} prepad (@var{x}, @var{l}, @var{c}) |
8545
faccdb98d953
postpad.m, prepad.m: doc fix
John W. Eaton <jwe@octave.org>
parents:
8436
diff
changeset
|
22 ## @deftypefnx {Function File} {} prepad (@var{x}, @var{l}, @var{c}, @var{dim}) |
12639
4d777e05d47c
doc: Review and update documentation for "Matrix Manipulation" chapter.
Rik <octave@nomad.inbox5.com>
parents:
11539
diff
changeset
|
23 ## Prepend the scalar value @var{c} to the vector @var{x} until it is of length |
4d777e05d47c
doc: Review and update documentation for "Matrix Manipulation" chapter.
Rik <octave@nomad.inbox5.com>
parents:
11539
diff
changeset
|
24 ## @var{l}. If @var{c} is not given, a value of 0 is used. |
2311 | 25 ## |
20368
7503499a252b
doc: Update docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
20038
diff
changeset
|
26 ## If @code{length (@var{x}) > @var{l}}, elements from the beginning of @var{x} |
7503499a252b
doc: Update docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
20038
diff
changeset
|
27 ## are removed until a vector of length @var{l} is obtained. |
3652 | 28 ## |
29 ## If @var{x} is a matrix, elements are prepended or removed from each row. | |
4870 | 30 ## |
20368
7503499a252b
doc: Update docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
20038
diff
changeset
|
31 ## If the optional argument @var{dim} is given, operate along this dimension. |
18581
d28c4c4547ef
Allow postpad and prepad to expand singleton dimensions (bug #41633)
Mike Miller <mtmiller@ieee.org>
parents:
17744
diff
changeset
|
32 ## |
20368
7503499a252b
doc: Update docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
20038
diff
changeset
|
33 ## If @var{dim} is larger than the dimensions of @var{x}, the result will have |
7503499a252b
doc: Update docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
20038
diff
changeset
|
34 ## @var{dim} dimensions. |
11539
6bac61388876
Add undocumented postpad function to documentation.
Rik <octave@nomad.inbox5.com>
parents:
11523
diff
changeset
|
35 ## @seealso{postpad, cat, resize} |
3428 | 36 ## @end deftypefn |
37 | |
38 ## Author: Tony Richardson <arichard@stark.cc.oh.us> | |
39 ## Created: June 1994 | |
1337 | 40 |
4862 | 41 function y = prepad (x, l, c, dim) |
42 | |
43 if (nargin < 2 || nargin > 4) | |
6046 | 44 print_usage (); |
4862 | 45 endif |
46 | |
47 if (nargin < 3 || isempty (c)) | |
48 c = 0; | |
49 else | |
50 if (! isscalar (c)) | |
14406
98cf8aa9f8ba
prepad.m: Improve input validation messages and tests.
Rik <octave@nomad.inbox5.com>
parents:
14363
diff
changeset
|
51 error ("prepad: pad value C must be empty or a scalar"); |
4862 | 52 endif |
53 endif | |
559 | 54 |
4862 | 55 nd = ndims (x); |
56 sz = size (x); | |
57 if (nargin < 4) | |
12674
9493880928c8
Use common idiom in m-files for finding first non-singleton dimension.
Rik <octave@nomad.inbox5.com>
parents:
12639
diff
changeset
|
58 ## Find the first non-singleton dimension. |
9493880928c8
Use common idiom in m-files for finding first non-singleton dimension.
Rik <octave@nomad.inbox5.com>
parents:
12639
diff
changeset
|
59 (dim = find (sz > 1, 1)) || (dim = 1); |
4862 | 60 else |
18581
d28c4c4547ef
Allow postpad and prepad to expand singleton dimensions (bug #41633)
Mike Miller <mtmiller@ieee.org>
parents:
17744
diff
changeset
|
61 if (!(isscalar (dim) && dim == fix (dim) && dim >= 1)) |
10690
35adf2a71f3f
Use common code block to find first non-singleton dimension.
Rik <octave@nomad.inbox5.com>
parents:
10635
diff
changeset
|
62 error ("prepad: DIM must be an integer and a valid dimension"); |
4862 | 63 endif |
559 | 64 endif |
65 | |
5459 | 66 if (! isscalar (l) || l < 0) |
14406
98cf8aa9f8ba
prepad.m: Improve input validation messages and tests.
Rik <octave@nomad.inbox5.com>
parents:
14363
diff
changeset
|
67 error ("prepad: length L must be a positive scalar"); |
559 | 68 endif |
69 | |
5459 | 70 if (dim > nd) |
71 sz(nd+1:dim) = 1; | |
72 endif | |
73 | |
20038
9fc020886ae9
maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents:
19898
diff
changeset
|
74 d = sz(dim); |
4862 | 75 |
76 if (d >= l) | |
12676
2783fa95cab7
Use common code idiom for creating cell array for indexing ND-arrays
Rik <octave@nomad.inbox5.com>
parents:
12674
diff
changeset
|
77 idx = repmat ({':'}, nd, 1); |
7208 | 78 idx{dim} = d-l+1:d; |
79 y = x(idx{:}); | |
559 | 80 else |
20038
9fc020886ae9
maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents:
19898
diff
changeset
|
81 sz(dim) = l - d; |
19855
107ec5195a47
prepad.m: Fix str-to-num warnings (bug #44162)
Mike Miller <mtmiller@ieee.org>
parents:
19047
diff
changeset
|
82 y = cat (dim, c(ones (sz)), x); |
559 | 83 endif |
84 | |
85 endfunction | |
12795
9e7ebbaf69ff
codesprint: new tests for files in scripts/general directory
John W. Eaton <jwe@octave.org>
parents:
12676
diff
changeset
|
86 |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
87 |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
88 %!assert (prepad ([1,2], 4), [0,0,1,2]) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
89 %!assert (prepad ([1;2], 4), [0;0;1;2]) |
12795
9e7ebbaf69ff
codesprint: new tests for files in scripts/general directory
John W. Eaton <jwe@octave.org>
parents:
12676
diff
changeset
|
90 |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
91 %!assert (prepad ([1,2], 4, 2), [2,2,1,2]) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
92 %!assert (prepad ([1;2], 4, 2), [2;2;1;2]) |
12795
9e7ebbaf69ff
codesprint: new tests for files in scripts/general directory
John W. Eaton <jwe@octave.org>
parents:
12676
diff
changeset
|
93 |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
94 %!assert (prepad ([1,2], 2, 2, 1), [2,2;1,2]) |
12795
9e7ebbaf69ff
codesprint: new tests for files in scripts/general directory
John W. Eaton <jwe@octave.org>
parents:
12676
diff
changeset
|
95 |
18581
d28c4c4547ef
Allow postpad and prepad to expand singleton dimensions (bug #41633)
Mike Miller <mtmiller@ieee.org>
parents:
17744
diff
changeset
|
96 %!assert (prepad ([1,2], 2, 2, 3), reshape ([2,2,1,2], 1, 2, 2)) |
d28c4c4547ef
Allow postpad and prepad to expand singleton dimensions (bug #41633)
Mike Miller <mtmiller@ieee.org>
parents:
17744
diff
changeset
|
97 %!assert (prepad ([1;2], 2, 2, 3), reshape ([2;2;1;2], 2, 1, 2)) |
d28c4c4547ef
Allow postpad and prepad to expand singleton dimensions (bug #41633)
Mike Miller <mtmiller@ieee.org>
parents:
17744
diff
changeset
|
98 |
19855
107ec5195a47
prepad.m: Fix str-to-num warnings (bug #44162)
Mike Miller <mtmiller@ieee.org>
parents:
19047
diff
changeset
|
99 %! ## Test with string concatenation (bug #44162) |
107ec5195a47
prepad.m: Fix str-to-num warnings (bug #44162)
Mike Miller <mtmiller@ieee.org>
parents:
19047
diff
changeset
|
100 %!assert (prepad ("Octave", 16, "x"), "xxxxxxxxxxOctave") |
107ec5195a47
prepad.m: Fix str-to-num warnings (bug #44162)
Mike Miller <mtmiller@ieee.org>
parents:
19047
diff
changeset
|
101 %!assert (prepad ("Octave", 4), "tave") |
107ec5195a47
prepad.m: Fix str-to-num warnings (bug #44162)
Mike Miller <mtmiller@ieee.org>
parents:
19047
diff
changeset
|
102 |
19047
7bbe3658c5ef
maint: Use "FIXME:" coding convention in m-files.
Rik <rik@octave.org>
parents:
18581
diff
changeset
|
103 ## FIXME: We need tests for multidimensional arrays. |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
104 |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
105 %!error prepad () |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
106 %!error prepad (1) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
107 %!error prepad (1,2,3,4,5) |
14406
98cf8aa9f8ba
prepad.m: Improve input validation messages and tests.
Rik <octave@nomad.inbox5.com>
parents:
14363
diff
changeset
|
108 %!error <C must be empty or a scalar> prepad ([1,2], 2, ones (2)) |
98cf8aa9f8ba
prepad.m: Improve input validation messages and tests.
Rik <octave@nomad.inbox5.com>
parents:
14363
diff
changeset
|
109 %!error <DIM must be an integer> prepad ([1,2], 2, 2, ones (3)) |
98cf8aa9f8ba
prepad.m: Improve input validation messages and tests.
Rik <octave@nomad.inbox5.com>
parents:
14363
diff
changeset
|
110 %!error <DIM must be an integer> prepad ([1,2], 2, 2, 1.1) |
98cf8aa9f8ba
prepad.m: Improve input validation messages and tests.
Rik <octave@nomad.inbox5.com>
parents:
14363
diff
changeset
|
111 %!error <L must be a positive scalar> prepad ([1,2], ones (2)) |
98cf8aa9f8ba
prepad.m: Improve input validation messages and tests.
Rik <octave@nomad.inbox5.com>
parents:
14363
diff
changeset
|
112 %!error <L must be a positive scalar> prepad ([1,2], -1) |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
113 |