Dear gsl team,
FFTs functions won't work along with "atof()". It report nan. if comment
the atof,the reusult is correct.
I attached the test codes and results as following,pls check the
attachments.
fft_right.c and fft_right.result.txt: with no atof function
fft_wrong.c and fft_wrong.result.txt: with atof function
GSL version :gsl-1.15
OS: windowx xp sp3
$ gcc -v
Using built-in specs.
Target: i686-pc-cygwin
Configured with:
/gnu/gcc/releases/respins/4.3.4-4/gcc4-4.3.4-4/src/gcc-4.3.4/configure
--srcdir=/gnu/gcc/releases/respins/4.3.4-4/gcc4-4.3.4-4/src/gcc-4.3.4
--prefix=/usr --exec-p
refix=/usr --bindir=/usr/bin --sbindir=/usr/sbin --libexecdir=/usr/lib
--datadir=/usr/share --localstatedir=/var --sysconfdir=/etc
--infodir=/usr/share/info --mandir=/usr/share/man
-C --datadir=/usr/share --infodir=/usr/share/info --mandir=/usr/share/man
-v --with-gmp=/usr --with-mpfr=/usr --enable-bootstrap
--enable-version-specific-runtime-libs --with-slib
dir=/usr/bin --libexecdir=/usr/lib --enable-static --enable-shared
--enable-shared-libgcc --disable-__cxa_atexit --with-gnu-ld --with-gnu-as
--with-dwarf2 --disable-sjlj-exceptions
--enable-languages=ada,c,c++,fortran,java,objc,obj-c++ --disable-symvers
--enable-libjava --program-suffix=-4 --enable-libgomp --enable-libssp
--enable-libada --enable-threads=pos
ix --with-arch=i686 --with-tune=generic --enable-libgcj-sublibs CC=gcc-4
CXX=g++-4 CC_FOR_TARGET=gcc-4 CXX_FOR_TARGET=g++-4
GNATMAKE_FOR_TARGET=gnatmake GNATBIND_FOR_TARGET=gnatbin
d --with-ecj-jar=/usr/share/java/ecj.jar
Thread model: posix
gcc version 4.3.4 20090804 (release) 1 (GCC)
Thanks a lot!
tangnengfu
(0,2.000000,0.000000) (512.000000,0.000000)
(1,0.000000,2.342592) (-0.000000,0.000000)
(2,0.000000,-1.402860) (0.000000,-0.000000)
(3,0.000000,-2.665310) (-0.000000,-0.000000)
(4,0.000000,-1.720267) (-0.000000,-0.000000)
(5,0.000000,-1.366667) (0.000000,-0.000000)
(6,0.000000,-0.179836) (0.000000,-0.000000)
(7,0.000000,-0.923201) (0.000000,-0.000000)
(8,0.000000,0.236348) (0.000000,0.000000)
(9,0.000000,1.806965) (-0.000000,0.000000)
(10,0.000000,-2.560757) (-0.000000,-0.000000)
(11,0.000000,-0.796709) (0.000000,-0.000000)
(12,0.000000,-0.026343) (0.000000,-0.000000)
(13,0.000000,-2.970171) (-0.000000,-0.000000)
(14,0.000000,-2.408778) (-0.000000,-0.000000)
(15,0.000000,-1.003780) (0.000000,-0.000000)
(16,0.000000,2.406060) (-0.000000,0.000000)
(17,0.000000,-0.326647) (0.000000,-0.000000)
(18,0.000000,2.953758) (-0.000000,0.000000)
(19,0.000000,0.526526) (0.000000,0.000000)
(20,0.000000,-2.748617) (-0.000000,-0.000000)
(21,0.000000,-1.045799) (0.000000,-0.000000)
(22,0.000000,-2.124860) (-0.000000,-0.000000)
(23,0.000000,0.820339) (0.000000,0.000000)
(24,0.000000,2.224442) (-0.000000,0.000000)
(25,0.000000,-1.825120) (-0.000000,-0.000000)
(26,0.000000,-0.928210) (0.000000,-0.000000)
(27,0.000000,1.480166) (0.000000,0.000000)
(28,0.000000,0.063368) (0.000000,0.000000)
(29,0.000000,2.224384) (-0.000000,0.000000)
(30,0.000000,2.899432) (-0.000000,0.000000)
(31,0.000000,2.487401) (-0.000000,0.000000)
(32,0.000000,1.366641) (0.000000,0.000000)
(33,0.000000,-1.620792) (-0.000000,-0.000000)
(34,0.000000,2.584878) (-0.000000,0.000000)
(35,0.000000,0.219265) (0.000000,0.000000)
(36,0.000000,2.217033) (-0.000000,0.000000)
(37,0.000000,-2.813820) (-0.000000,-0.000000)
(38,0.000000,-1.579736) (-0.000000,-0.000000)
(39,0.000000,2.878791) (-0.000000,0.000000)
(40,0.000000,0.202657) (0.000000,0.000000)
(41,0.000000,-3.136902) (-0.000000,-0.000000)
(42,0.000000,-1.953077) (-0.000000,-0.000000)
(43,0.000000,-0.405102) (0.000000,-0.000000)
(44,0.000000,1.923063) (-0.000000,0.000000)
(45,0.000000,2.325884) (-0.000000,0.000000)
(46,0.000000,-2.689597) (-0.000000,-0.000000)
(47,0.000000,0.433801) (0.000000,0.000000)
(48,0.000000,-3.099869) (-0.000000,-0.000000)
(49,0.000000,2.605426) (-0.000000,0.000000)
(50,3.000000,-0.523599) (332.553755,-192.000000)
(51,0.000000,-0.041907) (0.000000,-0.000000)
(52,0.000000,1.972211) (-0.000000,0.000000)
(53,0.000000,-1.415571) (0.000000,-0.000000)
(54,0.000000,-1.622563) (-0.000000,-0.000000)
(55,0.000000,0.388399) (0.000000,0.000000)
(56,0.000000,-1.295537) (0.000000,-0.000000)
(57,0.000000,0.232754) (0.000000,0.000000)
(58,0.000000,-0.119390) (0.000000,-0.000000)
(59,0.000000,-0.836448) (0.000000,-0.000000)
(60,0.000000,0.106039) (0.000000,0.000000)
(61,0.000000,0.335717) (0.000000,0.000000)
(62,0.000000,0.232767) (0.000000,0.000000)
(63,0.000000,1.169596) (0.000000,0.000000)
(64,0.000000,2.158799) (-0.000000,0.000000)
(65,0.000000,0.912857) (0.000000,0.000000)
(66,0.000000,2.306111) (-0.000000,0.000000)
(67,0.000000,-2.024289) (-0.000000,-0.000000)
(68,0.000000,1.815164) (-0.000000,0.000000)
(69,0.000000,-0.909173) (0.000000,-0.000000)
(70,0.000000,0.802963) (0.000000,0.000000)
(71,0.000000,-1.590531) (-0.000000,-0.000000)
(72,0.000000,-0.293299) (0.000000,-0.000000)
(73,0.000000,1.007679) (0.000000,0.000000)
(74,0.000000,3.090184) (-0.000000,0.000000)
(75,1.500000,1.570796) (0.000000,192.000000)
(76,0.000000,1.727560) (-0.000000,0.000000)
(77,0.000000,1.964301) (-0.000000,0.000000)
(78,0.000000,2.479745) (-0.000000,0.000000)
(79,0.000000,2.960755) (-0.000000,0.000000)
(80,0.000000,2.696683) (-0.000000,0.000000)
(81,0.000000,1.927000) (-0.000000,0.000000)
(82,0.000000,0.810314) (0.000000,0.000000)
(83,0.000000,2.709527) (-0.000000,0.000000)
(84,0.000000,2.912513) (-0.000000,0.000000)
(85,0.000000,-3.029752) (-0.000000,-0.000000)
(86,0.000000,-2.822517) (-0.000000,-0.000000)
(87,0.000000,-1.505556) (0.000000,-0.000000)
(88,0.000000,-1.770020) (-0.000000,-0.000000)
(89,0.000000,-1.860235) (-0.000000,-0.000000)
(90,0.000000,0.169014) (0.000000,0.000000)
(91,0.000000,-1.566537) (0.000000,-0.000000)
(92,0.000000,-2.641377) (-0.000000,-0.000000)
(93,0.000000,3.112215) (-0.000000,0.000000)
(94,0.000000,-1.911670) (-0.000000,-0.000000)
(95,0.000000,-2.581341) (-0.000000,-0.000000)
(96,0.000000,2.745210) (-0.000000,0.000000)
(97,0.000000,-1.963588) (-0.000000,-0.000000)
(98,0.000000,-2.119725) (-0.000000,-0.000000)
(99,0.000000,-0.962289) (0.000000,-0.000000)
(100,0.000000,-0.270373) (0.000000,-0.000000)
(101,0.000000,0.115352) (0.000000,0.000000)
(102,0.000000,-0.230430) (0.000000,-0.000000)
(103,0.000000,-0.119516) (0.000000,-0.000000)
(104,0.000000,0.179433) (0.000000,0.000000)
(105,0.000000,-0.128136) (0.000000,-0.000000)
(106,0.000000,1.168754) (0.000000,0.000000)
(107,0.000000,1.292020) (0.000000,0.000000)
(108,0.000000,0.886994) (0.000000,0.000000)
(109,0.000000,-2.641457) (-0.000000,-0.000000)
(110,0.000000,-2.483948) (-0.000000,-0.000000)
(111,0.000000,-1.778750) (-0.000000,-0.000000)
(112,0.000000,0.959296) (0.000000,0.000000)
(113,0.000000,-1.187197) (0.000000,-0.000000)
(114,0.000000,-0.244979) (0.000000,-0.000000)
(115,0.000000,1.680666) (-0.000000,0.000000)
(116,0.000000,-1.977667) (-0.000000,-0.000000)
(117,0.000000,-0.429199) (0.000000,-0.000000)
(118,0.000000,-0.019562) (0.000000,-0.000000)
(119,0.000000,1.302736) (0.000000,0.000000)
(120,0.000000,-0.474844) (0.000000,-0.000000)
(121,0.000000,1.125235) (0.000000,0.000000)
(122,0.000000,1.105999) (0.000000,0.000000)
(123,0.000000,2.713335) (-0.000000,0.000000)
(124,0.000000,2.619380) (-0.000000,0.000000)
(125,0.000000,2.205329) (-0.000000,0.000000)
(126,0.000000,1.785670) (-0.000000,0.000000)
(127,0.000000,2.965093) (-0.000000,0.000000)
(128,0.000000,3.141593) (-0.000000,0.000000)
(129,0.000000,-2.972852) (-0.000000,-0.000000)
(130,0.000000,-1.942269) (-0.000000,-0.000000)
(131,0.000000,-2.181582) (-0.000000,-0.000000)
(132,0.000000,-2.619380) (-0.000000,-0.000000)
(133,0.000000,-2.684619) (-0.000000,-0.000000)
(134,0.000000,-1.114298) (0.000000,-0.000000)
(135,0.000000,-1.119580) (0.000000,-0.000000)
(136,0.000000,0.474844) (0.000000,0.000000)
(137,0.000000,-1.303381) (0.000000,-0.000000)
(138,0.000000,0.013642) (0.000000,0.000000)
(139,0.000000,0.403394) (0.000000,0.000000)
(140,0.000000,1.977667) (-0.000000,0.000000)
(141,0.000000,-1.677028) (-0.000000,-0.000000)
(142,0.000000,0.187429) (0.000000,0.000000)
(143,0.000000,1.187038) (0.000000,0.000000)
(144,0.000000,-0.959296) (0.000000,-0.000000)
(145,0.000000,1.792000) (-0.000000,0.000000)
(146,0.000000,2.673193) (-0.000000,0.000000)
(147,0.000000,2.628886) (-0.000000,0.000000)
(148,0.000000,-0.886994) (0.000000,-0.000000)
(149,0.000000,-1.282714) (0.000000,-0.000000)
(150,0.000000,-1.062003) (0.000000,-0.000000)
(151,0.000000,0.135043) (0.000000,0.000000)
(152,0.000000,-0.179433) (0.000000,-0.000000)
(153,0.000000,0.138329) (0.000000,0.000000)
(154,0.000000,0.260717) (0.000000,0.000000)
(155,0.000000,-0.118840) (0.000000,-0.000000)
(156,0.000000,0.270373) (0.000000,0.000000)
(157,0.000000,0.958397) (0.000000,0.000000)
(158,0.000000,2.141444) (-0.000000,0.000000)
(159,0.000000,1.965901) (-0.000000,0.000000)
(160,0.000000,-2.745210) (-0.000000,-0.000000)
(161,0.000000,2.588956) (-0.000000,0.000000)
(162,0.000000,1.877786) (-0.000000,0.000000)
(163,0.000000,-3.089263) (-0.000000,-0.000000)
(164,0.000000,2.641377) (-0.000000,0.000000)
(165,0.000000,1.612590) (-0.000000,0.000000)
(166,0.000000,-0.131392) (0.000000,-0.000000)
(167,0.000000,1.840695) (-0.000000,0.000000)
(168,0.000000,1.770020) (-0.000000,0.000000)
(169,0.000000,1.505560) (0.000000,0.000000)
(170,0.000000,2.819771) (-0.000000,0.000000)
(171,0.000000,3.018499) (-0.000000,0.000000)
(172,0.000000,-2.912513) (-0.000000,-0.000000)
(173,0.000000,-2.703269) (-0.000000,-0.000000)
(174,0.000000,-0.767704) (0.000000,-0.000000)
(175,0.000000,-1.923675) (-0.000000,-0.000000)
(176,0.000000,-2.696683) (-0.000000,-0.000000)
(177,0.000000,-2.946122) (-0.000000,-0.000000)
(178,0.000000,-2.480549) (-0.000000,-0.000000)
(179,0.000000,-1.966050) (-0.000000,-0.000000)
(180,0.000000,-1.727560) (-0.000000,-0.000000)
(181,1.500000,-1.570796) (0.000000,-192.000000)
(182,0.000000,-3.095517) (-0.000000,-0.000000)
(183,0.000000,-1.008675) (0.000000,-0.000000)
(184,0.000000,0.293299) (0.000000,0.000000)
(185,0.000000,1.603428) (-0.000000,0.000000)
(186,0.000000,-0.783631) (0.000000,-0.000000)
(187,0.000000,0.909562) (0.000000,0.000000)
(188,0.000000,-1.815164) (-0.000000,-0.000000)
(189,0.000000,2.023715) (-0.000000,0.000000)
(190,0.000000,-2.251144) (-0.000000,-0.000000)
(191,0.000000,-0.904550) (0.000000,-0.000000)
(192,0.000000,-2.158799) (-0.000000,-0.000000)
(193,0.000000,-1.166312) (0.000000,-0.000000)
(194,0.000000,-0.188383) (0.000000,-0.000000)
(195,0.000000,-0.341397) (0.000000,-0.000000)
(196,0.000000,-0.106039) (0.000000,-0.000000)
(197,0.000000,0.852613) (0.000000,0.000000)
(198,0.000000,0.074254) (0.000000,0.000000)
(199,0.000000,-0.243233) (0.000000,-0.000000)
(200,0.000000,1.295537) (0.000000,0.000000)
(201,0.000000,-0.400668) (0.000000,-0.000000)
(202,0.000000,1.579965) (-0.000000,0.000000)
(203,0.000000,1.371667) (0.000000,0.000000)
(204,0.000000,-1.972211) (-0.000000,-0.000000)
(205,0.000000,0.040838) (0.000000,0.000000)
(206,3.000000,0.523599) (332.553755,192.000000)
(207,0.000000,-2.606836) (-0.000000,-0.000000)
(208,0.000000,3.099869) (-0.000000,0.000000)
(209,0.000000,-0.433291) (0.000000,-0.000000)
(210,0.000000,2.711553) (-0.000000,0.000000)
(211,0.000000,-2.330728) (-0.000000,-0.000000)
(212,0.000000,-1.923063) (-0.000000,-0.000000)
(213,0.000000,0.398873) (0.000000,0.000000)
(214,0.000000,1.970030) (-0.000000,0.000000)
(215,0.000000,-3.135765) (-0.000000,-0.000000)
(216,0.000000,-0.202657) (0.000000,-0.000000)
(217,0.000000,-2.893503) (-0.000000,-0.000000)
(218,0.000000,1.705998) (-0.000000,0.000000)
(219,0.000000,2.762704) (-0.000000,0.000000)
(220,0.000000,-2.217033) (-0.000000,-0.000000)
(221,0.000000,-0.185707) (0.000000,-0.000000)
(222,0.000000,-2.584274) (-0.000000,-0.000000)
(223,0.000000,1.624720) (-0.000000,0.000000)
(224,0.000000,-1.366641) (0.000000,-0.000000)
(225,0.000000,-2.499339) (-0.000000,-0.000000)
(226,0.000000,-2.795656) (-0.000000,-0.000000)
(227,0.000000,-2.222984) (-0.000000,-0.000000)
(228,0.000000,-0.063368) (0.000000,-0.000000)
(229,0.000000,-1.473438) (0.000000,-0.000000)
(230,0.000000,0.933669) (0.000000,0.000000)
(231,0.000000,1.815275) (-0.000000,0.000000)
(232,0.000000,-2.224442) (-0.000000,-0.000000)
(233,0.000000,-0.820274) (0.000000,-0.000000)
(234,0.000000,2.085325) (-0.000000,0.000000)
(235,0.000000,1.034485) (0.000000,0.000000)
(236,0.000000,2.748617) (-0.000000,0.000000)
(237,0.000000,-0.530435) (0.000000,-0.000000)
(238,0.000000,-2.910894) (-0.000000,-0.000000)
(239,0.000000,0.324915) (0.000000,0.000000)
(240,0.000000,-2.406060) (-0.000000,-0.000000)
(241,0.000000,1.006896) (0.000000,0.000000)
(242,0.000000,2.446854) (-0.000000,0.000000)
(243,0.000000,2.968481) (-0.000000,0.000000)
(244,0.000000,0.026343) (0.000000,0.000000)
(245,0.000000,0.875047) (0.000000,0.000000)
(246,0.000000,2.604181) (-0.000000,0.000000)
(247,0.000000,-1.808277) (-0.000000,-0.000000)
(248,0.000000,-0.236348) (0.000000,-0.000000)
(249,0.000000,0.921046) (0.000000,0.000000)
(250,0.000000,0.082198) (0.000000,0.000000)
(251,0.000000,1.325935) (0.000000,0.000000)
(252,0.000000,1.720267) (-0.000000,0.000000)
(253,0.000000,2.672567) (-0.000000,0.000000)
(254,0.000000,1.389691) (0.000000,0.000000)
(255,0.000000,-2.340182) (-0.000000,-0.000000)
#include <stdio.h>
#include <math.h>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_fft_complex.h>
#define REAL(z,i) ((z)[2*(i)])
#define IMAG(z,i) ((z)[2*(i)+1])
int main (void)
{
int i;
const int n = 256;
double data[2*n];
double omega1=50;
double omega2=75;
double phase1=-M_PI*30.0/180.0;
double phase2=M_PI*90.0/180.0;
double amplitude0=2.0;
double amplitude1=3.0;
double amplitude2=1.50;
double t=0.0;
gsl_fft_complex_wavetable * wavetable;
gsl_fft_complex_workspace * workspace;
for (i = 0; i < n; i++)
{
double fs =2.0*M_PI;
double fn = fs/n;
t= fn*i;
REAL(data,i) = amplitude0+amplitude1*cos(omega1*t+phase1)+amplitude2*cos(omega2*t+phase2);
IMAG(data,i) = 0.0;
char str[64]="123.45";
double test = atof(str);
}
wavetable = gsl_fft_complex_wavetable_alloc (n);
workspace = gsl_fft_complex_workspace_alloc (n);
gsl_fft_complex_forward (data, 1, n, wavetable, workspace);
for (i = 0; i < n; i++)
{
double R = sqrt(REAL(data,i)*REAL(data,i) + IMAG(data,i)*IMAG(data,i));
double amplitude =0;
double phase = 0.0;
phase = atan2(IMAG(data,i),REAL(data,i));
if(0 == i)
{
amplitude = R/n;
}
else
{
amplitude = R/n*2;
}
printf ("(%d,%f,%f) (%f,%f)\n", i, amplitude,phase,REAL(data,i),IMAG(data,i));
}
gsl_fft_complex_wavetable_free (wavetable);
gsl_fft_complex_workspace_free (workspace);
return 0;
}#include <stdio.h>
#include <math.h>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_fft_complex.h>
#define REAL(z,i) ((z)[2*(i)])
#define IMAG(z,i) ((z)[2*(i)+1])
int main (void)
{
int i;
const int n = 256;
double data[2*n];
double omega1=50;
double omega2=75;
double phase1=-M_PI*30.0/180.0;
double phase2=M_PI*90.0/180.0;
double amplitude0=2.0;
double amplitude1=3.0;
double amplitude2=1.50;
double t=0.0;
gsl_fft_complex_wavetable * wavetable;
gsl_fft_complex_workspace * workspace;
for (i = 0; i < n; i++)
{
double fs =2.0*M_PI;
double fn = fs/n;
t= fn*i;
REAL(data,i) = amplitude0+amplitude1*cos(omega1*t+phase1)+amplitude2*cos(omega2*t+phase2);
IMAG(data,i) = 0.0;
char str[64]="123.45";
//double test = atof(str);
}
wavetable = gsl_fft_complex_wavetable_alloc (n);
workspace = gsl_fft_complex_workspace_alloc (n);
gsl_fft_complex_forward (data, 1, n, wavetable, workspace);
for (i = 0; i < n; i++)
{
double R = sqrt(REAL(data,i)*REAL(data,i) + IMAG(data,i)*IMAG(data,i));
double amplitude =0;
double phase = 0.0;
phase = atan2(IMAG(data,i),REAL(data,i));
if(0 == i)
{
amplitude = R/n;
}
else
{
amplitude = R/n*2;
}
printf ("(%d,%f,%f) (%f,%f)\n", i, amplitude,phase,REAL(data,i),IMAG(data,i));
}
gsl_fft_complex_wavetable_free (wavetable);
gsl_fft_complex_workspace_free (workspace);
return 0;
}(0,nan,nan) (nan,nan)
(1,nan,nan) (nan,nan)
(2,nan,nan) (nan,nan)
(3,nan,nan) (nan,nan)
(4,nan,nan) (nan,nan)
(5,nan,nan) (nan,nan)
(6,nan,nan) (nan,nan)
(7,nan,nan) (nan,nan)
(8,nan,nan) (nan,nan)
(9,nan,nan) (nan,nan)
(10,nan,nan) (nan,nan)
(11,nan,nan) (nan,nan)
(12,nan,nan) (nan,nan)
(13,nan,nan) (nan,nan)
(14,nan,nan) (nan,nan)
(15,nan,nan) (nan,nan)
(16,nan,nan) (nan,nan)
(17,nan,nan) (nan,nan)
(18,nan,nan) (nan,nan)
(19,nan,nan) (nan,nan)
(20,nan,nan) (nan,nan)
(21,nan,nan) (nan,nan)
(22,nan,nan) (nan,nan)
(23,nan,nan) (nan,nan)
(24,nan,nan) (nan,nan)
(25,nan,nan) (nan,nan)
(26,nan,nan) (nan,nan)
(27,nan,nan) (nan,nan)
(28,nan,nan) (nan,nan)
(29,nan,nan) (nan,nan)
(30,nan,nan) (nan,nan)
(31,nan,nan) (nan,nan)
(32,nan,nan) (nan,nan)
(33,nan,nan) (nan,nan)
(34,nan,nan) (nan,nan)
(35,nan,nan) (nan,nan)
(36,nan,nan) (nan,nan)
(37,nan,nan) (nan,nan)
(38,nan,nan) (nan,nan)
(39,nan,nan) (nan,nan)
(40,nan,nan) (nan,nan)
(41,nan,nan) (nan,nan)
(42,nan,nan) (nan,nan)
(43,nan,nan) (nan,nan)
(44,nan,nan) (nan,nan)
(45,nan,nan) (nan,nan)
(46,nan,nan) (nan,nan)
(47,nan,nan) (nan,nan)
(48,nan,nan) (nan,nan)
(49,nan,nan) (nan,nan)
(50,nan,nan) (nan,nan)
(51,nan,nan) (nan,nan)
(52,nan,nan) (nan,nan)
(53,nan,nan) (nan,nan)
(54,nan,nan) (nan,nan)
(55,nan,nan) (nan,nan)
(56,nan,nan) (nan,nan)
(57,nan,nan) (nan,nan)
(58,nan,nan) (nan,nan)
(59,nan,nan) (nan,nan)
(60,nan,nan) (nan,nan)
(61,nan,nan) (nan,nan)
(62,nan,nan) (nan,nan)
(63,nan,nan) (nan,nan)
(64,nan,nan) (nan,nan)
(65,nan,nan) (nan,nan)
(66,nan,nan) (nan,nan)
(67,nan,nan) (nan,nan)
(68,nan,nan) (nan,nan)
(69,nan,nan) (nan,nan)
(70,nan,nan) (nan,nan)
(71,nan,nan) (nan,nan)
(72,nan,nan) (nan,nan)
(73,nan,nan) (nan,nan)
(74,nan,nan) (nan,nan)
(75,nan,nan) (nan,nan)
(76,nan,nan) (nan,nan)
(77,nan,nan) (nan,nan)
(78,nan,nan) (nan,nan)
(79,nan,nan) (nan,nan)
(80,nan,nan) (nan,nan)
(81,nan,nan) (nan,nan)
(82,nan,nan) (nan,nan)
(83,nan,nan) (nan,nan)
(84,nan,nan) (nan,nan)
(85,nan,nan) (nan,nan)
(86,nan,nan) (nan,nan)
(87,nan,nan) (nan,nan)
(88,nan,nan) (nan,nan)
(89,nan,nan) (nan,nan)
(90,nan,nan) (nan,nan)
(91,nan,nan) (nan,nan)
(92,nan,nan) (nan,nan)
(93,nan,nan) (nan,nan)
(94,nan,nan) (nan,nan)
(95,nan,nan) (nan,nan)
(96,nan,nan) (nan,nan)
(97,nan,nan) (nan,nan)
(98,nan,nan) (nan,nan)
(99,nan,nan) (nan,nan)
(100,nan,nan) (nan,nan)
(101,nan,nan) (nan,nan)
(102,nan,nan) (nan,nan)
(103,nan,nan) (nan,nan)
(104,nan,nan) (nan,nan)
(105,nan,nan) (nan,nan)
(106,nan,nan) (nan,nan)
(107,nan,nan) (nan,nan)
(108,nan,nan) (nan,nan)
(109,nan,nan) (nan,nan)
(110,nan,nan) (nan,nan)
(111,nan,nan) (nan,nan)
(112,nan,nan) (nan,nan)
(113,nan,nan) (nan,nan)
(114,nan,nan) (nan,nan)
(115,nan,nan) (nan,nan)
(116,nan,nan) (nan,nan)
(117,nan,nan) (nan,nan)
(118,nan,nan) (nan,nan)
(119,nan,nan) (nan,nan)
(120,nan,nan) (nan,nan)
(121,nan,nan) (nan,nan)
(122,nan,nan) (nan,nan)
(123,nan,nan) (nan,nan)
(124,nan,nan) (nan,nan)
(125,nan,nan) (nan,nan)
(126,nan,nan) (nan,nan)
(127,nan,nan) (nan,nan)
(128,nan,nan) (nan,nan)
(129,nan,nan) (nan,nan)
(130,nan,nan) (nan,nan)
(131,nan,nan) (nan,nan)
(132,nan,nan) (nan,nan)
(133,nan,nan) (nan,nan)
(134,nan,nan) (nan,nan)
(135,nan,nan) (nan,nan)
(136,nan,nan) (nan,nan)
(137,nan,nan) (nan,nan)
(138,nan,nan) (nan,nan)
(139,nan,nan) (nan,nan)
(140,nan,nan) (nan,nan)
(141,nan,nan) (nan,nan)
(142,nan,nan) (nan,nan)
(143,nan,nan) (nan,nan)
(144,nan,nan) (nan,nan)
(145,nan,nan) (nan,nan)
(146,nan,nan) (nan,nan)
(147,nan,nan) (nan,nan)
(148,nan,nan) (nan,nan)
(149,nan,nan) (nan,nan)
(150,nan,nan) (nan,nan)
(151,nan,nan) (nan,nan)
(152,nan,nan) (nan,nan)
(153,nan,nan) (nan,nan)
(154,nan,nan) (nan,nan)
(155,nan,nan) (nan,nan)
(156,nan,nan) (nan,nan)
(157,nan,nan) (nan,nan)
(158,nan,nan) (nan,nan)
(159,nan,nan) (nan,nan)
(160,nan,nan) (nan,nan)
(161,nan,nan) (nan,nan)
(162,nan,nan) (nan,nan)
(163,nan,nan) (nan,nan)
(164,nan,nan) (nan,nan)
(165,nan,nan) (nan,nan)
(166,nan,nan) (nan,nan)
(167,nan,nan) (nan,nan)
(168,nan,nan) (nan,nan)
(169,nan,nan) (nan,nan)
(170,nan,nan) (nan,nan)
(171,nan,nan) (nan,nan)
(172,nan,nan) (nan,nan)
(173,nan,nan) (nan,nan)
(174,nan,nan) (nan,nan)
(175,nan,nan) (nan,nan)
(176,nan,nan) (nan,nan)
(177,nan,nan) (nan,nan)
(178,nan,nan) (nan,nan)
(179,nan,nan) (nan,nan)
(180,nan,nan) (nan,nan)
(181,nan,nan) (nan,nan)
(182,nan,nan) (nan,nan)
(183,nan,nan) (nan,nan)
(184,nan,nan) (nan,nan)
(185,nan,nan) (nan,nan)
(186,nan,nan) (nan,nan)
(187,nan,nan) (nan,nan)
(188,nan,nan) (nan,nan)
(189,nan,nan) (nan,nan)
(190,nan,nan) (nan,nan)
(191,nan,nan) (nan,nan)
(192,nan,nan) (nan,nan)
(193,nan,nan) (nan,nan)
(194,nan,nan) (nan,nan)
(195,nan,nan) (nan,nan)
(196,nan,nan) (nan,nan)
(197,nan,nan) (nan,nan)
(198,nan,nan) (nan,nan)
(199,nan,nan) (nan,nan)
(200,nan,nan) (nan,nan)
(201,nan,nan) (nan,nan)
(202,nan,nan) (nan,nan)
(203,nan,nan) (nan,nan)
(204,nan,nan) (nan,nan)
(205,nan,nan) (nan,nan)
(206,nan,nan) (nan,nan)
(207,nan,nan) (nan,nan)
(208,nan,nan) (nan,nan)
(209,nan,nan) (nan,nan)
(210,nan,nan) (nan,nan)
(211,nan,nan) (nan,nan)
(212,nan,nan) (nan,nan)
(213,nan,nan) (nan,nan)
(214,nan,nan) (nan,nan)
(215,nan,nan) (nan,nan)
(216,nan,nan) (nan,nan)
(217,nan,nan) (nan,nan)
(218,nan,nan) (nan,nan)
(219,nan,nan) (nan,nan)
(220,nan,nan) (nan,nan)
(221,nan,nan) (nan,nan)
(222,nan,nan) (nan,nan)
(223,nan,nan) (nan,nan)
(224,nan,nan) (nan,nan)
(225,nan,nan) (nan,nan)
(226,nan,nan) (nan,nan)
(227,nan,nan) (nan,nan)
(228,nan,nan) (nan,nan)
(229,nan,nan) (nan,nan)
(230,nan,nan) (nan,nan)
(231,nan,nan) (nan,nan)
(232,nan,nan) (nan,nan)
(233,nan,nan) (nan,nan)
(234,nan,nan) (nan,nan)
(235,nan,nan) (nan,nan)
(236,nan,nan) (nan,nan)
(237,nan,nan) (nan,nan)
(238,nan,nan) (nan,nan)
(239,nan,nan) (nan,nan)
(240,nan,nan) (nan,nan)
(241,nan,nan) (nan,nan)
(242,nan,nan) (nan,nan)
(243,nan,nan) (nan,nan)
(244,nan,nan) (nan,nan)
(245,nan,nan) (nan,nan)
(246,nan,nan) (nan,nan)
(247,nan,nan) (nan,nan)
(248,nan,nan) (nan,nan)
(249,nan,nan) (nan,nan)
(250,nan,nan) (nan,nan)
(251,nan,nan) (nan,nan)
(252,nan,nan) (nan,nan)
(253,nan,nan) (nan,nan)
(254,nan,nan) (nan,nan)
(255,nan,nan) (nan,nan)
_______________________________________________
Bug-gsl mailing list
[email protected]
https://lists.gnu.org/mailman/listinfo/bug-gsl
