« previous next »
Pages: 1 [2]   Go Down

Author Topic: Particle engine  (Read 5650 times)

0 Members and 1 Guest are viewing this topic.

Offline STxAxTIC

  • Library Staff
  • Forum Resident
  • Posts: 1091
  • he lives
    • View Profile
Re: Particle engine
« Reply #15 on: August 30, 2020, 03:26:36 am »
Alright, this is the one you want to try...

Red objects are hot and make the particles rise.
Blue objects cool hot particles down quickly.
Green objects are one-way, upward-biased capillaries.

Draw red things with mouse1.
Draw blue things with mouse2.
Draw green things with mouse3.

I suggest drawing a shape in black space and then connecting it with green. See screenshot.

Code: QB64: [Select]
  2. SCREEN _NEWIMAGE(1024, 768, 32)
  4.  
  5. TYPE Vector
  6.     x AS DOUBLE
  7.     y AS DOUBLE
  9.  
  10. TYPE Pixel
  11.     active AS INTEGER
  12.     position AS Vector
  13.     velocity AS Vector
  14.     acceleration AS Vector
  15.     size AS DOUBLE
  16.     shadered AS DOUBLE
  17.     shadegreen AS DOUBLE
  18.     shadeblue AS DOUBLE
  19.     shadealpha AS DOUBLE
  21.  
  22. TYPE Object
  23.     position AS Vector
  24.     size AS DOUBLE
  25.     amplitude AS Vector
  26.     frequency AS DOUBLE
  27.     shadered AS DOUBLE
  28.     shadegreen AS DOUBLE
  29.     shadeblue AS DOUBLE
  30.     shadealpha AS DOUBLE
  32.  
  33. DIM SHARED PixelCloud(3000) AS Pixel
  34. DIM SHARED Level(5000) AS Object
  35. DIM SHARED ObjectCount AS INTEGER
  36. ObjectCount = 0
  37.  
  39. FOR k = 1 TO UBOUND(PixelCloud)
  40.     PixelCloud(k).active = 1
  41.     PixelCloud(k).size = 3
  42.     PixelCloud(k).acceleration.x = 0
  43.     PixelCloud(k).acceleration.y = 0
  44.     PixelCloud(k).velocity.x = 0
  45.     PixelCloud(k).velocity.y = 0
  46.     PixelCloud(k).position.x = (RND - .5) * _WIDTH * .9
  47.     PixelCloud(k).position.y = (RND - .5) * _HEIGHT * .9
  48.     CALL SetPixelShade(k, 0, 0, 255, 150)
  50.  
  51. CALL InitializeLevel
  52.  
  53. DIM SHARED whirlyparameter AS DOUBLE
  54. whirlyparameter = 0
  56.  
  57.     CALL UserInput
  58.     CLS
  59.     CALL DrawBorder
  60.     CALL DrawPixels
  61.     CALL DrawLevel
  62.     CALL Dynamics
  63.  
  65.     whirlyparameter = whirlyparameter + ATN(1) / 10
  66.     IF whirlyparameter > 16 * ATN(1) THEN whirlyparameter = whirlyparameter - 16 * ATN(1)
  67.     _DISPLAY
  68.     _LIMIT 30
  70.  
  72.  
  73. SUB Dynamics
  74.     DIM k AS INTEGER
  75.     FOR k = 1 TO UBOUND(PixelCloud)
  76.         CALL CalculateInfluence(k)
  77.         CALL UpdatePosition(k)
  78.     NEXT
  80.  
  81. SUB UpdatePosition (i AS INTEGER)
  82.     DIM dt AS DOUBLE
  83.     DIM damp AS DOUBLE
  84.     DIM brownian AS DOUBLE
  85.     dt = 2
  86.     damp = .25
  87.     brownian = .25
  88.     PixelCloud(i).velocity.x = damp * PixelCloud(i).velocity.x + dt * PixelCloud(i).acceleration.x
  89.     PixelCloud(i).velocity.y = damp * PixelCloud(i).velocity.y + dt * PixelCloud(i).acceleration.y
  90.     PixelCloud(i).position.x = PixelCloud(i).position.x + dt * PixelCloud(i).velocity.x + (RND - .5) * brownian
  91.     PixelCloud(i).position.y = PixelCloud(i).position.y + dt * PixelCloud(i).velocity.y + (RND - .5) * brownian
  92.     IF (PixelCloud(i).position.y <= -_HEIGHT / 2 + 2 * PixelCloud(i).size + 1) THEN
  93.         PixelCloud(i).position.y = _HEIGHT / 2 - 2 * PixelCloud(i).size
  94.     END IF
  96.  
  97. SUB DrawBorder
  98.     LINE (0, 0)-(_WIDTH, _HEIGHT), _RGB32(255, 0, 0, 255), BF
  99.     LINE (10, 10)-(_WIDTH - 10, _HEIGHT - 10), _RGB32(0, 0, 0, 255), BF
  101.  
  102. SUB DrawPixels
  103.     DIM k AS INTEGER
  104.     DIM x AS DOUBLE
  105.     DIM y AS DOUBLE
  106.     DIM s AS DOUBLE
  107.     FOR k = 1 TO UBOUND(PixelCloud)
  108.         IF (PixelCloud(k).active = 1) THEN
  109.             x = PixelCloud(k).position.x
  110.             y = PixelCloud(k).position.y
  111.             s = PixelCloud(k).size
  112.             CALL clinebf(x - s, y - s, x + s, y + s, _RGBA(PixelCloud(k).shadered, PixelCloud(k).shadegreen, PixelCloud(k).shadeblue, PixelCloud(k).shadealpha))
  113.         END IF
  114.     NEXT
  116.  
  117. SUB DrawLevel
  118.     DIM k AS INTEGER
  119.     DIM x AS DOUBLE
  120.     DIM y AS DOUBLE
  121.     DIM s AS DOUBLE
  122.     DIM t AS DOUBLE
  123.     t = whirlyparameter
  124.     FOR k = 1 TO ObjectCount
  125.         x = Level(k).position.x + Level(k).amplitude.x * SIN(Level(k).frequency * t)
  126.         y = Level(k).position.y + Level(k).amplitude.y * SIN(Level(k).frequency * t)
  127.         s = Level(k).size
  128.         CALL clinebf(x - s, y - s, x + s, y + s, _RGBA(Level(k).shadered, Level(k).shadegreen, Level(k).shadeblue, Level(k).shadealpha))
  129.     NEXT
  131.  
  132. FUNCTION cpoint& (x1 AS DOUBLE, y1 AS DOUBLE)
  133.     DIM TheReturn AS _UNSIGNED LONG
  134.     TheReturn = POINT(_WIDTH / 2 + x1, -y1 + _HEIGHT / 2)
  135.     cpoint = TheReturn
  137.  
  138. SUB clinebf (x1 AS DOUBLE, y1 AS DOUBLE, x2 AS DOUBLE, y2 AS DOUBLE, col AS _UNSIGNED LONG)
  139.     LINE (_WIDTH / 2 + x1, -y1 + _HEIGHT / 2)-(_WIDTH / 2 + x2 - 0, -y2 + _HEIGHT / 2 + 0), col, BF
  141.  
  142. SUB SetPixelShade (i AS INTEGER, r AS INTEGER, g AS INTEGER, b AS INTEGER, a AS INTEGER)
  143.     PixelCloud(i).shadered = r
  144.     PixelCloud(i).shadegreen = g
  145.     PixelCloud(i).shadeblue = b
  146.     PixelCloud(i).shadealpha = a
  148.  
  149. SUB SetObjectShade (i AS INTEGER, r AS INTEGER, g AS INTEGER, b AS INTEGER, a AS INTEGER)
  150.     Level(i).shadered = r
  151.     Level(i).shadegreen = g
  152.     Level(i).shadeblue = b
  153.     Level(i).shadealpha = a
  155.  
  156. SUB InitializeLevel
  157.     DIM xx AS DOUBLE
  158.     DIM yy AS DOUBLE
  159.     DIM mx AS DOUBLE
  160.     DIM my AS DOUBLE
  161.     FOR xx = -20 TO 20
  162.         FOR yy = -252 TO -250
  163.             mx = xx - xx MOD (3 * 2.5)
  164.             my = yy - yy MOD (3 * 2.5)
  165.             IF ((_RED32(cpoint(mx, my)) <> 255) AND (_GREEN32(cpoint(mx, my)) <> 255) AND (_BLUE32(cpoint(mx, my)) <> 255)) THEN
  166.                 ObjectCount = ObjectCount + 1
  167.                 Level(ObjectCount).size = 3
  168.                 Level(ObjectCount).position.x = mx
  169.                 Level(ObjectCount).position.y = my
  170.                 Level(ObjectCount).amplitude.x = 0
  171.                 Level(ObjectCount).amplitude.y = 70
  172.                 Level(ObjectCount).frequency = .5
  173.                 CALL SetObjectShade(ObjectCount, 0, 255, 0, 255)
  174.             END IF
  175.  
  176.         NEXT
  177.     NEXT
  178.     FOR xx = -27 TO -25
  179.         FOR yy = -250 - 70 TO -250 + 70
  180.             mx = xx - xx MOD (3 * 2.5)
  181.             my = yy - yy MOD (3 * 2.5)
  182.             IF ((_RED32(cpoint(mx, my)) <> 255) AND (_GREEN32(cpoint(mx, my)) <> 255) AND (_BLUE32(cpoint(mx, my)) <> 255)) THEN
  183.                 ObjectCount = ObjectCount + 1
  184.                 Level(ObjectCount).size = 3
  185.                 Level(ObjectCount).position.x = mx
  186.                 Level(ObjectCount).position.y = my
  187.                 Level(ObjectCount).amplitude.x = 0
  188.                 Level(ObjectCount).amplitude.y = 0
  189.                 Level(ObjectCount).frequency = 0
  190.                 CALL SetObjectShade(ObjectCount, 255, 0, 0, 255)
  191.             END IF
  192.  
  193.         NEXT
  194.     NEXT
  195.     FOR xx = 25 TO 27
  196.         FOR yy = -250 - 70 TO -250 + 70
  197.             mx = xx - xx MOD (3 * 2.5)
  198.             my = yy - yy MOD (3 * 2.5)
  199.             IF ((_RED32(cpoint(mx, my)) <> 255) AND (_GREEN32(cpoint(mx, my)) <> 255) AND (_BLUE32(cpoint(mx, my)) <> 255)) THEN
  200.                 ObjectCount = ObjectCount + 1
  201.                 Level(ObjectCount).size = 3
  202.                 Level(ObjectCount).position.x = mx
  203.                 Level(ObjectCount).position.y = my
  204.                 Level(ObjectCount).amplitude.x = 0
  205.                 Level(ObjectCount).amplitude.y = 0
  206.                 Level(ObjectCount).frequency = 0
  207.                 CALL SetObjectShade(ObjectCount, 255, 0, 0, 255)
  208.             END IF
  209.         NEXT
  210.     NEXT
  212.  
  213. SUB UserInput
  214.     DIM mb1 AS INTEGER
  215.     DIM mb2 AS INTEGER
  216.     DIM mb3 AS INTEGER
  217.     DIM mx AS DOUBLE
  218.     DIM my AS DOUBLE
  219.     DIM sz AS DOUBLE
  220.     IF (ObjectCount > 0) THEN
  221.         sz = Level(ObjectCount).size
  222.     ELSE
  223.         sz = 3
  224.     END IF
  225.     mb1 = 0
  226.     mb2 = 0
  227.     mb3 = 0
  229.         IF (_MOUSEBUTTON(1) = -1) AND (mb1 <> -1) THEN
  230.             mb1 = -1
  231.             mx = (_MOUSEX - _WIDTH / 2)
  232.             my = (-_MOUSEY + _HEIGHT / 2)
  233.             mx = mx - mx MOD (sz * 2.5)
  234.             my = my - my MOD (sz * 2.5)
  235.             IF ((_RED32(cpoint(mx, my)) <> 255) AND (_GREEN32(cpoint(mx, my)) <> 255) AND (_BLUE32(cpoint(mx, my)) <> 255)) THEN
  236.                 ObjectCount = ObjectCount + 1
  237.                 Level(ObjectCount).size = sz
  238.                 Level(ObjectCount).position.x = mx
  239.                 Level(ObjectCount).position.y = my
  240.                 CALL SetObjectShade(ObjectCount, 255, 0, 0, 255)
  241.             END IF
  242.         END IF
  243.         IF (_MOUSEBUTTON(2) = -1) AND (mb2 <> -1) THEN
  244.             mb2 = -1
  245.             mx = (_MOUSEX - _WIDTH / 2)
  246.             my = (-_MOUSEY + _HEIGHT / 2)
  247.             mx = mx - mx MOD (sz * 2.5)
  248.             my = my - my MOD (sz * 2.5)
  249.             IF ((_RED32(cpoint(mx, my)) <> 255) AND (_GREEN32(cpoint(mx, my)) <> 255) AND (_BLUE32(cpoint(mx, my)) <> 255)) THEN
  250.                 ObjectCount = ObjectCount + 1
  251.                 Level(ObjectCount).size = sz
  252.                 Level(ObjectCount).position.x = mx
  253.                 Level(ObjectCount).position.y = my
  254.                 CALL SetObjectShade(ObjectCount, 0, 0, 255, 255)
  255.             END IF
  256.         END IF
  257.         IF (_MOUSEBUTTON(3) = -1) AND (mb3 <> -1) THEN
  258.             mb3 = -1
  259.             mx = (_MOUSEX - _WIDTH / 2)
  260.             my = (-_MOUSEY + _HEIGHT / 2)
  261.             mx = mx - mx MOD (sz * 2.5)
  262.             my = my - my MOD (sz * 2.5)
  263.             IF ((_RED32(cpoint(mx, my)) <> 255) AND (_GREEN32(cpoint(mx, my)) <> 255) AND (_BLUE32(cpoint(mx, my)) <> 255)) THEN
  264.                 ObjectCount = ObjectCount + 1
  265.                 Level(ObjectCount).size = sz
  266.                 Level(ObjectCount).position.x = mx
  267.                 Level(ObjectCount).position.y = my
  268.                 Level(ObjectCount).frequency = 1
  269.                 CALL SetObjectShade(ObjectCount, 0, 255, 0, 255)
  270.             END IF
  271.         END IF
  272.     LOOP
  273.  
  274.     IF (INKEY$ = " ") THEN
  275.         ObjectCount = ObjectCount - 1
  276.         IF (ObjectCount < 0) THEN ObjectCount = 0
  277.     END IF
  278.  
  280.  
  281. SUB CalculateInfluence (i AS INTEGER)
  282.     DIM x AS DOUBLE
  283.     DIM y AS DOUBLE
  284.     DIM xr AS DOUBLE
  285.     DIM yr AS DOUBLE
  286.     DIM xg AS DOUBLE
  287.     DIM yg AS DOUBLE
  288.     DIM xb AS DOUBLE
  289.     DIM yb AS DOUBLE
  290.     DIM dx AS DOUBLE
  291.     DIM dy AS DOUBLE
  292.     DIM WPoint(9) AS _UNSIGNED LONG
  293.     DIM WRed(9) AS INTEGER
  294.     DIM WGreen(9) AS INTEGER
  295.     DIM WBlue(9) AS INTEGER
  296.     x = PixelCloud(i).position.x
  297.     y = PixelCloud(i).position.y
  298.     dx = 2 * PixelCloud(i).size
  299.     dy = 2 * PixelCloud(i).size
  300.     WPoint(7) = cpoint(x - dx, y + dy)
  301.     WPoint(8) = cpoint(x, y + dy)
  302.     WPoint(9) = cpoint(x + dx, y + dy)
  303.     WPoint(4) = cpoint(x - dx, y)
  304.     WPoint(6) = cpoint(x + dx, y)
  305.     WPoint(1) = cpoint(x - dx, y - dy)
  306.     WPoint(2) = cpoint(x, y - dy)
  307.     WPoint(3) = cpoint(x + dx, y - dy)
  308.     WRed(7) = _RED32(WPoint(7))
  309.     WRed(8) = _RED32(WPoint(8))
  310.     WRed(9) = _RED32(WPoint(9))
  311.     WRed(4) = _RED32(WPoint(4))
  312.     WRed(6) = _RED32(WPoint(6))
  313.     WRed(1) = _RED32(WPoint(1))
  314.     WRed(2) = _RED32(WPoint(2))
  315.     WRed(3) = _RED32(WPoint(3))
  316.     WGreen(7) = _GREEN32(WPoint(7))
  317.     WGreen(8) = _GREEN32(WPoint(8))
  318.     WGreen(9) = _GREEN32(WPoint(9))
  319.     WGreen(4) = _GREEN32(WPoint(4))
  320.     WGreen(6) = _GREEN32(WPoint(6))
  321.     WGreen(1) = _GREEN32(WPoint(1))
  322.     WGreen(2) = _GREEN32(WPoint(2))
  323.     WGreen(3) = _GREEN32(WPoint(3))
  324.     WBlue(7) = _BLUE32(WPoint(7))
  325.     WBlue(8) = _BLUE32(WPoint(8))
  326.     WBlue(9) = _BLUE32(WPoint(9))
  327.     WBlue(4) = _BLUE32(WPoint(4))
  328.     WBlue(6) = _BLUE32(WPoint(6))
  329.     WBlue(1) = _BLUE32(WPoint(1))
  330.     WBlue(2) = _BLUE32(WPoint(2))
  331.     WBlue(3) = _BLUE32(WPoint(3))
  332.  
  333.     xr = (WRed(6) - WRed(4) + (WRed(9) + WRed(3)) / SQR(2) - (WRed(7) + WRed(1)) / SQR(2))
  334.     yr = (WRed(8) - WRed(2) + (WRed(7) + WRed(9)) / SQR(2) - (WRed(1) + WRed(3)) / SQR(2))
  335.     'xg = (WGreen(6) - WGreen(4) + (WGreen(9) + WGreen(3)) / SQR(2) - (WGreen(7) + WGreen(1)) / SQR(2))
  336.     'yg = (WGreen(8) - WGreen(2) + (WGreen(7) + WGreen(9)) / SQR(2) - (WGreen(1) + WGreen(3)) / SQR(2))
  337.     xg = -(WGreen(6) - WGreen(4) + (WGreen(9) + WGreen(3)) / SQR(2) - (WGreen(7) + WGreen(1)) / SQR(2))
  338.     yg = -10 * (WGreen(8) - 0 + (WGreen(7) + WGreen(9)) / SQR(2) - (0 + 0) / SQR(2))
  339.     yg = 10 * (0 - WGreen(2) + (0 + 0) / SQR(2) - (WGreen(1) + WGreen(3)) / SQR(2))
  340.     xb = (WBlue(6) - WBlue(4) + (WBlue(9) + WBlue(3)) / SQR(2) - (WBlue(7) + WBlue(1)) / SQR(2))
  341.     yb = (WBlue(8) - WBlue(2) + (WBlue(7) + WBlue(9)) / SQR(2) - (WBlue(1) + WBlue(3)) / SQR(2))
  342.     x = xr + xg + xb
  343.     y = yr + yg + yb
  344.  
  345.     IF (xr ^ 2 + yr ^ 2 > xb ^ 2 + yb ^ 2) AND (xr ^ 2 + yr ^ 2 > 15) THEN
  346.         PixelCloud(i).shadered = PixelCloud(i).shadered + 128
  347.         IF (PixelCloud(i).shadered >= 255) THEN PixelCloud(i).shadered = 255
  348.         PixelCloud(i).shadeblue = PixelCloud(i).shadeblue - 128
  349.         IF (PixelCloud(i).shadeblue <= 0) THEN PixelCloud(i).shadeblue = 0
  350.     END IF
  351.  
  352.     IF (xb ^ 2 + yb ^ 2 > xr ^ 2 + yr ^ 2) THEN
  353.         PixelCloud(i).shadered = PixelCloud(i).shadered - 128
  354.         IF (PixelCloud(i).shadered <= 0) THEN PixelCloud(i).shadered = 0
  355.         PixelCloud(i).shadeblue = PixelCloud(i).shadeblue + 128
  356.         IF (PixelCloud(i).shadeblue >= 255) THEN PixelCloud(i).shadeblue = 255
  357.     END IF
  358.  
  359.     y = y - 2 * (PixelCloud(i).shadered - PixelCloud(i).shadeblue) / 255
  360.  
  361.     IF (ABS(x) < .001) THEN
  362.         PixelCloud(i).acceleration.x = 0
  363.     ELSE
  364.         PixelCloud(i).acceleration.x = -x / SQR(x * x + y * y)
  365.     END IF
  366.     IF (ABS(y) < .001) THEN
  367.         PixelCloud(i).acceleration.y = 0
  368.     ELSE
  369.         PixelCloud(i).acceleration.y = -y / SQR(x * x + y * y)
  370.     END IF
  371.  
  372.     PixelCloud(i).shadered = PixelCloud(i).shadered - 5
  373.     IF (PixelCloud(i).shadered <= 0) THEN PixelCloud(i).shadered = 0
  374.     PixelCloud(i).shadeblue = PixelCloud(i).shadeblue + 5
  375.     IF (PixelCloud(i).shadeblue >= 255) THEN PixelCloud(i).shadeblue = 255

ss.png
* ss.png (Filesize: 72.37 KB, Dimensions: 1024x768, Views: 202)
Logged
You're not done when it works, you're done when it's right.

Offline loudar

  • Newbie
  • Posts: 73
  • improve it bit by bit.
    • View Profile
Re: Particle engine
« Reply #16 on: September 04, 2020, 07:54:58 pm »
Built a little heat chamber which cooled down at the top :D

One thing I noticed is that particles cool down reaaally quickly, which is why I needed so much red to get it all the way to the top and keep it there.
Untitled-1.png
* Untitled-1.png (Filesize: 35.31 KB, Dimensions: 1024x794, Views: 187)
Logged
Check out what I do besides coding: http://loudar.myportfolio.com/

Offline STxAxTIC

  • Library Staff
  • Forum Resident
  • Posts: 1091
  • he lives
    • View Profile
Re: Particle engine
« Reply #17 on: September 04, 2020, 09:44:46 pm »
Thanks for trying it loudar!!!

About to post a new version and announce next changes over at https://www.qb64.org/forum/index.php?topic=2972.0

I love the screenshot too!
Logged
You're not done when it works, you're done when it's right.
Pages: 1 [2]   Go Up
« previous next »