which c complier and version you used to build binaries? On Sun, Jul 15, 2018, 7:02 AM Bill Heagy <[email protected]> wrote:
> JVERSION > Engine: j806/j32/linux > Beta: GPL3/2018-07-11T14:18:06 > Library: 8.05.14 > Platform: Linux 32 > Installer: unknown > InstallPath: /home/wheagy/tmp/jbld/j32 > Contact: [email protected] > > (current git) > But: it's been going on for awhile, showing up as a error in g410. I'm > pretty sure I tried it in one of the recent beta releases. It's in > intel 32 and 64nonavx. > > On 2018-07-14 06:54 PM, bill lam wrote: > > 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 > > > ---------------------------------------------------------------------- > For information about J forums see http://www.jsoftware.com/forums.htm ---------------------------------------------------------------------- For information about J forums see http://www.jsoftware.com/forums.htm
