Mercurial > octave
diff libinterp/corefcn/quadcc.cc @ 17787:175b392e91fe
Use GNU style coding conventions for code in libinterp/
* libinterp/corefcn/Cell.h, libinterp/corefcn/__contourc__.cc,
libinterp/corefcn/__dispatch__.cc, libinterp/corefcn/__lin_interpn__.cc,
libinterp/corefcn/__pchip_deriv__.cc, libinterp/corefcn/__qp__.cc,
libinterp/corefcn/balance.cc, libinterp/corefcn/besselj.cc,
libinterp/corefcn/betainc.cc, libinterp/corefcn/bitfcns.cc,
libinterp/corefcn/bsxfun.cc, libinterp/corefcn/c-file-ptr-stream.cc,
libinterp/corefcn/c-file-ptr-stream.h, libinterp/corefcn/cellfun.cc,
libinterp/corefcn/colloc.cc, libinterp/corefcn/comment-list.h,
libinterp/corefcn/conv2.cc, libinterp/corefcn/daspk.cc,
libinterp/corefcn/dasrt.cc, libinterp/corefcn/dassl.cc,
libinterp/corefcn/data.cc, libinterp/corefcn/debug.cc,
libinterp/corefcn/defaults.cc, libinterp/corefcn/defaults.in.h,
libinterp/corefcn/defun-int.h, libinterp/corefcn/defun.cc,
libinterp/corefcn/det.cc, libinterp/corefcn/dirfns.cc,
libinterp/corefcn/display.cc, libinterp/corefcn/dlmread.cc,
libinterp/corefcn/dot.cc, libinterp/corefcn/dynamic-ld.cc,
libinterp/corefcn/dynamic-ld.h, libinterp/corefcn/eig.cc,
libinterp/corefcn/ellipj.cc, libinterp/corefcn/error.cc,
libinterp/corefcn/error.h, libinterp/corefcn/event-queue.h,
libinterp/corefcn/fft.cc, libinterp/corefcn/fft2.cc, libinterp/corefcn/fftn.cc,
libinterp/corefcn/file-io.cc, libinterp/corefcn/filter.cc,
libinterp/corefcn/find.cc, libinterp/corefcn/gammainc.cc,
libinterp/corefcn/gcd.cc, libinterp/corefcn/getgrent.cc,
libinterp/corefcn/getpwent.cc, libinterp/corefcn/getrusage.cc,
libinterp/corefcn/givens.cc, libinterp/corefcn/gl-render.cc,
libinterp/corefcn/gl2ps-renderer.cc, libinterp/corefcn/gl2ps-renderer.h,
libinterp/corefcn/graphics.cc, libinterp/corefcn/graphics.in.h,
libinterp/corefcn/gripes.cc, libinterp/corefcn/gripes.h,
libinterp/corefcn/help.cc, libinterp/corefcn/hess.cc,
libinterp/corefcn/hex2num.cc, libinterp/corefcn/input.cc,
libinterp/corefcn/input.h, libinterp/corefcn/inv.cc,
libinterp/corefcn/jit-ir.h, libinterp/corefcn/jit-typeinfo.cc,
libinterp/corefcn/jit-typeinfo.h, libinterp/corefcn/jit-util.h,
libinterp/corefcn/kron.cc, libinterp/corefcn/load-path.cc,
libinterp/corefcn/load-path.h, libinterp/corefcn/load-save.cc,
libinterp/corefcn/load-save.h, libinterp/corefcn/lookup.cc,
libinterp/corefcn/ls-ascii-helper.cc, libinterp/corefcn/ls-hdf5.cc,
libinterp/corefcn/ls-hdf5.h, libinterp/corefcn/ls-mat-ascii.cc,
libinterp/corefcn/ls-mat-ascii.h, libinterp/corefcn/ls-mat4.cc,
libinterp/corefcn/ls-mat5.cc, libinterp/corefcn/ls-mat5.h,
libinterp/corefcn/ls-oct-ascii.cc, libinterp/corefcn/lsode.cc,
libinterp/corefcn/lu.cc, libinterp/corefcn/luinc.cc,
libinterp/corefcn/mappers.cc, libinterp/corefcn/matrix_type.cc,
libinterp/corefcn/max.cc, libinterp/corefcn/md5sum.cc,
libinterp/corefcn/mex.cc, libinterp/corefcn/mexproto.h,
libinterp/corefcn/mgorth.cc, libinterp/corefcn/mxarray.in.h,
libinterp/corefcn/nproc.cc, libinterp/corefcn/oct-hist.cc,
libinterp/corefcn/oct-lvalue.h, libinterp/corefcn/oct-map.cc,
libinterp/corefcn/oct-map.h, libinterp/corefcn/oct-obj.h,
libinterp/corefcn/oct-prcstrm.h, libinterp/corefcn/oct-stdstrm.h,
libinterp/corefcn/oct-stream.cc, libinterp/corefcn/oct-stream.h,
libinterp/corefcn/octave-link.cc, libinterp/corefcn/octave-link.h,
libinterp/corefcn/pager.cc, libinterp/corefcn/pinv.cc,
libinterp/corefcn/pr-output.cc, libinterp/corefcn/procstream.h,
libinterp/corefcn/profiler.cc, libinterp/corefcn/pt-jit.cc,
libinterp/corefcn/pt-jit.h, libinterp/corefcn/quad.cc,
libinterp/corefcn/quadcc.cc, libinterp/corefcn/qz.cc,
libinterp/corefcn/rand.cc, libinterp/corefcn/rcond.cc,
libinterp/corefcn/regexp.cc, libinterp/corefcn/schur.cc,
libinterp/corefcn/sighandlers.cc, libinterp/corefcn/sighandlers.h,
libinterp/corefcn/sparse-xdiv.cc, libinterp/corefcn/sparse-xdiv.h,
libinterp/corefcn/sparse-xpow.cc, libinterp/corefcn/sparse.cc,
libinterp/corefcn/spparms.cc, libinterp/corefcn/sqrtm.cc,
libinterp/corefcn/str2double.cc, libinterp/corefcn/strfind.cc,
libinterp/corefcn/strfns.cc, libinterp/corefcn/sub2ind.cc,
libinterp/corefcn/svd.cc, libinterp/corefcn/syl.cc,
libinterp/corefcn/symtab.cc, libinterp/corefcn/symtab.h,
libinterp/corefcn/syscalls.cc, libinterp/corefcn/sysdep.cc,
libinterp/corefcn/sysdep.h, libinterp/corefcn/time.cc,
libinterp/corefcn/toplev.cc, libinterp/corefcn/toplev.h,
libinterp/corefcn/tril.cc, libinterp/corefcn/txt-eng-ft.cc,
libinterp/corefcn/txt-eng-ft.h, libinterp/corefcn/txt-eng.h,
libinterp/corefcn/typecast.cc, libinterp/corefcn/urlwrite.cc,
libinterp/corefcn/utils.cc, libinterp/corefcn/variables.cc,
libinterp/corefcn/variables.h, libinterp/corefcn/xdiv.cc,
libinterp/corefcn/xdiv.h, libinterp/corefcn/xnorm.h, libinterp/corefcn/xpow.cc,
libinterp/corefcn/xpow.h, libinterp/corefcn/zfstream.cc,
libinterp/corefcn/zfstream.h, libinterp/dldfcn/__delaunayn__.cc,
libinterp/dldfcn/__dsearchn__.cc, libinterp/dldfcn/__eigs__.cc,
libinterp/dldfcn/__fltk_uigetfile__.cc, libinterp/dldfcn/__glpk__.cc,
libinterp/dldfcn/__init_fltk__.cc, libinterp/dldfcn/__init_gnuplot__.cc,
libinterp/dldfcn/__magick_read__.cc, libinterp/dldfcn/__voronoi__.cc,
libinterp/dldfcn/amd.cc, libinterp/dldfcn/ccolamd.cc, libinterp/dldfcn/chol.cc,
libinterp/dldfcn/colamd.cc, libinterp/dldfcn/convhulln.cc,
libinterp/dldfcn/dmperm.cc, libinterp/dldfcn/fftw.cc, libinterp/dldfcn/qr.cc,
libinterp/dldfcn/symbfact.cc, libinterp/dldfcn/symrcm.cc,
libinterp/dldfcn/tsearch.cc, libinterp/octave-value/ov-base-diag.cc,
libinterp/octave-value/ov-base-diag.h, libinterp/octave-value/ov-base-int.cc,
libinterp/octave-value/ov-base-int.h, libinterp/octave-value/ov-base-mat.h,
libinterp/octave-value/ov-base-scalar.cc,
libinterp/octave-value/ov-base-scalar.h,
libinterp/octave-value/ov-base-sparse.cc,
libinterp/octave-value/ov-base-sparse.h, libinterp/octave-value/ov-base.cc,
libinterp/octave-value/ov-base.h, libinterp/octave-value/ov-bool-mat.cc,
libinterp/octave-value/ov-bool-mat.h, libinterp/octave-value/ov-bool-sparse.cc,
libinterp/octave-value/ov-bool-sparse.h, libinterp/octave-value/ov-bool.cc,
libinterp/octave-value/ov-bool.h, libinterp/octave-value/ov-builtin.cc,
libinterp/octave-value/ov-builtin.h, libinterp/octave-value/ov-cell.cc,
libinterp/octave-value/ov-cell.h, libinterp/octave-value/ov-ch-mat.cc,
libinterp/octave-value/ov-ch-mat.h, libinterp/octave-value/ov-class.cc,
libinterp/octave-value/ov-class.h, libinterp/octave-value/ov-colon.h,
libinterp/octave-value/ov-complex.cc, libinterp/octave-value/ov-complex.h,
libinterp/octave-value/ov-cx-diag.cc, libinterp/octave-value/ov-cx-diag.h,
libinterp/octave-value/ov-cx-mat.cc, libinterp/octave-value/ov-cx-mat.h,
libinterp/octave-value/ov-cx-sparse.cc, libinterp/octave-value/ov-cx-sparse.h,
libinterp/octave-value/ov-dld-fcn.h, libinterp/octave-value/ov-fcn-handle.cc,
libinterp/octave-value/ov-fcn-handle.h,
libinterp/octave-value/ov-fcn-inline.cc,
libinterp/octave-value/ov-fcn-inline.h, libinterp/octave-value/ov-fcn.h,
libinterp/octave-value/ov-float.cc, libinterp/octave-value/ov-float.h,
libinterp/octave-value/ov-flt-complex.cc,
libinterp/octave-value/ov-flt-complex.h,
libinterp/octave-value/ov-flt-cx-diag.cc,
libinterp/octave-value/ov-flt-cx-diag.h,
libinterp/octave-value/ov-flt-cx-mat.cc,
libinterp/octave-value/ov-flt-cx-mat.h,
libinterp/octave-value/ov-flt-re-diag.cc,
libinterp/octave-value/ov-flt-re-diag.h,
libinterp/octave-value/ov-flt-re-mat.cc,
libinterp/octave-value/ov-flt-re-mat.h, libinterp/octave-value/ov-int16.cc,
libinterp/octave-value/ov-int32.cc, libinterp/octave-value/ov-int64.cc,
libinterp/octave-value/ov-int8.cc, libinterp/octave-value/ov-intx.h,
libinterp/octave-value/ov-java.cc, libinterp/octave-value/ov-lazy-idx.h,
libinterp/octave-value/ov-mex-fcn.cc, libinterp/octave-value/ov-mex-fcn.h,
libinterp/octave-value/ov-null-mat.cc, libinterp/octave-value/ov-null-mat.h,
libinterp/octave-value/ov-oncleanup.cc, libinterp/octave-value/ov-perm.cc,
libinterp/octave-value/ov-perm.h, libinterp/octave-value/ov-range.cc,
libinterp/octave-value/ov-range.h, libinterp/octave-value/ov-re-diag.cc,
libinterp/octave-value/ov-re-diag.h, libinterp/octave-value/ov-re-mat.cc,
libinterp/octave-value/ov-re-mat.h, libinterp/octave-value/ov-re-sparse.cc,
libinterp/octave-value/ov-re-sparse.h, libinterp/octave-value/ov-scalar.cc,
libinterp/octave-value/ov-scalar.h, libinterp/octave-value/ov-str-mat.cc,
libinterp/octave-value/ov-str-mat.h, libinterp/octave-value/ov-struct.cc,
libinterp/octave-value/ov-struct.h, libinterp/octave-value/ov-type-conv.h,
libinterp/octave-value/ov-typeinfo.cc, libinterp/octave-value/ov-typeinfo.h,
libinterp/octave-value/ov-uint16.cc, libinterp/octave-value/ov-uint32.cc,
libinterp/octave-value/ov-uint64.cc, libinterp/octave-value/ov-uint8.cc,
libinterp/octave-value/ov-usr-fcn.cc, libinterp/octave-value/ov-usr-fcn.h,
libinterp/octave-value/ov.cc, libinterp/octave-value/ov.h, libinterp/octave.cc,
libinterp/operators/op-b-bm.cc, libinterp/operators/op-b-sbm.cc,
libinterp/operators/op-bm-b.cc, libinterp/operators/op-bm-bm.cc,
libinterp/operators/op-cdm-cdm.cc, libinterp/operators/op-chm.cc,
libinterp/operators/op-class.cc, libinterp/operators/op-cm-cm.cc,
libinterp/operators/op-cm-cs.cc, libinterp/operators/op-cm-s.cc,
libinterp/operators/op-cm-scm.cc, libinterp/operators/op-cm-sm.cc,
libinterp/operators/op-cs-cm.cc, libinterp/operators/op-cs-cs.cc,
libinterp/operators/op-cs-scm.cc, libinterp/operators/op-cs-sm.cc,
libinterp/operators/op-dm-dm.cc, libinterp/operators/op-dm-scm.cc,
libinterp/operators/op-double-conv.cc, libinterp/operators/op-fcdm-fcdm.cc,
libinterp/operators/op-fcm-fcm.cc, libinterp/operators/op-fcm-fcs.cc,
libinterp/operators/op-fcm-fm.cc, libinterp/operators/op-fcm-fs.cc,
libinterp/operators/op-fcs-fcm.cc, libinterp/operators/op-fcs-fcs.cc,
libinterp/operators/op-fcs-fm.cc, libinterp/operators/op-fcs-fs.cc,
libinterp/operators/op-fdm-fdm.cc, libinterp/operators/op-float-conv.cc,
libinterp/operators/op-fm-fcm.cc, libinterp/operators/op-fm-fcs.cc,
libinterp/operators/op-fm-fm.cc, libinterp/operators/op-fm-fs.cc,
libinterp/operators/op-fs-fcm.cc, libinterp/operators/op-fs-fcs.cc,
libinterp/operators/op-fs-fm.cc, libinterp/operators/op-fs-fs.cc,
libinterp/operators/op-m-cm.cc, libinterp/operators/op-m-cs.cc,
libinterp/operators/op-m-m.cc, libinterp/operators/op-m-s.cc,
libinterp/operators/op-m-scm.cc, libinterp/operators/op-m-sm.cc,
libinterp/operators/op-pm-scm.cc, libinterp/operators/op-range.cc,
libinterp/operators/op-s-cm.cc, libinterp/operators/op-s-cs.cc,
libinterp/operators/op-s-scm.cc, libinterp/operators/op-sbm-b.cc,
libinterp/operators/op-sbm-bm.cc, libinterp/operators/op-sbm-sbm.cc,
libinterp/operators/op-scm-cm.cc, libinterp/operators/op-scm-cs.cc,
libinterp/operators/op-scm-m.cc, libinterp/operators/op-scm-s.cc,
libinterp/operators/op-scm-scm.cc, libinterp/operators/op-scm-sm.cc,
libinterp/operators/op-sm-cm.cc, libinterp/operators/op-sm-m.cc,
libinterp/operators/op-sm-s.cc, libinterp/operators/op-sm-scm.cc,
libinterp/operators/op-sm-sm.cc, libinterp/operators/op-str-m.cc,
libinterp/operators/op-str-s.cc, libinterp/operators/op-str-str.cc,
libinterp/operators/ops.h, libinterp/parse-tree/lex.h,
libinterp/parse-tree/parse.h, libinterp/parse-tree/pt-arg-list.cc,
libinterp/parse-tree/pt-arg-list.h, libinterp/parse-tree/pt-assign.cc,
libinterp/parse-tree/pt-assign.h, libinterp/parse-tree/pt-binop.cc,
libinterp/parse-tree/pt-binop.h, libinterp/parse-tree/pt-bp.h,
libinterp/parse-tree/pt-cbinop.cc, libinterp/parse-tree/pt-check.cc,
libinterp/parse-tree/pt-colon.cc, libinterp/parse-tree/pt-colon.h,
libinterp/parse-tree/pt-const.cc, libinterp/parse-tree/pt-decl.cc,
libinterp/parse-tree/pt-decl.h, libinterp/parse-tree/pt-eval.cc,
libinterp/parse-tree/pt-except.h, libinterp/parse-tree/pt-exp.h,
libinterp/parse-tree/pt-fcn-handle.cc, libinterp/parse-tree/pt-id.cc,
libinterp/parse-tree/pt-id.h, libinterp/parse-tree/pt-idx.cc,
libinterp/parse-tree/pt-idx.h, libinterp/parse-tree/pt-loop.h,
libinterp/parse-tree/pt-mat.cc, libinterp/parse-tree/pt-misc.cc,
libinterp/parse-tree/pt-misc.h, libinterp/parse-tree/pt-pr-code.cc,
libinterp/parse-tree/pt-select.h, libinterp/parse-tree/pt-stmt.h,
libinterp/parse-tree/token.h, libinterp/version.cc:
Use GNU style coding conventions for code in libinterp/
| author | Rik <rik@octave.org> |
|---|---|
| date | Mon, 28 Oct 2013 19:51:46 -0700 |
| parents | d63878346099 |
| children | 0850b5212619 446c46af4b42 |
line wrap: on
line diff
--- a/libinterp/corefcn/quadcc.cc Mon Oct 28 21:04:48 2013 -0400 +++ b/libinterp/corefcn/quadcc.cc Mon Oct 28 19:51:46 2013 -0700 @@ -38,7 +38,7 @@ #define DEBUG_QUADCC 0 /* Define the minimum size of the interval heap. */ -#define min_cquad_heapsize 200 +#define min_cquad_heapsize 200 /* Data of a single interval */ @@ -53,7 +53,8 @@ /* Some constants and matrices that we'll need. */ -static const double xi[33] = { +static const double xi[33] = +{ -1., -0.99518472667219688624, -0.98078528040323044912, -0.95694033573220886493, -0.92387953251128675612, -0.88192126434835502970, -0.83146961230254523708, @@ -69,7 +70,8 @@ 0.98078528040323044912, 0.99518472667219688624, 1. }; -static const double bee[68] = { +static const double bee[68] = +{ 0.00000000000000e+00, 2.28868854108532e-01, 0.00000000000000e+00, -8.15740215243451e-01, 0.00000000000000e+00, 5.31212715259731e-01, 0.00000000000000e+00, 1.38538036812454e-02, 0.00000000000000e+00, @@ -95,7 +97,8 @@ 0.00000000000000e+00, 5.40679930965238e-01 }; -static const double Lalpha[33] = { +static const double Lalpha[33] = +{ 5.77350269189626e-01, 5.16397779494322e-01, 5.07092552837110e-01, 5.03952630678970e-01, 5.02518907629606e-01, 5.01745206004255e-01, 5.01280411827603e-01, 5.00979432868120e-01, 5.00773395667191e-01, @@ -109,7 +112,8 @@ 5.00065049112355e-01, 5.00061046334395e-01, 5.00057401986298e-01 }; -static const double Lgamma[33] = { +static const double Lgamma[33] = +{ 0.0, 0.0, 5.16397779494322e-01, 5.07092552837110e-01, 5.03952630678970e-01, 5.02518907629606e-01, 5.01745206004255e-01, 5.01280411827603e-01, 5.00979432868120e-01, 5.00773395667191e-01, 5.00626174321759e-01, @@ -123,7 +127,8 @@ 5.00061046334395e-01 }; -static const double V1inv[5 * 5] = { +static const double V1inv[5 * 5] = +{ .47140452079103168293e-1, .37712361663282534635, .56568542494923801952, .37712361663282534635, .47140452079103168293e-1, -.81649658092772603273e-1, -.46188021535170061160, 0, @@ -135,7 +140,8 @@ .10774960475223581324 }; -static const double V2inv[9 * 9] = { +static const double V2inv[9 * 9] = +{ .11223917161691230546e-1, .10339219839658349826, .19754094204576565761, .25577315077753587922, .27835314560994251755, .25577315077753587922, .19754094204576565761, .10339219839658349826, .11223917161691230546e-1, @@ -169,7 +175,8 @@ .10916211417928767644, -.10916211417928767644, .54581057089643838221e-1 }; -static const double V3inv[17 * 17] = { +static const double V3inv[17 * 17] = +{ .27729677693590098996e-2, .26423663180333065153e-1, .53374068493933898312e-1, .77007854739523195947e-1, .98257061072911596869e-1, .11538049741786835604, .12832134344120884559, @@ -306,7 +313,8 @@ -.54971216929497681146e-1, .27485608464748840573e-1 }; -static const double V4inv[33 * 33] = { +static const double V4inv[33 * 33] = +{ .69120897476690862600e-3, .66419939766331555194e-2, .13600665164323186111e-1, .20122785860913684493e-1, .26583214101668429944e-1, .32712713318999268739e-1, @@ -882,7 +890,8 @@ -.27588282524939131481e-1, .13794141262469565740e-1 }; -static const double Tleft[33 * 33] = { +static const double Tleft[33 * 33] = +{ 1., -.86602540378443864678, 0., .33071891388307382381, 0., -.20728904939721249057, 0., .15128841196122722208, 0., -.11918864298744029244, 0., .98352013661686631224e-1, 0., @@ -1136,253 +1145,254 @@ 0., 0., .23283064365386962891e-9 }; -static const double Tright[33 * 33] = { +static const double Tright[33 * 33] = +{ 1., .86602540378443864678, 0., -.33071891388307382381, 0., - .20728904939721249057, 0., -.15128841196122722208, 0., - .11918864298744029244, 0., -.98352013661686631224e-1, 0., - .83727065404940845733e-1, 0., -.72893399403505841203e-1, 0., - .64544632643375022436e-1, 0., -.57913170372415565639e-1, 0., - .52518242575729562263e-1, 0., -.48043311993977520457e-1, 0., - .44271433659733990243e-1, 0., -.41048928022856771981e-1, 0., - .38263878662008271459e-1, 0., -.35832844026365304501e-1, 0., 0., - .50000000000000000000, .96824583655185422130, .57282196186948000082, - -.21650635094610966169, -.35903516540862679125, .97578093724974971969e-1, - .26203921611325660506, -.55792409597991015609e-1, -.20644078533943456204, - .36172381205961199479e-1, .17035068468874958194, - -.25371838001497225980e-1, -.14501953125000000000, - .18786835250972344757e-1, .12625507130328301066, - -.14473795929590520582e-1, -.11179458309419422675, - .11494434254897626155e-1, .10030855351241635862, - -.93498556820544479096e-2, -.90964264465390582629e-1, - .77546391824364392762e-2, .83213457337452292745e-1, - -.65358085945588638605e-2, -.76680372422574234569e-1, - .55835321940047427169e-2, .71098828931825789428e-1, - -.48253327982967591019e-2, -.66274981937248958553e-1, - .42118078245337801387e-2, .62064306433355646267e-1, - -.37083386598903548973e-2, 0., 0., .25000000000000000000, - .73950997288745200531, .83852549156242113615, .23175620272173946716, - -.37791833195149451496, -.25710129174850522325, .21608307321780204633, - .22844049245646009157, -.14009503000335388415, -.19897685605518413847, - .98264706042471226893e-1, .17445445004279014046, - -.72761100054958328401e-1, -.15463589893742108388, - .56056770591708784481e-1, .13855313872640495158, - -.44517752443294564781e-1, -.12534277657695128850, - .36211835346039665762e-1, .11434398255136139683, - -.30033588409423828125e-1, -.10506705408753910481, - .25313077840725783008e-1, .97149327637744872155e-1, - -.21624927200393328444e-1, -.90319582367202122625e-1, - .18688433567711780666e-1, .84372291635345108584e-1, - -.16312261561845420752e-1, -.79149526894804751586e-1, - .14362333871852474757e-1, 0., 0., 0., .12500000000000000000, - .49607837082461073572, .82265291131801144317, .59621200088559103072, - -.80054302859059362371e-1, -.42612156697795759420, - -.90098145270865592887e-1, .29769623255090078484, .13630307904779758221, - -.21638835185708931831, -.14600247270306082052, .16348801804014290453, - .14340708728599057249, -.12755243353979286190, -.13661523715071346961, - .10215585947881057394, .12864248070157166547, -.83592528025348693602e-1, - -.12066728689302565222, .69633728678718053052e-1, .11314245177331919532, - -.58882939251410088028e-1, -.10621835858758221487, - .50432266865187597572e-1, .99916834723527771581e-1, - -.43672094283057258509e-1, -.94206380251950852413e-1, - .38181356812697746418e-1, .89035739656537771225e-1, - -.33661934598216332678e-1, 0., 0., 0., 0., .62500000000000000000e-1, - .31093357409581873586, .67604086414949799246, .75644205980613611039, - .28990586430124175741, -.30648508196770360914, -.35801372616842500052, - .91326869828709014708e-1, .31127929687500000000, .90915752838698393094e-2, - -.25637381283965534330, -.57601077850322797594e-1, .21019685709225757945, - .81244992138514014256e-1, -.17375078516720988858, - -.92289437277967051125e-1, .14527351914265391374, - .96675340792832019889e-1, -.12289485697108543415, - -.97448175340011084006e-1, .10511755943298339844, - .96242247086378239657e-1, -.90822942272780513537e-1, - -.93966350452322132384e-1, .79189411876493712558e-1, - .91139307067989309325e-1, -.69613039934383197265e-1, - -.88062491671135767870e-1, .61646331729340817494e-1, 0., 0., 0., 0., 0., - .31250000000000000000e-1, .18684782411095934408, .50176689760410660236, - .74784031498626095398, .56472001151566251186, -.14842464993721351203e-1, - -.41162920273003120936, -.20243071230196532282, .23772054897172750436, - .24963810923972235950, -.12116179938394678936, -.24330535483519110663, - .47903849781124471359e-1, .22133299683101224293, - -.20542915138527200983e-2, -.19653465717678146728, - -.26818172626509178444e-1, .17319122357631210944, - .45065391411065545445e-1, -.15253391395444065941, - -.56543897711725408302e-1, .13469154928743585367, - .63632471400208840155e-1, -.11941684923913523817, - -.67828850207933293098e-1, .10636309084510652670, - .70095786922999181504e-1, -.95187373095150709082e-1, 0., 0., 0., 0., 0., - 0., .15625000000000000000e-1, .10909562534194485289, - .34842348626527747318, .64461114561628111443, .69382480527334683659, - .29551102358528827763, -.25527584713978439819, -.38878771718544715394, - -.82956185835347407489e-2, .31183177761966943912, .12831420840372374767, - -.22067618205599434368, -.17569196937129496961, .14598057000132284135, - .18864406621763419484, -.89921002550386645767e-1, -.18571835020187122114, - .48967672227195481777e-1, .17584685670380332798, - -.19267984545067426324e-1, -.16335437520503462738, - -.22598055455032407594e-2, .15032800884170631129, - .17883358353754640871e-1, -.13774837869432209951, - -.29227555960587143675e-1, .12604194747513151053, 0., 0., 0., 0., 0., 0., - 0., .78125000000000000000e-2, .62377810244809812496e-1, - .23080781467370883845, .50841310636012325368, .69834547012574056043, - .52572723156526459672, .11464215704954976471e-1, -.38698869011491210342, - -.26125646622255207507, .16951698812361607510, .29773875898928782269, - -.20130501202570367491e-1, -.26332493149159310198, - -.67734613690401207009e-1, .21207315477103762715, .11541543390889415193, - -.16249634759782417533, -.13885887405041735068, .11996491328010275427, - .14810432001630926895, -.85177658352556243411e-1, -.14918860659904380587, - .57317789510444151564e-1, .14569827645586660151, - -.35213090145965327390e-1, -.13975998126844578198, 0., 0., 0., 0., 0., 0., - 0., 0., .39062500000000000000e-2, .35101954600803571207e-1, - .14761284084133737720, .37655033076080192966, .62410290231517322776, - .64335622317683389875, .28188168266139524244, -.22488495672137010675, - -.39393811089283576186, -.75184777995770096714e-1, .28472023119398293003, - .20410910833705899572, -.15590046962908511750, -.23814567544617953125, - .54442805556829031204e-1, .22855930338589720954, .16303223615756629897e-1, - -.20172722433875559213, -.62723406421217419404e-1, .17012230831020922010, - .91754642766136561612e-1, -.13927644821381121197, -.10886600968068418181, - .11139075654373395292, .11797455976331702879, 0., 0., 0., 0., 0., 0., 0., - 0., 0., .19531250000000000000e-2, .19506820659607596598e-1, - .91865676095362231937e-1, .26604607809696493849, .51425874205091288223, - .66047561132505329292, .48660109511591303851, .17575661168678285615e-1, - -.36594333408055703366, -.29088854695378694533, .11318677346656537927, - .31110645235730182168, .60733219161008787341e-1, -.24333848233620420826, - -.15254312332655419708, .15995968483455388613, .19010344455215289289, - -.86040636766440260000e-1, -.19652589954665259945, - .27633388517205837713e-1, .18660848552712880387, .15942583868416775867e-1, - -.16902042462382064786, -.47278526495327740646e-1, 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., .97656250000000000000e-3, .10731084460857378207e-1, - .55939644713816406331e-1, .18118487371914493668, .39914857299829864263, - .60812322949933902435, .60011887183061967583, .26002695805835928795, - -.20883922404786010096, -.38988130966114638081, -.11797833550782589082, - .25231824756239520077, .24817859972953934712, -.90516417677868996417e-1, - -.26079073291293066798, -.30259468817169480161e-1, .22178195264114178432, - .10569877864302048175, -.16679648389266977455, -.14637718550245050850, - .11219272032739559870, .16359363640525750353, -.64358194509092101393e-1, - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., .48828125000000000000e-3, - .58542865274813470967e-2, .33461741635290096452e-1, .11979993155896201271, - .29580223766987206958, .51874761979436016742, .62861483498014306968, - .44868895761051453296, .12567502628371529386e-1, -.35040366183235474275, - -.30466868455569500886, .70903913601490112666e-1, .30822791893032512740, - .11969443264190207736, -.20764760317621313946, -.20629838355452128532, - .95269702915334718507e-1, .22432624768705133300, - -.33103381593477797101e-2, -.20570036048155716333, - -.62208282720094518964e-1, .17095309330441436348, 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., .24414062500000000000e-3, - .31714797501871532475e-2, .19721062526127334100e-1, - .77311181185536498246e-1, .21124871792841566575, .41777980401893650886, - .59401977834943551650, .56132417807488349048, .23433675061367565951, - -.20222775295220942126, -.38280372496506190127, -.14443804214023095767, - .22268950939178466797, .27211314150777981984, -.34184876506180717313e-1, - -.26006498895669734842, -.97650425186005090107e-1, .19024527660129101293, - .16789164198044635671, -.10875811641651905252, -.19276785058805921298, 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., .12207031250000000000e-3, - .17078941137247586143e-2, .11477733754843910060e-1, - .48887017020924625462e-1, .14634927241421789683, .32156282683019547854, - .52165811920227223937, .60001958466396926460, .41208501541480733755, - .11366945503190350975e-2, -.33968093962672089159, -.30955190935923386766, - .40657421856578262210e-1, .29873400409871531764, .16094481791768257440, - -.16876122436206497694, -.23650217045022161255, .33070260090574765012e-1, - .22985258456375907796, .68645651043827097771e-1, 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., .61035156250000000000e-4, - .91501857608428649078e-3, .66085179496951987952e-2, - .30383171695850355404e-1, .98840838845366876117e-1, .23855447246420318989, - .43322017468145613917, .58049033744876107191, .52533893203742699346, - .20681056202371946180, -.20180000924562504384, -.37503922291962681797, - -.15988102869837429062, .19823558102762374094, .28393023878803799622, - .11188133439357510403e-1, -.24730368377168229255, -.14731529061377942839, - .14878558042884266021, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., .30517578125000000000e-4, .48804277318479845551e-3, - .37696080990601968396e-2, .18603912108994738255e-1, - .65325006755649582964e-1, .17162960707938819795, .34411527956476971322, - .52289350347082497959, .57319653625674910592, .37662253421045430413, - -.14099055105384663902e-1, -.33265570610216904208, -.30921265572647566661, - .19911390594166455281e-1, .28738590811031797718, .18912130469738472647, - -.13235936203215819193, -.25076406142356675279, 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., .15258789062500000000e-4, - .25928719280954633249e-3, .21327398937568540428e-2, - .11244626133630732010e-1, .42375605740664331966e-1, .12031130345907846211, - .26352562258934426830, .44590628258512682078, .56682835613700749379, - .49116715128261660395, .17845943097110339078, -.20541650677432497477, - -.36739803642257458221, -.16776034069210108273, .17920950989905112908, - .28867732805385066532, .46473465543376206337e-1, 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., .76293945312500000000e-5, - .13727610943181290891e-3, .11979683091449349286e-2, - .67195313034570709806e-2, .27044920779931968175e-1, - .82472196498517457862e-1, .19570475044896150093, .36391620788543817693, - .52241392782736588032, .54727504974907879912, .34211551468813581183, - -.31580472732719957762e-1, -.32830006549176759667, -.30563797665254420769, - .64905014620683140120e-2, .27642986248995073032, 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., .38146972656250000000e-5, - .72454147007837596854e-4, .66859847582761390285e-3, - .39751311980366118437e-2, .17015198650201528366e-1, - .55443621868993855715e-1, .14157060481641692131, .28641242619559616836, - .45610665490966615415, .55262786406029265394, .45818352706035500108, - .14984403004611673047, -.21163807462970713245, -.36007252928843413718, - -.17030961385712954159, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., .19073486328125000000e-5, - .38135049864067468562e-4, .37101393638555730015e-3, - .23305339886279723213e-2, .10569913448297127219e-1, - .36640175162216897547e-1, .10010476414320235508, .21860074212675559892, - .38124757096345313719, .52020999209879669177, .52172632730659212045, - .30841620620308814614, -.50322546186721500184e-1, -.32577618885114899053, - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., .95367431640625000000e-6, .20021483206955925244e-4, - .20481807322420625431e-3, .13553476938058909882e-2, - .64919676350791905019e-2, .23848725425069251903e-1, - .69384632678886421292e-1, .16249711393618776934, .30736618106830314788, - .46399909601971539157, .53765031034002467225, .42598991476520183929, - .12130445348350215652, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., .47683715820312500000e-6, - .10487707828484902486e-4, .11254146162337528943e-3, - .78248929534271987118e-3, .39468337145306794566e-2, - .15313546659475671763e-1, .47249070825218564146e-1, .11804374107101480543, - .24031796927792491122, .39629215049166341285, .51629108968402548545, - .49622372075429782915, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., .23841857910156250000e-6, - .54823314130625337326e-5, .61575377321535518154e-4, - .44877834366497538134e-3, .23774612048621955857e-2, - .97136347645161687796e-2, .31671599547606636717e-1, - .84028665767000747480e-1, .18298487576742964949, .32647878537696945218, - .46970971486488895077, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., .11920928955078125000e-6, - .28604020001177375838e-5, .33559227978295551013e-4, - .25583821662860610560e-3, .14201552747787302339e-2, - .60938046986874414969e-2, .20930869247951926793e-1, - .58745021125678072911e-1, .13613725780285953720, .26083988356030237586, - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., .59604644775390625000e-7, - .14898180663526043291e-5, .18224991282807693921e-4, - .14504433444608833821e-3, .84184722720281809548e-3, - .37846965430000478789e-2, .13656355548211376864e-1, - .40409541997718853934e-1, .99226988101858325902e-1, 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., .29802322387695312500e-7, .77471708843445529468e-6, - .98649879372606876995e-5, .81814934772838523887e-4, - .49554483992403011328e-3, .23290922072351413938e-2, - .88068134250844034186e-2, .27393666952485719070e-1, 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., 0., .14901161193847656250e-7, .40226235946098233685e-6, - .53236418690561306700e-5, .45933829691164002269e-4, - .28982005232838857913e-3, .14212974043211018374e-2, - .56192363087488842264e-2, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - .74505805969238281250e-8, .20858299254133430408e-6, - .28648457300134381744e-5, .25677535898258910850e-4, - .16849420429491355445e-3, .86062824010315834002e-3, 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., .37252902984619140625e-8, .10801736017613096861e-6, - .15376606719887104015e-5, .14296523739727437959e-4, - .97419023656050887203e-4, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - .18626451492309570312e-8, .55871592916438890146e-7, - .82331193828137454068e-6, .79302250528382787666e-5, 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., .93132257461547851562e-9, - .28867244235852488244e-7, .43982811713864556957e-6, 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., .46566128730773925781e-9, - .14899342093408253335e-7, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., - 0., 0., .23283064365386962891e-9 + .20728904939721249057, 0., -.15128841196122722208, 0., + .11918864298744029244, 0., -.98352013661686631224e-1, 0., + .83727065404940845733e-1, 0., -.72893399403505841203e-1, 0., + .64544632643375022436e-1, 0., -.57913170372415565639e-1, 0., + .52518242575729562263e-1, 0., -.48043311993977520457e-1, 0., + .44271433659733990243e-1, 0., -.41048928022856771981e-1, 0., + .38263878662008271459e-1, 0., -.35832844026365304501e-1, 0., 0., + .50000000000000000000, .96824583655185422130, .57282196186948000082, + -.21650635094610966169, -.35903516540862679125, .97578093724974971969e-1, + .26203921611325660506, -.55792409597991015609e-1, -.20644078533943456204, + .36172381205961199479e-1, .17035068468874958194, + -.25371838001497225980e-1, -.14501953125000000000, + .18786835250972344757e-1, .12625507130328301066, + -.14473795929590520582e-1, -.11179458309419422675, + .11494434254897626155e-1, .10030855351241635862, + -.93498556820544479096e-2, -.90964264465390582629e-1, + .77546391824364392762e-2, .83213457337452292745e-1, + -.65358085945588638605e-2, -.76680372422574234569e-1, + .55835321940047427169e-2, .71098828931825789428e-1, + -.48253327982967591019e-2, -.66274981937248958553e-1, + .42118078245337801387e-2, .62064306433355646267e-1, + -.37083386598903548973e-2, 0., 0., .25000000000000000000, + .73950997288745200531, .83852549156242113615, .23175620272173946716, + -.37791833195149451496, -.25710129174850522325, .21608307321780204633, + .22844049245646009157, -.14009503000335388415, -.19897685605518413847, + .98264706042471226893e-1, .17445445004279014046, + -.72761100054958328401e-1, -.15463589893742108388, + .56056770591708784481e-1, .13855313872640495158, + -.44517752443294564781e-1, -.12534277657695128850, + .36211835346039665762e-1, .11434398255136139683, + -.30033588409423828125e-1, -.10506705408753910481, + .25313077840725783008e-1, .97149327637744872155e-1, + -.21624927200393328444e-1, -.90319582367202122625e-1, + .18688433567711780666e-1, .84372291635345108584e-1, + -.16312261561845420752e-1, -.79149526894804751586e-1, + .14362333871852474757e-1, 0., 0., 0., .12500000000000000000, + .49607837082461073572, .82265291131801144317, .59621200088559103072, + -.80054302859059362371e-1, -.42612156697795759420, + -.90098145270865592887e-1, .29769623255090078484, .13630307904779758221, + -.21638835185708931831, -.14600247270306082052, .16348801804014290453, + .14340708728599057249, -.12755243353979286190, -.13661523715071346961, + .10215585947881057394, .12864248070157166547, -.83592528025348693602e-1, + -.12066728689302565222, .69633728678718053052e-1, .11314245177331919532, + -.58882939251410088028e-1, -.10621835858758221487, + .50432266865187597572e-1, .99916834723527771581e-1, + -.43672094283057258509e-1, -.94206380251950852413e-1, + .38181356812697746418e-1, .89035739656537771225e-1, + -.33661934598216332678e-1, 0., 0., 0., 0., .62500000000000000000e-1, + .31093357409581873586, .67604086414949799246, .75644205980613611039, + .28990586430124175741, -.30648508196770360914, -.35801372616842500052, + .91326869828709014708e-1, .31127929687500000000, .90915752838698393094e-2, + -.25637381283965534330, -.57601077850322797594e-1, .21019685709225757945, + .81244992138514014256e-1, -.17375078516720988858, + -.92289437277967051125e-1, .14527351914265391374, + .96675340792832019889e-1, -.12289485697108543415, + -.97448175340011084006e-1, .10511755943298339844, + .96242247086378239657e-1, -.90822942272780513537e-1, + -.93966350452322132384e-1, .79189411876493712558e-1, + .91139307067989309325e-1, -.69613039934383197265e-1, + -.88062491671135767870e-1, .61646331729340817494e-1, 0., 0., 0., 0., 0., + .31250000000000000000e-1, .18684782411095934408, .50176689760410660236, + .74784031498626095398, .56472001151566251186, -.14842464993721351203e-1, + -.41162920273003120936, -.20243071230196532282, .23772054897172750436, + .24963810923972235950, -.12116179938394678936, -.24330535483519110663, + .47903849781124471359e-1, .22133299683101224293, + -.20542915138527200983e-2, -.19653465717678146728, + -.26818172626509178444e-1, .17319122357631210944, + .45065391411065545445e-1, -.15253391395444065941, + -.56543897711725408302e-1, .13469154928743585367, + .63632471400208840155e-1, -.11941684923913523817, + -.67828850207933293098e-1, .10636309084510652670, + .70095786922999181504e-1, -.95187373095150709082e-1, 0., 0., 0., 0., 0., + 0., .15625000000000000000e-1, .10909562534194485289, + .34842348626527747318, .64461114561628111443, .69382480527334683659, + .29551102358528827763, -.25527584713978439819, -.38878771718544715394, + -.82956185835347407489e-2, .31183177761966943912, .12831420840372374767, + -.22067618205599434368, -.17569196937129496961, .14598057000132284135, + .18864406621763419484, -.89921002550386645767e-1, -.18571835020187122114, + .48967672227195481777e-1, .17584685670380332798, + -.19267984545067426324e-1, -.16335437520503462738, + -.22598055455032407594e-2, .15032800884170631129, + .17883358353754640871e-1, -.13774837869432209951, + -.29227555960587143675e-1, .12604194747513151053, 0., 0., 0., 0., 0., 0., + 0., .78125000000000000000e-2, .62377810244809812496e-1, + .23080781467370883845, .50841310636012325368, .69834547012574056043, + .52572723156526459672, .11464215704954976471e-1, -.38698869011491210342, + -.26125646622255207507, .16951698812361607510, .29773875898928782269, + -.20130501202570367491e-1, -.26332493149159310198, + -.67734613690401207009e-1, .21207315477103762715, .11541543390889415193, + -.16249634759782417533, -.13885887405041735068, .11996491328010275427, + .14810432001630926895, -.85177658352556243411e-1, -.14918860659904380587, + .57317789510444151564e-1, .14569827645586660151, + -.35213090145965327390e-1, -.13975998126844578198, 0., 0., 0., 0., 0., 0., + 0., 0., .39062500000000000000e-2, .35101954600803571207e-1, + .14761284084133737720, .37655033076080192966, .62410290231517322776, + .64335622317683389875, .28188168266139524244, -.22488495672137010675, + -.39393811089283576186, -.75184777995770096714e-1, .28472023119398293003, + .20410910833705899572, -.15590046962908511750, -.23814567544617953125, + .54442805556829031204e-1, .22855930338589720954, .16303223615756629897e-1, + -.20172722433875559213, -.62723406421217419404e-1, .17012230831020922010, + .91754642766136561612e-1, -.13927644821381121197, -.10886600968068418181, + .11139075654373395292, .11797455976331702879, 0., 0., 0., 0., 0., 0., 0., + 0., 0., .19531250000000000000e-2, .19506820659607596598e-1, + .91865676095362231937e-1, .26604607809696493849, .51425874205091288223, + .66047561132505329292, .48660109511591303851, .17575661168678285615e-1, + -.36594333408055703366, -.29088854695378694533, .11318677346656537927, + .31110645235730182168, .60733219161008787341e-1, -.24333848233620420826, + -.15254312332655419708, .15995968483455388613, .19010344455215289289, + -.86040636766440260000e-1, -.19652589954665259945, + .27633388517205837713e-1, .18660848552712880387, .15942583868416775867e-1, + -.16902042462382064786, -.47278526495327740646e-1, 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., .97656250000000000000e-3, .10731084460857378207e-1, + .55939644713816406331e-1, .18118487371914493668, .39914857299829864263, + .60812322949933902435, .60011887183061967583, .26002695805835928795, + -.20883922404786010096, -.38988130966114638081, -.11797833550782589082, + .25231824756239520077, .24817859972953934712, -.90516417677868996417e-1, + -.26079073291293066798, -.30259468817169480161e-1, .22178195264114178432, + .10569877864302048175, -.16679648389266977455, -.14637718550245050850, + .11219272032739559870, .16359363640525750353, -.64358194509092101393e-1, + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., .48828125000000000000e-3, + .58542865274813470967e-2, .33461741635290096452e-1, .11979993155896201271, + .29580223766987206958, .51874761979436016742, .62861483498014306968, + .44868895761051453296, .12567502628371529386e-1, -.35040366183235474275, + -.30466868455569500886, .70903913601490112666e-1, .30822791893032512740, + .11969443264190207736, -.20764760317621313946, -.20629838355452128532, + .95269702915334718507e-1, .22432624768705133300, + -.33103381593477797101e-2, -.20570036048155716333, + -.62208282720094518964e-1, .17095309330441436348, 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., .24414062500000000000e-3, + .31714797501871532475e-2, .19721062526127334100e-1, + .77311181185536498246e-1, .21124871792841566575, .41777980401893650886, + .59401977834943551650, .56132417807488349048, .23433675061367565951, + -.20222775295220942126, -.38280372496506190127, -.14443804214023095767, + .22268950939178466797, .27211314150777981984, -.34184876506180717313e-1, + -.26006498895669734842, -.97650425186005090107e-1, .19024527660129101293, + .16789164198044635671, -.10875811641651905252, -.19276785058805921298, 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., .12207031250000000000e-3, + .17078941137247586143e-2, .11477733754843910060e-1, + .48887017020924625462e-1, .14634927241421789683, .32156282683019547854, + .52165811920227223937, .60001958466396926460, .41208501541480733755, + .11366945503190350975e-2, -.33968093962672089159, -.30955190935923386766, + .40657421856578262210e-1, .29873400409871531764, .16094481791768257440, + -.16876122436206497694, -.23650217045022161255, .33070260090574765012e-1, + .22985258456375907796, .68645651043827097771e-1, 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., .61035156250000000000e-4, + .91501857608428649078e-3, .66085179496951987952e-2, + .30383171695850355404e-1, .98840838845366876117e-1, .23855447246420318989, + .43322017468145613917, .58049033744876107191, .52533893203742699346, + .20681056202371946180, -.20180000924562504384, -.37503922291962681797, + -.15988102869837429062, .19823558102762374094, .28393023878803799622, + .11188133439357510403e-1, -.24730368377168229255, -.14731529061377942839, + .14878558042884266021, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., .30517578125000000000e-4, .48804277318479845551e-3, + .37696080990601968396e-2, .18603912108994738255e-1, + .65325006755649582964e-1, .17162960707938819795, .34411527956476971322, + .52289350347082497959, .57319653625674910592, .37662253421045430413, + -.14099055105384663902e-1, -.33265570610216904208, -.30921265572647566661, + .19911390594166455281e-1, .28738590811031797718, .18912130469738472647, + -.13235936203215819193, -.25076406142356675279, 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., .15258789062500000000e-4, + .25928719280954633249e-3, .21327398937568540428e-2, + .11244626133630732010e-1, .42375605740664331966e-1, .12031130345907846211, + .26352562258934426830, .44590628258512682078, .56682835613700749379, + .49116715128261660395, .17845943097110339078, -.20541650677432497477, + -.36739803642257458221, -.16776034069210108273, .17920950989905112908, + .28867732805385066532, .46473465543376206337e-1, 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., .76293945312500000000e-5, + .13727610943181290891e-3, .11979683091449349286e-2, + .67195313034570709806e-2, .27044920779931968175e-1, + .82472196498517457862e-1, .19570475044896150093, .36391620788543817693, + .52241392782736588032, .54727504974907879912, .34211551468813581183, + -.31580472732719957762e-1, -.32830006549176759667, -.30563797665254420769, + .64905014620683140120e-2, .27642986248995073032, 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., .38146972656250000000e-5, + .72454147007837596854e-4, .66859847582761390285e-3, + .39751311980366118437e-2, .17015198650201528366e-1, + .55443621868993855715e-1, .14157060481641692131, .28641242619559616836, + .45610665490966615415, .55262786406029265394, .45818352706035500108, + .14984403004611673047, -.21163807462970713245, -.36007252928843413718, + -.17030961385712954159, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., .19073486328125000000e-5, + .38135049864067468562e-4, .37101393638555730015e-3, + .23305339886279723213e-2, .10569913448297127219e-1, + .36640175162216897547e-1, .10010476414320235508, .21860074212675559892, + .38124757096345313719, .52020999209879669177, .52172632730659212045, + .30841620620308814614, -.50322546186721500184e-1, -.32577618885114899053, + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., .95367431640625000000e-6, .20021483206955925244e-4, + .20481807322420625431e-3, .13553476938058909882e-2, + .64919676350791905019e-2, .23848725425069251903e-1, + .69384632678886421292e-1, .16249711393618776934, .30736618106830314788, + .46399909601971539157, .53765031034002467225, .42598991476520183929, + .12130445348350215652, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., .47683715820312500000e-6, + .10487707828484902486e-4, .11254146162337528943e-3, + .78248929534271987118e-3, .39468337145306794566e-2, + .15313546659475671763e-1, .47249070825218564146e-1, .11804374107101480543, + .24031796927792491122, .39629215049166341285, .51629108968402548545, + .49622372075429782915, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., .23841857910156250000e-6, + .54823314130625337326e-5, .61575377321535518154e-4, + .44877834366497538134e-3, .23774612048621955857e-2, + .97136347645161687796e-2, .31671599547606636717e-1, + .84028665767000747480e-1, .18298487576742964949, .32647878537696945218, + .46970971486488895077, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., .11920928955078125000e-6, + .28604020001177375838e-5, .33559227978295551013e-4, + .25583821662860610560e-3, .14201552747787302339e-2, + .60938046986874414969e-2, .20930869247951926793e-1, + .58745021125678072911e-1, .13613725780285953720, .26083988356030237586, + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., .59604644775390625000e-7, + .14898180663526043291e-5, .18224991282807693921e-4, + .14504433444608833821e-3, .84184722720281809548e-3, + .37846965430000478789e-2, .13656355548211376864e-1, + .40409541997718853934e-1, .99226988101858325902e-1, 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., .29802322387695312500e-7, .77471708843445529468e-6, + .98649879372606876995e-5, .81814934772838523887e-4, + .49554483992403011328e-3, .23290922072351413938e-2, + .88068134250844034186e-2, .27393666952485719070e-1, 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., 0., .14901161193847656250e-7, .40226235946098233685e-6, + .53236418690561306700e-5, .45933829691164002269e-4, + .28982005232838857913e-3, .14212974043211018374e-2, + .56192363087488842264e-2, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + .74505805969238281250e-8, .20858299254133430408e-6, + .28648457300134381744e-5, .25677535898258910850e-4, + .16849420429491355445e-3, .86062824010315834002e-3, 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., .37252902984619140625e-8, .10801736017613096861e-6, + .15376606719887104015e-5, .14296523739727437959e-4, + .97419023656050887203e-4, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + .18626451492309570312e-8, .55871592916438890146e-7, + .82331193828137454068e-6, .79302250528382787666e-5, 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., .93132257461547851562e-9, + .28867244235852488244e-7, .43982811713864556957e-6, 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., .46566128730773925781e-9, + .14899342093408253335e-7, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., + 0., 0., .23283064365386962891e-9 }; /* Allocates a workspace for the given maximum number of intervals. @@ -1476,7 +1486,7 @@ /* The actual integration routine. */ DEFUN (quadcc, args, nargout, -"-*- texinfo -*-\n\ + "-*- texinfo -*-\n\ @deftypefn {Function File} {@var{q} =} quadcc (@var{f}, @var{a}, @var{b})\n\ @deftypefnx {Function File} {@var{q} =} quadcc (@var{f}, @var{a}, @var{b}, @var{tol})\n\ @deftypefnx {Function File} {@var{q} =} quadcc (@var{f}, @var{a}, @var{b}, @var{tol}, @var{sing})\n\ @@ -1585,12 +1595,12 @@ fcn = args(0).function_value (); else { - std::string fcn_name = unique_symbol_name ("__quadcc_fcn__"); - std::string fname = "function y = "; - fname.append (fcn_name); - fname.append ("(x) y = "); - fcn = extract_function (args(0), "quadcc", fcn_name, fname, - "; endfunction"); + std::string fcn_name = unique_symbol_name ("__quadcc_fcn__"); + std::string fname = "function y = "; + fname.append (fcn_name); + fname.append ("(x) y = "); + fcn = extract_function (args(0), "quadcc", fcn_name, fname, + "; endfunction"); } if (! args(1).is_real_scalar ()) @@ -1788,9 +1798,8 @@ h = (iv->b - iv->a) / 2; #if (DEBUG_QUADCC) - printf - ("quadcc: processing ival %i (of %i) with [%e,%e] int=%e, err=%e, depth=%i\n", - heap[0], nivals, iv->a, iv->b, iv->igral, iv->err, iv->depth); + printf ("quadcc: processing ival %i (of %i) with [%e,%e] int=%e, err=%e, depth=%i\n", + heap[0], nivals, iv->a, iv->b, iv->igral, iv->err, iv->depth); #endif /* Should we try to increase the degree? */ @@ -1891,14 +1900,13 @@ /* Should we drop this interval? */ if ((m + h * xi[0]) >= (m + h * xi[1]) || (m + h * xi[31]) >= (m + h * xi[32]) - || iv->err < fabs (iv->igral) * std::numeric_limits<double>::epsilon () * 10) + || iv->err < fabs (iv->igral) + * std::numeric_limits<double>::epsilon () * 10) { #if (DEBUG_QUADCC) - printf - ("quadcc: dropping ival %i (of %i) with [%e,%e] int=%e, err=%e, depth=%i\n", - heap[0], nivals, iv->a, iv->b, iv->igral, iv->err, - iv->depth); + printf ("quadcc: dropping ival %i (of %i) with [%e,%e] int=%e, err=%e, depth=%i\n", + heap[0], nivals, iv->a, iv->b, iv->igral, iv->err, iv->depth); #endif /* Keep this interval's contribution */ @@ -2207,10 +2215,8 @@ { iv = &(ivals[heap[nivals - 1]]); #if (DEBUG_QUADCC) - printf - ("quadcc: dropping ival %i (of %i) with [%e,%e] int=%e, err=%e, depth=%i\n", - heap[0], nivals, iv->a, iv->b, iv->igral, iv->err, - iv->depth); + printf ("quadcc: dropping ival %i (of %i) with [%e,%e] int=%e, err=%e, depth=%i\n", + heap[0], nivals, iv->a, iv->b, iv->igral, iv->err, iv->depth); #endif err_final += iv->err; igral_final += iv->igral; @@ -2233,10 +2239,9 @@ for (i = 0; i < nivals; i++) { iv = &(ivals[heap[i]]); - printf - ("quadcc: ival %i (%i) with [%e,%e], int=%e, err=%e, depth=%i, rdepth=%i, ndiv=%i\n", - i, heap[i], iv->a, iv->b, iv->igral, iv->err, iv->depth, - iv->rdepth, iv->ndiv); + printf ("quadcc: ival %i (%i) with [%e,%e], int=%e, err=%e, depth=%i, rdepth=%i, ndiv=%i\n", + i, heap[i], iv->a, iv->b, iv->igral, iv->err, iv->depth, + iv->rdepth, iv->ndiv); } #endif @@ -2264,16 +2269,16 @@ %!assert (quadcc (@(x) exp (-x .^ 2), -Inf, Inf), sqrt (pi), 1e-6) %!assert (quadcc (@(x) exp (-x .^ 2), -Inf, 0), sqrt (pi)/2, 1e-6) -## Test function with NaNs in interval +## Test function with NaNs in interval %!function y = __nansin (x) %! nan_locs = [-3*pi/4, -pi/4, 0, pi/3, pi/2, pi]; %! y = sin (x); -%! idx = min (abs (bsxfun (@minus, x(:), nan_locs)), [], 2); +%! idx = min (abs (bsxfun (@minus, x(:), nan_locs)), [], 2); %! y(idx < 1e-10) = NaN; -%!endfunction +%!endfunction %!test -%! [q, err, npoints] = quadcc ("__nansin", -pi, pi); +%! [q, err, npoints] = quadcc ("__nansin", -pi, pi); %! assert (q, 0, 1e-6); %! assert (err, 0, 15*eps);