There would be a non 0 probability of duplicates high enough to occur to
someone in next 50 years.
An algortithm that keeps the randomness generator of keeping the first
selection intact in a deal that is perhaps useful at least when there is low
probability of duplicates:
( ] (<@] (,~<) i.@#@] + {~) L:0 a:"_`(<@I.)@.(1 < +/)"1@:=(-.&a:@:)) 2 3 2 3 4
5 2
┌─────────────┬─────────┐
│┌─────┬─────┐│┌───┬───┐│
││2 3 4│0 2 6│││3 4│1 3││
│└─────┴─────┘│└───┴───┘│
└─────────────┴─────────┘
finds repeated indexes and increments replacements for duplicates, joins the
replacement and indexes in overall box, suitable for amendv verb below to scan/
replace all duplicates. Recursively applies (^:_) until no duplicates remain.
amendve =: 4 : 0 :: ] NB. x is boxed or itemed to be replace_values,indexes. y
is data to amend.
'a b' =. x
a b}y
)
amendv =: amendve :: ] NB. error no change useful for J corner cases in
building empty x
dealx =: 4 : '(<@] (,~ |.) ] (<@] (,~<) i.@#@] (y | +) {~) L:0 a:"_`(<@I.)@.(1
< +/)"1@:=(-.&a:@:))(amendv each/@:)(>@)^:_ x ?@:# y'
14 dealx 14
9 2 1 12 13 8 10 6 11 0 3 4 5 7
2 dealx 2980293480239480239480239480239480239482039x
2944004533308676397283493309982083700473234
612729966346708897309198185520497570771366
On Wednesday, June 14, 2023 at 04:52:12 p.m. EDT, 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
