Mercurial > octave
annotate liboctave/Quad.cc @ 11518:141b3fb5cef7
style fixes
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Thu, 13 Jan 2011 16:52:30 -0500 |
parents | 12884915a8e4 |
children | fd0a3ac60b0e |
rev | line source |
---|---|
3 | 1 /* |
2 | |
7017 | 3 Copyright (C) 1993, 1994, 1995, 1996, 1997, 1998, 2000, 2002, 2003, |
8920 | 4 2004, 2005, 2006, 2007, 2008 John W. Eaton |
3 | 5 |
6 This file is part of Octave. | |
7 | |
8 Octave is free software; you can redistribute it and/or modify it | |
9 under the terms of the GNU General Public License as published by the | |
7016 | 10 Free Software Foundation; either version 3 of the License, or (at your |
11 option) any later version. | |
3 | 12 |
13 Octave is distributed in the hope that it will be useful, but WITHOUT | |
14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
16 for more details. | |
17 | |
18 You should have received a copy of the GNU General Public License | |
7016 | 19 along with Octave; see the file COPYING. If not, see |
20 <http://www.gnu.org/licenses/>. | |
3 | 21 |
22 */ | |
23 | |
238 | 24 #ifdef HAVE_CONFIG_H |
1192 | 25 #include <config.h> |
3 | 26 #endif |
27 | |
28 #include "Quad.h" | |
1847 | 29 #include "f77-fcn.h" |
2292 | 30 #include "lo-error.h" |
4180 | 31 #include "quit.h" |
3 | 32 #include "sun-utils.h" |
33 | |
34 static integrand_fcn user_fcn; | |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
35 static float_integrand_fcn float_user_fcn; |
3 | 36 |
5775 | 37 // FIXME -- would be nice to not have to have this global |
260 | 38 // variable. |
39 // Nonzero means an error occurred in the calculation of the integrand | |
40 // function, and the user wants us to quit. | |
41 int quad_integration_error = 0; | |
42 | |
5275 | 43 typedef octave_idx_type (*quad_fcn_ptr) (double*, int&, double*); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
44 typedef octave_idx_type (*quad_float_fcn_ptr) (float*, int&, float*); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
45 |
3 | 46 extern "C" |
4552 | 47 { |
48 F77_RET_T | |
49 F77_FUNC (dqagp, DQAGP) (quad_fcn_ptr, const double&, const double&, | |
11518 | 50 const octave_idx_type&, const double*, |
51 const double&, const double&, double&, | |
52 double&, octave_idx_type&, octave_idx_type&, | |
53 const octave_idx_type&, const octave_idx_type&, | |
54 octave_idx_type&, octave_idx_type*, double*); | |
3 | 55 |
4552 | 56 F77_RET_T |
11518 | 57 F77_FUNC (dqagi, DQAGI) (quad_fcn_ptr, const double&, |
58 const octave_idx_type&, const double&, | |
59 const double&, double&, double&, | |
60 octave_idx_type&, octave_idx_type&, | |
61 const octave_idx_type&, const octave_idx_type&, | |
62 octave_idx_type&, octave_idx_type*, double*); | |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
63 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
64 F77_RET_T |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
65 F77_FUNC (qagp, QAGP) (quad_float_fcn_ptr, const float&, const float&, |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
66 const octave_idx_type&, const float*, const float&, |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
67 const float&, float&, float&, octave_idx_type&, |
11518 | 68 octave_idx_type&, const octave_idx_type&, |
69 const octave_idx_type&, octave_idx_type&, | |
70 octave_idx_type*, float*); | |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
71 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
72 F77_RET_T |
11518 | 73 F77_FUNC (qagi, QAGI) (quad_float_fcn_ptr, const float&, |
74 const octave_idx_type&, const float&, | |
75 const float&, float&, float&, octave_idx_type&, | |
76 octave_idx_type&, const octave_idx_type&, | |
77 const octave_idx_type&, octave_idx_type&, | |
78 octave_idx_type*, float*); | |
4552 | 79 } |
3 | 80 |
5275 | 81 static octave_idx_type |
3136 | 82 user_function (double *x, int& ierr, double *result) |
3 | 83 { |
4180 | 84 BEGIN_INTERRUPT_WITH_EXCEPTIONS; |
85 | |
9179
5be2e6696772
use access_double and assign_double on sparc only
Carsten Clark <tantumquantum+gnuoctave@gmail.com>
parents:
9169
diff
changeset
|
86 #if defined (__sparc) && defined (__GNUC__) |
3 | 87 double xx = access_double (x); |
88 #else | |
89 double xx = *x; | |
90 #endif | |
91 | |
260 | 92 quad_integration_error = 0; |
93 | |
3136 | 94 double xresult = (*user_fcn) (xx); |
95 | |
9179
5be2e6696772
use access_double and assign_double on sparc only
Carsten Clark <tantumquantum+gnuoctave@gmail.com>
parents:
9169
diff
changeset
|
96 #if defined (__sparc) && defined (__GNUC__) |
3136 | 97 assign_double (result, xresult); |
98 #else | |
99 *result = xresult; | |
100 #endif | |
260 | 101 |
102 if (quad_integration_error) | |
1251 | 103 ierr = -1; |
260 | 104 |
4180 | 105 END_INTERRUPT_WITH_EXCEPTIONS; |
106 | |
3136 | 107 return 0; |
3 | 108 } |
109 | |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
110 static octave_idx_type |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
111 float_user_function (float *x, int& ierr, float *result) |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
112 { |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
113 BEGIN_INTERRUPT_WITH_EXCEPTIONS; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
114 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
115 quad_integration_error = 0; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
116 |
9169
b1e82cc8a9f3
eliminate broken special case for copying floats on Sun systems
Carsten Clark <tantumquantum+gnuoctave@gmail.com>
parents:
8920
diff
changeset
|
117 *result = (*float_user_fcn) (*x); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
118 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
119 if (quad_integration_error) |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
120 ierr = -1; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
121 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
122 END_INTERRUPT_WITH_EXCEPTIONS; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
123 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
124 return 0; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
125 } |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
126 |
3 | 127 double |
5275 | 128 DefQuad::do_integrate (octave_idx_type& ier, octave_idx_type& neval, double& abserr) |
3 | 129 { |
5275 | 130 octave_idx_type npts = singularities.capacity () + 2; |
3 | 131 double *points = singularities.fortran_vec (); |
132 double result = 0.0; | |
1935 | 133 |
5275 | 134 octave_idx_type leniw = 183*npts - 122; |
10350
12884915a8e4
merge MArray classes & improve Array interface
Jaroslav Hajek <highegg@gmail.com>
parents:
10314
diff
changeset
|
135 Array<octave_idx_type> iwork (leniw, 1); |
5275 | 136 octave_idx_type *piwork = iwork.fortran_vec (); |
1935 | 137 |
5275 | 138 octave_idx_type lenw = 2*leniw - npts; |
10350
12884915a8e4
merge MArray classes & improve Array interface
Jaroslav Hajek <highegg@gmail.com>
parents:
10314
diff
changeset
|
139 Array<double> work (lenw, 1); |
1935 | 140 double *pwork = work.fortran_vec (); |
141 | |
3 | 142 user_fcn = f; |
5275 | 143 octave_idx_type last; |
3 | 144 |
289 | 145 double abs_tol = absolute_tolerance (); |
146 double rel_tol = relative_tolerance (); | |
147 | |
1935 | 148 F77_XFCN (dqagp, DQAGP, (user_function, lower_limit, upper_limit, |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
149 npts, points, abs_tol, rel_tol, result, |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
150 abserr, neval, ier, leniw, lenw, last, |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
151 piwork, pwork)); |
3 | 152 |
153 return result; | |
154 } | |
155 | |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
156 float |
7924 | 157 DefQuad::do_integrate (octave_idx_type&, octave_idx_type&, float&) |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
158 { |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
159 (*current_liboctave_error_handler) ("incorrect integration function called"); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
160 return 0.0; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
161 } |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
162 |
3 | 163 double |
5275 | 164 IndefQuad::do_integrate (octave_idx_type& ier, octave_idx_type& neval, double& abserr) |
3 | 165 { |
166 double result = 0.0; | |
1935 | 167 |
5275 | 168 octave_idx_type leniw = 128; |
10350
12884915a8e4
merge MArray classes & improve Array interface
Jaroslav Hajek <highegg@gmail.com>
parents:
10314
diff
changeset
|
169 Array<octave_idx_type> iwork (leniw, 1); |
5275 | 170 octave_idx_type *piwork = iwork.fortran_vec (); |
1935 | 171 |
5275 | 172 octave_idx_type lenw = 8*leniw; |
10350
12884915a8e4
merge MArray classes & improve Array interface
Jaroslav Hajek <highegg@gmail.com>
parents:
10314
diff
changeset
|
173 Array<double> work (lenw, 1); |
1935 | 174 double *pwork = work.fortran_vec (); |
175 | |
3 | 176 user_fcn = f; |
5275 | 177 octave_idx_type last; |
3 | 178 |
5275 | 179 octave_idx_type inf; |
3 | 180 switch (type) |
181 { | |
182 case bound_to_inf: | |
183 inf = 1; | |
184 break; | |
1360 | 185 |
3 | 186 case neg_inf_to_bound: |
187 inf = -1; | |
188 break; | |
1360 | 189 |
3 | 190 case doubly_infinite: |
191 inf = 2; | |
192 break; | |
1360 | 193 |
3 | 194 default: |
195 assert (0); | |
196 break; | |
197 } | |
198 | |
289 | 199 double abs_tol = absolute_tolerance (); |
200 double rel_tol = relative_tolerance (); | |
201 | |
1935 | 202 F77_XFCN (dqagi, DQAGI, (user_function, bound, inf, abs_tol, rel_tol, |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
203 result, abserr, neval, ier, leniw, lenw, |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
204 last, piwork, pwork)); |
3 | 205 |
206 return result; | |
207 } | |
208 | |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
209 float |
7924 | 210 IndefQuad::do_integrate (octave_idx_type&, octave_idx_type&, float&) |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
211 { |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
212 (*current_liboctave_error_handler) ("incorrect integration function called"); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
213 return 0.0; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
214 } |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
215 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
216 double |
7924 | 217 FloatDefQuad::do_integrate (octave_idx_type&, octave_idx_type&, double&) |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
218 { |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
219 (*current_liboctave_error_handler) ("incorrect integration function called"); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
220 return 0.0; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
221 } |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
222 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
223 float |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
224 FloatDefQuad::do_integrate (octave_idx_type& ier, octave_idx_type& neval, float& abserr) |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
225 { |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
226 octave_idx_type npts = singularities.capacity () + 2; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
227 float *points = singularities.fortran_vec (); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
228 float result = 0.0; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
229 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
230 octave_idx_type leniw = 183*npts - 122; |
10350
12884915a8e4
merge MArray classes & improve Array interface
Jaroslav Hajek <highegg@gmail.com>
parents:
10314
diff
changeset
|
231 Array<octave_idx_type> iwork (leniw, 1); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
232 octave_idx_type *piwork = iwork.fortran_vec (); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
233 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
234 octave_idx_type lenw = 2*leniw - npts; |
10350
12884915a8e4
merge MArray classes & improve Array interface
Jaroslav Hajek <highegg@gmail.com>
parents:
10314
diff
changeset
|
235 Array<float> work (lenw, 1); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
236 float *pwork = work.fortran_vec (); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
237 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
238 float_user_fcn = ff; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
239 octave_idx_type last; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
240 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
241 float abs_tol = single_precision_absolute_tolerance (); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
242 float rel_tol = single_precision_relative_tolerance (); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
243 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
244 F77_XFCN (qagp, QAGP, (float_user_function, lower_limit, upper_limit, |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
245 npts, points, abs_tol, rel_tol, result, |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
246 abserr, neval, ier, leniw, lenw, last, |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
247 piwork, pwork)); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
248 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
249 return result; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
250 } |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
251 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
252 double |
7924 | 253 FloatIndefQuad::do_integrate (octave_idx_type&, octave_idx_type&, double&) |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
254 { |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
255 (*current_liboctave_error_handler) ("incorrect integration function called"); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
256 return 0.0; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
257 } |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
258 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
259 float |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
260 FloatIndefQuad::do_integrate (octave_idx_type& ier, octave_idx_type& neval, float& abserr) |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
261 { |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
262 float result = 0.0; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
263 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
264 octave_idx_type leniw = 128; |
10350
12884915a8e4
merge MArray classes & improve Array interface
Jaroslav Hajek <highegg@gmail.com>
parents:
10314
diff
changeset
|
265 Array<octave_idx_type> iwork (leniw, 1); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
266 octave_idx_type *piwork = iwork.fortran_vec (); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
267 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
268 octave_idx_type lenw = 8*leniw; |
10350
12884915a8e4
merge MArray classes & improve Array interface
Jaroslav Hajek <highegg@gmail.com>
parents:
10314
diff
changeset
|
269 Array<float> work (lenw, 1); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
270 float *pwork = work.fortran_vec (); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
271 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
272 float_user_fcn = ff; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
273 octave_idx_type last; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
274 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
275 octave_idx_type inf; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
276 switch (type) |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
277 { |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
278 case bound_to_inf: |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
279 inf = 1; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
280 break; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
281 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
282 case neg_inf_to_bound: |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
283 inf = -1; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
284 break; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
285 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
286 case doubly_infinite: |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
287 inf = 2; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
288 break; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
289 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
290 default: |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
291 assert (0); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
292 break; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
293 } |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
294 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
295 float abs_tol = single_precision_absolute_tolerance (); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
296 float rel_tol = single_precision_relative_tolerance (); |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
297 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
298 F77_XFCN (qagi, QAGI, (float_user_function, bound, inf, abs_tol, rel_tol, |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
299 result, abserr, neval, ier, leniw, lenw, |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
300 last, piwork, pwork)); |
7805
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
301 |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
302 return result; |
62affb34e648
Make quad work with single precision
David Bateman <dbateman@free.fr>
parents:
7482
diff
changeset
|
303 } |