Small improvement: only select useful elements of t
ts 'nthHam 1e6'
0.22246 1314624
> nthHam=: 3 : 0
> 'ln2 ln3 ln5'=. ^. tdv=.2 3 5
> lo=. 0.01 -~ hi=. 1.693-~(*/6,y,ln2, ln3, ln5)^%3
> t=.,:0$0
> for_k. i. 1+ <. hi%ln5 do.
> for_j. i. 1+ <. ln3 %~ hi - p=. k*ln5 do.
> t=. t, j,k,(>.ln2%~lo-q),(<.ln2%~hi-q), q=. p + j * ln3
> end.
> end.
> c=. +/ 1+_2{"1 }. t
> z=.,:0$0
>
*t=. (#~(2&{<:3&{)"1) t**
*
> for_r. t do.
> for_i. ([+i.@>:@-~)/2{.2}.r do. z=.z, i,(2{.r), ({:r)+i*ln2 end.
> end.
> c=. c - y
> assert. 0 <: c
> assert. c <: # z
> __ q: inv tdv ,: x: 3 {. c { (\: {:"1) }. z
> )
>
>
>
--
Met vriendelijke groet,
=@@i
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
