bug-gsl
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Bug-gsl] questionable results for: gsl_sf_ellint_Kcomp --- redux


From: Gerald I. Evenden
Subject: [Bug-gsl] questionable results for: gsl_sf_ellint_Kcomp --- redux
Date: Sat, 24 Jun 2006 21:43:57 -0400
User-agent: KMail/1.8.2

Regarding my first submission I found that if I took the root of the arguments 
then the results agreed with the Abramowitz tables exactly:

#include <gsl_sf_ellint.h>
#include <stdio.h>
#include <math.h>
    int
main(int argc, char **argv) {
    double m;

    for (m = 0.01; m <=0.14; m+=0.01)
        printf("%.2f %.15f %.15f\n", m,
            gsl_sf_ellint_Kcomp(sqrt(m), GSL_PREC_DOUBLE),    <---------
                                            -------
            gsl_sf_ellint_Kcomp(sqrt(1. - m), GSL_PREC_DOUBLE));  <-------
                                            -----------
}

This does not make any sense to me!!

-- 
Jerry and the low-riders: Daisy Mae and Joshua
"Cogito cogito ergo cogito sum"
   Ambrose Bierce, The Devil's Dictionary




reply via email to

[Prev in Thread] Current Thread [Next in Thread]