Monte Carlo methods: Difference between revisions

fun dirvcts(): [2][30]inti32 =

fun grayCode(x: inti32): inti32 = (x >> 1) ^ x

fun testBit(n: inti32, ind: inti32): bool =

fun xorInds(n: inti32) (dir_vs: [num_bits]inti32): inti32 =

let reldv_vals = zipWith (fn \ dv i =->

fun sobolIndI (dir_vs: [m][num_bits]inti32, n: inti32): [m]inti32 =

fun sobolIndR(dir_vs: [m][num_bits]inti32) (n: inti32 ): [m]f32 =

in map (fn \ (x: inti32): f32 =-> f32(x) / divisor) arri

fun main(n: inti32): f32 =

let dists = map (fn \xy =->

let (x,y) = (xy[0],xy[1]) in sqrt32F32.sqrt(x*x + y*y))

let bs = map (fn \d =-> if d <= 1.0f32 then 1 else 0) dists