[Jprogramming] bug in latest beta and source

[Henry Rich]

Good shooting!  I introduced that bug in a recent change.  I will fix it 
for the next release.

Thank you for taking the time to run this down.

Henry Rich

On 7/14/2018 6:50 PM, Bill Heagy wrote:
> It looks like a problem with "^":
>    2^0.5
> 1.41421
>    2^0.5j3
> 1.41421
> For complex exponent with a real part of 0.5 (sqrt), it appears to 
> ignore the imaginary part.
>
> On 2018-07-02 04:52 PM, Bill Heagy wrote:
>> I haven't checked it all out, but I suspect the problem is that the 
>> Stirling formula is asymmptotic, and has a larger error near 1.
>>
>> On 2018-07-02 12:20 PM, Bill Heagy wrote:
>>> First four runs on 32bit linux:
>>>
>>> 1j7.26
>>>
>>> 1j8.53
>>>
>>> 1j_3.14
>>>
>>> 1j_2.67
>>>
>>> 64bit machine gives different imaginary parts, but real part=1.  
>>> Looks like 1j(whatever) may be causing the problem.
>>>
>>> On 2018-07-02 07:45 AM, Henry Rich wrote:
>>>> When this fails, what is the value in z that gives the error?
>>>>
>>>> BTW, Bill, that error with gcc that you reported a few weeks back 
>>>> turned out to be a program error that showed up only on that 
>>>> compiler.  If you switch back to gcc now it should be gone.
>>>>
>>>> Henry Rich
>>>>
>>>> On 6/29/2018 9:59 PM, Bill Heagy wrote:
>>>>> If I run test/g410.ijs often enough, I get:
>>>>>
>>>>> 100 RUN2 'g410'
>>>>> .......
>>>>> 1 1 1 1 1 1 1 1
>>>>>    (g = gamma) x=: _4.1 _4.9 4.1 4.9 j./_0.4 _0.5 0.4 0.5 _4 _5 4 5
>>>>> 1 1 1 1 1 1 1 1
>>>>> 1 1 1 1 1 1 1 1
>>>>> 1 1 1 1 1 1 1 1
>>>>> 1 1 1 1 1 1 1 1
>>>>>
>>>>>    (g =!.1e_12 gamma) 1.5
>>>>> 1
>>>>>    (g =!.1e_12 gamma) 1r3 1r2 2r3
>>>>> 1 1 1
>>>>>    (g =!.1e_12 gamma) x=:   1e_2*_7e2+?2 20$2e3
>>>>> 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
>>>>> 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
>>>>>    (g =!.1e_12 gamma) x=:   3e_2*_7e2+?2 20$2e3
>>>>> 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
>>>>> 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
>>>>>    (g =!.1e_12 gamma) z=:j./1e_2*_7e2+?2 20$2e3
>>>>> 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
>>>>> |assertion failure
>>>>> |[-103] /home/wheagy/git/jsource/test/g410.ijs
>>>>>
>>>>> Which specific value is wrong changes from failure to failure, but 
>>>>> it's always this line, and it appears to happen about every 50 to 
>>>>> 100 iterations.  This happens on debian linux 32 and 64nonavx, 
>>>>> with the latest beta release and when compiled from the latest 
>>>>> version on github.
>>>>>
>>>>> ---------------------------------------------------------------------- 
>>>>>
>>>>> For information about J forums see 
>>>>> http://www.jsoftware.com/forums.htm
>>>>
>>>>
>>>> ---
>>>> This email has been checked for viruses by AVG.
>>>> https://www.avg.com
>>>>
>>>> ----------------------------------------------------------------------
>>>> For information about J forums see http://www.jsoftware.com/forums.htm
>>>
>> ----------------------------------------------------------------------
>> For information about J forums see http://www.jsoftware.com/forums.htm
> ----------------------------------------------------------------------
> For information about J forums see http://www.jsoftware.com/forums.htm