« previous next »
Pages: [1]   Go Down

Author Topic: Molecules simulation (simple particle system, part II)  (Read 3857 times)

0 Members and 1 Guest are viewing this topic.

This topic contains a post which is marked as Best Answer. Press here if you would like to see it.

FellippeHeitor

  • Guest
Molecules simulation (simple particle system, part II)
« on: July 24, 2020, 12:32:31 am »


Similarly to the previous video, I create a simple particle system but this time I get them to notice each other in a different way, to simulate atoms coming together to form molecules.

Code:
Code: QB64: [Select]
  1. SCREEN _NEWIMAGE(800, 600, 32)
  2.  
  4.  
  5. TYPE object
  6.     x AS SINGLE
  7.     y AS SINGLE
  8.     xv AS SINGLE
  9.     yv AS SINGLE
  10.     size AS INTEGER
  11.     c AS _UNSIGNED LONG
  13.  
  14. DIM SHARED o(130) AS object
  16.  
  17. 'initialize array
  18. FOR i = 1 TO UBOUND(o)
  19.     o(i).x = RND * _WIDTH
  20.     o(i).y = RND * _HEIGHT
  21.     o(i).xv = 1 - RND * 2
  22.     o(i).yv = 1 - RND * 2
  23.     o(i).size = 10
  24.     o(i).c = _RGB32(RND * 255, RND * 255, RND * 255, 100)
  26.  
  27. DIM SHARED minDistance AS INTEGER
  28. minDistance = 50
  29.  
  32.     mx = _MOUSEX
  33.  
  34.     minDistance = map(mx, 0, _WIDTH, 0, 200)
  35.  
  36.     CLS
  37.  
  38.     LINE (0, _HEIGHT - 20)-(mx, _HEIGHT), _RGB32(255, 50), BF
  39.  
  40.     FOR i = 1 TO UBOUND(o)
  41.         move o(i)
  42.         connect o(i)
  43.         CircleFill o(i).x, o(i).y, o(i).size, o(i).c
  44.     NEXT
  45.  
  46.     _DISPLAY
  47.     _LIMIT 60
  50.  
  51. FUNCTION map! (value!, minRange!, maxRange!, newMinRange!, newMaxRange!)
  52.     map! = ((value! - minRange!) / (maxRange! - minRange!)) * (newMaxRange! - newMinRange!) + newMinRange!
  54.  
  55. SUB connect (this AS object)
  56.     DIM j AS INTEGER
  57.     FOR j = 1 TO UBOUND(o)
  58.         IF dist(o(j), this) <= minDistance THEN
  59.             LINE (this.x, this.y)-(o(j).x, o(j).y), this.c
  60.         END IF
  61.     NEXT
  63.  
  64. FUNCTION dist! (o1 AS object, o2 AS object)
  65.     x1! = o1.x
  66.     y1! = o1.y
  67.     x2! = o2.x
  68.     y2! = o2.y
  69.     dist! = _HYPOT((x2! - x1!), (y2! - y1!))
  71.  
  72. SUB move (this AS object)
  73.     this.x = this.x + this.xv
  74.     IF this.x < 0 OR this.x > _WIDTH THEN this.xv = this.xv * -1
  75.  
  76.     this.y = this.y + this.yv
  77.     IF this.y < 0 OR this.y > _HEIGHT THEN this.yv = this.yv * -1
  79.  
  80.  
  81. SUB CircleFill (CX AS INTEGER, CY AS INTEGER, R AS INTEGER, C AS _UNSIGNED LONG)
  82.     ' CX = center x coordinate
  83.     ' CY = center y coordinate
  84.     '  R = radius
  85.     '  C = fill color
  86.     DIM Radius AS INTEGER, RadiusError AS INTEGER
  87.     DIM X AS INTEGER, Y AS INTEGER
  88.     Radius = ABS(R)
  89.     RadiusError = -Radius
  90.     X = Radius
  91.     Y = 0
  92.     IF Radius = 0 THEN PSET (CX, CY), C: EXIT SUB
  93.     LINE (CX - X, CY)-(CX + X, CY), C, BF
  94.     WHILE X > Y
  95.         RadiusError = RadiusError + Y * 2 + 1
  96.         IF RadiusError >= 0 THEN
  97.             IF X <> Y + 1 THEN
  98.                 LINE (CX - Y, CY - X)-(CX + Y, CY - X), C, BF
  99.                 LINE (CX - Y, CY + X)-(CX + Y, CY + X), C, BF
  100.             END IF
  101.             X = X - 1
  102.             RadiusError = RadiusError - X * 2
  103.         END IF
  104.         Y = Y + 1
  105.         LINE (CX - X, CY - Y)-(CX + X, CY - Y), C, BF
  106.         LINE (CX - X, CY + Y)-(CX + X, CY + Y), C, BF
  107.     WEND
  109.  
Logged

Offline Ashish

  • Forum Resident
  • Posts: 630
  • Never Give Up!
    • View Profile
Re: Molecules simulation (simple particle system, part II)
« Reply #1 on: July 24, 2020, 04:58:20 am »
See, what I found :D . It does something similar to your code. -
Code: QB64: [Select]
  1. 'Coded by Ashish in Qb64
  2. 'inspire from Shiffman's video.
  3. TYPE vector
  4.     x AS _FLOAT
  5.     y AS _FLOAT
  6.     z AS _FLOAT
  8.  
  9. TYPE particle
  10.     pos AS vector
  11.     rr AS _FLOAT
  12.     i AS INTEGER
  13.     j AS INTEGER
  14.     r AS INTEGER
  15.     g AS INTEGER
  16.     b AS INTEGER
  19. _TITLE "Processing Network"
  20. SCREEN _NEWIMAGE(800, 500, 32)
  21.  
  22. CONST speed = .05
  23. DIM SHARED Particles(100) AS particle
  24. DIM SHARED glAllow AS _BYTE
  25. FOR i = 0 TO UBOUND(Particles)
  26.     Particles(i).pos.x = RND * _WIDTH
  27.     Particles(i).pos.y = RND * _HEIGHT
  28.     Particles(i).pos.z = RND * 100
  29.     Particles(i).j = RND * 7
  30.     Particles(i).i = 5
  31.     IF RND * 2 > 1 THEN Particles(i).i = -Particles(i).i
  32.     Particles(i).rr = RND * 3
  33.     SELECT CASE Particles(i).j
  34.         CASE 0: Particles(i).r = 5: Particles(i).g = 205: Particles(i).b = 229
  35.         CASE 1: Particles(i).r = 255: Particles(i).g = 184: Particles(i).b = 3
  36.         CASE 2: Particles(i).r = 91: Particles(i).g = 3: Particles(i).b = 255
  37.         CASE 3: Particles(i).r = 61: Particles(i).g = 62: Particles(i).b = 62
  38.         CASE 4: Particles(i).r = 255: Particles(i).g = 0: Particles(i).b = 255
  39.         CASE 5: Particles(i).r = 255: Particles(i).g = 0: Particles(i).b = 0
  40.         CASE 6: Particles(i).r = 0: Particles(i).g = 255: Particles(i).b = 0
  41.         CASE 7: Particles(i).r = 0: Particles(i).g = 0: Particles(i).b = 255
  42.     END SELECT
  43.     IF RND * 2 > 1 THEN Particles(i).j = -Particles(i).j
  45.  
  46. glAllow = -1
  48.     CLS , _RGB(255, 255, 255)
  49.     FOR i = 0 TO UBOUND(particles)
  50.         CIRCLE (Particles(i).pos.x, Particles(i).pos.y), Particles(i).rr, _RGB(Particles(i).r, Particles(i).g, Particles(i).b)
  51.         PAINT (Particles(i).pos.x, Particles(i).pos.y), _RGB(Particles(i).r, Particles(i).g, Particles(i).b), _RGB(Particles(i).r, Particles(i).g, Particles(i).b)
  52.         Particles(i).pos.x = Particles(i).pos.x + Particles(i).j * speed
  53.         Particles(i).pos.y = Particles(i).pos.y + Particles(i).i * speed
  54.         IF Particles(i).pos.y > _HEIGHT - Particles(i).rr THEN Particles(i).i = -Particles(i).i
  55.         IF Particles(i).pos.y < 0 + Particles(i).rr THEN Particles(i).i = -Particles(i).i
  56.         IF Particles(i).pos.x > _WIDTH - Particles(i).rr THEN Particles(i).j = -Particles(i).j
  57.         IF Particles(i).pos.x < 0 + Particles(i).rr THEN Particles(i).j = -Particles(i).j
  58.     NEXT
  59.     '   _LIMIT 60
  60.     _DISPLAY
  62.  
  63.  
  65.     IF NOT glAllow THEN EXIT SUB
  66.     '  _glClearColor 1, 1, 1, 1
  67.     ' _glClear _GL_COLOR_BUFFER_BIT OR _GL_DEPTH_BUFFER_BIT
  68.  
  69.     _glEnable _GL_BLEND
  70.     '   _glColor4f 1, 1, 1, .05
  71.     FOR i = 0 TO UBOUND(particles)
  72.         FOR j = i + 1 TO UBOUND(particles)
  73.             IF distance(Particles(i).pos.x, Particles(i).pos.y, Particles(j).pos.x, Particles(j).pos.y) < 80 THEN
  74.                 FOR n = j + 1 TO UBOUND(particles)
  75.                     IF distance(Particles(j).pos.x, Particles(j).pos.y, Particles(n).pos.x, Particles(n).pos.y) < 80 THEN
  76.                         _glColor4f Particles(i).r / 255, Particles(i).g / 255, Particles(i).b / 255, .1
  77.                         _glBegin _GL_TRIANGLES
  78.                         _glVertex3f toGlX(Particles(i).pos.x), toGlY(Particles(i).pos.y), 0
  79.                         _glVertex3f toGlX(Particles(j).pos.x), toGlY(Particles(j).pos.y), 0
  80.                         _glVertex3f toGlX(Particles(n).pos.x), toGlY(Particles(n).pos.y), 0
  81.                         _glEnd
  82.                     END IF
  83.                 NEXT
  84.             END IF
  85.         NEXT
  86.     NEXT
  87.     _glFlush
  89.  
  90. FUNCTION toGlX## (x#)
  91.     x# = (_WIDTH / 2) - x#
  92.     toGlX## = map(x#, -(_WIDTH / 2), _WIDTH / 2, 1, -1)
  94.  
  95. FUNCTION toGlY## (y#)
  96.     y# = (_HEIGHT / 2) - y#
  97.     toGlY## = map(y#, -(_HEIGHT / 2), _HEIGHT / 2, -1, 1)
  99. FUNCTION toGlZ## (z#)
  100.     toGlZ## = map(z#, 0, 100, 0, 1)
  102.  
  103. FUNCTION map## (value##, minRange##, maxRange##, newMinRange##, newMaxRange##)
  104.     map## = ((value## - minRange##) / (maxRange## - minRange##)) * (newMaxRange## - newMinRange##) + newMinRange##
  106.  
  107. FUNCTION distance## (x1##, y1##, x2##, y2##)
  108.     IF x1## > x2## THEN dx## = x1## - x2## ELSE dx## = x2## - x1##
  109.     IF y1## > y2## THEN dy## = y1## - y2## ELSE dy## = y2## - y1##
  110.     distance## = SQR(dx## * dx## + dy## * dy##)
  112.  
  113.  
  114.  
Logged
if (Me.success) {Me.improve()} else {Me.tryAgain()}


My Projects - https://github.com/AshishKingdom?tab=repositories
OpenGL tutorials - https://ashishkingdom.github.io/OpenGL-Tutorials

FellippeHeitor

  • Guest
Re: Molecules simulation (simple particle system, part II)
« Reply #2 on: July 24, 2020, 05:34:47 am »
Aha! That's why this was at the back of my mind! We'd written it before 😂

Thanks for sharing!
Logged
Marked as best answer by on October 06, 2024, 08:14:50 am

Offline bplus

  • Global Moderator
  • Forum Resident
  • Posts: 8053
  • b = b + ...
    • View Profile
Re: Molecules simulation (simple particle system, part II)
« Reply #3 on: July 24, 2020, 11:35:13 am »
  • Undo Best Answer
    • @FellippeHeitor Ah so that's that map function you are always using, very clear now that I've seen it applied. :)

    Also liked the mouse bar at bottom of screen.

    Next mod up is translucent triangular planes between 3 points.
    Logged

    FellippeHeitor

    • Guest
    Re: Molecules simulation (simple particle system, part II)
    « Reply #4 on: July 24, 2020, 12:10:54 pm »
    @bplus that's the one!

    Looking forward to your mod!
    Logged

    Offline bplus

    • Global Moderator
    • Forum Resident
    • Posts: 8053
    • b = b + ...
      • View Profile
    Re: Molecules simulation (simple particle system, part II)
    « Reply #5 on: July 24, 2020, 02:01:35 pm »
    Logged
    Pages: [1]   Go Up
    « previous next »