Mercurial > octave-nkf

comparison src/DLD-FUNCTIONS/gammainc.cc @ 11586:12df7854fa7c

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
strip trailing whitespace from source files
author John W. Eaton <jwe@octave.org>
date Thu, 20 Jan 2011 17:24:59 -0500
parents 01f703952eff
children 72c96de7a403
comparison
equal deleted inserted replaced
11585:1473d0cf86d2 11586:12df7854fa7c
113 if (a_arg.is_scalar_type ()) 113 if (a_arg.is_scalar_type ())
114 { 114 {
115 float a = a_arg.float_value (); 115 float a = a_arg.float_value ();
116 116
117 if (! error_state) 117 if (! error_state)
118 retval = lower ? gammainc (x, a) 118 retval = lower ? gammainc (x, a)
119 : static_cast<float>(1) - gammainc (x, a); 119 : static_cast<float>(1) - gammainc (x, a);
120 } 120 }
121 else 121 else
122 { 122 {
123 FloatNDArray a = a_arg.float_array_value (); 123 FloatNDArray a = a_arg.float_array_value ();
124 124
125 if (! error_state) 125 if (! error_state)
126 retval = lower ? gammainc (x, a) 126 retval = lower ? gammainc (x, a)
127 : static_cast<float>(1) - gammainc (x, a); 127 : static_cast<float>(1) - gammainc (x, a);
128 } 128 }
129 } 129 }
130 } 130 }
131 else 131 else
137 if (a_arg.is_scalar_type ()) 137 if (a_arg.is_scalar_type ())
138 { 138 {
139 float a = a_arg.float_value (); 139 float a = a_arg.float_value ();
140 140
141 if (! error_state) 141 if (! error_state)
142 retval = lower ? gammainc (x, a) 142 retval = lower ? gammainc (x, a)
143 : static_cast<float>(1) - gammainc (x, a); 143 : static_cast<float>(1) - gammainc (x, a);
144 } 144 }
145 else 145 else
146 { 146 {
147 FloatNDArray a = a_arg.float_array_value (); 147 FloatNDArray a = a_arg.float_array_value ();
148 148
149 if (! error_state) 149 if (! error_state)
150 retval = lower ? gammainc (x, a) 150 retval = lower ? gammainc (x, a)
151 : static_cast<float>(1) - gammainc (x, a); 151 : static_cast<float>(1) - gammainc (x, a);
152 } 152 }
153 } 153 }
154 } 154 }
155 } 155 }
214 %! x = [0 1 2 3 4]; 214 %! x = [0 1 2 3 4];
215 %! v1 = sqrt(pi)*erf(x)./gamma(a); 215 %! v1 = sqrt(pi)*erf(x)./gamma(a);
216 %! v3 = gammainc(x.*x,a); 216 %! v3 = gammainc(x.*x,a);
217 %! assert(v1, v3, sqrt(eps)); 217 %! assert(v1, v3, sqrt(eps));
218 218
219 %!assert (gammainc(0:4,0.5,"upper"), 1-gammainc(0:4,0.5),1e-10) 219 %!assert (gammainc(0:4,0.5,"upper"), 1-gammainc(0:4,0.5),1e-10)
220 220
221 %!test 221 %!test
222 %! a = single ([.5 .5 .5 .5 .5]); 222 %! a = single ([.5 .5 .5 .5 .5]);
223 %! x = single([0 1 2 3 4]); 223 %! x = single([0 1 2 3 4]);
224 %! v1 = sqrt(pi('single'))*erf(x)./gamma(a); 224 %! v1 = sqrt(pi('single'))*erf(x)./gamma(a);
225 %! v3 = gammainc(x.*x,a); 225 %! v3 = gammainc(x.*x,a);
226 %! assert(v1, v3, sqrt(eps('single'))); 226 %! assert(v1, v3, sqrt(eps('single')));
227 227
228 %!assert (gammainc(single(0:4),single(0.5),"upper"), single(1)-gammainc(single(0:4),single(0.5)),single(1e-7)) 228 %!assert (gammainc(single(0:4),single(0.5),"upper"), single(1)-gammainc(single(0:4),single(0.5)),single(1e-7))
229 229
230 */ 230 */