That's basically what jtroll() does for this case. Thanks,
-- Raul -- Raul On Wed, Jun 14, 2023 at 4:52 PM Henry Rich <henryhr...@gmail.com> wrote: > > 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 <henryhr...@gmail.com> 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 <rauldmil...@gmail.com> 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 > >> <b...@jsoftware.com> 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 ---------------------------------------------------------------------- For information about J forums see http://www.jsoftware.com/forums.htm
