[SciPy-dev] [SciPy-user] calculating factorial

[Travis Oliphant]

Pearu Peterson wrote:

>With the revision 1516 of scipy core arange seems to have ref.counting 
>problem:
>
>  
>
>>>>from scipy import *
>>>>prod(arange(2,101,dtype=object))
>>>>        
>>>>
>Illegal instruction
>
>  
>
>>>>from scipy import *
>>>>a=prod(arange(2,101,dtype=object))
>>>>a
>>>>        
>>>>
>Illegal instruction
>
>  
>
>>>>from scipy import *
>>>>a=arange(2,101,dtype=object)
>>>>b=prod(a)
>>>>b
>>>>        
>>>>
>9.3326215443944102e+157
>  
>
>
>Anyone can reproduce these segfaults?
>  
>
Yes, I can.

We need continued testing of the object arrays.  These can be tricky to 
get right.  Several object-array-based bugs are in Numeric, in fact, 
which led to some not using them.

-Travis