openscad/dotSCAD
1
0
Fork 0
mirror of https://github.com/JustinSDK/dotSCAD.git synced 2025-08-06 06:47:46 +02:00
Code Issues Releases Wiki Activity

add seed param

Browse Source
This commit is contained in:
Justin Lin
2022-06-30 10:54:22 +08:00
parent 964e6b5020
commit 183b3b74b6
2 changed files with 13 additions and 14 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
src/pp/_impl/_pp_poisson.scad
View File

@@ -1,6 +1,5 @@
use <../../matrix/m_replace.scad> use <../../matrix/m_replace.scad>
use <../../util/every.scad> use <../../util/every.scad>
use <../pp_disk.scad>
function sampling(size, r, start, k) = function sampling(size, r, start, k) =
let( let(
@@ -43,19 +42,19 @@ function sampling_noAdjacentNeighbor(s, sample, x, y) =
norm(sample - p) >= sampling_r(s)) norm(sample - p) >= sampling_r(s))
); );
function sampling_randomSample(s, pos) = function sampling_randomSample(s, pos, seed) =
let( let(
r = sampling_r(s), r = sampling_r(s),
a = rands(0, 360, 1)[0] a = rands(0, 360, 1, seed)[0]
) )
rands(r, 2 * r, 1)[0] * [cos(a), sin(a)] + pos; rands(r, 2 * r, 1, seed + 1)[0] * [cos(a), sin(a)] + pos;
function sampling_hasActive(s) = len(sampling_active(s)) > 0; function sampling_hasActive(s) = len(sampling_active(s)) > 0;
function sampling_kSamples(s, pos) = [ function sampling_kSamples(s, pos, seed) = [
for(n = [0:sampling_k(s) - 1]) for(n = [0:sampling_k(s) - 1])
let( let(
sample = sampling_randomSample(s, pos), sample = sampling_randomSample(s, pos, seed + n),
w = sampling_w(s), w = sampling_w(s),
x = floor(sample.x / w), x = floor(sample.x / w),
y = floor(sample.y / w) y = floor(sample.y / w)
@@ -77,12 +76,12 @@ function _sampling_minDistSample(s, pos, samples, min_s, min_dist, i = 1) =
) )
_sampling_minDistSample(s, pos, samples, params[0], params[1], i + 1); _sampling_minDistSample(s, pos, samples, params[0], params[1], i + 1);
function sampling_trySampleFromOneActive(s) = function sampling_trySampleFromOneActive(s, seed) =
let( let(
active = sampling_active(s), active = sampling_active(s),
i = floor(rands(0, len(active) - 1, 1)[0]), i = floor(rands(0, len(active) - 1, 1, seed)[0]),
pos = active[i], pos = active[i],
samples = sampling_kSamples(s, pos) samples = sampling_kSamples(s, pos, seed)
) )
len(samples) == 0 ? len(samples) == 0 ?
[ [
@@ -115,6 +114,6 @@ function sampling_trySampleFromOneActive(s) =
[each sampling_history(s), [pos, sample]] [each sampling_history(s), [pos, sample]]
]; ];
function _pp_poisson(s) = function _pp_poisson(s, seed, count = 0) =
!sampling_hasActive(s) ? s : _pp_poisson(sampling_trySampleFromOneActive(s)); !sampling_hasActive(s) ? s : _pp_poisson(sampling_trySampleFromOneActive(s, seed + count), seed, count + 1);

6
src/pp/pp_poisson.scad
View File

@@ -1,8 +1,8 @@
use <_impl/_pp_poisson.scad> use <_impl/_pp_poisson.scad>
function pp_poisson(size, r, start, k = 30, history = false) = function pp_poisson(size, r, start, k = 30, seed = undef, history = false) =
let( let(
s = _pp_poisson(sampling(size, r, start, k)), s = _pp_poisson(sampling(size, r, start, k), is_undef(seed) ? floor(rands(0, 1000, 1)[0]) : seed),
samples = [ samples = [
for(row = sampling_grid(s), sample = row) for(row = sampling_grid(s), sample = row)
if(!is_undef(sample)) sample if(!is_undef(sample)) sample
@@ -15,7 +15,7 @@ use <pp/pp_poisson.scad>
use <polyline_join.scad> use <polyline_join.scad>
pts_history = pp_poisson([100, 100], 5, history = true); pts_history = pp_poisson([100, 100], 5, seed = 1, history = true);
for(p = pts_history[0]) { for(p = pts_history[0]) {
translate(p) translate(p)
circle(1); circle(1);