Show Posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.


Messages - MasterGy

Pages: [1] 2 3 ... 22
1
Programs / Re: Garland of pearls
« on: April 05, 2022, 03:32:44 pm »
I added another curiosity. (but I haven't covered the point yet)
I recently had a question with _MAPTRIANGLE that calculates 2d points. I wanted to know so I could move with a mouse in 3d. This has now been achieved with a simple solution.
Stretched mesh in 3d. Walkable with WASD + mouse like in 3D games. If you press the Q key, you can grab the point with the mouse and move it. You can nail it with the SPACE key.



Code: QB64: [Select]
  1. mouse_sens = .5
  2. tightness = 1
  3. gravity = .005
  4. lass = .94
  5. power = -0.3
  6. adj = .8 'moving speed 3d point mouse
  7.  
  8.  
  10. CONST pip180 = 3.141592 / 180
  11. DIM SHARED monx, mony: monx = 800: mony = monx / _DESKTOPWIDTH * _DESKTOPHEIGHT: mon = _NEWIMAGE(monx, mony, 32): SCREEN mon: _FULLSCREEN: _MOUSEHIDE
  12.  
  13. DIM SHARED points_c_arr, conn_c_arr
  14. points_c_arr = 2000
  15. conn_c_arr = 2000
  16.  
  17. DIM SHARED points(points_c_arr, 19), conn(conn_c_arr, 9), cam(19), me(19), conn_c, points_c
  18.  
  19.  
  20. make_grid 1, -10, -10, 10, 10, -10, 0, -10, 10, 0, 10, 10, 0, 1.5
  21.  
  22.  
  23.  
  24.  
  25.  
  27.  
  28.     'control
  29.     mousex = 0: mousey = 0: mw = 0: WHILE _MOUSEINPUT: mousex = mousex + _MOUSEMOVEMENTX: mousey = mousey + _MOUSEMOVEMENTY: mw = mw + _MOUSEWHEEL: WEND
  30.     k_left = _KEYDOWN(19200) OR _KEYDOWN(ASC("a")): k_right = _KEYDOWN(19712) OR _KEYDOWN(ASC("d")): inkey2$ = INKEY$: IF inkey2$ = CHR$(27) THEN END
  31.     k_up = _KEYDOWN(18432) OR _KEYDOWN(ASC("w")): k_down = _KEYDOWN(20480) OR _KEYDOWN(ASC("s")): key_take = _KEYDOWN(ASC("q"))
  32.     me(2) = me(2) + mw: IF key_take THEN _MOUSESHOW ELSE _MOUSEHIDE
  33.     IF key_take = 0 THEN me(3) = me(3) - mousex * mouse_sens: me(4) = me(4) + mousey * mouse_sens: IF ABS(me(4)) > 80 THEN me(4) = 80 * SGN(me(4))
  34.     go = -.5 * ABS(k_left OR k_right OR k_up OR k_down): go_ang = -90 * (k_left - k_right) + (me(3)) - 90
  35.     ang1 = go_ang * pip180: ang2 = (90 - me(4)) * pip180: IF k_down THEN go = -go
  36.     me(0) = me(0) + SIN(ang2) * COS(ang1) * go: me(1) = me(1) + SIN(ang2) * SIN(ang1) * go: me(2) = me(2) + COS(ang2) * go
  37.     FOR t = 0 TO 9: cam(t) = me(t): NEXT t: cam(13) = cam(3) * pip180: cam(14) = cam(4) * pip180
  38.  
  39.  
  40.  
  41.     'calculate moving
  42.     FOR a = 0 TO points_c_arr - 1: FOR t = 0 TO 2: points(a, 11 + t) = points(a, t): NEXT t
  43.         IF (points(a, 14) = 0 OR points(a, 6) = 0) OR a = take_moving THEN _CONTINUE
  44.  
  45.         'connections vector
  46.         vec(0) = 0: vec(1) = 0: vec(2) = 0: vec_c = 0
  47.         FOR t = 0 TO conn_c - 1: IF conn(t, 3) = 0 THEN _CONTINUE
  48.             x = -1: FOR t2 = 0 TO 1: IF conn(t, t2) = a THEN x = conn(t, t2 XOR 1)
  49.             NEXT t2
  50.             IF x <> -1 THEN
  51.                 FOR av = 0 TO 2: disv = points(a, av) - points(x, av): vec(av) = vec(av) + ((ABS(disv) - conn(t, 4 + av))) * SGN(disv) * power
  52.                 NEXT av: vec_c = vec_c + 1
  53.             END IF
  54.         NEXT t
  55.  
  56.         FOR av = 0 TO 2: points(a, 3 + av) = (points(a, 3 + av) + vec(av) / vec_c) * lass: NEXT av: points(a, 4) = points(a, 4) - gravity
  57.         FOR av = 0 TO 2: points(a, 11 + av) = points(a, av) + points(a, 3 + av): NEXT av
  58.     NEXT a
  59.  
  60.     FOR a = 0 TO points_c_arr - 1: FOR t = 0 TO 2: points(a, t) = points(a, 11 + t): NEXT t, a
  61.  
  62.  
  63.  
  64.     'calculating points
  65.     FOR ap = 0 TO points_c_arr - 1: IF points(ap, 6) = 0 THEN _CONTINUE
  66.         points(ap, 8) = points(ap, 0): points(ap, 9) = points(ap, 1): points(ap, 10) = points(ap, 2)
  67.         rotate_display points(ap, 8), points(ap, 9), points(ap, 10)
  68.         FOR t = 0 TO 2: points(ap, 15 + t) = points(ap, 8 + t): NEXT t: point3d_to_pset points(ap, 15), points(ap, 16), points(ap, 17)
  69.     NEXT ap
  70.  
  71.     'draw connections
  72.     FOR ac = 0 TO conn_c_arr - 1: IF conn(ac, 3) = 0 THEN _CONTINUE
  73.         FOR t1 = 0 TO 1: FOR t2 = 0 TO 2: coo(t1, t2) = points(conn(ac, t1), 15 + t2): NEXT t2, t1
  74.         IF points(conn(ac, 0), 10) < -1 AND points(conn(ac, 1), 10) < -1 THEN
  75.             LINE (coo(0, 0), coo(0, 1))-(coo(1, 0), coo(1, 1)), _RGB32(coo(0, 2), coo(0, 2), coo(0, 2))
  76.         END IF
  77.     NEXT ac
  78.  
  79.     'fix points draw signal
  80.     FOR ap = 0 TO points_c_arr - 1: IF points(ap, 6) = 0 OR points(ap, 14) OR points(ap, 10) > -1 THEN _CONTINUE
  81.         CIRCLE (points(ap, 15), points(ap, 16)), 12000 / points(ap, 10), _RGB32(points(ap, 17), points(ap, 17), points(ap, 17))
  82.     NEXT ap
  83.  
  84.     'take point with mouse
  85.     IF key_take AND take_moving = -1 THEN
  86.         take_sens = monx * 0.01: mindis = 999999: take = -1
  87.         FOR ap = 0 TO points_c_arr - 1: IF points(ap, 6) = 0 OR points(ap, 10) > -1 THEN _CONTINUE
  88.             disx = ABS(_MOUSEX - points(ap, 15)): IF disx > take_sens THEN _CONTINUE
  89.             disy = ABS(_MOUSEY - points(ap, 16)): IF disy > take_sens THEN _CONTINUE
  90.             disv = dis(disx, disy, 0): IF disv > take_sens OR disv > mindis THEN _CONTINUE
  91.         take = ap: mindis = disv: NEXT ap: IF take <> -1 THEN CIRCLE (points(take, 15), points(take, 16)), 5000 / points(take, 10), _RGB32(255, 0, 0)
  92.     END IF
  93.     IF _MOUSEBUTTON(1) = 0 OR key_take = 0 THEN take_moving = -1
  94.     IF key_take AND _MOUSEBUTTON(1) AND take_moving = -1 AND take <> -1 THEN take_moving = take: dis_take = dis(points(take, 8), points(take, 9), points(take, 10))
  95.     IF take <> -1 AND inkey2$ = " " THEN points(take, 14) = points(take, 14) XOR 1
  96.  
  97.     'move with mouse 3d point
  98.     IF take_moving <> -1 THEN
  99.         IF dis(_MOUSEX - points(take_moving, 15), _MOUSEY - points(take_moving, 16), 0) > 20 THEN
  100.             mindis = 99999: FOR ang1 = 0 TO 360 * pip180 STEP 10 * pip180: FOR ang2 = 0 TO 180 * pip180 STEP 10 * pip180
  101.                     px = points(take_moving, 0) + adj * SIN(ang2) * COS(ang1): py = points(take_moving, 1) + adj * SIN(ang2) * SIN(ang1)
  102.                     pz = points(take_moving, 2) + adj * COS(ang2): vx = px: vy = py: vz = pz
  103.                     rotate_display vx, vy, vz: IF ABS(dis(vx, vy, vz) - dis_take) > 300 THEN _CONTINUE
  104.                     point3d_to_pset vx, vy, vz: disx = dis(_MOUSEX - vx, _MOUSEY - vy, 0)
  105.                     IF disx < mindis THEN mindis = disx: cx = px: cy = py: cz = pz
  106.             NEXT ang2, ang1: points(take_moving, 0) = cx: points(take_moving, 1) = cy: points(take_moving, 2) = cz
  107.         END IF
  108.     END IF
  109.  
  110.     mess$ = "moving WASD+mouse      Q-key fixing and moving point"
  111.     IF key_take THEN mess$ = "SPACE-fix/float point  MOUSEBUTTON - hold point"
  112.     IF take_moving <> -1 THEN mess$ = "moving the point with mouse"
  113.     COLOR _RGB32(200, 200, 200): LOCATE mony / 16 - 2: PRINT mess$;
  114.     _DISPLAY
  115.     CLS
  117.  
  118.  
  119. SUB point3d_to_pset (x, y, z): zlimit = -2: x = monx * .5 + (x + zlimit) / -z * 400: y = mony * .5 + (y + zlimit) / z * 400
  120.     z = 255 - ABS(z) * .1: IF z < 10 THEN grey = 10
  122.  
  123. SUB make_grid (acc, x0, y0, z0, x1, y1, z1, x2, y2, z2, x3, y3, z3, res): p_start = points_c
  124.     resA = INT((dis(x0 - x1, y0 - y1, z0 - z1) + dis(x2 - x3, y2 - y3, z2 - z3)) / 2 / res)
  125.     resB = INT((dis(x0 - x2, y0 - y2, z0 - z2) + dis(x1 - x3, y1 - y3, z1 - z3)) / 2 / res)
  126.     FOR a = 0 TO resA - 1: n = 1 / resA * a
  127.         tx1 = scale(x0, x1, n): ty1 = scale(y0, y1, n): tz1 = scale(z0, z1, n): tx2 = scale(x2, x3, n): ty2 = scale(y2, y3, n): tz2 = scale(z2, z3, n)
  128.         FOR b = 0 TO resB - 1: n = 1 / resB * b: l = p_start + b * resA + a
  129.             points(l, 0) = scale(tx1, tx2, n): points(l, 2) = scale(ty1, ty2, n): points(l, 1) = scale(tz1, tz2, n)
  130.             points(l, 6) = 1: points(l, 7) = acc: points(l, 14) = 1: IF (a = 0 OR a = resA - 1) AND (b = 0 OR b = resB - 1) THEN points(l, 14) = 0
  131.     NEXT b, a
  132.     FOR ty = 0 TO resB - 1: FOR tx = 0 TO resA - 2: p1 = tx + ty * resA: conn_add p1, p1 + 1, tightness: NEXT tx, ty
  133.     FOR tx = 0 TO resA - 1: FOR ty = 0 TO resB - 2: p1 = tx + ty * resA: conn_add p1, p1 + resA, tightness: NEXT ty, tx
  134.     points_c = points_c + resA * resB
  136.  
  137. FUNCTION dis (a, b, c): dis = SQR(a * a + b * b + c * c): END FUNCTION
  138. FUNCTION scale (a, b, n): scale = a + (b - a) * n: END FUNCTION
  139. SUB conn_add (p1, p2, tightness): IF p1 = p2 THEN EXIT SUB
  140.     FOR t = 0 TO conn_c_arr - 1: IF (conn(t, 0) = p1 AND conn(t, 1) = p2) OR (conn(t, 0) = p2 AND conn(t, 1) = p1) THEN EXIT SUB
  141.     NEXT t: REDIM disx(2): conn(conn_c, 0) = p1: conn(conn_c, 1) = p2
  142.     FOR t = 0 TO 2: conn(conn_c, 4 + t) = (ABS(points(p1, t) - points(p2, t))) * tightness: NEXT t
  143.     conn(conn_c, 2) = dis(conn(conn_c, 4), conn(conn_c, 5), conn(conn_c, 6)): conn(conn_c, 3) = 1: conn_c = conn_c + 1
  145. SUB rotate_2d (x, y, ang): x1 = x * COS(ang) - y * SIN(ang): y1 = x * SIN(ang) + y * COS(ang): x = x1: y = y1: END SUB
  146. SUB rotate_display (x, y, z)
  147.     x = (x - cam(0)) * 100: y = (y - cam(2)) * 100: z = -(z - cam(1)) * 100: rotate_2d x, z, cam(13): rotate_2d y, z, cam(14)
  149. FUNCTION degree (a, b): degreex = ATN(a / (b + .00001)) / pip180: degreex = degreex - 180 * ABS(0 > b): degreex = degreex - 360 * (degreex < 0): degree = degreex: END FUNCTION
  150.  
  151.  
  152.  
  153.  
  154.  
  155. 'POINTS array
  156. '0,1,2 : XYZ
  157. '3,4,5 : XYZ vector
  158. '6:active 7:identifier
  159. '8,9,10 :XYZ calculated to maptriangle
  160. '11,12,13 :XYZ moving temporary calculating
  161. '14: 0:FIX 1:MOVING
  162. '15,16,17 2D-X,Y, greyscale
  163.  
  164.  
  165. 'CONN array
  166. '0:point#1
  167. '1:point#2
  168. '2:distance
  169. '3:active
  170. '4:disX
  171. '5:disY
  172. '6:disZ
  173.  
  174.  
  175.  
  176.  
  177.  
  178.  

2
QB64 Discussion / Re: Variable v's Fixed Value
« on: April 04, 2022, 01:44:16 pm »
I actually use "CONST" if I want to see sub / function as well. In cases like "zoom" I leave. Because I might be thinking one and adding 1 line to the program that works with the "zoom" value. I don’t limit myself unnecessarily and I don’t put anything in it that I don’t have to.

3
Programs / Re: Garland of pearls
« on: April 03, 2022, 02:57:42 pm »
I'm very interested in this thing. Unfortunately I had little time, the season started. The above programs work in 2d. I started simulating the behavior of points in 3d. The goal would be collisions. I would like to see what I have described above. For now, you can see that it works nicely in 3D as well. Stretched mesh, gravity, WASD + mouse like in games. Use the space bar to move 1 point out of place. All other points will then respond and the system will attempt to return to the idle state.

Code: QB64: [Select]
  1. mouse_sens = .5
  2. tightness = 2
  3. gravity = .005
  4. lass = .94
  5. power = -0.3
  6.  
  7.  
  8.  
  10. CONST pip180 = 3.141592 / 180
  11. DIM SHARED monx, mony: monx = 800: mony = monx / _DESKTOPWIDTH * _DESKTOPHEIGHT: mon = _NEWIMAGE(monx, mony, 32): SCREEN mon: _FULLSCREEN: _MOUSEHIDE
  12.  
  13. DIM SHARED points_c_arr, conn_c_arr
  14. points_c_arr = 2000
  15. conn_c_arr = 2000
  16.  
  17. DIM SHARED points(points_c_arr, 19), conn(conn_c_arr, 9), cam(19), me(19), conn_c, points_c
  18.  
  19.  
  20. make_grid 1, -10, -10, 10, 10, -10, 0, -10, 10, 0, 10, 10, 0, 1.2
  21.  
  22.  
  23.  
  24.  
  25.  
  26. 'make_sphere 2, 0, 10, 0, 5, 8 'ident,X,Y,Z,RAD,POINTS
  28.  
  29.     'control
  30.     mousex = 0: mousey = 0: mw = 0: WHILE _MOUSEINPUT: mousex = mousex + _MOUSEMOVEMENTX: mousey = mousey + _MOUSEMOVEMENTY: mw = mw + _MOUSEWHEEL: WEND
  31.     k_left = _KEYDOWN(19200) OR _KEYDOWN(ASC("a")): k_right = _KEYDOWN(19712) OR _KEYDOWN(ASC("d"))
  32.     k_up = _KEYDOWN(18432) OR _KEYDOWN(ASC("w")): k_down = _KEYDOWN(20480) OR _KEYDOWN(ASC("s"))
  33.     me(2) = me(2) + mw: me(3) = me(3) - mousex * mouse_sens: me(4) = me(4) + mousey * mouse_sens: IF ABS(me(4)) > 80 THEN me(4) = 80 * SGN(me(4))
  34.     go = -.5 * ABS(k_left OR k_right OR k_up OR k_down): go_ang = -90 * (k_left - k_right) + (me(3)) - 90
  35.     ang1 = go_ang * pip180: ang2 = (90 - me(4)) * pip180: IF k_down THEN go = -go
  36.     me(0) = me(0) + SIN(ang2) * COS(ang1) * go: me(1) = me(1) + SIN(ang2) * SIN(ang1) * go: me(2) = me(2) + COS(ang2) * go
  37.     FOR t = 0 TO 9: cam(t) = me(t): NEXT t: cam(13) = cam(3) * pip180: cam(14) = cam(4) * pip180
  38.  
  39.  
  40.  
  41.     'calculate moving
  42.  
  43.     FOR a = 0 TO points_c_arr - 1: FOR t = 0 TO 2: points(a, 11 + t) = points(a, t): NEXT t: IF (points(a, 14) = 0 OR points(a, 6) = 0) THEN _CONTINUE
  44.  
  45.         'connections vector
  46.         vec(0) = 0: vec(1) = 0: vec(2) = 0: vec_c = 0
  47.         FOR t = 0 TO conn_c - 1: IF conn(t, 3) = 0 THEN _CONTINUE
  48.             x = -1: FOR t2 = 0 TO 1: IF conn(t, t2) = a THEN x = conn(t, t2 XOR 1)
  49.             NEXT t2
  50.             IF x <> -1 THEN
  51.                 FOR av = 0 TO 2: disv = points(a, av) - points(x, av): vec(av) = vec(av) + ((ABS(disv) - conn(t, 4 + av))) * SGN(disv) * power
  52.                 NEXT av: vec_c = vec_c + 1
  53.             END IF
  54.         NEXT t
  55.  
  56.         FOR av = 0 TO 2: points(a, 3 + av) = (points(a, 3 + av) + vec(av) / vec_c) * lass: NEXT av: points(a, 4) = points(a, 4) - gravity
  57.         FOR av = 0 TO 2: points(a, 11 + av) = points(a, av) + points(a, 3 + av): NEXT av
  58.     NEXT a
  59.  
  60.     FOR a = 0 TO points_c_arr - 1: FOR t = 0 TO 2: points(a, t) = points(a, 11 + t): NEXT t, a
  61.  
  62.     inkey2$ = INKEY$
  63.     IF inkey2$ = " " THEN
  64.         DO: a = INT(points_c_arr * RND(1)): LOOP UNTIL points(a, 14) AND points(a, 6)
  65.         rand = 30: FOR t = 0 TO 2: points(a, t) = points(a, t) + rand * RND(1) - rand / 2: NEXT t
  66.     END IF
  67.  
  68.     'calculating points
  69.     FOR ap = 0 TO points_c_arr - 1: IF points(ap, 6) = 0 THEN _CONTINUE
  70.         points(ap, 8) = points(ap, 0): points(ap, 9) = points(ap, 1): points(ap, 10) = points(ap, 2)
  71.     rotate_display points(ap, 8), points(ap, 9), points(ap, 10): NEXT ap
  72.  
  73.     'draw connections
  74.     FOR ac = 0 TO conn_c_arr - 1: IF conn(ac, 3) = 0 THEN _CONTINUE
  75.         FOR t1 = 0 TO 1: FOR t2 = 0 TO 2: coo(t1, t2) = points(conn(ac, t1), 8 + t2): NEXT t2, t1
  76.         line_3d coo(0, 0), coo(0, 1), coo(0, 2), coo(1, 0), coo(1, 1), coo(1, 2)
  77.     NEXT ac
  78.  
  79.     'fix points draw signal
  80.     FOR ap = 0 TO points_c_arr - 1: IF points(ap, 6) = 0 OR points(ap, 14) OR points(ap, 10) > -1 THEN _CONTINUE
  81.         FOR t2 = 0 TO 2: coo(0, t2) = points(ap, 8 + t2): NEXT t2
  82.         point3d_to_pset coo(0, 0), coo(0, 1), coo(0, 2): CIRCLE (coo(0, 0), coo(0, 1)), 12000 / coo(0, 2)
  83.     NEXT ap
  84.  
  85.     COLOR _RGB32(200, 200, 200): LOCATE mony / 16 - 2: PRINT "WASD+mouse      SPACE - impulse to a random point ";
  86.     _DISPLAY
  87.     CLS
  89.  
  90.  
  91. SUB line_3d (x0, y0, z0, x1, y1, z1)
  92.     IF z0 < -1 AND z1 < -1 THEN point3d_to_pset x0, y0, z0: point3d_to_pset x1, y1, z1: LINE (x0, y0)-(x1, y1)
  94.  
  95. SUB point3d_to_pset (x, y, z): multi = 400: zlimit = -2: x = monx / 2 + (x + zlimit) / -z * multi: y = mony / 2 + (y + zlimit) / z * multi
  96.     grey = 255 - ABS(z) * .1: IF grey < 10 THEN grey = 10
  97.     COLOR _RGB32(grey, grey, grey)
  99.  
  100. SUB make_grid (acc, x0, y0, z0, x1, y1, z1, x2, y2, z2, x3, y3, z3, res): p_start = points_c
  101.     resA = INT((dis(x0 - x1, y0 - y1, z0 - z1) + dis(x2 - x3, y2 - y3, z2 - z3)) / 2 / res)
  102.     resB = INT((dis(x0 - x2, y0 - y2, z0 - z2) + dis(x1 - x3, y1 - y3, z1 - z3)) / 2 / res)
  103.     FOR a = 0 TO resA - 1: n = 1 / resA * a
  104.         tx1 = scale(x0, x1, n): ty1 = scale(y0, y1, n): tz1 = scale(z0, z1, n): tx2 = scale(x2, x3, n): ty2 = scale(y2, y3, n): tz2 = scale(z2, z3, n)
  105.         FOR b = 0 TO resB - 1: n = 1 / resB * b: l = p_start + b * resA + a
  106.             points(l, 0) = scale(tx1, tx2, n): points(l, 2) = scale(ty1, ty2, n): points(l, 1) = scale(tz1, tz2, n)
  107.             points(l, 6) = 1: points(l, 7) = acc: points(l, 14) = 1: IF (a = 0 OR a = resA - 1) AND (b = 0 OR b = resB - 1) THEN points(l, 14) = 0
  108.     NEXT b, a
  109.     FOR ty = 0 TO resB - 1: FOR tx = 0 TO resA - 2: p1 = tx + ty * resA: conn_add p1, p1 + 1, tightness: NEXT tx, ty
  110.     FOR tx = 0 TO resA - 1: FOR ty = 0 TO resB - 2: p1 = tx + ty * resA: conn_add p1, p1 + resA, tightness: NEXT ty, tx
  111.     points_c = points_c + resA * resB
  113.  
  114. FUNCTION dis (a, b, c): dis = SQR(a * a + b * b + c * c): END FUNCTION
  115. FUNCTION scale (a, b, n): scale = a + (b - a) * n: END FUNCTION
  116. SUB conn_add (p1, p2, tightness): IF p1 = p2 THEN EXIT SUB
  117.     FOR t = 0 TO conn_c_arr - 1: IF (conn(t, 0) = p1 AND conn(t, 1) = p2) OR (conn(t, 0) = p2 AND conn(t, 1) = p1) THEN EXIT SUB
  118.     NEXT t: REDIM disx(2): conn(conn_c, 0) = p1: conn(conn_c, 1) = p2
  119.     FOR t = 0 TO 2: conn(conn_c, 4 + t) = (ABS(points(p1, t) - points(p2, t))) * tightness: NEXT t
  120.     conn(conn_c, 2) = dis(conn(conn_c, 4), conn(conn_c, 5), conn(conn_c, 6)): conn(conn_c, 3) = 1: conn_c = conn_c + 1
  122. SUB rotate_2d (x, y, ang): x1 = x * COS(ang) - y * SIN(ang): y1 = x * SIN(ang) + y * COS(ang): x = x1: y = y1: END SUB
  123. SUB rotate_display (x, y, z)
  124.     x = (x - cam(0)) * 100: y = (y - cam(2)) * 100: z = -(z - cam(1)) * 100: rotate_2d x, z, cam(13): rotate_2d y, z, cam(14)
  126. FUNCTION degree (a, b): degreex = ATN(a / (b + .00001)) / pip180: degreex = degreex - 180 * ABS(0 > b): degreex = degreex - 360 * (degreex < 0): degree = degreex: END FUNCTION
  127.  
  128.  
  129.  
  130.  
  131.  
  132. 'POINTS array
  133. '0,1,2 : X,Y,Z
  134. '3,4,5 : X,Y,Z vector
  135. '6:active 7:identifier
  136. '8,9,10 :X,Y,Z calculated to maptriangle
  137. '11,12,13 :X,Y,Z moving temporary calculating
  138. '14: 0:FIX 1:MOVING
  139.  
  140.  
  141. 'CONN array
  142. '0:point#1
  143. '1:point#2
  144. '2:distance
  145. '3:active
  146. '4:disX
  147. '5:disY
  148. '6:disZ
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  

4
QB64 Discussion / Re: Variable v's Fixed Value
« on: March 31, 2022, 09:58:55 am »
I'm sorry, I just saw the question, I was very busy for a few days.
In the program, the ZOOM variable only means how many times the coordinate points are multiplied to show well on the screen. The higher the magnification.
POWER is a repulsive or attractive force created by the connection of two points.
Fixed or fixed values? I prefer something to be associated with something. For example, I like to clarify the screen size with MONX, MONY at the beginning of the program, because when I place an object on the screen, I always specify it in proportion. If I want it to be roughly at the bottom left of the screen, but not in the corner, I specify that the image be MONX x .1 and MONY x .8 (example only) and MONX / 10 x MONX / 10- and size. So anytime, any time I resize the window, the position of the image will follow proportionally. That was just one example. I don't like fixed values.

5
Programs / Re: Garland of pearls
« on: March 23, 2022, 08:00:58 pm »
There is much more to it than interesting!
I asked the question on a Hungarian physics forum on facebook. If there is a meaningful answer, I will inform you! It needs to be addressed!

Hello !
It will be a little longer now. I want an answer to a question I haven't been able to sleep in days :)
Anyone who follows the group may know that I am involved in programming. It used to be machine control, more recently the math behind games. (Interestingly, behind the game programming are the most outdated applications of mathematics)
My recent creations:


I have only considered it important to note this now because I have been dealing with them feverishly lately. The basic development environment, in which I program, can only draw 1 triangle on the screen. That's why I wrote the graphic "engine", the camera movement, the physical background of the movements myself. This developer is so unpopular that I don't think QB64 is known in Hungary. I like to be active in a foreign forum because there are about a handful of people (a handful on a global scale) who struggle with this outdated developer and appreciate it.
Here a forum member wrote a simulation of a pendulum motion. This should be imagined as a fixed point + rigid connector + some weight + rigid connector + some weight .... What happens if we push? He returned the thing quite realistically on the screen. It’s some crazy math, a lot of computation, and limited because it works with 2 “weights,” and with 3, but 4 or more has so much computation that it crumbles all over and can’t be shown in real time.
I had an idea then. What if we created points in space / plane without giving it fixed coordinates. All you need to know about that point is that it is located relative to another point. There should be nothing else for that point other than to strive to maintain its relation to the other point specified.
Do this: "current situation" = "current situation" + ("where you should be" - "current situation") x "all"
This is important here. If "all" is 0, no change is made. If 1, an immediate change occurs. If all are between 0 and 1, it gives a proportionality to the dynamics / flexibility / speed of change.
Now I want to express the point somehow. Imagine a "space." It contains objects. Objects are defined by point clouds. In point clouds, we make connections between points. To an adjacent, or object center, or just a smooth random point, we connect the 1 given point within the point cloud that we just want to define, "smart." We determine for him how far away he should be from that other point ... We will try to make as many such "connections" as possible at all points in that point and that object so that he is as smart as possible and will know "where he is. ".
If this is the case, start the "simulation"! Suppose thousands of points are defined as above. Each point begins to move at such a speed / direction that the connections assigned to it can be fulfilled. All right, we're waiting ... after a while, all movement will stop. "system balanced". The shapes are nicely outlined in front of us.
Now move only 1 out of position !!! that is, ONE point !!!!!!!! This will affect the points you are connected to ... and also the points you are connected to ... and those .... so everyone is running that 1 pulse. A mechanism is starting. Impulses give and take ...
The video shows how I programmed this principle differently. Chain, an image, and a cobweb. In each case, I gave relationships to the points as above. The thing works. And here comes the point. Smart people have already written the mathematical model for every event in Newtonian physics. Rigid-body, wave curves, collision of bodies, pendulum motion, centrifugal force .... what else to say ???????????????
This idea, which I did, completely bypasses the knowledge and use of these. The screen is fully functional. You don't need any physical knowledge or complications. We need a point and the strengths of its connections and connections (or how long we tie it to another point) (perhaps a metaphor for mass?)
Using the principle, we can in principle simulate all Newtonian physics by completely circumventing their use. Example..something .... if we create a brick wall based on the above principle, we can shoot a brick out of it so that the wall collapses in reality..etc ... that's just 1 of the many options.
Here's my question. Unfortunately, I did not study mathematics or physics at higher levels. I'm sure this thing isn't some kind of invention. This can only be new to me. This may be a known, existing theory that is also taught in good places. I would like to ask those who were so lucky to learn things like what it is? what is it used for where can i find writing about this
Thank you!

(just a small philosophical remark: if we hit the wall! the above principle would show the tiny vibration of every part of the wall. what if the above principle simulates the connection between molecules? since reality "renders" so fast that we can't even measure it with instruments, we're just amazed at the end result, because it happens so fast that we can't see, so we don't understand)
To download programs running under video (running under windows) with source code:
https://drive.google.com/file/d/1hkIqbFx_JGzhYpGdQjZ1h0lZYFtiw3Ce/view?usp=sharing


6
Programs / Re: Garland of pearls
« on: March 22, 2022, 07:09:43 pm »
enter a picture in the program "temp = _LOADIMAGE("pic.jpg", 32)"
crush the picture with your mouse!

Code: QB64: [Select]
  2. CONST pip180 = 3.141592 / 180
  3. monx = 800: mony = monx / _DESKTOPWIDTH * _DESKTOPHEIGHT: mon = _NEWIMAGE(monx, mony, 32): SCREEN mon: _DISPLAYORDER _HARDWARE , _SOFTWARE
  5.  
  6. min_weight = 6
  7. gravity = 1
  8. obj_c = 1000 'number objects
  9. DIM obj(obj_c - 1, 9) '0,1 posX,Y 3fix/mover
  10. cent_x = monx / 2: cent_y = mony / 2: zoom = 2.4
  11.  
  12.  
  13.  
  14. conn_c = 1000: DIM conn(conn_c - 1, 4)
  15.  
  16.  
  17. temp = _LOADIMAGE("pic.jpg", 32): _DEST temp: _SETALPHA 100: pic = _COPYIMAGE(temp, 33): _FREEIMAGE temp
  18.  
  19. sizex = 300: sizey = INT(sizex / _WIDTH(pic) * _HEIGHT(pic))
  20. felbx = 15: felby = INT(felbx / _WIDTH(pic) * _HEIGHT(pic))
  21. unitx = INT(_WIDTH(pic) / (felbx - 1))
  22. unity = INT(_HEIGHT(pic) / (felby - 1))
  23.  
  24.  
  25. FOR ty = 0 TO felby - 1
  26.     FOR tx = 0 TO felbx - 1
  27.         aobj = tx + ty * felbx
  28.         obj(aobj, 0) = sizex / felbx * tx - sizex / 2
  29.         obj(aobj, 1) = sizey / felby * ty - sizey / 2
  30.         obj(aobj, 2) = 5
  31.         obj(aobj, 3) = 1
  32.         IF (ty = 0 OR ty = felby - 1) AND (tx = 0 OR tx = felbx - 1) THEN obj(aobj, 3) = 0
  33.         obj(aobj, 4) = 1
  34.     NEXT tx
  35. NEXT ty
  36.  
  37.  
  38. 'connections
  39. FOR ty = 0 TO felby - 1: FOR tx = 0 TO felbx - 2: conn(ac, 0) = tx + ty * felbx: conn(ac, 1) = conn(ac, 0) + 1
  40.         conn(ac, 2) = SQR((obj(conn(ac, 0), 0) - obj(conn(ac, 1), 0)) ^ 2 + (obj(conn(ac, 0), 1) - obj(conn(ac, 1), 1)) ^ 2) * .8
  41. conn(ac, 3) = 1: ac = ac + 1: NEXT tx, ty
  42.  
  43. FOR tx = 0 TO felbx - 1: FOR ty = 0 TO felby - 2: conn(ac, 0) = tx + ty * felbx: conn(ac, 1) = conn(ac, 0) + felbx
  44.         conn(ac, 2) = SQR((obj(conn(ac, 0), 0) - obj(conn(ac, 1), 0)) ^ 2 + (obj(conn(ac, 0), 1) - obj(conn(ac, 1), 1)) ^ 2) * .8
  45. conn(ac, 3) = 1: ac = ac + 1: NEXT ty, tx
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56. stiff = 1
  57. COLOR , 0
  58. DO: _LIMIT 30 '100
  59.  
  60.     'show info
  61.     wch = -1: FOR a = 0 TO obj_c - 1
  62.         IF _MOUSEX < 200 AND INT(_MOUSEY / 16) = a THEN wch = a: COLOR _RGB32(255, 255, 255) ELSE COLOR _RGB32(100, 100, 100)
  63.         ' LOCATE a + 1, 1: PRINT a; "# weight:"; INT(obj(a, 2) * 10) / 10;: IF obj(a, 3) = 0 THEN PRINT "FIX";
  64.         'IF wch = a THEN PRINT "<--- using mousewheel !" 'change weights with mousewheel
  65.         IF wch <> -1 THEN obj(wch, 2) = obj(wch, 2) + mw / 10: IF obj(wch, 2) < min_weight THEN obj(wch, 2) = min_weight
  66.     NEXT a
  67.     LOCATE mony / 16 - 2, 1
  68.     IF moving THEN PRINT "press SPACE to FIX/FLY" ELSE PRINT "grab the ball with the mouse and move it!"
  69.     '    PRINT "S-button - rigid / flexible connections";
  70.     'draw objects
  71.     FOR a = 0 TO obj_c - 1: IF obj(a, 4) = 0 THEN _CONTINUE
  72.         x = cent_x + obj(a, 0) * zoom: y = cent_y + obj(a, 1) * zoom
  73.         obj(a, 8) = x
  74.         obj(a, 9) = y
  75.         IF mouse_near = a THEN COLOR _RGB32(255, 0, 0) ELSE COLOR _RGB32(200, 200, 200)
  76.         '        CIRCLE (x, y), obj(a, 2) * zoom: _PRINTSTRING (x, y), LTRIM$(STR$(a))
  77.  
  78.     NEXT a
  79.  
  80.     'draw connections
  81.     COLOR _RGB32(200, 200, 200)
  82.     FOR a = 0 TO conn_c - 1: IF conn(a, 3) = 0 THEN _CONTINUE
  83.         IF obj(conn(a, 0), 4) = 0 OR obj(conn(a, 1), 4) = 0 THEN _CONTINUE
  84.         x1 = obj(conn(a, 0), 0) * zoom + cent_x: y1 = obj(conn(a, 0), 1) * zoom + cent_y
  85.         x2 = obj(conn(a, 1), 0) * zoom + cent_x: y2 = obj(conn(a, 1), 1) * zoom + cent_y
  86.         '    LINE (x1, y1)-(x2, y2), _RGBA32(0, 0, 0, 128)
  87.     NEXT a
  88.  
  89.  
  90.  
  91.     'calculate moving
  92.     FOR a = 0 TO obj_c - 1: obj(a, 6) = obj(a, 0): obj(a, 7) = obj(a, 1)
  93.         IF obj(a, 3) = 0 OR (mouse_near = a AND _MOUSEBUTTON(1)) OR obj(a, 4) = 0 THEN _CONTINUE
  94.         'connections vector
  95.         vec_x = 0: vec_y = 0: vec_c = 0
  96.         FOR t = 0 TO conn_c - 1: IF conn(t, 3) = 0 THEN _CONTINUE
  97.             x = -1: FOR t2 = 0 TO 1: IF conn(t, t2) = a THEN x = conn(t, t2 XOR 1)
  98.             NEXT t2
  99.             IF x <> -1 THEN
  100.                 IF obj(x, 4) THEN
  101.                     disx = obj(a, 0) - obj(x, 0)
  102.                     disy = obj(a, 1) - obj(x, 1)
  103.                     ang = (-degree(disx, disy) + 0) * pip180
  104.                     dis = SQR(disx * disx + disy * disy)
  105.                     power = (dis - conn(t, 2)) * 2
  106.                     IF stiff THEN power = power / obj(a, 2) ^ 2 * 4
  107.                     vec_x = vec_x + SIN(ang) * power
  108.                     vec_y = vec_y - COS(ang) * power
  109.                     vec_c = vec_c + 1
  110.                 END IF
  111.             END IF
  112.         NEXT t
  113.  
  114.         gravity = 0 '.1
  115.         lass = .93
  116.         obj(a, 4) = (obj(a, 4) + vec_x / vec_c) * lass
  117.         obj(a, 5) = (obj(a, 5) + vec_y / vec_c) * lass + gravity
  118.         obj(a, 6) = obj(a, 0) + obj(a, 4)
  119.         obj(a, 7) = obj(a, 1) + obj(a, 5)
  120.     NEXT a
  121.     FOR a = 0 TO obj_c - 1: FOR t = 0 TO 1: obj(a, t) = obj(a, 6 + t): NEXT t, a
  122.  
  123.  
  124.  
  125.  
  126.     'mouse moving activing
  127.     mw = 0: WHILE _MOUSEINPUT: mw = mw + _MOUSEWHEEL: WEND
  128.     mindis = 99999999
  129.     IF moving = 0 THEN
  130.         mouse_near = -1
  131.         FOR a = 0 TO obj_c - 1
  132.             disx = (cent_x + obj(a, 0) * zoom) - _MOUSEX
  133.             disy = (cent_y + obj(a, 1) * zoom) - _MOUSEY
  134.             dis = SQR(disx * disx + disy * disy)
  135.             IF dis < 60 AND mindis > dis THEN mouse_near = a: mindis = dis
  136.         NEXT a
  137.     END IF
  138.  
  139.     IF (_MOUSEBUTTON(1) AND mouse_near <> -1) OR moving THEN
  140.         moving = 1
  141.         obj(mouse_near, 0) = (_MOUSEX - cent_x) / zoom
  142.         obj(mouse_near, 1) = (_MOUSEY - cent_y) / zoom
  143.     END IF
  144.  
  145.     moving = moving AND _MOUSEBUTTON(1)
  146.  
  147.     stiff = stiff XOR ABS(inkey2$ = "s")
  148.  
  149.     'picture draw
  150.     FOR ty = 0 TO felby - 2: sy1 = INT(unity * ty): sy2 = sy1 + unity - 1
  151.         FOR tx = 0 TO felbx - 2: sx1 = INT(unitx * tx): sx2 = sx1 + unitx - 1: aobj = tx + ty * felbx
  152.             FOR t = 0 TO 7: p(t) = obj(aobj + (t AND 1) + SGN(t AND 2) * felbx, 8 + SGN(t AND 4)): NEXT t
  153.             _MAPTRIANGLE (sx1, sy1)-(sx2, sy1)-(sx1, sy2), pic TO(p(0), p(4))-(p(1), p(5))-(p(2), p(6))
  154.             _MAPTRIANGLE (sx2, sy2)-(sx2, sy1)-(sx1, sy2), pic TO(p(3), p(7))-(p(1), p(5))-(p(2), p(6))
  155.         NEXT tx
  156.     NEXT ty
  157.  
  158.     _DISPLAY: CLS
  159.  
  160.     'commands
  161.     inkey2$ = INKEY$: IF inkey2$ = CHR$(27) THEN END
  162.     IF inkey2$ = " " AND _MOUSEBUTTON(1) AND moving THEN obj(mouse_near, 3) = obj(mouse_near, 3) XOR 1
  163.  
  164.  
  166.  
  167.  
  168. FUNCTION degree (a, b): degreex = ATN(a / (b + .00001)) / pip180: degreex = degreex - 180 * ABS(0 > b): degreex = degreex - 360 * (degreex < 0): degree = degreex: END FUNCTION
  169. FUNCTION randx (x): randx = x * RND(1) - x / 2: END FUNCTION
  170.  

7
Programs / Re: Garland of pearls
« on: March 22, 2022, 03:36:20 pm »
grab the cobweb!

Code: QB64: [Select]
  2. CONST pip180 = 3.141592 / 180
  3. monx = 800: mony = 800: mon = _NEWIMAGE(monx, mony, 32): SCREEN mon
  4.  
  5. min_weight = 6
  6. gravity = 1
  7. obj_c = 1000 'number objects
  8. DIM obj(obj_c - 1, 9) '0,1 posX,Y 3fix/mover
  9. cent_x = monx / 2: cent_y = mony / 2: zoom = 1.8
  10.  
  11.  
  12.  
  13. 'cobweb
  14. ag0 = 320 'size
  15. ag1 = 18 'felbontas ag
  16. ag2 = 14 'felbontas net
  17. ag3 = ag0 * .18 / 2 'rad1
  18. ag4 = ag0 * .9 / 2 'rad2
  19.  
  20.  
  21. FOR t = 0 TO ag1 - 1
  22.     ang = (360 / ag1 * t + randx(360 / ag1 * .3)) * pip180
  23.     dis = ABS(ag3 - ag4)
  24.     x1 = SIN(ang) * ag3 + randx(ag3 * .1)
  25.     y1 = COS(ang) * ag3 + randx(ag3 * .1)
  26.     x2 = SIN(ang) * ag4 + randx(ag4 * .2)
  27.     y2 = COS(ang) * ag4 + randx(ag4 * .2)
  28.  
  29.     FOR t2 = 0 TO ag2 - 1
  30.         aobj = t2 + ag2 * t
  31.         obj(aobj, 0) = x1 + (x2 - x1) / (ag2 - 1) * t2
  32.         obj(aobj, 1) = y1 + (y2 - y1) / (ag2 - 1) * t2
  33.         obj(aobj, 2) = 1.9 'min_weight
  34.         obj(aobj, 3) = 1
  35.         '        IF t2 = ag2 - 1 THEN obj(aobj, 3) = 0
  36.  
  37.         obj(aobj, 4) = 1
  38.     NEXT t2
  40.  
  41.  
  42. 'connections
  43. conn_c = 1000: DIM conn(conn_c - 1, 4): ac = 0
  44.  
  45. FOR t = 0 TO ag2 - 2: FOR t2 = 0 TO ag1 - 1: conn(ac, 0) = t2 * ag2 + t: conn(ac, 1) = t2 * ag2 + ((t + 1) MOD ag2)
  46.         conn(ac, 2) = SQR((obj(conn(ac, 0), 0) - obj(conn(ac, 1), 0)) ^ 2 + (obj(conn(ac, 0), 1) - obj(conn(ac, 1), 1)) ^ 2)
  47.         conn(ac, 2) = conn(ac, 2) + randx(conn(ac, 2) * .1)
  48. conn(ac, 3) = 1: ac = ac + 1: NEXT t2, t
  49.  
  50. FOR t = 0 TO ag2 - 1: FOR t2 = 0 TO ag1 - 1: conn(ac, 0) = t + ag2 * t2: conn(ac, 1) = t + ag2 * ((t2 + 1) MOD ag1)
  51.         conn(ac, 2) = SQR((obj(conn(ac, 0), 0) - obj(conn(ac, 1), 0)) ^ 2 + (obj(conn(ac, 0), 1) - obj(conn(ac, 1), 1)) ^ 2)
  52.         conn(ac, 2) = conn(ac, 2) + randx(conn(ac, 2) * .1)
  53. conn(ac, 3) = 1: ac = ac + 1: NEXT t2, t
  54.  
  55. polar_p = ag1 * ag2
  56. FOR t = 0 TO 3
  57.     obj(polar_p + t, 0) = ((t AND 1) * 2 - 1) * ag0 / 2
  58.     obj(polar_p + t, 1) = (SGN(t AND 2) * 2 - 1) * ag0 / 2
  59.     obj(polar_p + t, 2) = 6 'min_weight
  60.     obj(polar_p + t, 3) = 0
  61.     obj(polar_p + t, 4) = 1
  62.     mindis = 999999999
  63.  
  64.     FOR a = 0 TO polar_p - 1
  65.         dis = SQR((obj(polar_p + t, 0) - obj(a, 0)) ^ 2 + (obj(polar_p + t, 1) - obj(a, 1)) ^ 2) * .8
  66.         IF dis < mindis THEN mindis = dis: conn_n = a
  67.     NEXT a
  68.     conn(ac, 0) = polar_p + t
  69.     conn(ac, 1) = conn_n
  70.     conn(ac, 2) = mindis
  71.     conn(ac, 3) = 1
  72.     ac = ac + 1
  73.  
  75.  
  76.  
  77.  
  78.  
  79.  
  80. stiff = 1
  81.  
  82. DO: _LIMIT 100
  83.  
  84.     'show info
  85.     wch = -1: FOR a = 0 TO obj_c - 1
  86.         IF _MOUSEX < 200 AND INT(_MOUSEY / 16) = a THEN wch = a: COLOR _RGB32(255, 255, 255) ELSE COLOR _RGB32(100, 100, 100)
  87.         ' LOCATE a + 1, 1: PRINT a; "# weight:"; INT(obj(a, 2) * 10) / 10;: IF obj(a, 3) = 0 THEN PRINT "FIX";
  88.         'IF wch = a THEN PRINT "<--- using mousewheel !" 'change weights with mousewheel
  89.         IF wch <> -1 THEN obj(wch, 2) = obj(wch, 2) + mw / 10: IF obj(wch, 2) < min_weight THEN obj(wch, 2) = min_weight
  90.     NEXT a
  91.     LOCATE mony / 16 - 2, 1
  92.     IF moving THEN PRINT "press SPACE to FIX/FLY" ELSE PRINT "grab the ball with the mouse and move it!"
  93.     '    PRINT "S-button - rigid / flexible connections";
  94.     'draw objects
  95.     FOR a = 0 TO obj_c - 1: IF obj(a, 4) = 0 THEN _CONTINUE
  96.         x = cent_x + obj(a, 0) * zoom: y = cent_y + obj(a, 1) * zoom
  97.         IF mouse_near = a THEN COLOR _RGB32(255, 0, 0) ELSE COLOR _RGB32(200, 200, 200)
  98.         '        CIRCLE (x, y), obj(a, 2) * zoom: _PRINTSTRING (x, y), LTRIM$(STR$(a))
  99.     NEXT a
  100.  
  101.     'draw connections
  102.  
  103.     COLOR _RGB32(200, 200, 200)
  104.     FOR a = 0 TO conn_c - 1: IF conn(a, 3) = 0 THEN _CONTINUE
  105.         IF obj(conn(a, 0), 4) = 0 OR obj(conn(a, 1), 4) = 0 THEN _CONTINUE
  106.         x1 = obj(conn(a, 0), 0) * zoom + cent_x: y1 = obj(conn(a, 0), 1) * zoom + cent_y
  107.         x2 = obj(conn(a, 1), 0) * zoom + cent_x: y2 = obj(conn(a, 1), 1) * zoom + cent_y
  108.     LINE (x1, y1)-(x2, y2): NEXT a
  109.  
  110.  
  111.  
  112.     'calculate moving
  113.     FOR a = 0 TO obj_c - 1: obj(a, 6) = obj(a, 0): obj(a, 7) = obj(a, 1)
  114.         IF obj(a, 3) = 0 OR (mouse_near = a AND _MOUSEBUTTON(1)) OR obj(a, 4) = 0 THEN _CONTINUE
  115.         'connections vector
  116.         vec_x = 0: vec_y = 0: vec_c = 0
  117.         FOR t = 0 TO conn_c - 1: IF conn(t, 3) = 0 THEN _CONTINUE
  118.             x = -1: FOR t2 = 0 TO 1: IF conn(t, t2) = a THEN x = conn(t, t2 XOR 1)
  119.             NEXT t2
  120.             IF x <> -1 THEN
  121.                 IF obj(x, 4) THEN
  122.                     disx = obj(a, 0) - obj(x, 0)
  123.                     disy = obj(a, 1) - obj(x, 1)
  124.                     ang = (-degree(disx, disy) + 0) * pip180
  125.                     dis = SQR(disx * disx + disy * disy)
  126.                     power = (dis - conn(t, 2)) * .6 ' obj(a, 2) ^ 2 '/ 20
  127.                     IF stiff THEN power = power / obj(a, 2) ^ 2 * 4
  128.                     vec_x = vec_x + SIN(ang) * power
  129.                     vec_y = vec_y - COS(ang) * power
  130.                     vec_c = vec_c + 1
  131.                 END IF
  132.             END IF
  133.         NEXT t
  134.  
  135.         gravity = 0 '.1
  136.         lass = .90
  137.         obj(a, 4) = (obj(a, 4) + vec_x / vec_c) * lass
  138.         obj(a, 5) = (obj(a, 5) + vec_y / vec_c) * lass + gravity
  139.         obj(a, 6) = obj(a, 0) + obj(a, 4)
  140.         obj(a, 7) = obj(a, 1) + obj(a, 5)
  141.     NEXT a
  142.     FOR a = 0 TO obj_c - 1: FOR t = 0 TO 1: obj(a, t) = obj(a, 6 + t): NEXT t, a
  143.  
  144.  
  145.  
  146.  
  147.     'mouse moving activing
  148.     mw = 0: WHILE _MOUSEINPUT: mw = mw + _MOUSEWHEEL: WEND
  149.     mindis = 99999999
  150.     IF moving = 0 THEN
  151.         mouse_near = -1
  152.         FOR a = 0 TO obj_c - 1
  153.             disx = (cent_x + obj(a, 0) * zoom) - _MOUSEX
  154.             disy = (cent_y + obj(a, 1) * zoom) - _MOUSEY
  155.             dis = SQR(disx * disx + disy * disy)
  156.             IF dis < 60 AND mindis > dis THEN mouse_near = a:mindis = dis
  157.         NEXT a
  158.     END IF
  159.  
  160.     IF (_MOUSEBUTTON(1) AND mouse_near <> -1) OR moving THEN
  161.         moving = 1
  162.         obj(mouse_near, 0) = (_MOUSEX - cent_x) / zoom
  163.         obj(mouse_near, 1) = (_MOUSEY - cent_y) / zoom
  164.     END IF
  165.  
  166.     moving = moving AND _MOUSEBUTTON(1)
  167.  
  168.     stiff = stiff XOR ABS(inkey2$ = "s")
  169.  
  170.  
  171.     _DISPLAY: CLS
  172.  
  173.     'commands
  174.     inkey2$ = INKEY$: IF inkey2$ = CHR$(27) THEN END
  175.     IF inkey2$ = " " AND _MOUSEBUTTON(1) AND moving THEN obj(mouse_near, 3) = obj(mouse_near, 3) XOR 1
  176.  
  177.  
  179.  
  180.  
  181.  
  182.  
  183. FUNCTION degree (a, b): degreex = ATN(a / (b + .00001)) / pip180: degreex = degreex - 180 * ABS(0 > b): degreex = degreex - 360 * (degreex < 0): degree = degreex: END FUNCTION
  184.  
  185. FUNCTION randx (x): randx = x * RND(1) - x / 2: END FUNCTION

8
Programs / Re: Garland of pearls
« on: March 22, 2022, 10:33:05 am »
The trouble is, because of the translation, I can’t express myself the way I want to.
It doesn't matter what we measure. Distance, volume, size, speed, rotation, anything.
Quantity "A". This quantity wants to be "B".
Then use: B = B + (A-B) x "all"
This offers many opportunities.
I hope this thought, which I really want to go through, has gone through somewhat.
I wrote a simple illustration on it.

Code: QB64: [Select]
  1. SCREEN _NEWIMAGE(800, 600, 32)
  2.  
  3.  
  4. CHANGING = .05
  5. 'if CHANGING 0:no move
  6. 'if CHANGING 1:immediate displacement
  7. 'if CHANGING between 0,1 elastic displacement
  8.  
  9.  
  11.     _LIMIT 50
  13.     x = x + (_MOUSEX - x) * CHANGING
  14.     y = y + (_MOUSEY - y) * CHANGING
  15.     _PRINTSTRING (x, y), "X"
  16.     _DISPLAY
  17.     CLS
  19.  

It shows movement. But the principle can be distance, volume, size, speed, rotation, anything ...

9
Programs / Re: Garland of pearls
« on: March 22, 2022, 10:11:44 am »
Thanks for the comments. I honestly admit, I myself am amazed that this works. In games, I long ago figured it out to display text, animate, move when I want change to be flexible. Let A and B. Let "A" be the current value. Let B be the new value. By definition B = A. But this makes it inflexible in timing, animation, positioning work.
Therefore, I usually do it so that not B = A, but B = B + (A-B) x "all", for example. "All" must be between 0 and 1! This will be a constant value. If 0, no offset occurs. If 1, an inelastic, immediate displacement occurs. On the other hand, if I write e.g.
 The more times we perform it, the closer it gets to B = A. The larger the difference between A and B, the faster it makes the change, and when it is almost reached, it slows down. In fact, it never reaches.
Implementing STxAxTIC with the pendulum program returns the cleanest physics. Then came my idea of ​​what it would be like to forget that much complicated physics and use the method described above to move the ball. Ball "B" wants to be in place of "A" ... B = B + (A-B) x "all"
In principle, since this seems to work, this method could also be used to move, collide without the use of rigid-body physics or any more serious mathematics. In theory. Pendulum, brick wall collapses by bricks, car overturns, precise collision of two bodies with rotations and directions. It seems as if this simple method bypasses the part of physics where pulse retention, angular velocity, vectors should be taken into account. The system simply extracts from itself.
I’m very curious as to how far this can go, and is the latter suggestion true? does it work
I described these so that if anyone cares and wants to deal with it, sees opportunities in it like I do, then see how simple it is. try it !
I set up the program to eliminate flexibility. From here, the things described above should work in principle.
The program starts and you can turn flexibility on and off with the "S" key. If we turn it off, we get approximately the pendulum motion of STxAxTIC. When turned on, it will be flexible as if they were connected with rubber.

Code: QB64: [Select]
  2. CONST pip180 = 3.141592 / 180
  3. monx = 800: mony = 800: mon = _NEWIMAGE(monx, mony, 32): SCREEN mon
  4.  
  5. min_weight = 6
  6. gravity = 1
  7. obj_c = 10 'number objects
  8. DIM obj(obj_c - 1, 9) '0,1 posX,Y 3fix/mover
  9. cent_x = monx / 2: cent_y = mony * .2: zoom = 2
  10.  
  11. 'garland
  12. FOR a = 0 TO obj_c - 1
  13.     obj(a, 0) = (RND(1) - .5) * 200
  14.     obj(a, 1) = (RND(1) - .5) * 200
  15.     obj(a, 2) = min_weight + min_weight / 2 * RND(1)
  16.     obj(a, 3) = 1 '<---- 0 if fix, 1 if fly
  17.     obj(a, 4) = 1 'active point
  18.  
  20. obj(0, 0) = -64.5: obj(0, 1) = 5: obj(0, 3) = 0
  21. 'obj(obj_c - 1, 0) = 183.5: obj(obj_c - 1, 1) = -.5: obj(obj_c - 1, 3) = 0
  22.  
  23. 'connections
  24. conn_c = 100: DIM conn(conn_c - 1, 4)
  25. FOR a = 0 TO obj_c - 2
  26.     conn(a, 0) = a
  27.     conn(a, 1) = a + 1
  28.     conn(a, 2) = 15 'rope width
  29.     conn(a, 3) = 1 'active connect
  31.  
  32.  
  33.  
  34. stiff = 1
  35.  
  36. DO: _LIMIT 100
  37.  
  38.     'show info
  39.     wch = -1: FOR a = 0 TO obj_c - 1
  40.         IF _MOUSEX < 200 AND INT(_MOUSEY / 16) = a THEN wch = a: COLOR _RGB32(255, 255, 255) ELSE COLOR _RGB32(100, 100, 100)
  41.         LOCATE a + 1, 1: PRINT a; "# weight:"; INT(obj(a, 2) * 10) / 10;: IF obj(a, 3) = 0 THEN PRINT "FIX";
  42.         IF wch = a THEN PRINT "<--- using mousewheel !" 'change weights with mousewheel
  43.         IF wch <> -1 THEN obj(wch, 2) = obj(wch, 2) + mw / 10: IF obj(wch, 2) < min_weight THEN obj(wch, 2) = min_weight
  44.     NEXT a
  45.     LOCATE mony / 16 - 2, 1
  46.     IF moving THEN PRINT "press SPACE to FIX/FLY" ELSE PRINT "grab the ball with the mouse and move it!"
  47.     PRINT "S-button - rigid / flexible connections";
  48.     'draw objects
  49.     FOR a = 0 TO obj_c - 1: IF obj(a, 4) = 0 THEN _CONTINUE
  50.         x = cent_x + obj(a, 0) * zoom: y = cent_y + obj(a, 1) * zoom
  51.         IF mouse_near = a THEN COLOR _RGB32(255, 0, 0) ELSE COLOR _RGB32(200, 200, 200)
  52.         CIRCLE (x, y), obj(a, 2) * zoom: _PRINTSTRING (x, y), LTRIM$(STR$(a))
  53.     NEXT a
  54.  
  55.     'draw connections
  56.  
  57.     COLOR _RGB32(150, 150, 150)
  58.     FOR a = 0 TO conn_c - 1: IF conn(a, 3) = 0 THEN _CONTINUE
  59.         IF obj(conn(a, 0), 4) = 0 OR obj(conn(a, 1), 4) = 0 THEN _CONTINUE
  60.         x1 = obj(conn(a, 0), 0) * zoom + cent_x: y1 = obj(conn(a, 0), 1) * zoom + cent_y
  61.         x2 = obj(conn(a, 1), 0) * zoom + cent_x: y2 = obj(conn(a, 1), 1) * zoom + cent_y
  62.     LINE (x1, y1)-(x2, y2): NEXT a
  63.  
  64.  
  65.  
  66.     'calculate moving
  67.     FOR a = 0 TO obj_c - 1: obj(a, 6) = obj(a, 0): obj(a, 7) = obj(a, 1)
  68.         IF obj(a, 3) = 0 OR (mouse_near = a AND _MOUSEBUTTON(1)) OR obj(a, 4) = 0 THEN _CONTINUE
  69.         'connections vector
  70.         vec_x = 0: vec_y = 0: vec_c = 0
  71.         FOR t = 0 TO conn_c - 1: IF conn(t, 3) = 0 THEN _CONTINUE
  72.             x = -1: FOR t2 = 0 TO 1: IF conn(t, t2) = a THEN x = conn(t, t2 XOR 1)
  73.             NEXT t2
  74.             IF x <> -1 THEN
  75.                 IF obj(x, 4) THEN
  76.                     disx = obj(a, 0) - obj(x, 0)
  77.                     disy = obj(a, 1) - obj(x, 1)
  78.                     ang = (-degree(disx, disy) + 0) * pip180
  79.                     dis = SQR(disx * disx + disy * disy)
  80.                     power = (dis - (obj(a, 2) + obj(x, 2) + conn(t, 2))) * 1 ' obj(a, 2) ^ 2 '/ 20
  81.                     IF stiff THEN power = power / obj(a, 2) ^ 2 * 4
  82.                     vec_x = vec_x + SIN(ang) * power
  83.                     vec_y = vec_y - COS(ang) * power
  84.                     vec_c = vec_c + 1
  85.                 END IF
  86.             END IF
  87.         NEXT t
  88.  
  89.         gravity = .1: lass = .98
  90.         obj(a, 4) = (obj(a, 4) + vec_x / vec_c) * lass: obj(a, 5) = (obj(a, 5) + vec_y / vec_c) * lass + gravity
  91.         obj(a, 6) = obj(a, 0) + obj(a, 4): obj(a, 7) = obj(a, 1) + obj(a, 5)
  92.     NEXT a
  93.     FOR a = 0 TO obj_c - 1: FOR t = 0 TO 1: obj(a, t) = obj(a, 6 + t): NEXT t, a
  94.  
  95.  
  96.  
  97.  
  98.     'mouse moving activing
  99.     mw = 0: WHILE _MOUSEINPUT: mw = mw + _MOUSEWHEEL: WEND
  100.     IF moving = 0 THEN
  101.         mouse_near = -1
  102.         FOR a = 0 TO obj_c - 1
  103.             disx = (cent_x + obj(a, 0) * zoom) - _MOUSEX
  104.             disy = (cent_y + obj(a, 1) * zoom) - _MOUSEY
  105.             dis = SQR(disx * disx + disy * disy)
  106.             IF dis < 10 THEN mouse_near = a
  107.         NEXT a
  108.     END IF
  109.  
  110.     IF (_MOUSEBUTTON(1) AND mouse_near <> -1) OR moving THEN
  111.         moving = 1
  112.         obj(mouse_near, 0) = (_MOUSEX - cent_x) / zoom
  113.         obj(mouse_near, 1) = (_MOUSEY - cent_y) / zoom
  114.     END IF
  115.  
  116.     moving = moving AND _MOUSEBUTTON(1)
  117.  
  118.     stiff = stiff XOR ABS(inkey2$ = "s")
  119.  
  120.     _DISPLAY: CLS
  121.  
  122.     'commands
  123.     inkey2$ = INKEY$: IF inkey2$ = CHR$(27) THEN END
  124.     IF inkey2$ = " " AND _MOUSEBUTTON(1) AND moving THEN obj(mouse_near, 3) = obj(mouse_near, 3) XOR 1
  125.  
  126.  
  128.  
  129.  
  130.  
  131.  
  132. FUNCTION degree (a, b): degreex = ATN(a / (b + .00001)) / pip180: degreex = degreex - 180 * ABS(0 > b): degreex = degreex - 360 * (degreex < 0): degree = degreex: END FUNCTION
  133.  
  134. FUNCTION randx (x): randx = x * RND(1) - x / 2: END FUNCTION

10
Programs / Garland of pearls
« on: March 21, 2022, 08:12:03 pm »
in principle, if I create a cobweb of points and connections, it will be a real cobweb in 3d that takes over all the vibrations of the spider


Code: QB64: [Select]
  1.  
  3. CONST pip180 = 3.141592 / 180
  4. monx = 800: mony = 800: mon = _NEWIMAGE(monx, mony, 32): SCREEN mon
  5.  
  6. gravity = 1
  7. obj_c = 15 'number objects
  8. DIM obj(obj_c - 1, 9) '0,1 posX,Y 3fix/mover
  9.  
  10.  
  11.  
  12. cent_x = monx / 2: cent_y = 20: zoom = 2
  13.  
  14. 'install "garland of pearls"
  15.  
  16.  
  17. 'garland
  18. min_weight = 6
  19. FOR a = 0 TO obj_c - 1
  20.     obj(a, 0) = (RND(1) - .5) * 100
  21.     obj(a, 1) = (RND(1) - .5) * 100
  22.     obj(a, 2) = min_weight + min_weight / 2 * RND(1)
  23.     obj(a, 3) = 1 '<---- 0 if fix, 1 if fly
  24.  
  26. obj(0, 0) = -64.5: obj(0, 1) = 5: obj(0, 3) = 0
  27. obj(obj_c - 1, 0) = 183.5: obj(obj_c - 1, 1) = -.5: obj(obj_c - 1, 3) = 0
  28.  
  29. 'connections
  30. conn_c = 100: DIM conn(conn_c - 1, 4)
  31. FOR a = 0 TO obj_c - 2
  32.     conn(a, 0) = a
  33.     conn(a, 1) = a + 1
  34.     conn(a, 2) = 20 'rope width
  35.     conn(a, 3) = 1
  37.  
  38.  
  39.  
  40.  
  41.  
  42. DO: _LIMIT 100
  43.  
  44.     'show info
  45.     wch = -1: FOR a = 0 TO obj_c - 1
  46.         IF _MOUSEX < 200 AND INT(_MOUSEY / 16) = a THEN wch = a: COLOR _RGB32(255, 255, 255) ELSE COLOR _RGB32(100, 100, 100)
  47.         LOCATE a + 1, 1: PRINT a; "# weight:"; INT(obj(a, 2) * 10) / 10;: IF obj(a, 3) = 0 THEN PRINT "FIX";
  48.         IF wch = a THEN PRINT "<--- using mousewheel !"
  49.         IF wch <> -1 THEN obj(wch, 2) = obj(wch, 2) + mw / 10: IF obj(wch, 2) < min_weight THEN obj(wch, 2) = min_weight
  50.     NEXT a
  51.     LOCATE mony / 16 - 1, 1: IF moving THEN PRINT "press SPACE to FIX/FLY"; ELSE PRINT "grab the ball with the mouse and move it!";
  52.  
  53.     'draw objects
  54.     FOR a = 0 TO obj_c - 1
  55.         x = cent_x + obj(a, 0) * zoom
  56.         y = cent_y + obj(a, 1) * zoom
  57.         IF mouse_near = a THEN COLOR _RGB32(255, 0, 0) ELSE COLOR _RGB32(200, 200, 200)
  58.         CIRCLE (x, y), obj(a, 2) * zoom
  59.         _PRINTSTRING (x, y), LTRIM$(STR$(a))
  60.     NEXT a
  61.  
  62.     'draw connections
  63.  
  64.     COLOR _RGB32(150, 150, 150)
  65.     FOR a = 0 TO conn_c - 1: IF conn(a, 3) = 0 THEN _CONTINUE
  66.         x1 = obj(conn(a, 0), 0) * zoom + cent_x
  67.         y1 = obj(conn(a, 0), 1) * zoom + cent_y
  68.         x2 = obj(conn(a, 1), 0) * zoom + cent_x
  69.         y2 = obj(conn(a, 1), 1) * zoom + cent_y
  70.     LINE (x1, y1)-(x2, y2): NEXT a
  71.  
  72.     'calculate moving
  73.  
  74.  
  75.     FOR a = 0 TO obj_c - 1: IF obj(a, 3) = 0 OR (mouse_near = a AND _MOUSEBUTTON(1)) THEN _CONTINUE
  76.  
  77.         'connections vector
  78.         vec_x = 0: vec_y = 0: vec_c = 0
  79.         FOR t = 0 TO conn_c - 1: IF conn(t, 3) = 0 THEN _CONTINUE
  80.             x = -1: FOR t2 = 0 TO 1: IF conn(t, t2) = a THEN x = conn(t, t2 XOR 1)
  81.             NEXT t2
  82.             IF x <> -1 THEN
  83.                 disx = obj(a, 0) - obj(x, 0)
  84.                 disy = obj(a, 1) - obj(x, 1)
  85.                 ang = (-degree(disx, disy) + 0) * pip180
  86.                 dis = SQR(disx * disx + disy * disy)
  87.                 power = dis / obj(a, 2) / 5
  88.                 '  IF dis < (obj(a, 2) + obj(x, 2)) THEN power = -power / 10000
  89.                 vec_x = vec_x + SIN(ang) * power
  90.                 vec_y = vec_y - COS(ang) * power
  91.                 vec_c = vec_c + 1
  92.             END IF
  93.         NEXT t
  94.  
  95.         gravity = .1
  96.         lass = .98
  97.  
  98.         obj(a, 4) = (obj(a, 4) + vec_x / vec_c) * lass
  99.         obj(a, 5) = (obj(a, 5) + vec_y / vec_c) * lass + gravity
  100.  
  101.         lx = obj(a, 0) + obj(a, 4)
  102.         ly = obj(a, 1) + obj(a, 5)
  103.  
  104.         obj(a, 0) = obj(a, 0) + obj(a, 4)
  105.         obj(a, 1) = obj(a, 1) + obj(a, 5)
  106.  
  107.  
  108.  
  109.     NEXT a
  110.  
  111.     mw = 0: WHILE _MOUSEINPUT: mw = mw + _MOUSEWHEEL: WEND
  112.  
  113.     'mouse moving activing
  114.     IF moving = 0 THEN
  115.         mouse_near = -1
  116.         FOR a = 0 TO obj_c - 1
  117.             disx = (cent_x + obj(a, 0) * zoom) - _MOUSEX
  118.             disy = (cent_y + obj(a, 1) * zoom) - _MOUSEY
  119.             dis = SQR(disx * disx + disy * disy)
  120.             IF dis < 10 THEN mouse_near = a
  121.         NEXT a
  122.     END IF
  123.  
  124.     IF (_MOUSEBUTTON(1) AND mouse_near <> -1) OR moving THEN
  125.         moving = 1
  126.         obj(mouse_near, 0) = (_MOUSEX - cent_x) / zoom
  127.         obj(mouse_near, 1) = (_MOUSEY - cent_y) / zoom
  128.     END IF
  129.  
  130.     moving = moving AND _MOUSEBUTTON(1)
  131.  
  132.  
  133.     'change weights with mousewheel
  134.  
  135.  
  136.  
  137.     _DISPLAY: CLS
  138.  
  139.     'commands
  140.     inkey2$ = INKEY$: IF inkey2$ = CHR$(27) THEN END
  141.     IF inkey2$ = " " AND _MOUSEBUTTON(1) AND moving THEN obj(mouse_near, 3) = obj(mouse_near, 3) XOR 1
  142.  
  143.  
  145.  
  146.  
  147.  
  148.  
  149. FUNCTION degree (a, b): degreex = ATN(a / (b + .00001)) / pip180: degreex = degreex - 180 * ABS(0 > b): degreex = degreex - 360 * (degreex < 0): degree = degreex: END FUNCTION
  150.  

11
Programs / Re: Double Pendulum Studies
« on: March 20, 2022, 06:02:56 pm »
Bplus, the "drawpully" is a very excellent SUB !!! :)

12
Programs / Re: Double Pendulum Studies
« on: March 20, 2022, 04:10:07 pm »
I would like to present the work of STxAxTIC. Man starts and something moves relative to something. the clock arms can be described in a few lines.

mon = _NEWIMAGE(600, 600, 32): SCREEN mon
r(0) = 12: r(1) = 60: r(2) = 60: r(3) = 100
DO
    t$ = TIME$ + ":" + RIGHT$("00" + LTRIM$(STR$(TIMER)), 2)
    x1 = _WIDTH(mon) / 2: y1 = _HEIGHT(mon) / 2
    FOR t = 0 TO 3
        ang = 360 / r(t) * VAL(MID$(t$, t * 3 + 1, 2)) * (3.141592 / 180)
        rad = (4 - t) ^ 1 * (_WIDTH(mon) * .05)
        x2 = SIN(ang) * rad + x1
        y2 = -COS(ang) * rad + y1
        LINE (x1, y1)-(x2, y2)
        x1 = x2
        y1 = y2
    NEXT t

    _DISPLAY
    CLS
LOOP





  Or something like this. Could be clarified, it is irrelevant now. I am writing an example to Bplus that there are two things here!
STxAxTIC was brave and climbed into the most cruel pots of physics. I was wondering if it could be simplified. I tried, but those results didn't come out. This cannot be achieved by "cheap" means.
  Congratulations to STxAxTIC for his patience and fanaticism for showing a great connection without neglecting the rigor of physics! Congratulations!

13
Programs / Re: Double Pendulum Studies
« on: March 20, 2022, 08:20:38 am »
I was thinking of a different approach. The way you calculate it is too complicated for me. I was wondering what if the masses only had a vector. No other. Strength and direction. An inertia. How that would change would determine what it connects to. The distances need to be somewhat flexible, so you don’t necessarily have to calculate with just one huge calculation who is exactly where. For now, it works with these values. I try to take the elasticity as low as possible, the length of the ropes will always be fixed (apparently). This will have the effect of being like a pendulum like yours.

Obj_c is the number of bullets. The obj array contains the position, vector.
The conn array contains the connections. For example.
conn (x, 0) index obj connected to conn (x, 1) obj, rope of conn (x, 2) length

Code: QB64: [Select]
  2. CONST pip180 = 3.141592 / 180
  3. monx = 800: mony = 600: mon = _NEWIMAGE(monx, mony, 32): SCREEN mon
  4.  
  5. gravity = 1
  6. obj_c = 10 'number objects
  7. DIM obj(obj_c - 1, 9) '0,1 posX,Y 3fix/mover
  8. rad = 10 'circle rad
  9.  
  10.  
  11. cent_x = monx / 2: cent_y = 20: zoom = .05
  12.  
  13.  
  14. 'objects location
  15. FOR a = 1 TO obj_c - 1
  16.     obj(a, 0) = (RND(1) - .5) * 10000
  17.     obj(a, 1) = (RND(1) - .5) * 5000
  18.     obj(a, 3) = 1 '<---- 0 if fix, 1 if fly
  20.  
  21. 'connections
  22. conn_c = obj_c - 1: DIM conn(conn_c - 1, 2)
  23. FOR a = 0 TO conn_c - 1
  24.     conn(a, 0) = a: conn(a, 1) = a + 1
  25.     conn(a, 2) = 30 'rope width
  27.  
  28.  
  29. DO: _LIMIT 100: IF INKEY$ = CHR$(27) THEN END
  30.     'draw objects
  31.     FOR a = 0 TO obj_c - 1
  32.         x = cent_x + obj(a, 0) * zoom
  33.         y = cent_y + obj(a, 1) * zoom
  34.         IF mouse_near = a THEN COLOR _RGB32(255, 0, 0) ELSE COLOR _RGB32(200, 200, 200)
  35.         CIRCLE (x, y), rad
  36.         _PRINTSTRING (x, y), LTRIM$(STR$(a))
  37.     NEXT a
  38.  
  39.     'draw connections
  40.     COLOR _RGB32(150, 150, 150)
  41.     FOR a = 0 TO conn_c - 1:
  42.         x1 = obj(conn(a, 0), 0) * zoom + cent_x
  43.         y1 = obj(conn(a, 0), 1) * zoom + cent_y
  44.         x2 = obj(conn(a, 1), 0) * zoom + cent_x
  45.         y2 = obj(conn(a, 1), 1) * zoom + cent_y
  46.     LINE (x1, y1)-(x2, y2): NEXT a
  47.  
  48.     'calculate moving
  49.     FOR a = 0 TO obj_c - 1: IF obj(a, 3) = 0 OR (mouse_near = a AND _MOUSEBUTTON(1)) THEN _CONTINUE
  50.  
  51.         'connections vector
  52.         vec_x = 0: vec_y = 0: vec_c = 0
  53.         FOR t = 0 TO conn_c - 1: x = -1: FOR t2 = 0 TO 1: IF conn(t, t2) = a THEN x = conn(t, t2 XOR 1)
  54.             NEXT t2
  55.             IF x <> -1 THEN
  56.                 disx = obj(a, 0) - obj(x, 0)
  57.                 disy = obj(a, 1) - obj(x, 1)
  58.                 ang = (-degree(disx, disy) + 0) * pip180
  59.                 dis = SQR(disx * disx + disy * disy)
  60.                 power = ABS(dis - conn(t, 2)) / 70
  61.                 vec_x = vec_x + SIN(ang) * power
  62.                 vec_y = vec_y - COS(ang) * power
  63.                 vec_c = vec_c + 1
  64.             END IF
  65.         NEXT t
  66.  
  67.         obj(a, 4) = (obj(a, 4) + vec_x / vec_c) * .99
  68.         obj(a, 5) = (obj(a, 5) + vec_y / vec_c + gravity) * .99
  69.         obj(a, 0) = obj(a, 0) + obj(a, 4)
  70.         obj(a, 1) = obj(a, 1) + obj(a, 5)
  71.  
  72.     NEXT a
  73.  
  75.  
  76.     'mouse moving activing
  77.     IF moving = 0 THEN
  78.         mouse_near = -1
  79.         FOR a = 0 TO obj_c - 1
  80.             disx = (cent_x + obj(a, 0) * zoom) - _MOUSEX
  81.             disy = (cent_y + obj(a, 1) * zoom) - _MOUSEY
  82.             dis = SQR(disx * disx + disy * disy)
  83.             IF dis < 10 THEN mouse_near = a
  84.         NEXT a
  85.     END IF
  86.  
  87.     IF (_MOUSEBUTTON(1) AND mouse_near <> -1) OR moving THEN
  88.         moving = 1
  89.         obj(mouse_near, 0) = (_MOUSEX - cent_x) / zoom
  90.         obj(mouse_near, 1) = (_MOUSEY - cent_y) / zoom
  91.     END IF
  92.  
  93.     moving = moving AND _MOUSEBUTTON(1)
  94.  
  95.     LOCATE 1, 1: PRINT "grab the ball with the mouse and move it!"
  96.     _DISPLAY: CLS
  97.  
  99.  
  100.  
  101.  
  102.  
  103. FUNCTION degree (a, b): degreex = ATN(a / (b + .00001)) / pip180: degreex = degreex - 180 * ABS(0 > b): degreex = degreex - 360 * (degreex < 0): degree = degreex: END FUNCTION
  104.  

14
Programs / Re: Double Pendulum Studies
« on: March 19, 2022, 03:21:24 pm »
I still didn’t delve into the math part. works fantastic with 3 wrists! Will I need to write a control for N joints separately? one for each?

15
Programs / Re: EVIL SPIDERS <--- only for the nervous system!
« on: March 19, 2022, 07:21:46 am »
I think I'm done with it. I have no more ideas. The walls collapse when energy is depleted, which is accompanied by a little claustrophobia in addition to arachnophobia. I stated on the difficulties, I think now everyone can find the level to be able to fulfill their mission.




Pages: [1] 2 3 ... 22