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

Author Topic: U.S. Flag  (Read 4634 times)

0 Members and 1 Guest are viewing this topic.

Offline SierraKen

  • Forum Resident
  • Posts: 1454
    • View Profile
U.S. Flag
« on: May 12, 2020, 12:30:28 am »
This is the U.S. Flag I made using numbers and commands. I was able to only use the DATA lines for one single star and copied it enough times to be 50 stars on the flag. I did similar loops with Line commands for the stripes. Feel free to use any or all of this code in your own applications or games. I included a picture of it below.

Code: QB64: [Select]
  1. 'Made to honor the U.S. Flag.
  2. 'Feel free to use any or all of this code in your own applications or games.
  3. 'By Sierraken on May 11, 2020.
  4.  
  5. _LIMIT 300
  6. _TITLE "U.S. Flag"
  7. SCREEN _NEWIMAGE(800, 600, 32)
  8. x = 150
  9. y = 100
  10.  
  11. 'Sky
  12. PAINT (2, 2), _RGB32(0, 205, 255)
  13.  
  14. 'Stars
  15. LINE (x, y)-(x + 185, y + 130), _RGB32(0, 0, 255), BF
  16. FOR xx = 155 TO 345 STEP 32
  17.     FOR yy = 105 TO 220 STEP 28
  18.         FOR stary = 1 TO 10
  19.             FOR starx = 1 TO 12
  20.                 READ c
  21.                 PSET (starx + xx, stary + yy), _RGB32(c, c, 255)
  22.             NEXT starx
  23.         NEXT stary
  24.         RESTORE Stars
  25.     NEXT yy
  26. NEXT xx
  27. RESTORE Stars
  28. FOR xx = 172 TO 329 STEP 32
  29.     FOR yy = 118.9 TO 213.05 STEP 28
  30.         FOR stary = 1 TO 10
  31.             FOR starx = 1 TO 12
  32.                 READ c
  33.                 PSET (starx + xx, stary + yy), _RGB32(c, c, 255)
  34.             NEXT starx
  35.         NEXT stary
  36.         RESTORE Stars
  37.     NEXT yy
  38. NEXT xx
  39.  
  40. 'Stripes
  41. FOR rs = 100 TO 230 STEP 37.2
  42.     w = w + 1
  43.     LINE (335, rs)-(612.5, rs + 18.6), _RGB32(255, 0, 0), BF
  44.     IF w > 3 THEN GOTO nex:
  45.     LINE (335, rs + 18.6)-(612.5, rs + 37.2), _RGB32(255, 255, 255), BF
  46. NEXT rs
  47. nex:
  48. w = 0
  49. FOR rs = 230 TO 341.6 STEP 37.2
  50.     r = r + 1
  51.     LINE (150, rs)-(612.5, rs + 18.6), _RGB32(255, 255, 255), BF
  52.     IF r > 3 THEN GOTO nex2:
  53.     LINE (150, rs + 18.6)-(612.5, rs + 37.2), _RGB32(255, 0, 0), BF
  54. NEXT rs
  55. nex2:
  56. r = 0
  57.  
  58. 'Flag Pole
  59. FOR sz = .25 TO 10 STEP .25
  60.     CIRCLE (145, 80), sz, _RGB32(122, 128, 166)
  61. NEXT sz
  62. LINE (142, 80)-(147, 600), _RGB32(122, 128, 166), BF
  63.  
  64. 'Wait for Esc key.
  66.     a$ = INKEY$
  67. LOOP UNTIL a$ = CHR$(27)
  68.  
  69. Stars:
  70. DATA 0,0,0,0,0,255,255,0,0,0,0,0
  71. DATA 0,0,0,0,255,255,255,0,0,0,0,0
  72. DATA 255,255,255,255,255,255,255,255,255,255,255,255
  73. DATA 255,255,255,255,255,255,255,255,255,255,255,255
  74. DATA 0,0,0,255,255,255,255,255,255,0,0,0
  75. DATA 0,0,255,255,255,255,255,255,255,255,0,0
  76. DATA 0,255,255,255,255,255,255,255,255,255,255,0
  77. DATA 255,255,255,0,0,0,0,0,0,0,255,255,255
  78. DATA 255,255,0,0,0,0,0,0,0,0,255,255
  79. DATA 255,0,0,0,0,0,0,0,0,0,0,255
  80.  

Ken's US Flag Program.jpg
* Ken's US Flag Program.jpg (Filesize: 359.2 KB, Dimensions: 1909x988, Views: 335)
« Last Edit: May 12, 2020, 12:55:48 am by SierraKen »
Logged

Offline Ashish

  • Forum Resident
  • Posts: 630
  • Never Give Up!
    • View Profile
Re: U.S. Flag
« Reply #1 on: May 12, 2020, 04:45:23 am »
Nice! Now, make it wave around.
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

Offline Pete

  • Forum Resident
  • Posts: 2361
  • Cuz I sez so, varmint!
    • View Profile
Re: U.S. Flag
« Reply #2 on: May 12, 2020, 05:19:00 am »
Quote from: Ashish on May 12, 2020, 04:45:23 am
Nice! Now, make it wave around.

You guys just make it too easy for me. :D

Pete
Logged
Want to learn how to write code on cave walls? https://www.tapatalk.com/groups/qbasic/qbasic-f1/

Offline TerryRitchie

  • Seasoned Forum Regular
  • Posts: 495
  • Semper Fidelis
    • View Profile
Re: U.S. Flag
« Reply #3 on: May 12, 2020, 08:25:16 am »
Logged
In order to understand recursion, one must first understand recursion.

Offline SierraKen

  • Forum Resident
  • Posts: 1454
    • View Profile
Re: U.S. Flag
« Reply #4 on: May 12, 2020, 12:08:36 pm »
Thanks guys.
Logged

Offline Pete

  • Forum Resident
  • Posts: 2361
  • Cuz I sez so, varmint!
    • View Profile
Re: U.S. Flag
« Reply #5 on: May 12, 2020, 01:30:38 pm »
Nice flag design Ken, and cool tutorial, Terryl!

When TheBOB had to do all this stuff with data fields, like Ken's approach, well, that was impressive, but it was also very time consuming, especially when it got into doing sprites for animation. The added ability of QB64 to make use of pre-made images, as in Terry's waiving flag tutorial, shown below, could almost get a stubborn SCREEN 0 Hero, like me, interested in graphics programming.

Looking forward to task 15b, Chinese flag on fire!

 
task15figure12.gif


Pete
Logged
Want to learn how to write code on cave walls? https://www.tapatalk.com/groups/qbasic/qbasic-f1/

Offline SierraKen

  • Forum Resident
  • Posts: 1454
    • View Profile
Re: U.S. Flag
« Reply #6 on: May 12, 2020, 02:46:08 pm »
That's awesome Pete! Yeah I thought about animating mine but it wouldn't come out good I think. But that's OK, it still can be used in programs that people want little or no added files to theirs.
Logged

Offline SierraKen

  • Forum Resident
  • Posts: 1454
    • View Profile
Re: U.S. Flag
« Reply #7 on: May 12, 2020, 03:21:02 pm »
Stxaxtic had a great wavy flag awhile back. I just tried the link on this post of Petr's flag, but the link no long works on there.
Here is the original post, but I will also put Stxaxtic's code on here, its worth checking out.


https://www.qb64.org/forum/index.php?topic=44.msg276#msg276


Code: QB64: [Select]
  1. '#lang "qb" '...for freebasic compiler. Keep disabled for QB64.
  2.  
  3. ' *** Video settings. ***
  4.  
  5. screenwidth = 640
  6. screenheight = 480
  8.  
  9. ' *** Performance settings. ***
  10.  
  11. bignumber = 500000 ' Maximum objects per array (determined by memory).
  12. globaldelayinit = 1000 ' Loop damping factor (fine adjustment below).
  14.  
  15. ' *** Initialize counters and array sizes. ***
  16.  
  17. numparticleorig = bignumber
  18. numparticlevisible = bignumber
  19.  
  20. ' *** Define basis arrays and structures. ***
  21.  
  22. ' Screen vectors in three-space.
  23. ' These vectors define the camera angle.
  24. DIM uhat(1 TO 3), vhat(1 TO 3), nhat(1 TO 3)
  25.  
  26. ' Basis vectors defined in three-space.
  27. DIM xhat(1 TO 4), yhat(1 TO 4), zhat(1 TO 4)
  28. xhat(1) = 1: xhat(2) = 0: xhat(3) = 0: xhat(4) = 4
  29. yhat(1) = 0: yhat(2) = 1: yhat(3) = 0: yhat(4) = 2
  30. zhat(1) = 0: zhat(2) = 0: zhat(3) = 1: zhat(4) = 1
  31.  
  32. ' Basis vectors projected into uv two-space.
  33. DIM xhatp(1 TO 2), yhatp(1 TO 2), zhatp(1 TO 2)
  34. DIM xhatp.old(1 TO 2), yhatp.old(1 TO 2), zhatp.old(1 TO 2)
  35.  
  36. ' *** Define particle arrays and structures. ***
  37.  
  38. ' Particle vectors defined in three-space.
  39. DIM vec(numparticleorig, 4)
  40. DIM vecdotnhat(numparticleorig)
  41. DIM vecdotnhatunit.old(numparticleorig)
  42. DIM vecvisible(numparticlevisible, 4)
  43. DIM vecvisibledotnhat(numparticlevisible)
  44. DIM vecvisibledotnhatunit.old(numparticlevisible)
  45.  
  46. ' Particle vectors projected onto infinite uv two-space.
  47. DIM vecpuv(numparticleorig, 1 TO 2)
  48. DIM vecpuv.old(numparticleorig, 1 TO 2)
  49. DIM vecvisiblepuv(numparticlevisible, 1 TO 2)
  50. DIM vecvisiblepuv.old(numparticlevisible, 1 TO 2)
  51.  
  52. ' Particle projections adjusted for screen uv two-space.
  53. DIM vecpuvs(numparticleorig, 1 TO 2)
  54. DIM vecpuvs.old(numparticleorig, 1 TO 2)
  55. DIM vecvisiblepuvs(numparticlevisible, 1 TO 2)
  56. DIM vecvisiblepuvs.old(numparticlevisible, 1 TO 2)
  57.  
  58. ' *** Define specialized arrays and structures. ***
  59.  
  60. ' Arrays for tech 2 mesh plotting.
  61. DIM plotflag(numparticleorig)
  62. DIM plotflag.old(numparticleorig)
  63.  
  64. ' *** Set mathematical constants. ***
  65.  
  66. pi = 3.1415926536
  67. ee = 2.7182818285
  68.  
  69. mainstart:
  70.  
  71. ' *** Initialize user input variables. ***
  72. KEY$ = ""
  73. 'mousekey$ = ""
  74.  
  75. ' *** Process first user input. ***
  76.  
  77. a$ = "13"
  78.  
  80.     CASE "13": genscheme$ = "wave2d": plotmode$ = "simplemesh"
  82.  
  83. substart:
  84.  
  85. ' *** Zero counters and array sizes. ***
  86.  
  87. numparticleorig = 0
  88. numparticlevisible = 0
  89.  
  90. ' *** Constants and switch control. ***
  91.  
  92. ' Perspective and animation switches/defaults.
  93. centerx = screenwidth / 2
  94. centery = screenheight / 2
  95. speedconst = 50
  96. falsedepth = .01
  97. zoom = 30
  98. timevar = 0
  99. T = 0
  100. togglehud = 1
  101. toggleatomnumbers = -1
  102. toggletimeanimate = -1
  103. toggletimealert = 0
  104. camx = 0
  105. camy = 0
  106. camz = 0
  107. uhat(1) = -3: uhat(2) = 5: uhat(3) = 1 / 4
  108. vhat(1) = -1: vhat(2) = -1: vhat(3) = 8
  109.  
  111.  
  112. SELECT CASE LCASE$(genscheme$)
  113.     CASE "wave2d":
  114.         globaldelay = globaldelayinit * 1500
  115.         uhat(1) = .7802773: uhat(2) = -.4759201: uhat(3) = .4058135
  116.         vhat(1) = .2502121: vhat(2) = .8321912: vhat(3) = .4948249
  117.         togglehud = 1
  118.         toggletimealert = 1
  119.         GOSUB genscheme.wave2d.init
  120.         numparticleorig = pcountparticleorig
  121.         REDIM vec2dz(xrange, yrange)
  122.         REDIM vec2dztemp(xrange, yrange)
  123.         REDIM vec2dzprev(xrange, yrange)
  124.         GOSUB genschemeUSAcolors
  125.         GOSUB genscheme.wave2d.gridinit
  127.  
  128. ' Move objects to accomodate initial camera position.
  129. IF camx <> 0 AND camy <> 0 AND camz <> 0 THEN
  130.     FOR i = 1 TO numparticleorig
  131.         vec(i, 1) = vec(i, 1) + camx
  132.         vec(i, 2) = vec(i, 2) + camy
  133.         vec(i, 3) = vec(i, 3) + camz
  134.     NEXT
  136.  
  137. GOSUB redraw
  138.  
  139. ' *** Begin main loop. ***
  141.     IF toggletimeanimate = 1 THEN
  142.         GOSUB timeanimate
  143.         flagredraw = 1
  144.     END IF
  145.     IF flagredraw = 1 THEN
  146.         GOSUB redraw
  147.         flagredraw = -1
  148.     END IF
  149.     GOSUB keyprocess
  150.     IF toggletimeanimate = 1 THEN
  151.         FOR delaycount = 1 TO globaldelay: NEXT
  152.     END IF
  154. ' *** End main loop. ***
  155.  
  156. ' *** Begin function definitions. ***
  157.  
  158. ' Comment out the conents of this gosub for non-QB64 compiler.
  159. mouseprocess:
  161.  
  162. keyprocess:
  164. IF KEY$ <> "" THEN
  165.     flagredraw = 1
  168.     CASE "8":
  169.         GOSUB rotate.vhat.plus
  170.     CASE "2":
  171.         GOSUB rotate.vhat.minus
  172.     CASE "4":
  173.         GOSUB rotate.uhat.minus
  174.     CASE "6":
  175.         GOSUB rotate.uhat.plus
  176.     CASE "7":
  177.         GOSUB rotate.clockwise
  178.     CASE "9":
  179.         GOSUB rotate.counterclockwise
  180.     CASE "1":
  181.         GOSUB rotate.uhat.minus: GOSUB normalize.SCREEN.vectors: GOSUB rotate.clockwise
  182.     CASE "3":
  183.         GOSUB rotate.uhat.plus: GOSUB normalize.SCREEN.vectors: GOSUB rotate.counterclockwise
  184.     CASE "q"
  185.         GOSUB strafe.objects.vhat.plus
  186.         GOSUB strafe.camera.vhat.plus
  187.     CASE "e"
  188.         GOSUB strafe.objects.vhat.minus
  189.         GOSUB strafe.camera.vhat.minus
  190.     CASE "x"
  191.         uhat(1) = 0: uhat(2) = 1: uhat(3) = 0
  192.         vhat(1) = 0: vhat(2) = 0: vhat(3) = 1
  193.     CASE "y"
  194.         uhat(1) = 0: uhat(2) = 0: uhat(3) = 1
  195.         vhat(1) = 1: vhat(2) = 0: vhat(3) = 0
  196.     CASE "z"
  197.         uhat(1) = 1: uhat(2) = 0: uhat(3) = 0
  198.         vhat(1) = 0: vhat(2) = 1: vhat(3) = 0
  199.     CASE "t"
  200.         toggletimeanimate = -toggletimeanimate
  201.     CASE "r"
  202.         timevar = 0
  203.     CASE "f"
  204.         GOTO substart
  205.     CASE "g"
  206.         GOTO mainstart
  207.     CASE "-"
  208.         globaldelay = globaldelay * 1.1
  209.     CASE "="
  210.         globaldelay = globaldelay * 0.9
  211.     CASE "`"
  212.         globaldelayinit = globaldelay
  213.     CASE " "
  214.         togglehud = -togglehud
  215.         CLS
  216.     CASE CHR$(27)
  217.         END
  220.  
  221. convert:
  222. ' Convert graphics from uv-cartesian coordinates to monitor coordinates.
  223. x0 = x: y0 = y
  224. x = x0 + centerx
  225. y = -y0 + centery
  227.  
  228. ' *** Define functions for view translation and rotation. ***
  229.  
  230. rotate.uhat.plus:
  231. uhat(1) = nhat(1) + speedconst * uhat(1)
  232. uhat(2) = nhat(2) + speedconst * uhat(2)
  233. uhat(3) = nhat(3) + speedconst * uhat(3)
  235.  
  236. rotate.uhat.minus:
  237. uhat(1) = -nhat(1) + speedconst * uhat(1)
  238. uhat(2) = -nhat(2) + speedconst * uhat(2)
  239. uhat(3) = -nhat(3) + speedconst * uhat(3)
  241.  
  242. rotate.vhat.plus:
  243. vhat(1) = nhat(1) + speedconst * vhat(1)
  244. vhat(2) = nhat(2) + speedconst * vhat(2)
  245. vhat(3) = nhat(3) + speedconst * vhat(3)
  247.  
  248. rotate.vhat.minus:
  249. vhat(1) = -nhat(1) + speedconst * vhat(1)
  250. vhat(2) = -nhat(2) + speedconst * vhat(2)
  251. vhat(3) = -nhat(3) + speedconst * vhat(3)
  253.  
  254. rotate.counterclockwise:
  255. v1 = vhat(1)
  256. v2 = vhat(2)
  257. v3 = vhat(3)
  258. vhat(1) = uhat(1) + speedconst * vhat(1)
  259. vhat(2) = uhat(2) + speedconst * vhat(2)
  260. vhat(3) = uhat(3) + speedconst * vhat(3)
  261. uhat(1) = -v1 + speedconst * uhat(1)
  262. uhat(2) = -v2 + speedconst * uhat(2)
  263. uhat(3) = -v3 + speedconst * uhat(3)
  265.  
  266. rotate.clockwise:
  267. v1 = vhat(1)
  268. v2 = vhat(2)
  269. v3 = vhat(3)
  270. vhat(1) = -uhat(1) + speedconst * vhat(1)
  271. vhat(2) = -uhat(2) + speedconst * vhat(2)
  272. vhat(3) = -uhat(3) + speedconst * vhat(3)
  273. uhat(1) = v1 + speedconst * uhat(1)
  274. uhat(2) = v2 + speedconst * uhat(2)
  275. uhat(3) = v3 + speedconst * uhat(3)
  278.  
  279. strafe.objects.uhat.plus:
  280. FOR i = 1 TO numparticleorig
  281.     vec(i, 1) = vec(i, 1) + uhat(1) * 1 / zoom
  282.     vec(i, 2) = vec(i, 2) + uhat(2) * 1 / zoom
  283.     vec(i, 3) = vec(i, 3) + uhat(3) * 1 / zoom
  286.  
  287. strafe.objects.uhat.minus:
  288. FOR i = 1 TO numparticleorig
  289.     vec(i, 1) = vec(i, 1) - uhat(1) * 1 / zoom
  290.     vec(i, 2) = vec(i, 2) - uhat(2) * 1 / zoom
  291.     vec(i, 3) = vec(i, 3) - uhat(3) * 1 / zoom
  294.  
  295. strafe.objects.vhat.plus:
  296. FOR i = 1 TO numparticleorig
  297.     vec(i, 1) = vec(i, 1) + vhat(1) * 1 / zoom
  298.     vec(i, 2) = vec(i, 2) + vhat(2) * 1 / zoom
  299.     vec(i, 3) = vec(i, 3) + vhat(3) * 1 / zoom
  302.  
  303. strafe.objects.vhat.minus:
  304. FOR i = 1 TO numparticleorig
  305.     vec(i, 1) = vec(i, 1) - vhat(1) * 1 / zoom
  306.     vec(i, 2) = vec(i, 2) - vhat(2) * 1 / zoom
  307.     vec(i, 3) = vec(i, 3) - vhat(3) * 1 / zoom
  310.  
  311. strafe.objects.nhat.plus:
  312. FOR i = 1 TO numparticleorig
  313.     vec(i, 1) = vec(i, 1) + nhat(1) * 1 / zoom
  314.     vec(i, 2) = vec(i, 2) + nhat(2) * 1 / zoom
  315.     vec(i, 3) = vec(i, 3) + nhat(3) * 1 / zoom
  318.  
  319. strafe.objects.nhat.minus:
  320. FOR i = 1 TO numparticleorig
  321.     vec(i, 1) = vec(i, 1) - nhat(1) * 1 / zoom
  322.     vec(i, 2) = vec(i, 2) - nhat(2) * 1 / zoom
  323.     vec(i, 3) = vec(i, 3) - nhat(3) * 1 / zoom
  326.  
  327. strafe.camera.uhat.plus:
  328. camx = camx + uhat(1) * 1 / zoom
  329. camy = camy + uhat(2) * 1 / zoom
  330. camz = camz + uhat(3) * 1 / zoom
  332.  
  333. strafe.camera.uhat.minus:
  334. camx = camx - uhat(1) * 1 / zoom
  335. camy = camy - uhat(2) * 1 / zoom
  336. camz = camz - uhat(3) * 1 / zoom
  338.  
  339. strafe.camera.vhat.plus:
  340. camx = camx + vhat(1) * 1 / zoom
  341. camy = camy + vhat(2) * 1 / zoom
  342. camz = camz + vhat(3) * 1 / zoom
  344.  
  345. strafe.camera.vhat.minus:
  346. camx = camx - vhat(1) * 1 / zoom
  347. camy = camy - vhat(2) * 1 / zoom
  348. camz = camz - vhat(3) * 1 / zoom
  350.  
  351. strafe.camera.nhat.plus:
  352. camx = camx + nhat(1) * 1 / zoom
  353. camy = camy + nhat(2) * 1 / zoom
  354. camz = camz + nhat(3) * 1 / zoom
  356.  
  357. strafe.camera.nhat.minus:
  358. camx = camx - nhat(1) * 1 / zoom
  359. camy = camy - nhat(2) * 1 / zoom
  360. camz = camz - nhat(3) * 1 / zoom
  362.  
  363. ' *** Define core functions. ***
  364.  
  365. timeanimate:
  366. timevar = timevar + 1
  367. IF timevar > 10 ^ 6 THEN timevar = 1
  368. SELECT CASE genscheme$
  369.     CASE "wave2d": GOSUB genscheme.wave2d.timeanimate
  372.  
  373. normalize.SCREEN.vectors:
  374. 'normalize the two vectors that define the screen orientation
  375. uhatmag = SQR(uhat(1) ^ 2 + uhat(2) ^ 2 + uhat(3) ^ 2)
  376. uhat(1) = uhat(1) / uhatmag: uhat(2) = uhat(2) / uhatmag: uhat(3) = uhat(3) / uhatmag
  377. vhatmag = SQR(vhat(1) ^ 2 + vhat(2) ^ 2 + vhat(3) ^ 2)
  378. vhat(1) = vhat(1) / vhatmag: vhat(2) = vhat(2) / vhatmag: vhat(3) = vhat(3) / vhatmag
  379. uhatdotvhat = uhat(1) * vhat(1) + uhat(2) * vhat(2) + uhat(3) * vhat(3)
  380. IF SQR(uhatdotvhat ^ 2) > .0005 THEN
  381.     CLS: COLOR 15: LOCATE 5, 5: PRINT "Screen vectors are not perpendicular. Press ESC to quit."
  382.     'DO: LOOP UNTIL INKEY$ = CHR$(27): END
  384. ' Compute the normal vector to the view plane.
  385. ' The normal vector points toward the eye, away from view frustum.
  386. nhat(1) = uhat(2) * vhat(3) - uhat(3) * vhat(2)
  387. nhat(2) = uhat(3) * vhat(1) - uhat(1) * vhat(3)
  388. nhat(3) = uhat(1) * vhat(2) - uhat(2) * vhat(1)
  389. nhatmag = SQR(nhat(1) ^ 2 + nhat(2) ^ 2 + nhat(3) ^ 2)
  390. nhat(1) = nhat(1) / nhatmag: nhat(2) = nhat(2) / nhatmag: nhat(3) = nhat(3) / nhatmag
  392.  
  393. redraw:
  394. GOSUB normalize.SCREEN.vectors
  395. ' Project the three-space basis vectors onto the screen plane.
  396. xhatp(1) = xhat(1) * uhat(1) + xhat(2) * uhat(2) + xhat(3) * uhat(3)
  397. xhatp(2) = xhat(1) * vhat(1) + xhat(2) * vhat(2) + xhat(3) * vhat(3)
  398. yhatp(1) = yhat(1) * uhat(1) + yhat(2) * uhat(2) + yhat(3) * uhat(3)
  399. yhatp(2) = yhat(1) * vhat(1) + yhat(2) * vhat(2) + yhat(3) * vhat(3)
  400. zhatp(1) = zhat(1) * uhat(1) + zhat(2) * uhat(2) + zhat(3) * uhat(3)
  401. zhatp(2) = zhat(1) * vhat(1) + zhat(2) * vhat(2) + zhat(3) * vhat(3)
  402. IF numparticleorig > 0 THEN
  403.     GOSUB compute.visible.particles
  404.     GOSUB project.particles
  405.     GOSUB depth.adjust.particles
  407. GOSUB DRAW.all.objects
  408. GOSUB store.SCREEN.projections
  410.  
  411. reverse.uvnhat:
  412. uhat(1) = -uhat(1)
  413. uhat(2) = -uhat(2)
  414. uhat(3) = -uhat(3)
  415. GOSUB normalize.SCREEN.vectors
  417.  
  418. compute.visible.particles:
  419. numparticlevisible = 0
  420. FOR i = 1 TO numparticleorig
  421.     IF falsedepth = 0 THEN
  422.     ELSE
  423.         numparticlevisible = numparticlevisible + 1
  424.         vecvisible(numparticlevisible, 1) = vec(i, 1)
  425.         vecvisible(numparticlevisible, 2) = vec(i, 2)
  426.         vecvisible(numparticlevisible, 3) = vec(i, 3)
  427.         vecvisible(numparticlevisible, 4) = vec(i, 4)
  428.     END IF
  431.  
  432. project.particles:
  433. ' Project object vectors onto the screen plane.
  434. FOR i = 1 TO numparticlevisible
  435.     vecvisibledotnhat(i) = vecvisible(i, 1) * nhat(1) + vecvisible(i, 2) * nhat(2) + vecvisible(i, 3) * nhat(3)
  436.     vecvisiblepuv(i, 1) = (vecvisible(i, 1) * uhat(1) + vecvisible(i, 2) * uhat(2) + vecvisible(i, 3) * uhat(3))
  437.     vecvisiblepuv(i, 2) = (vecvisible(i, 1) * vhat(1) + vecvisible(i, 2) * vhat(2) + vecvisible(i, 3) * vhat(3))
  440.  
  441. depth.adjust.particles:
  442. IF falsedepth = 0 THEN
  443.     FOR i = 1 TO numparticlevisible
  444.         vecvisiblepuvs(i, 1) = vecvisiblepuv(i, 1) * fovd / vecvisibledotnhat(i)
  445.         vecvisiblepuvs(i, 2) = vecvisiblepuv(i, 2) * fovd / vecvisibledotnhat(i)
  446.     NEXT
  448.     FOR i = 1 TO numparticlevisible
  449.         vecvisiblepuvs(i, 1) = vecvisiblepuv(i, 1) * (1 + falsedepth * vecvisibledotnhat(i))
  450.         vecvisiblepuvs(i, 2) = vecvisiblepuv(i, 2) * (1 + falsedepth * vecvisibledotnhat(i))
  451.     NEXT
  454.  
  455. DRAW.all.objects:
  456. SELECT CASE plotmode$
  457.     CASE "simplemesh": GOSUB plotmode.simplemesh
  460. LOCATE 28, 23: PRINT "SPACE = toggle HUD,  ESC = quit."
  461. IF togglehud = 1 THEN
  462.     ' Replace basis vector triad.
  463.     x = 50 * xhatp.old(1): y = 50 * xhatp.old(2): GOSUB convert
  464.     LINE (centerx, centery)-(x, y), 0
  465.     x = 50 * yhatp.old(1): y = 50 * yhatp.old(2): GOSUB convert
  466.     LINE (centerx, centery)-(x, y), 0
  467.     x = 50 * zhatp.old(1): y = 50 * zhatp.old(2): GOSUB convert
  468.     LINE (centerx, centery)-(x, y), 0
  469.     x = 50 * xhatp(1): y = 50 * xhatp(2): GOSUB convert
  470.     LINE (centerx, centery)-(x, y), xhat(4)
  471.     x = 50 * yhatp(1): y = 50 * yhatp(2): GOSUB convert
  472.     LINE (centerx, centery)-(x, y), yhat(4)
  473.     x = 50 * zhatp(1): y = 50 * zhatp(2): GOSUB convert
  474.     LINE (centerx, centery)-(x, y), zhat(4)
  475.     COLOR 14
  476.     LOCATE 25, 68: PRINT "-   VIEW   -"
  477.     COLOR 15
  478.     LOCATE 26, 68: PRINT "  8  "
  479.     LOCATE 27, 68: PRINT "4   6"
  480.     LOCATE 28, 68: PRINT "  2  "
  481.     COLOR 7
  482.     LOCATE 26, 75: PRINT "7   9"
  483.     LOCATE 27, 75: PRINT "     "
  484.     LOCATE 28, 75: PRINT "1   3"
  486. IF toggletimealert = 1 THEN
  487.     COLOR 7
  488.     LOCATE 1, 25: PRINT "Press 'T' to toggle animation."
  489.     LOCATE 2, 25: PRINT "Use numeric keypad to rotate."
  490.     LOCATE 3, 25: PRINT "      Press G to reset."
  493.  
  494. store.SCREEN.projections:
  495. xhatp.old(1) = xhatp(1): xhatp.old(2) = xhatp(2)
  496. yhatp.old(1) = yhatp(1): yhatp.old(2) = yhatp(2)
  497. zhatp.old(1) = zhatp(1): zhatp.old(2) = zhatp(2)
  498. FOR i = 1 TO numparticlevisible
  499.     vecvisiblepuvs.old(i, 1) = vecvisiblepuvs(i, 1)
  500.     vecvisiblepuvs.old(i, 2) = vecvisiblepuvs(i, 2)
  502. numparticlevisible.old = numparticlevisible
  504.  
  505. plotmode.simplemesh:
  506. FOR i = 1 TO numparticlevisible - 1
  507.     IF i MOD yrange <> 0 THEN 'point obeys normal neighbor-connect scheme
  508.         ' Erase old graphics.
  509.         x = zoom * vecvisiblepuvs.old(i, 1): y = zoom * vecvisiblepuvs.old(i, 2): GOSUB convert: x1 = x: y1 = y
  510.         x = zoom * vecvisiblepuvs.old(i + 1, 1): y = zoom * vecvisiblepuvs.old(i + 1, 2): GOSUB convert: x2 = x: y2 = y
  511.         x = zoom * vecvisiblepuvs.old(i + yrange, 1): y = zoom * vecvisiblepuvs.old(i + yrange, 2): GOSUB convert: x3 = x: y3 = y
  512.         LINE (x1, y1)-(x2, y2), 0
  513.         IF i < (numparticlevisible - yrange) THEN LINE (x1, y1)-(x3, y3), 0
  514.         ' Draw new graphics.
  515.         x = zoom * vecvisiblepuvs(i, 1): y = zoom * vecvisiblepuvs(i, 2): GOSUB convert: x1 = x: y1 = y
  516.         x = zoom * vecvisiblepuvs(i + 1, 1): y = zoom * vecvisiblepuvs(i + 1, 2): GOSUB convert: x2 = x: y2 = y
  517.         x = zoom * vecvisiblepuvs(i + yrange, 1): y = zoom * vecvisiblepuvs(i + yrange, 2): GOSUB convert: x3 = x: y3 = y
  518.         LINE (x1, y1)-(x2, y2), vecvisible(i, 4)
  519.         IF i < (numparticlevisible - yrange) THEN LINE (x1, y1)-(x3, y3), vecvisible(i, 4)
  520.     ELSE 'point does not obey normal neighbor-connect scheme
  521.         ' Erase old graphics.
  522.         x = zoom * vecvisiblepuvs.old(i, 1): y = zoom * vecvisiblepuvs.old(i, 2): GOSUB convert: x1 = x: y1 = y
  523.         x = zoom * vecvisiblepuvs.old(i + yrange, 1): y = zoom * vecvisiblepuvs.old(i + yrange, 2): GOSUB convert: x3 = x: y3 = y
  524.         IF i < (numparticlevisible - yrange + 1) THEN LINE (x1, y1)-(x3, y3), 0
  525.         ' Draw new graphics.
  526.         x = zoom * vecvisiblepuvs(i, 1): y = zoom * vecvisiblepuvs(i, 2): GOSUB convert: x1 = x: y1 = y
  527.         x = zoom * vecvisiblepuvs(i + yrange, 1): y = zoom * vecvisiblepuvs(i + yrange, 2): GOSUB convert: x3 = x: y3 = y
  528.         IF i < (numparticlevisible - yrange + 1) THEN LINE (x1, y1)-(x3, y3), vecvisible(i, 4)
  529.     END IF
  532.  
  533. genschemeUSAcolors:
  534. 'delinearize
  535. pcountparticleorig = 0
  536. FOR i = 1 TO xrange
  537.     FOR j = 1 TO yrange
  538.         pcountparticleorig = pcountparticleorig + 1
  539.         vec2dztemp(i, j) = vec(pcountparticleorig, 4)
  540.     NEXT
  542. FOR i = 1 TO xrange
  543.     FOR jj = 0 TO 12 STEP 2
  544.         FOR j = jj * yrange / 13 + 1 TO (jj + 1) * yrange / 13
  545.             vec2dztemp(i, j) = 4
  546.         NEXT
  547.     NEXT
  548.     FOR jj = 1 TO 11 STEP 2
  549.         FOR j = jj * yrange / 13 + 1 TO (jj + 1) * yrange / 13
  550.             vec2dztemp(i, j) = 7
  551.         NEXT
  552.     NEXT
  554. starflag = -1
  555. FOR i = 1 TO xrange * .76 / 1.9
  556.     FOR j = yrange TO yrange * .5385 STEP -1
  557.         IF starflag = 1 THEN
  558.             vec2dztemp(i, j) = 15
  559.         ELSE
  560.             vec2dztemp(i, j) = 1
  561.         END IF
  562.         starflag = -starflag
  563.     NEXT
  565. 'relinearize
  566. pcountparticleorig = 0
  567. FOR i = 1 TO xrange
  568.     FOR j = 1 TO yrange
  569.         pcountparticleorig = pcountparticleorig + 1
  570.         vec(pcountparticleorig, 4) = vec2dztemp(i, j)
  571.     NEXT
  574.  
  575. genscheme.wave2d.init:
  576. xl = -1.9: xr = 1.9
  577. yl = -1: yr = 1
  578. xl = xl * 4: xr = xr * 4: yl = yl * 4: yr = yr * 4
  579. Dx = .32
  580. Dy = .32
  581. xrange = 1 + INT((-xl + xr) / Dx)
  582. yrange = 1 + INT((-yl + yr) / Dy)
  583. alpha = .25
  584. pcountparticleorig = 0
  585. FOR i = xl TO xr STEP Dx
  586.     FOR j = yl TO yr STEP Dy
  587.         pcountparticleorig = pcountparticleorig + 1
  588.         vec(pcountparticleorig, 1) = i
  589.         vec(pcountparticleorig, 2) = j
  590.         vec(pcountparticleorig, 3) = 0
  591.         '*' vec(pcountparticleorig, 4) = 14 'use special color scheme
  592.     NEXT
  595.  
  596. genscheme.wave2d.gridinit:
  597. 'delinearize
  598. pcountparticleorig = 0
  599. FOR i = 1 TO xrange
  600.     FOR j = 1 TO yrange
  601.         pcountparticleorig = pcountparticleorig + 1
  602.         vec2dztemp(i, j) = vec(pcountparticleorig, 3)
  603.     NEXT
  605. 'initial position condition
  606. 'FOR i = 1 TO xrange 'random high points
  607. '  FOR j = 1 TO yrange
  608. '      nrand = RND * 1000
  609. '      IF nrand < 10 THEN
  610. '          vec2dz(i, j) = 5
  611. '      END IF
  612. '  NEXT
  613. 'NEXT
  614. 'vec2dz(xrange * .8, yrange * .2) = -2.5 'single plucked point
  615. 'FOR i = 1 TO xrange 'cross arm
  616. '   vec2dz(i, yrange / 3) = 2
  617. 'NEXT
  618. 'FOR j = 1 TO yrange 'cross arm
  619. '   vec2dz(xrange / 2, j) = 1
  620. 'NEXT
  621. 'sync
  622. pcountparticleorig = 0
  623. FOR i = 1 TO xrange
  624.     FOR j = 1 TO yrange
  625.         pcountparticleorig = pcountparticleorig + 1
  626.         vec2dzprev(i, j) = vec2dz(i, j)
  627.         vec2dztemp(i, j) = vec2dz(i, j)
  628.     NEXT
  630. 'initial velocity condition
  631. vec2dzprev(xrange * .8, yrange * .8) = 1.5 'single struck point
  632. 'relinearize
  633. pcountparticleorig = 0
  634. FOR i = 1 TO xrange
  635.     FOR j = 1 TO yrange
  636.         pcountparticleorig = pcountparticleorig + 1
  637.         vec(pcountparticleorig, 3) = vec2dz(i, j)
  638.     NEXT
  641.  
  642. genscheme.wave2d.timeanimate:
  643. 'delinearize
  644. pcountparticleorig = 0
  645. FOR i = 1 TO xrange
  646.     FOR j = 1 TO yrange
  647.         pcountparticleorig = pcountparticleorig + 1
  648.         vec2dztemp(i, j) = vec(pcountparticleorig, 3)
  649.     NEXT
  651. 'begin propagation process
  652. FOR i = 2 TO xrange - 1 'boDy, no edges
  653.     FOR j = 2 TO yrange - 1
  654.         wp1 = alpha * (vec2dztemp(i + 1, j) + vec2dztemp(i - 1, j)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  655.         wp2 = alpha * (vec2dztemp(i, j + 1) + vec2dztemp(i, j - 1)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  656.         vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  657.     NEXT
  659. 'comment out this section for fixed edges (or pieces of this section)
  660. i = 1 'left edge
  661. FOR j = 2 TO yrange - 1
  662.     wfp = vec2dztemp(i, j)
  663.     wp1 = alpha * (vec2dztemp(i + 1, j) + wfp) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  664.     wp2 = alpha * (vec2dztemp(i, j + 1) + vec2dztemp(i, j - 1)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  665.     'vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  667. i = xrange 'right edge
  668. FOR j = 2 TO yrange - 1
  669.     wfp = vec2dztemp(i, j)
  670.     wp1 = alpha * (wfp + vec2dztemp(i - 1, j)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  671.     wp2 = alpha * (vec2dztemp(i, j + 1) + vec2dztemp(i, j - 1)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  672.     vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  674. j = 1 'bottom edge
  675. FOR i = 2 TO xrange - 1
  676.     wfp = vec2dztemp(i, j)
  677.     wp2 = alpha * (vec2dztemp(i, j + 1) + wfp) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  678.     wp1 = alpha * (vec2dztemp(i + 1, j) + vec2dztemp(i - 1, j)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  679.     vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  681. j = yrange 'top edge
  682. FOR i = 2 TO xrange - 1
  683.     wfp = vec2dztemp(i, j)
  684.     wp2 = alpha * (wfp + vec2dztemp(i, j - 1)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  685.     wp1 = alpha * (vec2dztemp(i + 1, j) + vec2dztemp(i - 1, j)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  686.     vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  688. 'bottom left corner
  689. i = 1: j = 1
  690. wfp1 = vec2dztemp(i, j)
  691. wfp2 = vec2dztemp(i, j)
  692. wp1 = alpha * (vec2dztemp(i + 1, j) + wfp1) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  693. wp2 = alpha * (vec2dztemp(i, j + 1) + wfp2) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  694. 'vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  695. 'bottom right corner
  696. i = xrange: j = 1
  697. wfp1 = vec2dztemp(i, j)
  698. wfp2 = vec2dztemp(i, j)
  699. wp1 = alpha * (wfp1 + vec2dztemp(i - 1, j)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  700. wp2 = alpha * (vec2dztemp(i, j + 1) + wfp2) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  701. vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  702. 'top left corner
  703. i = 1: j = yrange
  704. wfp1 = vec2dztemp(i, j)
  705. wfp2 = vec2dztemp(i, j)
  706. wp1 = alpha * (vec2dztemp(i + 1, j) + wfp1) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  707. wp2 = alpha * (wfp2 + vec2dztemp(i, j - 1)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  708. 'vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  709. 'top right corner
  710. i = xrange: j = yrange
  711. wfp1 = vec2dztemp(i, j)
  712. wfp2 = vec2dztemp(i, j)
  713. wp1 = alpha * (wfp1 + vec2dztemp(i - 1, j)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  714. wp2 = alpha * (wfp2 + vec2dztemp(i, j - 1)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  715. vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  716. 'start special movements
  717. T = timevar / 5
  718. IF T < pi THEN 'wave left edge just once
  719.     FOR j = 1 TO yrange
  720.         i = 1
  721.         vec2dz(i, j) = (j / yrange) * 1.5 * SIN(T)
  722.     NEXT
  724. IF T < pi THEN 'wave bottom edge just once
  725.     FOR i = 1 TO xrange
  726.         j = 1
  727.         vec2dz(i, j) = (i / xrange) * .45 * SIN(T)
  728.     NEXT
  730. 'relinearize
  731. pcountparticleorig = 0
  732. FOR i = 1 TO xrange
  733.     FOR j = 1 TO yrange
  734.         pcountparticleorig = pcountparticleorig + 1
  735.         vec2dzprev(i, j) = vec2dztemp(i, j)
  736.         vec2dztemp(i, j) = vec2dz(i, j)
  737.         vec(pcountparticleorig, 3) = vec2dz(i, j)
  738.     NEXT
  741.  
  742.  
« Last Edit: May 12, 2020, 03:27:59 pm by SierraKen »
Logged

Offline bplus

  • Global Moderator
  • Forum Resident
  • Posts: 8053
  • b = b + ...
    • View Profile
Re: U.S. Flag
« Reply #8 on: May 12, 2020, 03:28:11 pm »
Quote from: SierraKen on May 12, 2020, 03:21:02 pm
I just found Static's wavy flag I have on my computer, which you can find on a reply to Petr's other flag post here:


https://www.qb64.org/forum/index.php?topic=44.msg276#msg276


That link leads to file not found? Do you have the code? Update: Why yes you do! :)

Here is nice one by Ashish mod of one I started:
https://www.qb64.org/forum/index.php?topic=1384.msg105761#msg105761
« Last Edit: May 12, 2020, 03:29:29 pm by bplus »
Logged

Offline SierraKen

  • Forum Resident
  • Posts: 1454
    • View Profile
Re: U.S. Flag
« Reply #9 on: May 12, 2020, 05:44:53 pm »
Incredible bplus and Ashish! Good job.
Logged

Offline _vince

  • Seasoned Forum Regular
  • Posts: 422
    • View Profile
Re: U.S. Flag
« Reply #10 on: May 13, 2020, 10:31:52 am »
Quote from: Pete on May 12, 2020, 01:30:38 pm
task15figure12.gif


Pete

I'd like to see this one but generated from a single flat 2D image in QB64, would be more impressive without GL.  Then you wouldn't even need the image because the stars and stripes are fairly easy to draw.
Logged

Offline _vince

  • Seasoned Forum Regular
  • Posts: 422
    • View Profile
Re: U.S. Flag
« Reply #11 on: May 15, 2020, 01:16:26 am »
Ok, I went for it.  I used these instructions for drawing the flag:  https://en.wikipedia.org/wiki/Flag_of_the_United_States#Specifications

Code: QB64: [Select]
  1. deflng a-z
  2.  
  3. sw = 640
  4. sh = 480
  5.  
  7. pi = 4*atn(1)
  8.  
  9. screen _newimage(sw, sh*2, 32)
  10.  
  11. h = 300
  12. w = 1.9*h
  13. a = h/7
  14.  
  15. img = _newimage(w, h, 32)
  16. _dest img
  17. x0 = 0
  18. y0 = 0
  19.  
  20. line (0, 0)-step(w, h),_rgb(255,255,255),bf
  21. for i=0 to 6
  22.         line (0, i*h*2/13)-step(w, h/13),_rgb(255*0.698,255*0.132,255*0.203),bf
  24. line (0, 0)-step(w*2/5, h*7/13),_rgb(255*0.234,255*0.233,255*0.430),bf
  25.  
  26. for i=0 to 4
  27. for j=0 to 5
  28.         starf (j*2 + 1)*w*2/(5*12), (i*2 + 1)*h*7/130, h*4/(13*5*2), _rgb(255,255,255)
  31.  
  32. for i=1 to 4
  33. for j=1 to 5
  34.         starf (j*2)*w*2/(5*12), (i*2)*h*7/130, h*4/(13*5*2), _rgb(255,255,255)
  37.  
  39. _putimage (sw/2 - w/2, sh/2 - h/2), img
  41.  
  42. x0 = sw/2 - w/2
  43. y0 = sh/2 - h/2 + sh
  44.  
  47.  
  51.         t = t + 0.2
  52.  
  53.         line (0,sh)-step(sw, sh),_rgb(0,0,0),bf
  54.  
  55.         for y=0 to h + a*0.707 step 1
  56.         for x=0 to w + a*0.707 step 1
  57.                 z = (0.1 + 0.4*(x/w))*a*sin(x/35 - y/70 - t) + 0.5*a
  58.                 dz = 50*a*cos(x/35 - y/70 - t)/35
  59.  
  60.                 xx = x + z*0.707 - a*0.707
  61.                 yy = y - z*0.707
  62.  
  63.                 if (int(xx) >=0 and int(xx) < w - 1 and int(yy) >= 0 and int(yy) < h - 1) then
  64.                         tl = point(int(xx), int(yy))
  65.                         tr = point(int(xx) + 1, int(yy))
  66.                         bl = point(int(xx), int(yy) + 1)
  67.                         br = point(int(xx) + 1, int(yy) + 1)
  68.  
  69.                         dx = xx - int(xx)
  70.                         dy = yy - int(yy)
  71.  
  72.                         r =_round((1 - dy)*((1 - dx)*  _red(tl) + dx*  _red(tr)) + dy*((1 - dx)*  _red(bl) + dx*  _red(br)))
  73.                         g = _round((1 - dy)*((1 - dx)*_green(tl) + dx*_green(tr)) + dy*((1 - dx)*_green(bl) + dx*_green(br)))
  74.                         b = _round((1 - dy)*((1 - dx)* _blue(tl) + dx* _blue(tr)) + dy*((1 - dx)* _blue(bl) + dx* _blue(br)))
  75.  
  76.                         r = r + dz
  77.                         g = g + dz
  78.                         b = b + dz
  79.  
  80.                         if r<0 then r = 0
  81.                         if r>255 then r = 255
  82.                         if g<0 then g = 0
  83.                         if g>255 then g = 255
  84.                         if b<0 then b = 0
  85.                         if b>255 then b = 255
  86.  
  87.                         pset (x0 + x, y0 - a*0.707 + y), _rgb(r,g,b)
  88.                 end if
  89.         next
  90.         next
  91.  
  92.         _display
  93.         _limit 50
  95.  
  98.  
  99. sub starf(x, y, r, c)
  100.         pset (x + r*cos(pi/2), y - r*sin(pi/2)),c
  101.         for i = 0 to 5
  102.                 xx = r*cos(i*4*pi/5 + pi/2)
  103.                 yy = r*sin(i*4*pi/5 + pi/2)
  104.                 line -(x + xx, y - yy),c
  105.         next
  106.         paint (x, y),c
  107.         for i = 0 to 5
  108.                 xx = r*cos(i*4*pi/5 + pi/2)/2
  109.                 yy = r*sin(i*4*pi/5 + pi/2)/2
  110.                 paint (x + xx, y - yy),c
  111.         next
  113.  
qb64flag.png
* qb64flag.png (Filesize: 70.23 KB, Dimensions: 640x960, Views: 245)
Logged

Offline Ashish

  • Forum Resident
  • Posts: 630
  • Never Give Up!
    • View Profile
Re: U.S. Flag
« Reply #12 on: May 15, 2020, 02:24:48 am »
@vince Brilliant!
I think it could also be done if one :-
1.  Create the US Flag image
2. Divide the image into many smaller rectangles and store its coordinates
3. Use _maptriangle 3d command and oscillate the z value with sin() or cos().
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

Offline STxAxTIC

  • Library Staff
  • Forum Resident
  • Posts: 1091
  • he lives
    • View Profile
Re: U.S. Flag
« Reply #13 on: May 15, 2020, 04:18:55 am »
Nice work boys.

Here's one that obeys the wave equation and plots in 3D:

Code: QB64: [Select]
  1. '#lang "qb" '...for freebasic compiler. Keep disabled for QB64.
  2.  
  3. ' *** Video settings. ***
  4.  
  5. screenwidth = 640
  6. screenheight = 480
  8.  
  9. ' *** Performance settings. ***
  10.  
  11. bignumber = 500000 ' Maximum objects per array (determined by memory).
  12. globaldelayinit = 1000 ' Loop damping factor (fine adjustment below).
  14.  
  15. ' *** Initialize counters and array sizes. ***
  16.  
  17. numparticleorig = bignumber
  18. numparticlevisible = bignumber
  19.  
  20. ' *** Define basis arrays and structures. ***
  21.  
  22. ' Screen vectors in three-space.
  23. ' These vectors define the camera angle.
  24. DIM uhat(1 TO 3), vhat(1 TO 3), nhat(1 TO 3)
  25.  
  26. ' Basis vectors defined in three-space.
  27. DIM xhat(1 TO 4), yhat(1 TO 4), zhat(1 TO 4)
  28. xhat(1) = 1: xhat(2) = 0: xhat(3) = 0: xhat(4) = 4
  29. yhat(1) = 0: yhat(2) = 1: yhat(3) = 0: yhat(4) = 2
  30. zhat(1) = 0: zhat(2) = 0: zhat(3) = 1: zhat(4) = 1
  31.  
  32. ' Basis vectors projected into uv two-space.
  33. DIM xhatp(1 TO 2), yhatp(1 TO 2), zhatp(1 TO 2)
  34. DIM xhatp.old(1 TO 2), yhatp.old(1 TO 2), zhatp.old(1 TO 2)
  35.  
  36. ' *** Define particle arrays and structures. ***
  37.  
  38. ' Particle vectors defined in three-space.
  39. DIM vec(numparticleorig, 4)
  40. DIM vecdotnhat(numparticleorig)
  41. DIM vecdotnhatunit.old(numparticleorig)
  42. DIM vecvisible(numparticlevisible, 4)
  43. DIM vecvisibledotnhat(numparticlevisible)
  44. DIM vecvisibledotnhatunit.old(numparticlevisible)
  45.  
  46. ' Particle vectors projected onto infinite uv two-space.
  47. DIM vecpuv(numparticleorig, 1 TO 2)
  48. DIM vecpuv.old(numparticleorig, 1 TO 2)
  49. DIM vecvisiblepuv(numparticlevisible, 1 TO 2)
  50. DIM vecvisiblepuv.old(numparticlevisible, 1 TO 2)
  51.  
  52. ' Particle projections adjusted for screen uv two-space.
  53. DIM vecpuvs(numparticleorig, 1 TO 2)
  54. DIM vecpuvs.old(numparticleorig, 1 TO 2)
  55. DIM vecvisiblepuvs(numparticlevisible, 1 TO 2)
  56. DIM vecvisiblepuvs.old(numparticlevisible, 1 TO 2)
  57.  
  58. ' *** Define specialized arrays and structures. ***
  59.  
  60. ' Arrays for tech 2 mesh plotting.
  61. DIM plotflag(numparticleorig)
  62. DIM plotflag.old(numparticleorig)
  63.  
  64. ' *** Set mathematical constants. ***
  65.  
  66. pi = 3.1415926536
  67. ee = 2.7182818285
  68.  
  69. mainstart:
  70.  
  71. ' *** Initialize user input variables. ***
  72. key$ = ""
  73. 'mousekey$ = ""
  74.  
  75. ' *** Process first user input. ***
  76.  
  77. a$ = "13"
  78.  
  80.     CASE "13": genscheme$ = "wave2d": plotmode$ = "simplemesh"
  82.  
  83. substart:
  84.  
  85. ' *** Zero counters and array sizes. ***
  86.  
  87. numparticleorig = 0
  88. numparticlevisible = 0
  89.  
  90. ' *** Constants and switch control. ***
  91.  
  92. ' Perspective and animation switches/defaults.
  93. centerx = screenwidth / 2
  94. centery = screenheight / 2
  95. speedconst = 50
  96. falsedepth = .01
  97. zoom = 30
  98. timevar = 0
  99. T = 0
  100. togglehud = 1
  101. toggleatomnumbers = -1
  102. toggletimeanimate = -1
  103. toggletimealert = 0
  104. camx = 0
  105. camy = 0
  106. camz = 0
  107. uhat(1) = -3: uhat(2) = 5: uhat(3) = 1 / 4
  108. vhat(1) = -1: vhat(2) = -1: vhat(3) = 8
  109.  
  111.  
  112. SELECT CASE LCASE$(genscheme$)
  113.     CASE "wave2d":
  114.         globaldelay = globaldelayinit * 1500
  115.         uhat(1) = .7802773: uhat(2) = -.4759201: uhat(3) = .4058135
  116.         vhat(1) = .2502121: vhat(2) = .8321912: vhat(3) = .4948249
  117.         togglehud = 1
  118.         toggletimealert = 1
  119.         GOSUB genscheme.wave2d.init
  120.         numparticleorig = pcountparticleorig
  121.         REDIM vec2dz(xrange, yrange)
  122.         REDIM vec2dztemp(xrange, yrange)
  123.         REDIM vec2dzprev(xrange, yrange)
  124.         GOSUB genschemeUSAcolors
  125.         GOSUB genscheme.wave2d.gridinit
  127.  
  128. ' Move objects to accomodate initial camera position.
  129. IF camx <> 0 AND camy <> 0 AND camz <> 0 THEN
  130.     FOR i = 1 TO numparticleorig
  131.         vec(i, 1) = vec(i, 1) + camx
  132.         vec(i, 2) = vec(i, 2) + camy
  133.         vec(i, 3) = vec(i, 3) + camz
  134.     NEXT
  136.  
  137. GOSUB redraw
  138.  
  139. ' *** Begin main loop. ***
  141.     IF toggletimeanimate = 1 THEN
  142.         GOSUB timeanimate
  143.         flagredraw = 1
  144.     END IF
  145.     IF flagredraw = 1 THEN
  146.         GOSUB redraw
  147.         flagredraw = -1
  148.     END IF
  149.     GOSUB keyprocess
  150.     IF toggletimeanimate = 1 THEN
  151.         FOR delaycount = 1 TO globaldelay: NEXT
  152.     END IF
  154. ' *** End main loop. ***
  155.  
  156. ' *** Begin function definitions. ***
  157.  
  158. ' Comment out the conents of this gosub for non-QB64 compiler.
  159. mouseprocess:
  161.  
  162. keyprocess:
  164. IF key$ <> "" THEN
  165.     flagredraw = 1
  168.     CASE "8":
  169.         GOSUB rotate.vhat.plus
  170.     CASE "2":
  171.         GOSUB rotate.vhat.minus
  172.     CASE "4":
  173.         GOSUB rotate.uhat.minus
  174.     CASE "6":
  175.         GOSUB rotate.uhat.plus
  176.     CASE "7":
  177.         GOSUB rotate.clockwise
  178.     CASE "9":
  179.         GOSUB rotate.counterclockwise
  180.     CASE "1":
  181.         GOSUB rotate.uhat.minus: GOSUB normalize.screen.vectors: GOSUB rotate.clockwise
  182.     CASE "3":
  183.         GOSUB rotate.uhat.plus: GOSUB normalize.screen.vectors: GOSUB rotate.counterclockwise
  184.     CASE "q"
  185.         GOSUB strafe.objects.vhat.plus
  186.         GOSUB strafe.camera.vhat.plus
  187.     CASE "e"
  188.         GOSUB strafe.objects.vhat.minus
  189.         GOSUB strafe.camera.vhat.minus
  190.     CASE "x"
  191.         uhat(1) = 0: uhat(2) = 1: uhat(3) = 0
  192.         vhat(1) = 0: vhat(2) = 0: vhat(3) = 1
  193.     CASE "y"
  194.         uhat(1) = 0: uhat(2) = 0: uhat(3) = 1
  195.         vhat(1) = 1: vhat(2) = 0: vhat(3) = 0
  196.     CASE "z"
  197.         uhat(1) = 1: uhat(2) = 0: uhat(3) = 0
  198.         vhat(1) = 0: vhat(2) = 1: vhat(3) = 0
  199.     CASE "t"
  200.         toggletimeanimate = -toggletimeanimate
  201.     CASE "r"
  202.         timevar = 0
  203.     CASE "f"
  204.         GOTO substart
  205.     CASE "g"
  206.         GOTO mainstart
  207.     CASE "-"
  208.         globaldelay = globaldelay * 1.1
  209.     CASE "="
  210.         globaldelay = globaldelay * 0.9
  211.     CASE "`"
  212.         globaldelayinit = globaldelay
  213.     CASE " "
  214.         togglehud = -togglehud
  215.         CLS
  216.     CASE CHR$(27)
  217.         END
  220.  
  221. convert:
  222. ' Convert graphics from uv-cartesian coordinates to monitor coordinates.
  223. x0 = x: y0 = y
  224. x = x0 + centerx
  225. y = -y0 + centery
  227.  
  228. ' *** Define functions for view translation and rotation. ***
  229.  
  230. rotate.uhat.plus:
  231. uhat(1) = nhat(1) + speedconst * uhat(1)
  232. uhat(2) = nhat(2) + speedconst * uhat(2)
  233. uhat(3) = nhat(3) + speedconst * uhat(3)
  235.  
  236. rotate.uhat.minus:
  237. uhat(1) = -nhat(1) + speedconst * uhat(1)
  238. uhat(2) = -nhat(2) + speedconst * uhat(2)
  239. uhat(3) = -nhat(3) + speedconst * uhat(3)
  241.  
  242. rotate.vhat.plus:
  243. vhat(1) = nhat(1) + speedconst * vhat(1)
  244. vhat(2) = nhat(2) + speedconst * vhat(2)
  245. vhat(3) = nhat(3) + speedconst * vhat(3)
  247.  
  248. rotate.vhat.minus:
  249. vhat(1) = -nhat(1) + speedconst * vhat(1)
  250. vhat(2) = -nhat(2) + speedconst * vhat(2)
  251. vhat(3) = -nhat(3) + speedconst * vhat(3)
  253.  
  254. rotate.counterclockwise:
  255. v1 = vhat(1)
  256. v2 = vhat(2)
  257. v3 = vhat(3)
  258. vhat(1) = uhat(1) + speedconst * vhat(1)
  259. vhat(2) = uhat(2) + speedconst * vhat(2)
  260. vhat(3) = uhat(3) + speedconst * vhat(3)
  261. uhat(1) = -v1 + speedconst * uhat(1)
  262. uhat(2) = -v2 + speedconst * uhat(2)
  263. uhat(3) = -v3 + speedconst * uhat(3)
  265.  
  266. rotate.clockwise:
  267. v1 = vhat(1)
  268. v2 = vhat(2)
  269. v3 = vhat(3)
  270. vhat(1) = -uhat(1) + speedconst * vhat(1)
  271. vhat(2) = -uhat(2) + speedconst * vhat(2)
  272. vhat(3) = -uhat(3) + speedconst * vhat(3)
  273. uhat(1) = v1 + speedconst * uhat(1)
  274. uhat(2) = v2 + speedconst * uhat(2)
  275. uhat(3) = v3 + speedconst * uhat(3)
  278.  
  279. strafe.objects.uhat.plus:
  280. FOR i = 1 TO numparticleorig
  281.     vec(i, 1) = vec(i, 1) + uhat(1) * 1 / zoom
  282.     vec(i, 2) = vec(i, 2) + uhat(2) * 1 / zoom
  283.     vec(i, 3) = vec(i, 3) + uhat(3) * 1 / zoom
  286.  
  287. strafe.objects.uhat.minus:
  288. FOR i = 1 TO numparticleorig
  289.     vec(i, 1) = vec(i, 1) - uhat(1) * 1 / zoom
  290.     vec(i, 2) = vec(i, 2) - uhat(2) * 1 / zoom
  291.     vec(i, 3) = vec(i, 3) - uhat(3) * 1 / zoom
  294.  
  295. strafe.objects.vhat.plus:
  296. FOR i = 1 TO numparticleorig
  297.     vec(i, 1) = vec(i, 1) + vhat(1) * 1 / zoom
  298.     vec(i, 2) = vec(i, 2) + vhat(2) * 1 / zoom
  299.     vec(i, 3) = vec(i, 3) + vhat(3) * 1 / zoom
  302.  
  303. strafe.objects.vhat.minus:
  304. FOR i = 1 TO numparticleorig
  305.     vec(i, 1) = vec(i, 1) - vhat(1) * 1 / zoom
  306.     vec(i, 2) = vec(i, 2) - vhat(2) * 1 / zoom
  307.     vec(i, 3) = vec(i, 3) - vhat(3) * 1 / zoom
  310.  
  311. strafe.objects.nhat.plus:
  312. FOR i = 1 TO numparticleorig
  313.     vec(i, 1) = vec(i, 1) + nhat(1) * 1 / zoom
  314.     vec(i, 2) = vec(i, 2) + nhat(2) * 1 / zoom
  315.     vec(i, 3) = vec(i, 3) + nhat(3) * 1 / zoom
  318.  
  319. strafe.objects.nhat.minus:
  320. FOR i = 1 TO numparticleorig
  321.     vec(i, 1) = vec(i, 1) - nhat(1) * 1 / zoom
  322.     vec(i, 2) = vec(i, 2) - nhat(2) * 1 / zoom
  323.     vec(i, 3) = vec(i, 3) - nhat(3) * 1 / zoom
  326.  
  327. strafe.camera.uhat.plus:
  328. camx = camx + uhat(1) * 1 / zoom
  329. camy = camy + uhat(2) * 1 / zoom
  330. camz = camz + uhat(3) * 1 / zoom
  332.  
  333. strafe.camera.uhat.minus:
  334. camx = camx - uhat(1) * 1 / zoom
  335. camy = camy - uhat(2) * 1 / zoom
  336. camz = camz - uhat(3) * 1 / zoom
  338.  
  339. strafe.camera.vhat.plus:
  340. camx = camx + vhat(1) * 1 / zoom
  341. camy = camy + vhat(2) * 1 / zoom
  342. camz = camz + vhat(3) * 1 / zoom
  344.  
  345. strafe.camera.vhat.minus:
  346. camx = camx - vhat(1) * 1 / zoom
  347. camy = camy - vhat(2) * 1 / zoom
  348. camz = camz - vhat(3) * 1 / zoom
  350.  
  351. strafe.camera.nhat.plus:
  352. camx = camx + nhat(1) * 1 / zoom
  353. camy = camy + nhat(2) * 1 / zoom
  354. camz = camz + nhat(3) * 1 / zoom
  356.  
  357. strafe.camera.nhat.minus:
  358. camx = camx - nhat(1) * 1 / zoom
  359. camy = camy - nhat(2) * 1 / zoom
  360. camz = camz - nhat(3) * 1 / zoom
  362.  
  363. ' *** Define core functions. ***
  364.  
  365. timeanimate:
  366. timevar = timevar + 1
  367. IF timevar > 10 ^ 6 THEN timevar = 1
  368. SELECT CASE genscheme$
  369.     CASE "wave2d": GOSUB genscheme.wave2d.timeanimate
  372.  
  373. normalize.screen.vectors:
  374. 'normalize the two vectors that define the screen orientation
  375. uhatmag = SQR(uhat(1) ^ 2 + uhat(2) ^ 2 + uhat(3) ^ 2)
  376. uhat(1) = uhat(1) / uhatmag: uhat(2) = uhat(2) / uhatmag: uhat(3) = uhat(3) / uhatmag
  377. vhatmag = SQR(vhat(1) ^ 2 + vhat(2) ^ 2 + vhat(3) ^ 2)
  378. vhat(1) = vhat(1) / vhatmag: vhat(2) = vhat(2) / vhatmag: vhat(3) = vhat(3) / vhatmag
  379. uhatdotvhat = uhat(1) * vhat(1) + uhat(2) * vhat(2) + uhat(3) * vhat(3)
  380. IF SQR(uhatdotvhat ^ 2) > .0005 THEN
  381.     CLS: COLOR 15: LOCATE 5, 5: PRINT "Screen vectors are not perpendicular. Press ESC to quit."
  382.     'DO: LOOP UNTIL INKEY$ = CHR$(27): END
  384. ' Compute the normal vector to the view plane.
  385. ' The normal vector points toward the eye, away from view frustum.
  386. nhat(1) = uhat(2) * vhat(3) - uhat(3) * vhat(2)
  387. nhat(2) = uhat(3) * vhat(1) - uhat(1) * vhat(3)
  388. nhat(3) = uhat(1) * vhat(2) - uhat(2) * vhat(1)
  389. nhatmag = SQR(nhat(1) ^ 2 + nhat(2) ^ 2 + nhat(3) ^ 2)
  390. nhat(1) = nhat(1) / nhatmag: nhat(2) = nhat(2) / nhatmag: nhat(3) = nhat(3) / nhatmag
  392.  
  393. redraw:
  394. GOSUB normalize.screen.vectors
  395. ' Project the three-space basis vectors onto the screen plane.
  396. xhatp(1) = xhat(1) * uhat(1) + xhat(2) * uhat(2) + xhat(3) * uhat(3)
  397. xhatp(2) = xhat(1) * vhat(1) + xhat(2) * vhat(2) + xhat(3) * vhat(3)
  398. yhatp(1) = yhat(1) * uhat(1) + yhat(2) * uhat(2) + yhat(3) * uhat(3)
  399. yhatp(2) = yhat(1) * vhat(1) + yhat(2) * vhat(2) + yhat(3) * vhat(3)
  400. zhatp(1) = zhat(1) * uhat(1) + zhat(2) * uhat(2) + zhat(3) * uhat(3)
  401. zhatp(2) = zhat(1) * vhat(1) + zhat(2) * vhat(2) + zhat(3) * vhat(3)
  402. IF numparticleorig > 0 THEN
  403.     GOSUB compute.visible.particles
  404.     GOSUB project.particles
  405.     GOSUB depth.adjust.particles
  407. GOSUB draw.all.objects
  408. GOSUB store.screen.projections
  410.  
  411. reverse.uvnhat:
  412. uhat(1) = -uhat(1)
  413. uhat(2) = -uhat(2)
  414. uhat(3) = -uhat(3)
  415. GOSUB normalize.screen.vectors
  417.  
  418. compute.visible.particles:
  419. numparticlevisible = 0
  420. FOR i = 1 TO numparticleorig
  421.     IF falsedepth = 0 THEN
  422.     ELSE
  423.         numparticlevisible = numparticlevisible + 1
  424.         vecvisible(numparticlevisible, 1) = vec(i, 1)
  425.         vecvisible(numparticlevisible, 2) = vec(i, 2)
  426.         vecvisible(numparticlevisible, 3) = vec(i, 3)
  427.         vecvisible(numparticlevisible, 4) = vec(i, 4)
  428.     END IF
  431.  
  432. project.particles:
  433. ' Project object vectors onto the screen plane.
  434. FOR i = 1 TO numparticlevisible
  435.     vecvisibledotnhat(i) = vecvisible(i, 1) * nhat(1) + vecvisible(i, 2) * nhat(2) + vecvisible(i, 3) * nhat(3)
  436.     vecvisiblepuv(i, 1) = (vecvisible(i, 1) * uhat(1) + vecvisible(i, 2) * uhat(2) + vecvisible(i, 3) * uhat(3))
  437.     vecvisiblepuv(i, 2) = (vecvisible(i, 1) * vhat(1) + vecvisible(i, 2) * vhat(2) + vecvisible(i, 3) * vhat(3))
  440.  
  441. depth.adjust.particles:
  442. IF falsedepth = 0 THEN
  443.     FOR i = 1 TO numparticlevisible
  444.         vecvisiblepuvs(i, 1) = vecvisiblepuv(i, 1) * fovd / vecvisibledotnhat(i)
  445.         vecvisiblepuvs(i, 2) = vecvisiblepuv(i, 2) * fovd / vecvisibledotnhat(i)
  446.     NEXT
  448.     FOR i = 1 TO numparticlevisible
  449.         vecvisiblepuvs(i, 1) = vecvisiblepuv(i, 1) * (1 + falsedepth * vecvisibledotnhat(i))
  450.         vecvisiblepuvs(i, 2) = vecvisiblepuv(i, 2) * (1 + falsedepth * vecvisibledotnhat(i))
  451.     NEXT
  454.  
  455. draw.all.objects:
  456. SELECT CASE plotmode$
  457.     CASE "simplemesh": GOSUB plotmode.simplemesh
  460. LOCATE 28, 23: PRINT "SPACE = toggle HUD,  ESC = quit."
  461. IF togglehud = 1 THEN
  462.     ' Replace basis vector triad.
  463.     x = 50 * xhatp.old(1): y = 50 * xhatp.old(2): GOSUB convert
  464.     LINE (centerx, centery)-(x, y), 0
  465.     x = 50 * yhatp.old(1): y = 50 * yhatp.old(2): GOSUB convert
  466.     LINE (centerx, centery)-(x, y), 0
  467.     x = 50 * zhatp.old(1): y = 50 * zhatp.old(2): GOSUB convert
  468.     LINE (centerx, centery)-(x, y), 0
  469.     x = 50 * xhatp(1): y = 50 * xhatp(2): GOSUB convert
  470.     LINE (centerx, centery)-(x, y), xhat(4)
  471.     x = 50 * yhatp(1): y = 50 * yhatp(2): GOSUB convert
  472.     LINE (centerx, centery)-(x, y), yhat(4)
  473.     x = 50 * zhatp(1): y = 50 * zhatp(2): GOSUB convert
  474.     LINE (centerx, centery)-(x, y), zhat(4)
  475.     COLOR 14
  476.     LOCATE 25, 68: PRINT "-   VIEW   -"
  477.     COLOR 15
  478.     LOCATE 26, 68: PRINT "  8  "
  479.     LOCATE 27, 68: PRINT "4   6"
  480.     LOCATE 28, 68: PRINT "  2  "
  481.     COLOR 7
  482.     LOCATE 26, 75: PRINT "7   9"
  483.     LOCATE 27, 75: PRINT "     "
  484.     LOCATE 28, 75: PRINT "1   3"
  486. IF toggletimealert = 1 THEN
  487.     COLOR 7
  488.     LOCATE 1, 25: PRINT "Press 'T' to toggle animation."
  489.     LOCATE 2, 25: PRINT "Use numeric keypad to rotate."
  490.     LOCATE 3, 25: PRINT "      Press G to reset."
  493.  
  494. store.screen.projections:
  495. xhatp.old(1) = xhatp(1): xhatp.old(2) = xhatp(2)
  496. yhatp.old(1) = yhatp(1): yhatp.old(2) = yhatp(2)
  497. zhatp.old(1) = zhatp(1): zhatp.old(2) = zhatp(2)
  498. FOR i = 1 TO numparticlevisible
  499.     vecvisiblepuvs.old(i, 1) = vecvisiblepuvs(i, 1)
  500.     vecvisiblepuvs.old(i, 2) = vecvisiblepuvs(i, 2)
  502. numparticlevisible.old = numparticlevisible
  504.  
  505. plotmode.simplemesh:
  506. FOR i = 1 TO numparticlevisible - 1
  507.     IF i MOD yrange <> 0 THEN 'point obeys normal neighbor-connect scheme
  508.         ' Erase old graphics.
  509.         x = zoom * vecvisiblepuvs.old(i, 1): y = zoom * vecvisiblepuvs.old(i, 2): GOSUB convert: x1 = x: y1 = y
  510.         x = zoom * vecvisiblepuvs.old(i + 1, 1): y = zoom * vecvisiblepuvs.old(i + 1, 2): GOSUB convert: x2 = x: y2 = y
  511.         x = zoom * vecvisiblepuvs.old(i + yrange, 1): y = zoom * vecvisiblepuvs.old(i + yrange, 2): GOSUB convert: x3 = x: y3 = y
  512.         LINE (x1, y1)-(x2, y2), 0
  513.         IF i < (numparticlevisible - yrange) THEN LINE (x1, y1)-(x3, y3), 0
  514.         ' Draw new graphics.
  515.         x = zoom * vecvisiblepuvs(i, 1): y = zoom * vecvisiblepuvs(i, 2): GOSUB convert: x1 = x: y1 = y
  516.         x = zoom * vecvisiblepuvs(i + 1, 1): y = zoom * vecvisiblepuvs(i + 1, 2): GOSUB convert: x2 = x: y2 = y
  517.         x = zoom * vecvisiblepuvs(i + yrange, 1): y = zoom * vecvisiblepuvs(i + yrange, 2): GOSUB convert: x3 = x: y3 = y
  518.         LINE (x1, y1)-(x2, y2), vecvisible(i, 4)
  519.         IF i < (numparticlevisible - yrange) THEN LINE (x1, y1)-(x3, y3), vecvisible(i, 4)
  520.     ELSE 'point does not obey normal neighbor-connect scheme
  521.         ' Erase old graphics.
  522.         x = zoom * vecvisiblepuvs.old(i, 1): y = zoom * vecvisiblepuvs.old(i, 2): GOSUB convert: x1 = x: y1 = y
  523.         x = zoom * vecvisiblepuvs.old(i + yrange, 1): y = zoom * vecvisiblepuvs.old(i + yrange, 2): GOSUB convert: x3 = x: y3 = y
  524.         IF i < (numparticlevisible - yrange + 1) THEN LINE (x1, y1)-(x3, y3), 0
  525.         ' Draw new graphics.
  526.         x = zoom * vecvisiblepuvs(i, 1): y = zoom * vecvisiblepuvs(i, 2): GOSUB convert: x1 = x: y1 = y
  527.         x = zoom * vecvisiblepuvs(i + yrange, 1): y = zoom * vecvisiblepuvs(i + yrange, 2): GOSUB convert: x3 = x: y3 = y
  528.         IF i < (numparticlevisible - yrange + 1) THEN LINE (x1, y1)-(x3, y3), vecvisible(i, 4)
  529.     END IF
  532.  
  533. genschemeUSAcolors:
  534. 'delinearize
  535. pcountparticleorig = 0
  536. FOR i = 1 TO xrange
  537.     FOR j = 1 TO yrange
  538.         pcountparticleorig = pcountparticleorig + 1
  539.         vec2dztemp(i, j) = vec(pcountparticleorig, 4)
  540.     NEXT
  542. FOR i = 1 TO xrange
  543.     FOR jj = 0 TO 12 STEP 2
  544.         FOR j = jj * yrange / 13 + 1 TO (jj + 1) * yrange / 13
  545.             vec2dztemp(i, j) = 4
  546.         NEXT
  547.     NEXT
  548.     FOR jj = 1 TO 11 STEP 2
  549.         FOR j = jj * yrange / 13 + 1 TO (jj + 1) * yrange / 13
  550.             vec2dztemp(i, j) = 7
  551.         NEXT
  552.     NEXT
  554. starflag = -1
  555. FOR i = 1 TO xrange * .76 / 1.9
  556.     FOR j = yrange TO yrange * .5385 STEP -1
  557.         IF starflag = 1 THEN
  558.             vec2dztemp(i, j) = 15
  559.         ELSE
  560.             vec2dztemp(i, j) = 1
  561.         END IF
  562.         starflag = -starflag
  563.     NEXT
  565. 'relinearize
  566. pcountparticleorig = 0
  567. FOR i = 1 TO xrange
  568.     FOR j = 1 TO yrange
  569.         pcountparticleorig = pcountparticleorig + 1
  570.         vec(pcountparticleorig, 4) = vec2dztemp(i, j)
  571.     NEXT
  574.  
  575. genscheme.wave2d.init:
  576. xl = -1.9: xr = 1.9
  577. yl = -1: yr = 1
  578. xl = xl * 4: xr = xr * 4: yl = yl * 4: yr = yr * 4
  579. Dx = .32
  580. Dy = .32
  581. xrange = 1 + INT((-xl + xr) / Dx)
  582. yrange = 1 + INT((-yl + yr) / Dy)
  583. alpha = .25
  584. pcountparticleorig = 0
  585. FOR i = xl TO xr STEP Dx
  586.     FOR j = yl TO yr STEP Dy
  587.         pcountparticleorig = pcountparticleorig + 1
  588.         vec(pcountparticleorig, 1) = i
  589.         vec(pcountparticleorig, 2) = j
  590.         vec(pcountparticleorig, 3) = 0
  591.         '*' vec(pcountparticleorig, 4) = 14 'use special color scheme
  592.     NEXT
  595.  
  596. genscheme.wave2d.gridinit:
  597. 'delinearize
  598. pcountparticleorig = 0
  599. FOR i = 1 TO xrange
  600.     FOR j = 1 TO yrange
  601.         pcountparticleorig = pcountparticleorig + 1
  602.         vec2dztemp(i, j) = vec(pcountparticleorig, 3)
  603.     NEXT
  605. 'initial position condition
  606. 'FOR i = 1 TO xrange 'random high points
  607. '  FOR j = 1 TO yrange
  608. '      nrand = RND * 1000
  609. '      IF nrand < 10 THEN
  610. '          vec2dz(i, j) = 5
  611. '      END IF
  612. '  NEXT
  613. 'NEXT
  614. 'vec2dz(xrange * .8, yrange * .2) = -2.5 'single plucked point
  615. 'FOR i = 1 TO xrange 'cross arm
  616. '   vec2dz(i, yrange / 3) = 2
  617. 'NEXT
  618. 'FOR j = 1 TO yrange 'cross arm
  619. '   vec2dz(xrange / 2, j) = 1
  620. 'NEXT
  621. 'sync
  622. pcountparticleorig = 0
  623. FOR i = 1 TO xrange
  624.     FOR j = 1 TO yrange
  625.         pcountparticleorig = pcountparticleorig + 1
  626.         vec2dzprev(i, j) = vec2dz(i, j)
  627.         vec2dztemp(i, j) = vec2dz(i, j)
  628.     NEXT
  630. 'initial velocity condition
  631. vec2dzprev(xrange * .8, yrange * .8) = 1.5 'single struck point
  632. 'relinearize
  633. pcountparticleorig = 0
  634. FOR i = 1 TO xrange
  635.     FOR j = 1 TO yrange
  636.         pcountparticleorig = pcountparticleorig + 1
  637.         vec(pcountparticleorig, 3) = vec2dz(i, j)
  638.     NEXT
  641.  
  642. genscheme.wave2d.timeanimate:
  643. 'delinearize
  644. pcountparticleorig = 0
  645. FOR i = 1 TO xrange
  646.     FOR j = 1 TO yrange
  647.         pcountparticleorig = pcountparticleorig + 1
  648.         vec2dztemp(i, j) = vec(pcountparticleorig, 3)
  649.     NEXT
  651. 'begin propagation process
  652. FOR i = 2 TO xrange - 1 'boDy, no edges
  653.     FOR j = 2 TO yrange - 1
  654.         wp1 = alpha * (vec2dztemp(i + 1, j) + vec2dztemp(i - 1, j)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  655.         wp2 = alpha * (vec2dztemp(i, j + 1) + vec2dztemp(i, j - 1)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  656.         vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  657.     NEXT
  659. 'comment out this section for fixed edges (or pieces of this section)
  660. i = 1 'left edge
  661. FOR j = 2 TO yrange - 1
  662.     wfp = vec2dztemp(i, j)
  663.     wp1 = alpha * (vec2dztemp(i + 1, j) + wfp) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  664.     wp2 = alpha * (vec2dztemp(i, j + 1) + vec2dztemp(i, j - 1)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  665.     'vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  667. i = xrange 'right edge
  668. FOR j = 2 TO yrange - 1
  669.     wfp = vec2dztemp(i, j)
  670.     wp1 = alpha * (wfp + vec2dztemp(i - 1, j)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  671.     wp2 = alpha * (vec2dztemp(i, j + 1) + vec2dztemp(i, j - 1)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  672.     vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  674. j = 1 'bottom edge
  675. FOR i = 2 TO xrange - 1
  676.     wfp = vec2dztemp(i, j)
  677.     wp2 = alpha * (vec2dztemp(i, j + 1) + wfp) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  678.     wp1 = alpha * (vec2dztemp(i + 1, j) + vec2dztemp(i - 1, j)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  679.     vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  681. j = yrange 'top edge
  682. FOR i = 2 TO xrange - 1
  683.     wfp = vec2dztemp(i, j)
  684.     wp2 = alpha * (wfp + vec2dztemp(i, j - 1)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  685.     wp1 = alpha * (vec2dztemp(i + 1, j) + vec2dztemp(i - 1, j)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  686.     vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  688. 'bottom left corner
  689. i = 1: j = 1
  690. wfp1 = vec2dztemp(i, j)
  691. wfp2 = vec2dztemp(i, j)
  692. wp1 = alpha * (vec2dztemp(i + 1, j) + wfp1) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  693. wp2 = alpha * (vec2dztemp(i, j + 1) + wfp2) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  694. 'vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  695. 'bottom right corner
  696. i = xrange: j = 1
  697. wfp1 = vec2dztemp(i, j)
  698. wfp2 = vec2dztemp(i, j)
  699. wp1 = alpha * (wfp1 + vec2dztemp(i - 1, j)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  700. wp2 = alpha * (vec2dztemp(i, j + 1) + wfp2) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  701. vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  702. 'top left corner
  703. i = 1: j = yrange
  704. wfp1 = vec2dztemp(i, j)
  705. wfp2 = vec2dztemp(i, j)
  706. wp1 = alpha * (vec2dztemp(i + 1, j) + wfp1) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  707. wp2 = alpha * (wfp2 + vec2dztemp(i, j - 1)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  708. 'vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  709. 'top right corner
  710. i = xrange: j = yrange
  711. wfp1 = vec2dztemp(i, j)
  712. wfp2 = vec2dztemp(i, j)
  713. wp1 = alpha * (wfp1 + vec2dztemp(i - 1, j)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  714. wp2 = alpha * (wfp2 + vec2dztemp(i, j - 1)) + 2 * (1 - alpha) * vec2dztemp(i, j) - vec2dzprev(i, j)
  715. vec2dz(i, j) = (1 / 2) * (wp1 + wp2)
  716. 'start special movements
  717. T = timevar / 5
  718. IF T < pi THEN 'wave left edge just once
  719.     FOR j = 1 TO yrange
  720.         i = 1
  721.         vec2dz(i, j) = (j / yrange) * 1.5 * SIN(T)
  722.     NEXT
  724. IF T < pi THEN 'wave bottom edge just once
  725.     FOR i = 1 TO xrange
  726.         j = 1
  727.         vec2dz(i, j) = (i / xrange) * .45 * SIN(T)
  728.     NEXT
  730. 'relinearize
  731. pcountparticleorig = 0
  732. FOR i = 1 TO xrange
  733.     FOR j = 1 TO yrange
  734.         pcountparticleorig = pcountparticleorig + 1
  735.         vec2dzprev(i, j) = vec2dztemp(i, j)
  736.         vec2dztemp(i, j) = vec2dz(i, j)
  737.         vec(pcountparticleorig, 3) = vec2dz(i, j)
  738.     NEXT
  741.  
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: U.S. Flag
« Reply #14 on: May 15, 2020, 11:08:28 am »
@_vince

Very nice coloring with wave, beautiful!
Logged
Pages: [1] 2   Go Up
« previous next »