Mercurial > octave-nkf
annotate scripts/specfun/factor.m @ 20521:4bb41929286b
Deprecate bitmax.
* NEWS: Announce deprecation and replacement functions.
* scripts/deprecated/bitmax.m: New function. Add deprecation message to
docstring. Add deprecation warning to m-file.
* scripts/deprecated/module.mk: Add deprecated versions to build system.
* libinterp/corefcn/bitfcns.cc: Remove build-in function. Remove function
from docstrings.
* libinterp/corefcn/data.cc: Remove function from docstrings.
* scripts/general/bitcmp.m: Replace bitmax function calls with flintmax.
Remove function from docstrings.
* scripts/general/bitget.m: Replace bitmax function calls with flintmax.
Remove function from docstrings.
* scripts/general/bitset.m: Replace bitmax function calls with flintmax.
Remove function from docstrings.
* scripts/specfun/factor.m: Replace bitmax function calls with flintmax.
Remove function from docstrings.
* scripts/statistics/distributions/unidpdf.m: Remove function from docstrings.
* scripts/general/num2str.m: Replaced function in comment.
* doc/interpreter/doccheck/aspell-octave.en.pws: Remove functions from list.
* doc/interpreter/numbers.txi: Remove function from manual.
| author | Kai T. Ohlhus <k.ohlhus@gmail.com> |
|---|---|
| date | Wed, 09 Sep 2015 18:43:52 +0200 |
| parents | 9fc020886ae9 |
| children |
| rev | line source |
|---|---|
|
19731
4197fc428c7d
maint: Update copyright notices for 2015.
John W. Eaton <jwe@octave.org>
parents:
19630
diff
changeset
|
1 ## Copyright (C) 2000-2015 Paul Kienzle |
| 5827 | 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. | |
| 5827 | 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/>. | |
| 5827 | 18 |
| 19 ## -*- texinfo -*- | |
| 19131 | 20 ## @deftypefn {Function File} {@var{pf} =} factor (@var{q}) |
| 21 ## @deftypefnx {Function File} {[@var{pf}, @var{n}] =} factor (@var{q}) | |
| 22 ## Return the prime factorization of @var{q}. | |
| 5827 | 23 ## |
| 19131 | 24 ## The prime factorization is defined as @code{prod (@var{pf}) == @var{q}} |
| 25 ## where every element of @var{pf} is a prime number. If @code{@var{q} == 1}, | |
| 26 ## return 1. | |
|
19627
446c46af4b42
strip trailing whitespace from most source files
John W. Eaton <jwe@octave.org>
parents:
17744
diff
changeset
|
27 ## |
| 19131 | 28 ## With two output arguments, return the unique prime factors @var{pf} and |
| 29 ## their multiplicities. That is, @code{prod (@var{pf} .^ @var{n}) == @var{q}}. | |
|
19630
0e1f5a750d00
maint: Periodic merge of gui-release to default.
John W. Eaton <jwe@octave.org>
diff
changeset
|
30 ## |
|
19186
0f9c5a15c8fa
doc: Periodic grammarcheck of documentation.
Rik <rik@octave.org>
parents:
19131
diff
changeset
|
31 ## Implementation Note: The input @var{q} must be less than |
| 20521 | 32 ## @code{flintmax} (9.0072e+15) in order to factor correctly. |
| 19131 | 33 ## @seealso{gcd, lcm, isprime, primes} |
| 5827 | 34 ## @end deftypefn |
| 35 | |
| 36 ## Author: Paul Kienzle | |
| 37 | |
| 38 ## 2002-01-28 Paul Kienzle | |
| 39 ## * remove recursion; only check existing primes for multiplicity > 1 | |
| 40 ## * return multiplicity as suggested by Dirk Laurie | |
| 41 ## * add error handling | |
| 42 | |
| 19131 | 43 function [pf, n] = factor (q) |
| 5827 | 44 |
| 19131 | 45 if (nargin != 1) |
| 5827 | 46 print_usage (); |
| 47 endif | |
| 48 | |
| 19131 | 49 if (! isreal (q) || ! isscalar (q) || q != fix (q)) |
| 50 error ("factor: Q must be a real integer"); | |
| 5827 | 51 endif |
| 52 | |
|
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
53 ## Special case of no primes less than sqrt(q). |
|
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
54 if (q < 4) |
| 19131 | 55 pf = q; |
|
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
56 n = 1; |
| 5827 | 57 return; |
|
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
58 endif |
| 5827 | 59 |
|
17716
e48f5a52e838
factor.m: Warn when the input is too large to calculate reliable answer.
Doug Stewart <doug.dastew@gmail.com>
parents:
14868
diff
changeset
|
60 q = double (q); # For the time being, calcs rely on double precision var. |
|
e48f5a52e838
factor.m: Warn when the input is too large to calculate reliable answer.
Doug Stewart <doug.dastew@gmail.com>
parents:
14868
diff
changeset
|
61 qorig = q; |
| 19131 | 62 pf = []; |
|
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
63 ## There is at most one prime greater than sqrt(q), and if it exists, |
| 5827 | 64 ## it has multiplicity 1, so no need to consider any factors greater |
| 19131 | 65 ## than sqrt(q) directly. [If there were two factors p1, p2 > sqrt(q), |
|
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
66 ## then q >= p1*p2 > sqrt(q)*sqrt(q) == q. Contradiction.] |
|
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
67 p = primes (sqrt (q)); |
|
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
68 while (q > 1) |
|
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
69 ## Find prime factors in remaining q. |
|
17716
e48f5a52e838
factor.m: Warn when the input is too large to calculate reliable answer.
Doug Stewart <doug.dastew@gmail.com>
parents:
14868
diff
changeset
|
70 p = p(rem (q, p) == 0); |
| 5827 | 71 if (isempty (p)) |
|
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
72 ## Can't be reduced further, so q must itself be a prime. |
|
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
73 p = q; |
| 5827 | 74 endif |
| 19131 | 75 pf = [pf, p]; |
|
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
76 ## Reduce q. |
|
17716
e48f5a52e838
factor.m: Warn when the input is too large to calculate reliable answer.
Doug Stewart <doug.dastew@gmail.com>
parents:
14868
diff
changeset
|
77 q /= prod (p); |
| 5827 | 78 endwhile |
| 19131 | 79 pf = sort (pf); |
| 5827 | 80 |
|
17716
e48f5a52e838
factor.m: Warn when the input is too large to calculate reliable answer.
Doug Stewart <doug.dastew@gmail.com>
parents:
14868
diff
changeset
|
81 ## Verify algorithm was succesful |
| 19131 | 82 q = prod (pf); |
|
17716
e48f5a52e838
factor.m: Warn when the input is too large to calculate reliable answer.
Doug Stewart <doug.dastew@gmail.com>
parents:
14868
diff
changeset
|
83 if (q != qorig) |
| 19131 | 84 error ("factor: Q too large to factor"); |
| 20521 | 85 elseif (q >= flintmax ()) |
| 19131 | 86 warning ("factor: Q too large. Answer is unreliable"); |
|
17716
e48f5a52e838
factor.m: Warn when the input is too large to calculate reliable answer.
Doug Stewart <doug.dastew@gmail.com>
parents:
14868
diff
changeset
|
87 endif |
|
e48f5a52e838
factor.m: Warn when the input is too large to calculate reliable answer.
Doug Stewart <doug.dastew@gmail.com>
parents:
14868
diff
changeset
|
88 |
| 8506 | 89 ## Determine muliplicity. |
| 5827 | 90 if (nargout > 1) |
| 19131 | 91 idx = find ([0, pf] != [pf, 0]); |
| 92 pf = pf(idx(1:length (idx)-1)); | |
|
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
93 n = diff (idx); |
| 5827 | 94 endif |
| 95 | |
| 96 endfunction | |
| 97 | |
|
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
98 |
|
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
99 %!assert (factor (1), 1) |
|
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
100 %!test |
|
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
101 %! for i = 2:20 |
| 19131 | 102 %! pf = factor (i); |
| 103 %! assert (prod (pf), i); | |
| 104 %! assert (all (isprime (pf))); | |
| 105 %! [pf, n] = factor (i); | |
| 106 %! assert (prod (pf.^n), i); | |
| 107 %! assert (all ([0,pf] != [pf,0])); | |
|
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
108 %! endfor |
|
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10476
diff
changeset
|
109 |
|
19867
9fc020886ae9
maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents:
19731
diff
changeset
|
110 ## Test input validation |
|
17716
e48f5a52e838
factor.m: Warn when the input is too large to calculate reliable answer.
Doug Stewart <doug.dastew@gmail.com>
parents:
14868
diff
changeset
|
111 %!error factor () |
| 19131 | 112 %!error factor (1,2) |
| 113 %!error <Q must be a real integer> factor (6i) | |
| 114 %!error <Q must be a real integer> factor ([1,2]) | |
| 115 %!error <Q must be a real integer> factor (1.5) | |
|
17716
e48f5a52e838
factor.m: Warn when the input is too large to calculate reliable answer.
Doug Stewart <doug.dastew@gmail.com>
parents:
14868
diff
changeset
|
116 |