diff options
Diffstat (limited to 'source/blender/blenlib/intern/jitter_2d.c')
-rw-r--r-- | source/blender/blenlib/intern/jitter_2d.c | 277 |
1 files changed, 143 insertions, 134 deletions
diff --git a/source/blender/blenlib/intern/jitter_2d.c b/source/blender/blenlib/intern/jitter_2d.c index bbf1948be4d..aed43657b3f 100644 --- a/source/blender/blenlib/intern/jitter_2d.c +++ b/source/blender/blenlib/intern/jitter_2d.c @@ -31,148 +31,157 @@ #include "BLI_strict_flags.h" - void BLI_jitterate1(float (*jit1)[2], float (*jit2)[2], int num, float rad1) { - int i, j, k; - float vecx, vecy, dvecx, dvecy, x, y, len; - - for (i = num - 1; i >= 0; i--) { - dvecx = dvecy = 0.0; - x = jit1[i][0]; - y = jit1[i][1]; - for (j = num - 1; j >= 0; j--) { - if (i != j) { - vecx = jit1[j][0] - x - 1.0f; - vecy = jit1[j][1] - y - 1.0f; - for (k = 3; k > 0; k--) { - if (fabsf(vecx) < rad1 && fabsf(vecy) < rad1) { - len = sqrtf(vecx * vecx + vecy * vecy); - if (len > 0 && len < rad1) { - len = len / rad1; - dvecx += vecx / len; - dvecy += vecy / len; - } - } - vecx += 1.0f; - - if (fabsf(vecx) < rad1 && fabsf(vecy) < rad1) { - len = sqrtf(vecx * vecx + vecy * vecy); - if (len > 0 && len < rad1) { - len = len / rad1; - dvecx += vecx / len; - dvecy += vecy / len; - } - } - vecx += 1.0f; - - if (fabsf(vecx) < rad1 && fabsf(vecy) < rad1) { - len = sqrtf(vecx * vecx + vecy * vecy); - if (len > 0 && len < rad1) { - len = len / rad1; - dvecx += vecx / len; - dvecy += vecy / len; - } - } - vecx -= 2.0f; - vecy += 1.0f; - } - } - } - - x -= dvecx / 18.0f; - y -= dvecy / 18.0f; - x -= floorf(x); - y -= floorf(y); - jit2[i][0] = x; - jit2[i][1] = y; - } - memcpy(jit1, jit2, 2 * (unsigned int)num * sizeof(float)); + int i, j, k; + float vecx, vecy, dvecx, dvecy, x, y, len; + + for (i = num - 1; i >= 0; i--) { + dvecx = dvecy = 0.0; + x = jit1[i][0]; + y = jit1[i][1]; + for (j = num - 1; j >= 0; j--) { + if (i != j) { + vecx = jit1[j][0] - x - 1.0f; + vecy = jit1[j][1] - y - 1.0f; + for (k = 3; k > 0; k--) { + if (fabsf(vecx) < rad1 && fabsf(vecy) < rad1) { + len = sqrtf(vecx * vecx + vecy * vecy); + if (len > 0 && len < rad1) { + len = len / rad1; + dvecx += vecx / len; + dvecy += vecy / len; + } + } + vecx += 1.0f; + + if (fabsf(vecx) < rad1 && fabsf(vecy) < rad1) { + len = sqrtf(vecx * vecx + vecy * vecy); + if (len > 0 && len < rad1) { + len = len / rad1; + dvecx += vecx / len; + dvecy += vecy / len; + } + } + vecx += 1.0f; + + if (fabsf(vecx) < rad1 && fabsf(vecy) < rad1) { + len = sqrtf(vecx * vecx + vecy * vecy); + if (len > 0 && len < rad1) { + len = len / rad1; + dvecx += vecx / len; + dvecy += vecy / len; + } + } + vecx -= 2.0f; + vecy += 1.0f; + } + } + } + + x -= dvecx / 18.0f; + y -= dvecy / 18.0f; + x -= floorf(x); + y -= floorf(y); + jit2[i][0] = x; + jit2[i][1] = y; + } + memcpy(jit1, jit2, 2 * (unsigned int)num * sizeof(float)); } void BLI_jitterate2(float (*jit1)[2], float (*jit2)[2], int num, float rad2) { - int i, j; - float vecx, vecy, dvecx, dvecy, x, y; - - for (i = num - 1; i >= 0; i--) { - dvecx = dvecy = 0.0; - x = jit1[i][0]; - y = jit1[i][1]; - for (j = num - 1; j >= 0; j--) { - if (i != j) { - vecx = jit1[j][0] - x - 1.0f; - vecy = jit1[j][1] - y - 1.0f; - - if (fabsf(vecx) < rad2) { dvecx += vecx * rad2; } - vecx += 1.0f; - if (fabsf(vecx) < rad2) { dvecx += vecx * rad2; } - vecx += 1.0f; - if (fabsf(vecx) < rad2) { dvecx += vecx * rad2; } - - if (fabsf(vecy) < rad2) { dvecy += vecy * rad2; } - vecy += 1.0f; - if (fabsf(vecy) < rad2) { dvecy += vecy * rad2; } - vecy += 1.0f; - if (fabsf(vecy) < rad2) { dvecy += vecy * rad2; } - - } - } - - x -= dvecx / 2.0f; - y -= dvecy / 2.0f; - x -= floorf(x); - y -= floorf(y); - jit2[i][0] = x; - jit2[i][1] = y; - } - memcpy(jit1, jit2, (unsigned int)num * sizeof(float[2])); + int i, j; + float vecx, vecy, dvecx, dvecy, x, y; + + for (i = num - 1; i >= 0; i--) { + dvecx = dvecy = 0.0; + x = jit1[i][0]; + y = jit1[i][1]; + for (j = num - 1; j >= 0; j--) { + if (i != j) { + vecx = jit1[j][0] - x - 1.0f; + vecy = jit1[j][1] - y - 1.0f; + + if (fabsf(vecx) < rad2) { + dvecx += vecx * rad2; + } + vecx += 1.0f; + if (fabsf(vecx) < rad2) { + dvecx += vecx * rad2; + } + vecx += 1.0f; + if (fabsf(vecx) < rad2) { + dvecx += vecx * rad2; + } + + if (fabsf(vecy) < rad2) { + dvecy += vecy * rad2; + } + vecy += 1.0f; + if (fabsf(vecy) < rad2) { + dvecy += vecy * rad2; + } + vecy += 1.0f; + if (fabsf(vecy) < rad2) { + dvecy += vecy * rad2; + } + } + } + + x -= dvecx / 2.0f; + y -= dvecy / 2.0f; + x -= floorf(x); + y -= floorf(y); + jit2[i][0] = x; + jit2[i][1] = y; + } + memcpy(jit1, jit2, (unsigned int)num * sizeof(float[2])); } - void BLI_jitter_init(float (*jitarr)[2], int num) { - float (*jit2)[2]; - float num_fl, num_fl_sqrt; - float x, rad1, rad2, rad3; - RNG *rng; - int i; - - if (num == 0) { - return; - } - - num_fl = (float)num; - num_fl_sqrt = sqrtf(num_fl); - - jit2 = MEM_mallocN(12 + (unsigned int)num * sizeof(float[2]), "initjit"); - rad1 = 1.0f / num_fl_sqrt; - rad2 = 1.0f / num_fl; - rad3 = num_fl_sqrt / num_fl; - - rng = BLI_rng_new(31415926 + (unsigned int)num); - - x = 0; - for (i = 0; i < num; i++) { - jitarr[i][0] = x + rad1 * (float)(0.5 - BLI_rng_get_double(rng)); - jitarr[i][1] = (float)i / num_fl + rad1 * (float)(0.5 - BLI_rng_get_double(rng)); - x += rad3; - x -= floorf(x); - } - - BLI_rng_free(rng); - - for (i = 0; i < 24; i++) { - BLI_jitterate1(jitarr, jit2, num, rad1); - BLI_jitterate1(jitarr, jit2, num, rad1); - BLI_jitterate2(jitarr, jit2, num, rad2); - } - - MEM_freeN(jit2); - - /* finally, move jittertab to be centered around (0, 0) */ - for (i = 0; i < num; i++) { - jitarr[i][0] -= 0.5f; - jitarr[i][1] -= 0.5f; - } + float(*jit2)[2]; + float num_fl, num_fl_sqrt; + float x, rad1, rad2, rad3; + RNG *rng; + int i; + + if (num == 0) { + return; + } + + num_fl = (float)num; + num_fl_sqrt = sqrtf(num_fl); + + jit2 = MEM_mallocN(12 + (unsigned int)num * sizeof(float[2]), "initjit"); + rad1 = 1.0f / num_fl_sqrt; + rad2 = 1.0f / num_fl; + rad3 = num_fl_sqrt / num_fl; + + rng = BLI_rng_new(31415926 + (unsigned int)num); + + x = 0; + for (i = 0; i < num; i++) { + jitarr[i][0] = x + rad1 * (float)(0.5 - BLI_rng_get_double(rng)); + jitarr[i][1] = (float)i / num_fl + rad1 * (float)(0.5 - BLI_rng_get_double(rng)); + x += rad3; + x -= floorf(x); + } + + BLI_rng_free(rng); + + for (i = 0; i < 24; i++) { + BLI_jitterate1(jitarr, jit2, num, rad1); + BLI_jitterate1(jitarr, jit2, num, rad1); + BLI_jitterate2(jitarr, jit2, num, rad2); + } + + MEM_freeN(jit2); + + /* finally, move jittertab to be centered around (0, 0) */ + for (i = 0; i < num; i++) { + jitarr[i][0] -= 0.5f; + jitarr[i][1] -= 0.5f; + } } |