gammai in 2.0.8

classic Classic list List threaded Threaded
3 messages Options
Reply | Threaded
Open this post in threaded view
|

gammai in 2.0.8

Ted.Harding
Hi Folks,

I have found that Kurt Hornik's gammai(a,x) function seg-faults
when run in my octave-2.0.8.

Have other people noticed this or should I have a look at my hardware?

(for gammai(3,.5), octave-core gives a 109x1 vector A with values from
0.125 down to 4.68191016873419e-215).

Best wishes,
Ted.                                    [hidden email]


Reply | Threaded
Open this post in threaded view
|

Re: gammai in 2.0.8

Mario Storti-4

> Hi Folks,
>
> I have found that Kurt Hornik's gammai(a,x) function seg-faults
> when run in my octave-2.0.8.
>
> Have other people noticed this or should I have a look at my hardware?
>
> (for gammai(3,.5), octave-core gives a 109x1 vector A with values from
> 0.125 down to 4.68191016873419e-215).
>
> Best wishes,
> Ted.                                    [hidden email]

It seems to work well in 2.0.9 (Linux  1.2.13  i586):

> Octave, version 2.0.9 (i586-pc-linux-gnu).
> Copyright (C) 1996, 1997 John W. Eaton.
> This is free software with ABSOLUTELY NO WARRANTY.
> For details, type `warranty'.
>
> octave> gammai(3,.5)
> ans = 0.014388
> octave>

Mario

%%%%%%<>%%%%%%<>%%%%%%<>%%%%%%<>%%%%%%<>%%%%%%<>%%%%%%<>%%%%%
Mario Alberto Storti               | Fax: (54)(42) 55.09.44 |
Grupo de Tecnologia Mecanica       | Tel: (54)(42) 55.91.75 |
INTEC, Guemes 3450 - 3000 Santa Fe | Home: Gob. Vera 3161   |
Argentina                          |       (54)(42) 55.00.23|
Reply: [hidden email]  |                        |
                http://venus.unl.edu.ar/gtm-eng.html        |


Reply | Threaded
Open this post in threaded view
|

Re: gammai in 2.0.8

Ted.Harding
> > I have found that Kurt Hornik's gammai(a,x) function seg-faults
> > when run in my octave-2.0.8.
>
> It seems to work well in 2.0.9 (Linux  1.2.13  i586):
>
> > Octave, version 2.0.9 (i586-pc-linux-gnu).

Thanks, Mario. If it's any clue, I have since found that bessel
gives a FORTRAN exception:

    octave-2:1> bessel(1.3,5)
     ***MESSAGE FROM ROUTINE DGAMMA IN LIBRARY SLATEC.
     ***FATAL ERROR, PROGRAM ABORTED.
     *  X IS 0
     *  ERROR NUMBER = 4
     *
     ***END OF MESSAGE

     ***JOB ABORT DUE TO FATAL ERROR.
    error: exception encountered in Fortran subroutine xdgamma

However, my mileage varies:

    octave-2:1> bessel(.5,5)
    ans = -0.34217
    octave-2:2> bessel(3.5,2.5)
    error: Segmentation fault -- stopping myself...
    attempting to save variables to `octave-core'...
    save to `octave-core' complete

Ted.                                    [hidden email]