Welcome to mirror list, hosted at ThFree Co, Russian Federation.

git.blender.org/blender.git - Unnamed repository; edit this file 'description' to name the repository.
summaryrefslogtreecommitdiff
path: root/source/blender/blenlib/intern/jitter_2d.c
diff options
context:
space:
mode:
Diffstat (limited to 'source/blender/blenlib/intern/jitter_2d.c')
-rw-r--r--source/blender/blenlib/intern/jitter_2d.c277
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;
+ }
}
generated by cgit v1.2.3 (git 2.25.1) at 2024-07-04 07:59:55 +0300