As a workaround, if you just roll 5 numbers from that domain you will get the right answer. That might be an acceptable 'fix' for the source code too.
Henry Rich On Wed, Jun 14, 2023, 10:46 PM Henry Rich <[email protected]> wrote: > The i0 gives domain error if the argument exceeds IMAX, which is what you > are seeing. This looks like a bug to me. I think the algorithm for deal > would work if you just called roll with the extended arguments. Raul, would > you like to take that on? > > Henry Rich > > On Wed, Jun 14, 2023, 6:53 PM Raul Miller <[email protected]> wrote: > >> This error is coming from >> >> #0 jtjsignal (jt=0x7ffff7fc0200, e=3) at ../../../../jsrc/d.c:361 >> #1 0x00007ffff34b131a in jtcvt () from >> /mnt/c/cygwin64/home/15712/jsoftware/jsource/bin/linux/j64avx2/libj.so >> #2 0x00007ffff350173e in jtvi (jt=0x7ffff7fc0200, w=0x4b4000) at >> ../../../../jsrc/u.c:570 >> #3 0x00007ffff3501562 in jti0 (jt=0x7ffff7fc0200, w=0x4b4000) at >> ../../../../jsrc/u.c:238 >> #4 0x00007ffff36fa6a1 in jtdeal (jt=0x7ffff7fc0200, a=0x7ffff3a3a680 >> <Bnum+1088>, w=0x4b4000, >> self=0x7ffff3a3b8c0 <primtab+2176>) at ../../../../jsrc/vrand.c:794 >> >> The issue seems to be the second i0 here >> >> RE(m=i0(a)); RE(c=n=i0(w)); // c starts as max#+1 >> >> In u.c, the comment on i0 is: >> >> // Extract the integer value from w, return it. Set error if >> non-integral or non-atomic. Values whose abs > IMAX are converted to >> IMAX/-IMAX >> >> I haven't yet figured out what's wrong with jti0 in this instance. >> >> FYI, >> >> -- >> Raul >> >> On Wed, Jun 14, 2023 at 11:54 AM 'Veter Kamenev' via Beta >> <[email protected]> wrote: >> > >> > Hi! >> > >> > 5 ? 2980293480239480239480239480239480239482039x >> > >> > |domain error, executing dyad ? >> > |y must be a positive integer >> > >> > I read in NuVoc: >> > >> > y range of random number >> > 0 floating-point value in interval (0,1) >> > 1 0 (always) >> > 2 Boolean >> > >1 integer or extended integer from list i.y >> > >> > https://code.jsoftware.com/wiki/Vocabulary/query#dyadic >> > >> > Best wishes, >> > Veter >> > ---------------------------------------------------------------------- >> > 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
