« previous next »
Pages: [1]   Go Down

Author Topic: Trivial 2D particle explosions I can't stop watching...  (Read 4564 times)

0 Members and 1 Guest are viewing this topic.

Offline STxAxTIC

  • Library Staff
  • Forum Resident
  • Posts: 1091
  • he lives
    • View Profile
Trivial 2D particle explosions I can't stop watching...
« on: February 02, 2019, 11:36:59 pm »
Code: QB64: [Select]
  1. '#lang "qb"              ' freebasic edit 2011-06
  2. delayconst = 900000 ' freebasic edit 2011-06
  3. 'CLEAR                  ' freebasic edit 2011-06
  4.  
  8.  
  9. LOCATE 1, 2: INPUT "Enter number of particles (default is 80): ", num
  10. IF num = 0 THEN num = 80
  11. DIM x(num), y(num), xold(num), yold(num), v0x(num), v0y(num), col(num)
  12.  
  13. start:
  15. iterations = 0
  16. 'g = RND * 10 + 20
  17. g = RND * 15 + 25
  18. xdamp = RND * .15 + .55
  19. ydamp = RND * .15 + .55
  20. exploderadius = 200 '75
  21. numobstacles = 0
  22. iterationmax = 1200
  23.  
  24. choosecol:
  25. bgcol = INT(RND * 14)
  26. wallcol = 0 'INT(RND * 14)'change to zero for spider mode
  27. IF bgcol = wallcol THEN GOTO choosecol
  28.  
  29. LINE (1, 1)-(639, 479), bgcol, BF
  30. LINE (1, 1)-(639, 479), wallcol, B
  31.  
  32. 'Draw obstacles randomly.
  33. FOR i = 1 TO numobstacles
  34.     LINE (RND * 640, RND * 480)-(RND * 640, RND * 480), wallcol, B
  36.  
  37. 'Make predetermined obstacles.
  38. 'LINE (50, 75)-(600, 125), wallcol, B
  39.  
  40. 'Toggle for random starting position.
  41. xshift = RND * 640
  42. yshift = RND * 480
  43. 'Toggle for fixed starting position
  44. 'xshift = 100
  45. 'yshift = 100
  46.  
  47. FOR i = 1 TO num
  48.     speed = RND * 90
  49.     ang1 = RND * 2 * 3.141592653589793#
  50.     ang2 = RND * 2 * 3.141592653589793#
  51.     x(i) = xshift + RND * exploderadius * COS(ang1)
  52.     y(i) = yshift + RND * exploderadius * SIN(ang1)
  53.     v0x(i) = 1.5 * speed * COS(ang2)
  54.     v0y(i) = speed * SIN(ang2)
  55.     dotcol:
  56.     col(i) = INT(RND * 13 + 1)
  57.     IF col(i) = bgcol OR col(i) = wallcol THEN GOTO dotcol
  58.     IF POINT(x(i), y(i)) = wallcol OR x(i) < 1 OR x(i) > 639 OR y(i) < 1 OR y(i) > 479 THEN i = i - 1
  59.     dv = SQR((v0x(i)) ^ 2 + (v0y(i)) ^ 2)
  60.     IF dv > vmax THEN vmax = dv
  61.     PSET (x(i), y(i)), col(i)
  63.  
  64. dt = .995 / vmax
  65. 'PRINT dt
  66.  
  68.  
  70.     idel = 0: DO: idel = idel + 1: LOOP UNTIL idel > delayconst ' freebasic edit 2011-06
  71.  
  72.     iterations = iterations + 1
  73.     smax = 0
  74.     FOR i = 1 TO num
  75.         xold(i) = x(i)
  76.         yold(i) = y(i)
  77.         v0x(i) = v0x(i) + 0 * dt
  78.         v0y(i) = v0y(i) + g * dt
  79.         xtmp = x(i) + v0x(i) * dt
  80.         ytmp = y(i) + v0y(i) * dt
  81.         IF POINT(xtmp, yold(i)) = wallcol THEN v0x(i) = v0x(i) * -1 * xdamp
  82.         IF POINT(xold(i), ytmp) = wallcol THEN v0y(i) = v0y(i) * -1 * ydamp
  83.         x(i) = x(i) + v0x(i) * dt
  84.         y(i) = y(i) + v0y(i) * dt
  85.         'Recolor stagnant particles.
  86.         xx = x(i) - xold(i)
  87.         yy = y(i) - yold(i)
  88.         IF SQR(xx ^ 2 + yy ^ 2) < .05 THEN col(i) = bgcol
  89.         PSET (xold(i), yold(i)), 0 'bgcol
  90.         PSET (x(i), y(i)), col(i)
  91.         ds = SQR((y(i) - yold(i)) ^ 2 + (x(i) - xold(i)) ^ 2)
  92.         IF ds > smax THEN smax = ds
  93.     NEXT
  94.     IF smax > .95 THEN dt = dt * (1 - .01)
  95.     IF smax < .9 THEN dt = dt * (1 + .01)
  96.     IF iterations > iterationmax THEN
  97.         SLEEP 2
  98.         GOTO start
  99.     END IF
  100.     LINE (19, 459)-(151, 471), wallcol, B
  101.     LINE (20, 460)-(20 + 130 * (iterations / iterationmax), 470), 15, BF
  104.  
Logged
You're not done when it works, you're done when it's right.

Offline bplus

  • Global Moderator
  • Forum Resident
  • Posts: 8053
  • b = b + ...
    • View Profile
Re: Trivial 2D particle explosions I can't stop watching...
« Reply #1 on: February 03, 2019, 03:45:16 am »
Here is a mod
Code: QB64: [Select]
  1. _TITLE "trivial 2D explosions" 'STxAxTIC mod B+ 2019-02-03
  2. CONST xmax = 1200
  3. CONST ymax = 700
  4. SCREEN _NEWIMAGE(xmax, ymax, 32)
  7.  
  8. np = 150 'npber of particles
  9. DIM x(np), y(np), xold(np), yold(np), v0x(np), v0y(np), col(np) AS _UNSIGNED LONG
  10.  
  11. nr = 1500 'number of rocks
  12. DIM rx(nr), ry(nr), rw(nr), rh(nr), rc(nr) AS _UNSIGNED LONG
  13. FOR i = 0 TO nr
  14.     rx(i) = -50 + RND * xmax
  15.     ry(i) = .66 * ymax + RND * .5 * ymax
  16.     rw(i) = 30 + RND * 150
  17.     rh(i) = 10 + RND * 40
  18.     r = 200 * RND
  19.     rc(i) = _RGB32(r, .45 * r, .2 * r)
  21.  
  22. no = 40
  23. DIM ox(no), oy(no), ow(no), oh(no)
  24. FOR i = 0 TO no
  25.     ox(i) = RND * xmax
  26.     oy(i) = .66 * ymax + RND * .5 * ymax
  27.     ow(i) = 30 + RND * 150
  28.     oh(i) = 10 + RND * 40
  30.  
  31.  
  33. wallcol = _RGB32(200, 100, 50)
  34. g = 55
  35. xdamp = .1
  36. ydamp = .1
  37. exploderadius = 10
  38.  
  39.  
  40. start:
  41. iterations = 0
  42. 'Toggle for random starting position.
  43. xshift = RND * xmax
  44. yshift = RND * ymax
  45.  
  46. 'sky
  47. FOR y = 0 TO ymax
  48.     LINE (0, y)-(xmax, y), _RGB32(.05 * i, .05 * i, .25 * i), BF
  50. 'rocks
  51. FOR i = 0 TO nr
  52.     LINE (rx(i), ry(i))-STEP(rw(i), rh(i)), rc(i), BF
  54.  
  55. 'Draw obstacles randomly
  56. FOR i = o TO no
  57.     LINE (ox(i), oy(i))-STEP(ow(i), oh(i)), wallcol, BF
  59.  
  60. FOR i = 1 TO np
  61.     speed = RND * 85 + 5
  62.     ang1 = RND * 2 * 3.141592653589793#
  63.     ang2 = RND * 2 * 3.141592653589793#
  64.     x(i) = xshift + RND * exploderadius * COS(ang1)
  65.     y(i) = yshift + RND * exploderadius * SIN(ang1)
  66.     v0x(i) = speed * COS(ang2)
  67.     v0y(i) = speed * SIN(ang2)
  68.     dotcol:
  69.     col(i) = _RGB32(RND * 255, RND * 255, RND * 255)
  70.     IF col(i) = bgcol OR col(i) = wallcol THEN GOTO dotcol
  71.     IF POINT(x(i), y(i)) = wallcol OR x(i) < 0 OR x(i) > xmax OR y(i) < 0 OR y(i) > ymax THEN i = i - 1
  72.     dv = SQR((v0x(i)) ^ 2 + (v0y(i)) ^ 2)
  73.     IF dv > vmax THEN vmax = dv
  74.     fcirc x(i), y(i), 3, col(i)
  77. dt = .995 / vmax
  78.  
  80.  
  81.     'sky
  82.     FOR y = 0 TO ymax
  83.         LINE (0, y)-(xmax, y), _RGB32(.05 * i, .05 * i, .25 * i), BF
  84.     NEXT
  85.     'rocks
  86.     FOR i = 0 TO nr
  87.         LINE (rx(i), ry(i))-STEP(rw(i), rh(i)), rc(i), BF
  88.     NEXT i
  89.  
  90.     'Draw obstacles randomly
  91.     FOR i = o TO no
  92.         LINE (ox(i), oy(i))-STEP(ow(i), oh(i)), wallcol, BF
  93.     NEXT i
  94.  
  95.  
  96.     iterations = iterations + 1
  97.     smax = 0
  98.     FOR i = 1 TO np
  99.         xold(i) = x(i)
  100.         yold(i) = y(i)
  101.         v0x(i) = v0x(i) + 0 * dt
  102.         v0y(i) = v0y(i) + g * dt
  103.         xtmp = x(i) + v0x(i) * dt
  104.         ytmp = y(i) + v0y(i) * dt
  105.         IF POINT(xtmp, yold(i)) = wallcol THEN v0x(i) = v0x(i) * -1 * xdamp
  106.         IF POINT(xold(i), ytmp) = wallcol THEN v0y(i) = v0y(i) * -1 * ydamp
  107.         x(i) = x(i) + v0x(i) * dt
  108.         y(i) = y(i) + v0y(i) * dt
  109.         fcirc x(i), y(i), 3, col(i)
  110.         ds = SQR((y(i) - yold(i)) ^ 2 + (x(i) - xold(i)) ^ 2)
  111.         IF ds > smax THEN smax = ds
  112.     NEXT
  113.     IF smax > .95 THEN dt = dt * (1 - .01)
  114.     IF smax < .9 THEN dt = dt * (1 + .01)
  115.     _DISPLAY
  116.     _LIMIT 200
  117.     IF iterations > 1500  THEN GOTO start
  119.  
  120.  
  121. SUB fcirc (CX AS LONG, CY AS LONG, R AS LONG, k AS _UNSIGNED LONG)
  122.     DIM subRadius AS LONG, RadiusError AS LONG
  123.     DIM X AS LONG, Y AS LONG
  124.  
  125.     subRadius = ABS(R)
  126.     RadiusError = -subRadius
  127.     X = subRadius
  128.     Y = 0
  129.  
  130.     IF subRadius = 0 THEN PSET (CX, CY), k: EXIT SUB
  131.  
  132.     ' Draw the middle span here so we don't draw it twice in the main loop,
  133.     ' which would be a problem with blending turned on.
  134.     LINE (CX - X, CY)-(CX + X, CY), k, BF
  135.  
  136.     WHILE X > Y
  137.         RadiusError = RadiusError + Y * 2 + 1
  138.         IF RadiusError >= 0 THEN
  139.             IF X <> Y + 1 THEN
  140.                 LINE (CX - Y, CY - X)-(CX + Y, CY - X), k, BF
  141.                 LINE (CX - Y, CY + X)-(CX + Y, CY + X), k, BF
  142.             END IF
  143.             X = X - 1
  144.             RadiusError = RadiusError - X * 2
  145.         END IF
  146.         Y = Y + 1
  147.         LINE (CX - X, CY - Y)-(CX + X, CY - Y), k, BF
  148.         LINE (CX - X, CY + Y)-(CX + X, CY + Y), k, BF
  149.     WEND
  151.  

hmm... this runs nicely for awhile and then freezes up?
« Last Edit: February 03, 2019, 04:01:43 am by bplus »
Logged

Offline STxAxTIC

  • Library Staff
  • Forum Resident
  • Posts: 1091
  • he lives
    • View Profile
Re: Trivial 2D particle explosions I can't stop watching...
« Reply #2 on: February 03, 2019, 08:31:20 am »
Glad to have caught your attention for a moment bplus. Always keeping it real with mods!

Bplus is not just a man, but an operator! The formula for making a fun mod is evidently:

bplus(program) = program++
Logged
You're not done when it works, you're done when it's right.

Offline bplus

  • Global Moderator
  • Forum Resident
  • Posts: 8053
  • b = b + ...
    • View Profile
Re: Trivial 2D particle explosions I can't stop watching...
« Reply #3 on: February 03, 2019, 12:07:33 pm »
Another mod :-))

Code: QB64: [Select]
  1. _TITLE "trivial 2D explosions B+ mod 2" 'STxAxTIC mod B+ 2019-02-03
  2. ' 2019-02-03 use a bg& for background instead of redraws each loop, that should cool down CPU
  3. '  rounder rocks, more fiddle with numbers
  4.  
  5. CONST xmax = 1200
  6. CONST ymax = 700
  7. SCREEN _NEWIMAGE(xmax, ymax, 32)
  8. _SCREENMOVE 100, 20
  10.  
  11. np = 500 'number of particles
  12. DIM x(np), y(np), xold(np), yold(np), v0x(np), v0y(np), col(np) AS _UNSIGNED LONG
  13.  
  14. nr = 3500 'number of rocks
  15. DIM rx(nr), ry(nr), rw(nr), rh(nr), rc(nr) AS _UNSIGNED LONG
  16. FOR i = 0 TO nr \ 2
  17.     rx(i) = RND * xmax
  18.     ry(i) = .5 * ymax + RND * .25 * ymax + rx(i) * .09
  19.     rw(i) = ry(i) * ry(i) * .00015
  20.     rh(i) = .3 * rw(i)
  21.     r = 200 * RND
  22.     rc(i) = _RGB32(r, .45 * r, .2 * r)
  24. FOR i = nr \ 2 + 1 TO nr
  25.     rx(i) = RND * xmax
  26.     ry(i) = .5 * ymax + RND * .75 * ymax + rx(i) * .09
  27.     rw(i) = ry(i) * ry(i) * .00015
  28.     rh(i) = .3 * rw(i)
  29.     r = 200 * RND
  30.     rc(i) = _RGB32(r, .45 * r, .2 * r)
  32.  
  33. no = 80 'number of rock bounce reflectors
  34. DIM ox(no), oy(no), ow(no), oh(no)
  35. FOR i = 0 TO no
  36.     ox(i) = RND * xmax
  37.     oy(i) = .5 * ymax + RND * .75 * ymax + ox(i) * .09
  38.     ow(i) = oy(i) * oy(i) * .00015
  39.     oh(i) = .3 * ow(i)
  41.  
  42.  
  44. wallcol = _RGB32(200, 100, 50)
  45. g = 95
  46. xdamp = .07
  47. ydamp = .07
  48. exploderadius = 10
  49.  
  50. 'draw background
  51. bgrd& = _NEWIMAGE(xmax, ymax, 32)
  52. _DEST bgrd&
  53. 'sky
  54. FOR y = 0 TO ymax
  55.     LINE (0, y)-(xmax, y), _RGB32(.1 * y, .1 * y, .15 * y), BF
  57. 'rocks
  58. FOR i = 0 TO nr
  59.     EllipseFill rx(i), ry(i), rw(i), rh(i), rc(i)
  60.     'LINE (rx(i), ry(i))-STEP(rw(i), rh(i)), rc(i), BF
  62.  
  63. 'Draw obstacles randomly
  64. FOR i = o TO no
  65.     EllipseFill ox(i), oy(i), ow(i), oh(i), wallcol
  66.     'LINE (ox(i), oy(i))-STEP(ow(i), oh(i)), wallcol, BF
  68.  
  70.  
  71. start:
  72. iterations = 0
  73. 'Toggle for random starting position.
  74. xshift = RND * xmax
  75. yshift = RND * ymax * .6
  76.  
  77.  
  78. FOR i = 1 TO np
  79.     speed = RND * 150 + 1
  80.     ang1 = RND * 2 * 3.141592653589793#
  81.     ang2 = RND * 2 * 3.141592653589793#
  82.     x(i) = xshift + RND * exploderadius * COS(ang1)
  83.     y(i) = yshift + RND * exploderadius * SIN(ang1)
  84.     v0x(i) = speed * COS(ang2)
  85.     v0y(i) = speed * SIN(ang2)
  86.     dotcol:
  87.     col(i) = _RGB32(RND * 255, RND * 255, RND * 255)
  88.     IF col(i) = bgcol OR col(i) = wallcol THEN GOTO dotcol
  89.     IF POINT(x(i), y(i)) = wallcol OR x(i) < 0 OR x(i) > xmax OR y(i) < 0 OR y(i) > ymax THEN i = i - 1
  90.     dv = SQR((v0x(i)) ^ 2 + (v0y(i)) ^ 2)
  91.     IF dv > vmax THEN vmax = dv
  93. dt = .995 / vmax
  94.  
  96.     _PUTIMAGE , bgrd&, 0
  97.     iterations = iterations + 1
  98.     smax = 0
  99.     FOR i = 1 TO np
  100.         xold(i) = x(i)
  101.         yold(i) = y(i)
  102.         v0x(i) = v0x(i) + .1 * dt
  103.         v0y(i) = v0y(i) + g * dt + .2 'more gravity
  104.         xtmp = x(i) + v0x(i) * dt
  105.         ytmp = y(i) + v0y(i) * dt
  106.         IF POINT(xtmp, yold(i)) = wallcol THEN v0x(i) = v0x(i) * -1 * xdamp
  107.         IF POINT(xold(i), ytmp) = wallcol THEN v0y(i) = v0y(i) * -1 * ydamp
  108.         x(i) = x(i) + v0x(i) * dt
  109.         y(i) = y(i) + v0y(i) * dt
  110.         EllipseFill x(i), y(i), 3, 3, col(i)
  111.         ds = SQR((y(i) - yold(i)) ^ 2 + (x(i) - xold(i)) ^ 2)
  112.         IF ds > smax THEN smax = ds
  113.     NEXT
  114.     IF smax > .95 THEN dt = dt * (1 - .01)
  115.     IF smax < .9 THEN dt = dt * (1 + .01)
  116.     _DISPLAY
  117.     _LIMIT 200
  118.     IF iterations > 1500 THEN GOTO start
  120.  
  121. ' with Steve's EllipseFill, who needs CircleFill?
  122. SUB EllipseFill (cx AS INTEGER, cy AS INTEGER, rx AS INTEGER, ry AS INTEGER, c AS _UNSIGNED LONG)
  123.     DIM a AS LONG, b AS LONG
  124.     DIM x AS LONG, y AS LONG
  125.     DIM xx AS LONG, yy AS LONG
  126.     DIM sx AS LONG, sy AS LONG
  127.     DIM e AS LONG
  128.  
  129.     a = 2 * rx * rx
  130.     b = 2 * ry * ry
  131.     x = rx
  132.     xx = ry * ry * (1 - rx - rx)
  133.     yy = rx * rx
  134.     sx = b * rx
  135.  
  136.     DO WHILE sx >= sy
  137.         LINE (cx - x, cy - y)-(cx + x, cy - y), c, BF
  138.         IF y <> 0 THEN LINE (cx - x, cy + y)-(cx + x, cy + y), c, BF
  139.  
  140.         y = y + 1
  141.         sy = sy + a
  142.         e = e + yy
  143.         yy = yy + a
  144.  
  145.         IF (e + e + xx) > 0 THEN
  146.             x = x - 1
  147.             sx = sx - b
  148.             e = e + xx
  149.             xx = xx + b
  150.         END IF
  151.     LOOP
  152.  
  153.     x = 0
  154.     y = ry
  155.     xx = rx * ry
  156.     yy = rx * rx * (1 - ry - ry)
  157.     e = 0
  158.     sx = 0
  159.     sy = a * ry
  160.  
  161.     DO WHILE sx <= sy
  162.         LINE (cx - x, cy - y)-(cx + x, cy - y), c, BF
  163.         LINE (cx - x, cy + y)-(cx + x, cy + y), c, BF
  164.  
  165.         DO
  166.             x = x + 1
  167.             sx = sx + b
  168.             e = e + xx
  169.             xx = xx + b
  170.         LOOP UNTIL (e + e + yy) > 0
  171.  
  172.         y = y - 1
  173.         sy = sy - a
  174.         e = e + yy
  175.         yy = yy + a
  176.  
  177.     LOOP
  178.  
  180.  
  181.  
« Last Edit: February 03, 2019, 12:16:33 pm by bplus »
Logged

Offline STxAxTIC

  • Library Staff
  • Forum Resident
  • Posts: 1091
  • he lives
    • View Profile
Re: Trivial 2D particle explosions I can't stop watching...
« Reply #4 on: February 03, 2019, 02:22:31 pm »
Man, that one's beautiful. Makes this thread earn its title.
Logged
You're not done when it works, you're done when it's right.

Offline Pete

  • Forum Resident
  • Posts: 2361
  • Cuz I sez so, varmint!
    • View Profile
Re: Trivial 2D particle explosions I can't stop watching...
« Reply #5 on: February 03, 2019, 03:51:56 pm »
Mark has too much 4th dimension on his hands... at least I hope that's time. Anyway, back to livin' in a 2D world! Now if someone can just invent a time machine and travel this stuff back to when Atari was born, they'd make a fortune. It's kind of a shame this kind of math is not my kind of thing, but it isn't. I mean I'm always amazed at how much happens on the screen with so few statements to make things happen; that's quite something.

Anyway, nice work guys...

Pete
Logged
Want to learn how to write code on cave walls? https://www.tapatalk.com/groups/qbasic/qbasic-f1/
Pages: [1]   Go Up
« previous next »