Re: [Bug-gsl] gsl_ran_discrete_preproc "bug"

[Srećko Jurić-Kavelj]

Hi,

I've just tested the patch on Versalogic's Cobra (Pioneer 3-DX onboard
comp.) and it works.

As for the sampling results in the test case, they seem to be sound,
but I haven't tested thoroughly.

$ ./randist_test
0 52
1 265
2 152
3 20
4 63
5 145
6 25
7 142
8 94
9 221

Thanks.

2009/5/14 Brian Gough <address@hidden>:
> At Mon, 4 May 2009 15:52:49 +0200,
> Srećko Jurić-Kavelj wrote:
>> The problem is with bigs and smalls stack (randist/discrete.c).
>> Although the same code is used when calculating the size of stacks and
>> when inserting (that's when the problem manifests) in them it seems it
>> works differently (FPU unit?). I've solved this temporarily with
>> allocating stacks of fixed size (Kevents).
>
> Could you try the attached patch and see if it fixes the problem (I
> think I see how the probelm occurs but I wasn't able to reproduce it
> locally on x86 or x86_64).  Thanks.
>
> --
> Brian Gough
> (GSL Maintainer)
>
> Support freedom by joining the FSF!
> http://www.fsf.org/associate/support_freedom/join_fsf?referrer=37
>
>



-- 
Srećko Jurić-Kavelj
+385913444837

Sanctus Eusebius Hieronymus: "Parce mihi, Domine, quia dalmata sum!"