Mercurial > octave-nkf
annotate libinterp/corefcn/gammainc.cc @ 20616:fd0efcdb3718
use new string_value method to handle value extraction errors
* dirfns.cc, file-io.cc, gammainc.cc, help.cc, load-path.cc,
octave-link.cc, qz.cc, regexp.cc, strfns.cc, syscalls.cc, time.cc,
variables.cc: Use new string_value method.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Thu, 08 Oct 2015 19:00:51 -0400 |
parents | b10432a40432 |
children |
rev | line source |
---|---|
3254 | 1 /* |
2 | |
19731
4197fc428c7d
maint: Update copyright notices for 2015.
John W. Eaton <jwe@octave.org>
parents:
19185
diff
changeset
|
3 Copyright (C) 1997-2015 John W. Eaton |
3254 | 4 |
5 This file is part of Octave. | |
6 | |
7 Octave is free software; you can redistribute it and/or modify it | |
8 under the terms of the GNU General Public License as published by the | |
7016 | 9 Free Software Foundation; either version 3 of the License, or (at your |
10 option) any later version. | |
3254 | 11 |
12 Octave is distributed in the hope that it will be useful, but WITHOUT | |
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
15 for more details. | |
16 | |
17 You should have received a copy of the GNU General Public License | |
7016 | 18 along with Octave; see the file COPYING. If not, see |
19 <http://www.gnu.org/licenses/>. | |
3254 | 20 |
21 */ | |
22 | |
23 #ifdef HAVE_CONFIG_H | |
24 #include <config.h> | |
25 #endif | |
26 | |
27 #include "lo-specfun.h" | |
28 | |
15039
e753177cde93
maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents:
14501
diff
changeset
|
29 #include "defun.h" |
3254 | 30 #include "error.h" |
31 #include "gripes.h" | |
32 #include "oct-obj.h" | |
33 #include "utils.h" | |
34 | |
15039
e753177cde93
maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents:
14501
diff
changeset
|
35 DEFUN (gammainc, args, , |
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
36 "-*- texinfo -*-\n\ |
10840 | 37 @deftypefn {Mapping Function} {} gammainc (@var{x}, @var{a})\n\ |
10698
fa00ccf7b1f9
Implementary complementary incomplete gamma function (bug #30088)
David Bateman <dbateman@free.fr>
parents:
10155
diff
changeset
|
38 @deftypefnx {Mapping Function} {} gammainc (@var{x}, @var{a}, \"lower\")\n\ |
fa00ccf7b1f9
Implementary complementary incomplete gamma function (bug #30088)
David Bateman <dbateman@free.fr>
parents:
10155
diff
changeset
|
39 @deftypefnx {Mapping Function} {} gammainc (@var{x}, @var{a}, \"upper\")\n\ |
20207
4f45eaf83908
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19773
diff
changeset
|
40 Compute the normalized incomplete gamma function.\n\ |
4f45eaf83908
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19773
diff
changeset
|
41 \n\ |
4f45eaf83908
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19773
diff
changeset
|
42 This is defined as\n\ |
3444 | 43 @tex\n\ |
44 $$\n\ | |
11415
e7ed20f87e82
Add @tex blocks to gammainc and legendre docstrings.
Michael Godfrey <godfrey@isl.stanford.edu>
parents:
10840
diff
changeset
|
45 \\gamma (x, a) = {1 \\over {\\Gamma (a)}}\\displaystyle{\\int_0^x t^{a-1} e^{-t} dt}\n\ |
3444 | 46 $$\n\ |
47 @end tex\n\ | |
8517
81d6ab3ac93c
Allow documentation tobe built for other formats than tex and info
sh@sh-laptop
parents:
7815
diff
changeset
|
48 @ifnottex\n\ |
3254 | 49 \n\ |
11415
e7ed20f87e82
Add @tex blocks to gammainc and legendre docstrings.
Michael Godfrey <godfrey@isl.stanford.edu>
parents:
10840
diff
changeset
|
50 @example\n\ |
10840 | 51 @group\n\ |
14437
62a2ef6c09a4
doc: Correct spacing of integrals in Info version of docstrings.
Rik <octave@nomad.inbox5.com>
parents:
14360
diff
changeset
|
52 x\n\ |
62a2ef6c09a4
doc: Correct spacing of integrals in Info version of docstrings.
Rik <octave@nomad.inbox5.com>
parents:
14360
diff
changeset
|
53 1 /\n\ |
3683 | 54 gammainc (x, a) = --------- | exp (-t) t^(a-1) dt\n\ |
14360
97883071e8e4
doc: Correct off-by-1 spacings in all .cc docstrings
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
55 gamma (a) /\n\ |
97883071e8e4
doc: Correct off-by-1 spacings in all .cc docstrings
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
56 t=0\n\ |
10840 | 57 @end group\n\ |
11415
e7ed20f87e82
Add @tex blocks to gammainc and legendre docstrings.
Michael Godfrey <godfrey@isl.stanford.edu>
parents:
10840
diff
changeset
|
58 @end example\n\ |
5517 | 59 \n\ |
8517
81d6ab3ac93c
Allow documentation tobe built for other formats than tex and info
sh@sh-laptop
parents:
7815
diff
changeset
|
60 @end ifnottex\n\ |
5517 | 61 with the limiting value of 1 as @var{x} approaches infinity.\n\ |
19074
0850b5212619
doc: Add @nospell macro around proper names in documentation.
Rik <rik@octave.org>
parents:
17861
diff
changeset
|
62 The standard notation is @math{P(a,x)}, e.g., @nospell{Abramowitz} and\n\ |
0850b5212619
doc: Add @nospell macro around proper names in documentation.
Rik <rik@octave.org>
parents:
17861
diff
changeset
|
63 @nospell{Stegun} (6.5.1).\n\ |
3254 | 64 \n\ |
3683 | 65 If @var{a} is scalar, then @code{gammainc (@var{x}, @var{a})} is returned\n\ |
3444 | 66 for each element of @var{x} and vice versa.\n\ |
67 \n\ | |
3683 | 68 If neither @var{x} nor @var{a} is scalar, the sizes of @var{x} and\n\ |
11415
e7ed20f87e82
Add @tex blocks to gammainc and legendre docstrings.
Michael Godfrey <godfrey@isl.stanford.edu>
parents:
10840
diff
changeset
|
69 @var{a} must agree, and @code{gammainc} is applied element-by-element.\n\ |
10698
fa00ccf7b1f9
Implementary complementary incomplete gamma function (bug #30088)
David Bateman <dbateman@free.fr>
parents:
10155
diff
changeset
|
70 \n\ |
fa00ccf7b1f9
Implementary complementary incomplete gamma function (bug #30088)
David Bateman <dbateman@free.fr>
parents:
10155
diff
changeset
|
71 By default the incomplete gamma function integrated from 0 to @var{x} is\n\ |
17281
bc924baa2c4e
doc: Add new @qcode macro for code samples which are quoted.
Rik <rik@octave.org>
parents:
15195
diff
changeset
|
72 computed. If @qcode{\"upper\"} is given then the complementary function\n\ |
bc924baa2c4e
doc: Add new @qcode macro for code samples which are quoted.
Rik <rik@octave.org>
parents:
15195
diff
changeset
|
73 integrated from @var{x} to infinity is calculated. It should be noted that\n\ |
10698
fa00ccf7b1f9
Implementary complementary incomplete gamma function (bug #30088)
David Bateman <dbateman@free.fr>
parents:
10155
diff
changeset
|
74 \n\ |
fa00ccf7b1f9
Implementary complementary incomplete gamma function (bug #30088)
David Bateman <dbateman@free.fr>
parents:
10155
diff
changeset
|
75 @example\n\ |
11415
e7ed20f87e82
Add @tex blocks to gammainc and legendre docstrings.
Michael Godfrey <godfrey@isl.stanford.edu>
parents:
10840
diff
changeset
|
76 gammainc (@var{x}, @var{a}) @equiv{} 1 - gammainc (@var{x}, @var{a}, \"upper\")\n\ |
10698
fa00ccf7b1f9
Implementary complementary incomplete gamma function (bug #30088)
David Bateman <dbateman@free.fr>
parents:
10155
diff
changeset
|
77 @end example\n\ |
19185
e09318118fd6
doc: Add note to gamma docstring about using gammaln when input is large.
Rik <rik@octave.org>
parents:
19074
diff
changeset
|
78 @seealso{gamma, gammaln}\n\ |
5642 | 79 @end deftypefn") |
3254 | 80 { |
81 octave_value retval; | |
10698
fa00ccf7b1f9
Implementary complementary incomplete gamma function (bug #30088)
David Bateman <dbateman@free.fr>
parents:
10155
diff
changeset
|
82 bool lower = true; |
3254 | 83 |
84 int nargin = args.length (); | |
85 | |
10698
fa00ccf7b1f9
Implementary complementary incomplete gamma function (bug #30088)
David Bateman <dbateman@free.fr>
parents:
10155
diff
changeset
|
86 if (nargin == 3) |
fa00ccf7b1f9
Implementary complementary incomplete gamma function (bug #30088)
David Bateman <dbateman@free.fr>
parents:
10155
diff
changeset
|
87 { |
20616
fd0efcdb3718
use new string_value method to handle value extraction errors
John W. Eaton <jwe@octave.org>
parents:
20589
diff
changeset
|
88 std::string s = args(2).string_value ("gammainc: third argument must be \"lower\" or \"upper\""); |
fd0efcdb3718
use new string_value method to handle value extraction errors
John W. Eaton <jwe@octave.org>
parents:
20589
diff
changeset
|
89 |
fd0efcdb3718
use new string_value method to handle value extraction errors
John W. Eaton <jwe@octave.org>
parents:
20589
diff
changeset
|
90 std::transform (s.begin (), s.end (), s.begin (), tolower); |
fd0efcdb3718
use new string_value method to handle value extraction errors
John W. Eaton <jwe@octave.org>
parents:
20589
diff
changeset
|
91 |
fd0efcdb3718
use new string_value method to handle value extraction errors
John W. Eaton <jwe@octave.org>
parents:
20589
diff
changeset
|
92 if (s == "upper") |
fd0efcdb3718
use new string_value method to handle value extraction errors
John W. Eaton <jwe@octave.org>
parents:
20589
diff
changeset
|
93 lower = false; |
fd0efcdb3718
use new string_value method to handle value extraction errors
John W. Eaton <jwe@octave.org>
parents:
20589
diff
changeset
|
94 else if (s != "lower") |
11553
01f703952eff
Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents:
11523
diff
changeset
|
95 error ("gammainc: third argument must be \"lower\" or \"upper\""); |
10698
fa00ccf7b1f9
Implementary complementary incomplete gamma function (bug #30088)
David Bateman <dbateman@free.fr>
parents:
10155
diff
changeset
|
96 } |
fa00ccf7b1f9
Implementary complementary incomplete gamma function (bug #30088)
David Bateman <dbateman@free.fr>
parents:
10155
diff
changeset
|
97 |
20589
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
98 if (nargin < 2 || nargin > 3) |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
99 print_usage (); |
3254 | 100 |
20589
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
101 octave_value x_arg = args(0); |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
102 octave_value a_arg = args(1); |
3254 | 103 |
20589
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
104 // FIXME: Can we make a template version of the duplicated code below |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
105 if (x_arg.is_single_type () || a_arg.is_single_type ()) |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
106 { |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
107 if (x_arg.is_scalar_type ()) |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
108 { |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
109 float x = x_arg.float_value (); |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
110 |
20589
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
111 if (a_arg.is_scalar_type ()) |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
112 { |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
113 float a = a_arg.float_value (); |
3254 | 114 |
20589
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
115 retval = lower ? gammainc (x, a) : 1.0f - gammainc (x, a); |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
116 } |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
117 else |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
118 { |
20589
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
119 FloatNDArray a = a_arg.float_array_value (); |
3254 | 120 |
20589
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
121 retval = lower ? gammainc (x, a) : 1.0f - gammainc (x, a); |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
122 } |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
123 } |
3254 | 124 else |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
125 { |
20589
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
126 FloatNDArray x = x_arg.float_array_value (); |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
127 |
20589
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
128 if (a_arg.is_scalar_type ()) |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
129 { |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
130 float a = a_arg.float_value (); |
3254 | 131 |
20589
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
132 retval = lower ? gammainc (x, a) : 1.0f - gammainc (x, a); |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
133 } |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
134 else |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
135 { |
20589
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
136 FloatNDArray a = a_arg.float_array_value (); |
3254 | 137 |
20589
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
138 retval = lower ? gammainc (x, a) : 1.0f - gammainc (x, a); |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
139 } |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
140 } |
3254 | 141 } |
142 else | |
20589
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
143 { |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
144 if (x_arg.is_scalar_type ()) |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
145 { |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
146 double x = x_arg.double_value (); |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
147 |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
148 if (a_arg.is_scalar_type ()) |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
149 { |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
150 double a = a_arg.double_value (); |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
151 |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
152 retval = lower ? gammainc (x, a) : 1. - gammainc (x, a); |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
153 } |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
154 else |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
155 { |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
156 NDArray a = a_arg.array_value (); |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
157 |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
158 retval = lower ? gammainc (x, a) : 1. - gammainc (x, a); |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
159 } |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
160 } |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
161 else |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
162 { |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
163 NDArray x = x_arg.array_value (); |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
164 |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
165 if (a_arg.is_scalar_type ()) |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
166 { |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
167 double a = a_arg.double_value (); |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
168 |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
169 retval = lower ? gammainc (x, a) : 1. - gammainc (x, a); |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
170 } |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
171 else |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
172 { |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
173 NDArray a = a_arg.array_value (); |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
174 |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
175 retval = lower ? gammainc (x, a) : 1. - gammainc (x, a); |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
176 } |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
177 } |
b10432a40432
eliminate more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents:
20207
diff
changeset
|
178 } |
3254 | 179 |
180 return retval; | |
181 } | |
182 | |
183 /* | |
7815
a41df65f3f00
Add some single precision test code and fix resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
184 %!test |
a41df65f3f00
Add some single precision test code and fix resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
185 %! a = [.5 .5 .5 .5 .5]; |
a41df65f3f00
Add some single precision test code and fix resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
186 %! x = [0 1 2 3 4]; |
14501
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14437
diff
changeset
|
187 %! v1 = sqrt (pi)*erf (x)./gamma (a); |
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14437
diff
changeset
|
188 %! v3 = gammainc (x.*x, a); |
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14437
diff
changeset
|
189 %! assert (v1, v3, sqrt (eps)); |
7815
a41df65f3f00
Add some single precision test code and fix resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
190 |
14501
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14437
diff
changeset
|
191 %!assert (gammainc (0:4,0.5, "upper"), 1-gammainc (0:4,0.5), 1e-10) |
10698
fa00ccf7b1f9
Implementary complementary incomplete gamma function (bug #30088)
David Bateman <dbateman@free.fr>
parents:
10155
diff
changeset
|
192 |
7815
a41df65f3f00
Add some single precision test code and fix resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
193 %!test |
a41df65f3f00
Add some single precision test code and fix resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
194 %! a = single ([.5 .5 .5 .5 .5]); |
14501
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14437
diff
changeset
|
195 %! x = single ([0 1 2 3 4]); |
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14437
diff
changeset
|
196 %! v1 = sqrt (pi ("single"))*erf (x)./gamma (a); |
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14437
diff
changeset
|
197 %! v3 = gammainc (x.*x, a); |
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14437
diff
changeset
|
198 %! assert (v1, v3, sqrt (eps ("single"))); |
7815
a41df65f3f00
Add some single precision test code and fix resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
199 |
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
200 %!assert (gammainc (single (0:4), single (0.5), "upper"), |
14501
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14437
diff
changeset
|
201 %! single (1)-gammainc (single (0:4), single (0.5)), |
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14437
diff
changeset
|
202 %! single (1e-7)) |
7815
a41df65f3f00
Add some single precision test code and fix resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
203 */ |