This is the mail archive of the
gsl-discuss@sourceware.org
mailing list for the GSL project.
Re: Improved Mathieu equation eigenvalues and testing
On Wednesday 14 June 2006 12:06 pm, Brian Gough wrote:
> Lowell Johnson writes:
> > I've attached a patch for the Mathieu function eigenvalue accuracy
> > and testing in the specfunc directory. Based on input from Brian
> > Gladman, I've increased the number of terms in the eigenvalue
> > matrix to ensure double precision accuracy.
>
> Thanks, I've added it and updated the documentation for that
> function.
Ah, yes. I overlooked updating the documentation. Thanks Brian.
I've attached an additional change to the documentation describing the
new q argument to gsl_sf_mathieu_workspace().
--
Lowell Johnson
Index: specfunc-mathieu.texi
===================================================================
RCS file: /cvs/gsl/gsl/doc/specfunc-mathieu.texi,v
retrieving revision 1.3
diff -u -r1.3 specfunc-mathieu.texi
--- specfunc-mathieu.texi 14 Jun 2006 17:06:16 -0000 1.3
+++ specfunc-mathieu.texi 15 Jun 2006 01:08:32 -0000
@@ -186,5 +186,12 @@
@deftypefun gsl_sf_mathieu_workspace *gsl_sf_mathieu_alloc (size_t @var{n}, double @var{q});
@deftypefunx void gsl_sf_mathieu_free (gsl_sf_mathieu_workspace @var{*workspace});
-These routines handle memory allocation and freeing for workspace arrays required by the ``_array'' versions of the Mathieu function routines. The variable @var{n} specifies the maximum order of the functions to be computed, not the number of functions to be computed.
+These routines handle memory allocation and freeing for workspace
+arrays required by the ``_array'' versions of the Mathieu function
+routines. The variable @var{n} specifies the maximum order of the
+functions to be computed, not the number of functions to be computed.
+The variable @var{q} specifies the maximum value of @math{q} for which
+Mathieu functions will be calculated. This is required in order to
+properly terminate the infinite eigenvalue matrix for high precision
+solutions.
@end deftypefun