- Fix use of id in *formula* sketch (change to index).

 - Change to `+=` syntax in *Number_003_A_Ngan*'s `circle.py`.
 - First works-as-expected commit for *Number_003_B_Ngan*.

Author: tildebyte

Whitney Artport - (Software) Structures/#003/A/Interpretation/Number_003_A_Ngan/circle.py

@@ -96,8 +96,8 @@ def intersect(self, other):
         return True
 
     def repel(self, angle):
-        self.x = self.x + cos(angle) / 4
-        self.y = self.y + sin(angle) / 4
+        self.x += cos(angle) / 4
+        self.y += sin(angle) / 4
 
     def renderHair(self):
         for hair in self.hairs:

Whitney Artport - (Software) Structures/#003/B/Interpretation/Number_003_B_Ngan/Number_003_B_Ngan.pde

@@ -14,7 +14,7 @@
 
 
 float[][] field;
-int[][] fieldShade;
+// int[][] fieldShade;
 int gaph, gapv;
 int marginh, marginv; // margin
 float cnt = 0;
@@ -34,7 +34,7 @@ int cTimer = 0;
 void setup()
 {
   size( 500, 500 );
-  frameRate( 30 );
+  // frameRate( 30 );
 
   gaph = 3;
   gapv = 3;
@@ -52,19 +52,19 @@ void setup()
 
   // field
   field = new float[(width-marginh*2)/gaph][(height-marginv*2)/gapv];
-  fieldShade = new int[field.length][field[0].length];
+  // fieldShade = new int[field.length][field[0].length];
 
   for (int i=0; i<field.length; i++) {
     for (int k=0; k<field[0].length; k++) {
       field[i][k] = PI2-PI/3;
-      fieldShade[i][k] = 1;
+      // fieldShade[i][k] = 1;
     }
   }
 
   circles = new Circle[100];
 
-  ellipseMode( CENTER );
-  noFill();
+  // ellipseMode( CENTER );
+  // noFill();
 }
 
 void draw()
@@ -90,11 +90,11 @@ void draw()
     }
   }
 
-  noStroke();
-  noFill();
+  // noStroke();
+  // noFill();
 
   for (int i =0; i<cCounter; i++) {
-    circles[i].render();
+    circles[i].move();
     circles[i].getGrid();
   }
 
@@ -158,7 +158,7 @@ class Circle
   }
 
   void move() {
-    float angle = sin( sp1 ) - cos(sp2);
+    float angle = getSin( sp1 ) - getCos(sp2);
 
     sp1+=ac1;
     sp2+=ac2;
@@ -182,8 +182,8 @@ class Circle
   }
 
   void repel( float angle ) {
-    x = x + getCos(angle)/10;
-    y = y + getSin(angle)/10;
+    x = x + getCos(angle)*0.1;
+    y = y + getSin(angle)*0.1;
   }
 
   void setState( float px, float py ) {
@@ -199,7 +199,7 @@ class Circle
         float dx = circles[i].x - x;
         float dy = circles[i].y - y;
         float drr = dx*dx + dy*dy;
-        float brr = circles[i].rr + 2*circles[i].r*r + rr;
+        // float brr = circles[i].rr + 2*circles[i].r*r + rr;
         float d = sqrt( drr );
 
         if ( d>r+circles[i].r || d<abs(r-circles[i].r) ) continue; // no solution
@@ -228,22 +228,22 @@ class Circle
     int numy = (int)floor(d/gapv);
 
     int[] pos;
-    float dx, dy, ang;
+    // float dx, dy, ang;
 
     for (int i=sx; i<sx+numx; i++) {
       if (i>=0 && i<field.length) {
         for (int k=sy; k<sy+numy; k++) {
           if (k>=0 && k<field[0].length) {
             if (over) {
               pos = getLocation(i, k);
-              dx = pos[0]-x;
-              dy = pos[1]-y;
+              // dx = pos[0]-x;
+              // dy = pos[1]-y;
 
               if (dist( x, y, pos[0], pos[1]) < r) {
                 float da = atan2( pos[1]-iny, pos[0]-inx );
                 if (field[i][k] < da) field[i][k] += PI/20;
                 else if (field[i][k] > da) field[i][k] -= PI/20;
-                fieldShade[i][k] = 2;
+                // fieldShade[i][k] = 2;
               }
             }
           }

Whitney Artport - (Software) Structures/#003/B/Interpretation/Number_003_B_Ngan/Number_003_B_Ngan.pyde

@@ -11,203 +11,34 @@ Processing v.68 <http://processing.org>
 
 Port to Processing.py/Processing 2.0 by Ben Alkov 5 September 2014
 '''
+from circle import Circle
+from grid import Grid
 
-field
-fieldShade
-gaph
-gapv
-marginh, marginv # margin
-cnt = 0
-sintable = 628
-costable = 628
-PI2 = 2 * PI
-tiltAngle = 0
-circles = Circle[]
+
+grid = None
+NumCircles = 10
+circles = []
 counter = 0
-cCounter = 0
-cTimer = 0
+CirclePosition = 250
+Radius = 40
 
 
 def setup():
+    global grid
     size(500, 500)
-    frameRate(30)
-    gaph = 3
-    gapv = 3
-    marginh = 20
-    marginv = 20
-
-    # lookup table
-    for i = 0 in range( i < sintable.length): # i++
-        sintable[i] = sin(i / 100.0)
-    for i = 0 in range( i < costable.length): # i++
-        costable[i] = cos(i / 100.0)
-
-    # field
-    field = (width - marginh * 2) / gaph][(height - marginv * 2) / gapv]
-    fieldShade = field.length][field[0].length]
-    for i = 0 in range( i < field.length): # i++
-        for k = 0 in range( k < field[0].length): # k++
-            field[i][k] = PI2 - PI / 3
-            fieldShade[i][k] = 1
-    circles = Circle[100]
-    ellipseMode(CENTER)
-    noFill()
+    # Init here because of sketch `width` & `height` being undefined prior to
+    # this point.
+    grid = Grid(TAU - PI / 3, 20, 3)
 
 
 def draw():
+    global counter
     background(50)
     stroke(255, 255, 255, 50)
-    ax, ay
-    cTimer++
-    if cTimer > 5 and cCounter < circles.length - 1:
-        circles[cCounter] = Circle(250, 250, 40, cCounter)
-        cCounter++
-        cTimer = 0
-    len = 10
-    for i = 0 in range( i < field.length): # i++
-        for k = 0 in range( k < field[0].length): # k++
-            ax = marginh + i * gaph
-            ay = marginv + k * gapv
-            line(ax, ay, ax + len * getCos(field[i][k]), ay + len * getSin(field[i][k]))
-    noStroke()
-    noFill()
-    for i =0 in range( i < cCounter): # i++
-        circles[i].render()
-        circles[i].getGrid()
-
-
-def getLocation(i, k):
-    return marginh + i*gaph, marginv + k*gapv
-
-
-def getSin(val):
-    if val < 0:
-        val = 6.28 + val
-    if val >= 6.28:
-        val -= 6.28
-    val = min(6.27, max(0, val))
-    return sintable[floor(val * 100)]
-
-
-def getCos(val):
-    if val < 0:
-        val = 6.28 + val
-    if val >= 6.28:
-        val -= 6.28
-    val = min(6.27, max(0, val))
-    return costable[floor(val * 100)]
-
-
-class Circle(object):
-    x
-    y
-    r
-    d
-    rr
-    ac1
-    ac2
-    ac3
-    sp1
-    sp2
-    sp3
-    id
-    inx = 0
-    iny = 0
-    over = True
-
-    def __init__(self, px, py, pr, id):
-        x = px
-        y = py
-        r = pr
-        d = r * 2
-        rr = r * r
-        this.id = id
-        sp1 = random(2)
-        sp2 = random(2)
-        sp3 = random(2)
-        ac1 = random(0.5) - random(0.5)
-        ac2 = random(0.5) - random(0.5)
-        ac3 = random(0.5) - random(0.5)
-
-
-    def render(self):
-        move()
-        ellipse(x, y, d, d)
-
-
-    def move(self):
-        angle = sin(sp1) - cos(sp2)
-        sp1 += ac1
-        sp2 += ac2
-        sp3 += ac3
-        angle = (angle < 0) ? angle + PI2: ((angle >= PI2) ? angle - PI2: angle)
-        x = x + getSin(angle)
-        y = y + getCos(angle)
-        checkBounds()
-        checkOverlap()
-
-    def checkBounds(self):
-        if x > width:
-            x = 0
-        if x < 0:
-            x = width
-        if y > height:
-            y = 0
-        if y < 0:
-            y = height
-
-
-    def repel(self, angle):
-        x = x + getCos(angle) / 10
-        y = y + getSin(angle) / 10
-
-
-    def setState(self, px, py):
-        inx = px
-        iny = py
-        over = True
-
-
-    def checkOverlap(self):
-        for i = id + 1 in range( i < cCounter): # i++
-            if i != id:
-                dx = circles[i].x - x
-                dy = circles[i].y - y
-                drr = dx * dx + dy * dy
-                brr = circles[i].rr + 2 * circles[i].r * r + rr
-                d = sqrt(drr)
-                if d > r + circles[i].r or d < abs(r - circles[i].r):
-                    continue # no solution
-                ang = atan2(dy, dx)
-                repel(ang + PI)
-                circles[i].repel(ang)
-                a = (rr - circles[i].rr + drr) / (2 * d)
-                h = sqrt(rr - a * a)
-                x2 = x + a * (circles[i].x - x) / d
-                y2 = y + a * (circles[i].y - y) / d
-                setState(x2, y2)
-                circles[i].setState(x2, y2)
-
-
-    def getGrid(self):
-        sx = ceil((x - r - marginh) / gaph)
-        sy = ceil((y - r - marginv) / gapv)
-        numx = floor(d / gaph)
-        numy = floor(d / gapv)
-        pos
-        dx, dy, ang
-        for i = sx in range( i < sx + numx): # i++
-            if i >= 0 and i < field.length:
-                for k = sy in range( k < sy + numy): # k++
-                    if k >= 0 and k < field[0].length:
-                        if over:
-                            pos = getLocation(i, k)
-                            dx = pos[0] - x
-                            dy = pos[1] - y
-                            if dist(x, y, pos[0], pos[1]) < r:
-                                da = atan2(pos[1] - iny, pos[0] - inx)
-                                if field[i][k] < da:
-                                    field[i][k] += PI / 20
-                                elif field[i][k] > da:
-                                    fieldShade[i][k] = 2
-        over = False
+    if frameCount % 5 == 0 and counter < NumCircles:
+        circles.append(Circle(CirclePosition, CirclePosition, Radius, counter))
+        counter += 1
+    grid.update()
+    for circle in circles:
+        circle.move(circles)
+        circle.getGrid(grid)