what's the JVERSION On Sun, Jul 15, 2018, 6:50 AM Bill Heagy <[email protected]> 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 ---------------------------------------------------------------------- For information about J forums see http://www.jsoftware.com/forums.htm
