« previous next »
Pages: [1]   Go Down

Author Topic: Sanctum 3D World by STxAxTIC  (Read 9882 times)

0 Members and 1 Guest are viewing this topic.

Offline The Librarian

  • Moderator
  • Newbie
  • Posts: 39
    • View Profile
Sanctum 3D World by STxAxTIC
« on: August 13, 2019, 09:53:46 pm »
Sanctum 3D World

Author: STxAxTIC
Source: qb64.org Forum
URL: https://www.qb64.org/forum/index.php?topic=1600.0
Version: QB64
Tags: [3d] [graphics] [open world]

Description:
Sanctum is an open world 3D engine that makes everything from CIRCLE and LINE. Suggested left hand on "WSAD", and right hand on numeric keypad (hybrid gamer setup)... Discover heaven, discover hell, find the pyramid, find the megalith, blow up the moon, play with blocks.

Source Code:
Code: QB64: [Select]
  1. '$EXEICON:'sanctum.ico'
  2. 'REM $include:'Color32.BI'
  3.  
  4. CONST Aquamarine = _RGB32(127, 255, 212)
  5. CONST Blue = _RGB32(0, 0, 255)
  6. CONST BlueViolet = _RGB32(138, 43, 226)
  7. CONST Chocolate = _RGB32(210, 105, 30)
  8. CONST Cyan = _RGB32(0, 255, 255)
  9. CONST DarkBlue = _RGB32(0, 0, 139)
  10. CONST DarkGoldenRod = _RGB32(184, 134, 11)
  11. CONST DarkGray = _RGB32(169, 169, 169)
  12. CONST DarkKhaki = _RGB32(189, 183, 107)
  13. CONST DeepPink = _RGB32(255, 20, 147)
  14. CONST DodgerBlue = _RGB32(30, 144, 255)
  15. CONST ForestGreen = _RGB32(34, 139, 34)
  16. CONST Gray = _RGB32(128, 128, 128)
  17. CONST Green = _RGB32(0, 128, 0)
  18. CONST Indigo = _RGB32(75, 0, 130)
  19. CONST Ivory = _RGB32(255, 255, 240)
  20. CONST LightSeaGreen = _RGB32(32, 178, 170)
  21. CONST Lime = _RGB32(0, 255, 0)
  22. CONST LimeGreen = _RGB32(50, 205, 50)
  23. CONST Magenta = _RGB32(255, 0, 255)
  24. CONST PaleGoldenRod = _RGB32(238, 232, 170)
  25. CONST Purple = _RGB32(128, 0, 128)
  26. CONST Red = _RGB32(255, 0, 0)
  27. CONST RoyalBlue = _RGB32(65, 105, 225)
  28. CONST SaddleBrown = _RGB32(139, 69, 19)
  29. CONST Sienna = _RGB32(160, 82, 45)
  30. CONST SlateGray = _RGB32(112, 128, 144)
  31. CONST Snow = _RGB32(255, 250, 250)
  32. CONST Sunglow = _RGB32(255, 207, 72)
  33. CONST SunsetOrange = _RGB32(253, 94, 83)
  34. CONST Teal = _RGB32(0, 128, 128)
  35. CONST White = _RGB32(255, 255, 255)
  36. CONST Yellow = _RGB32(255, 255, 0)
  37.  
  38. ' Constants.
  39. pi = 3.1415926536
  40. ee = 2.7182818285
  41.  
  42. ' Scale.
  43. DIM bignumber AS LONG
  44. bignumber = 3000000
  45.  
  46. ' Video.
  47. 'SCREEN _NEWIMAGE(640, 480, 32)
  48. SCREEN _NEWIMAGE(800, 600, 32)
  49. 'SCREEN _NEWIMAGE(1024, 768, 32)
  52. screenwidth = _WIDTH
  53. screenheight = _HEIGHT
  54.  
  55. ' Camera orientation vectors.
  56. DIM uhat(3), vhat(3), nhat(3)
  57.  
  58. ' Basis vectors defined in three-space.
  59. DIM xhat(3), yhat(3), zhat(3)
  60. xhat(1) = 1: xhat(2) = 0: xhat(3) = 0
  61. yhat(1) = 0: yhat(2) = 1: yhat(3) = 0
  62. zhat(1) = 0: zhat(2) = 0: zhat(3) = 1
  63.  
  64. ' Group structure.
  65. TYPE VectorGroupElement
  66.     Identity AS LONG
  67.     Pointer AS LONG
  68.     Lagger AS LONG
  69.     FirstVector AS LONG
  70.     LastVector AS LONG
  71.     GroupName AS STRING * 50
  72.     Visible AS INTEGER
  73.     ForceAnimate AS INTEGER
  74.     COMFixed AS INTEGER
  75.     COMx AS SINGLE ' Center of mass
  76.     COMy AS SINGLE
  77.     COMz AS SINGLE
  78.     ROTx AS SINGLE ' Center of rotation
  79.     ROTy AS SINGLE
  80.     ROTz AS SINGLE
  81.     REVx AS SINGLE ' Revolution speed
  82.     REVy AS SINGLE
  83.     REVz AS SINGLE
  84.     DIMx AS SINGLE ' Maximum volume
  85.     DIMy AS SINGLE
  86.     DIMz AS SINGLE
  88. DIM VectorGroup(bignumber) AS VectorGroupElement
  89.  
  90. ' World vectors.
  91. DIM vec(bignumber, 3) ' Relative Position
  92. DIM vec3Dpos(bignumber, 3) ' Position
  93. DIM vec3Dvel(bignumber, 3) ' Linear velocity
  94. DIM vec3Dacc(bignumber, 3) ' Linear acceleration
  95. DIM vec3Danv(bignumber, 3) ' Angular velocity
  96. DIM vec3Dvis(bignumber) ' Visible toggle
  97. DIM vec2D(bignumber, 2) ' Projection onto 2D plane
  98. DIM vec3Dcolor(bignumber) AS LONG ' Original color
  99. DIM vec2Dcolor(bignumber) AS LONG ' Projected color
  100.  
  101. ' Clipping planes.
  102. DIM nearplane(4), farplane(4), rightplane(4), leftplane(4), topplane(4), bottomplane(4)
  103.  
  104. ' State.
  106. nearplane(4) = 1
  107. farplane(4) = -100
  108. rightplane(4) = 0 '*' fovd * (nhat(1) * rightplane(1) + nhat(2) * rightplane(2) + nhat(3) * rightplane(3))
  109. leftplane(4) = 0
  110. topplane(4) = 0
  111. bottomplane(4) = 0
  112. midscreenx = screenwidth / 2
  113. midscreeny = screenheight / 2
  114. fovd = -256
  115. numgroupvisible = 0
  116. numvectorvisible = 0
  117. groupidticker = 0
  118. vecgroupid = 0
  119. vectorindex = 0
  120. rotspeed = 1 / 33
  121. linspeed = 3 / 2
  122. timestep = .001
  123. camx = -40
  124. camy = 30
  125. camz = 40
  126. uhat(1) = -.2078192: uhat(2) = -.9781672: uhat(3) = 0
  127. vhat(1) = 0: vhat(2) = 0: vhat(3) = 1
  128. toggletimeanimate = 1
  129. toggleinvertmouse = -1
  130. togglehud = 1
  131.  
  132. ' Prime main loop.
  133. GOSUB initialize.objects
  134. GOSUB redraw
  135.  
  136. ' Begin main loop.
  137. fpstimer = INT(TIMER)
  138. fps = 0
  140.     GOSUB redraw
  141.     GOSUB mouseprocess
  142.     GOSUB keyprocess
  143.  
  144.     fps = fps + 1
  145.     tt = INT(TIMER)
  146.     IF tt = fpstimer + 1 THEN
  147.         fpstimer = tt
  148.         fpsreport = fps
  149.         fps = 0
  150.     END IF
  151.  
  152.     _DISPLAY
  153.     _KEYCLEAR
  154.     _LIMIT 30
  156.  
  158.  
  159. ' Gosubs.
  160.  
  161. redraw:
  162. GOSUB normalize.screen.vectors
  163. GOSUB calculate.clippingplanes
  164. GOSUB compute.visible.groups
  165. GOSUB plot.visible.vectors
  167.  
  168. mouseprocess:
  169. 'mx = 0
  170. 'my = 0
  171. 'DO WHILE _MOUSEINPUT
  172. '    mx = mx + _MOUSEMOVEMENTX
  173. '    my = my + _MOUSEMOVEMENTY
  174. '    IF _MOUSEWHEEL > 0 THEN GOSUB rotate.clockwise
  175. '    IF _MOUSEWHEEL < 0 THEN GOSUB rotate.counterclockwise
  176. '    IF mx > 0 THEN
  177. '        GOSUB rotate.uhat.plus: GOSUB normalize.screen.vectors
  178. '    END IF
  179. '    IF mx < 0 THEN
  180. '        GOSUB rotate.uhat.minus: GOSUB normalize.screen.vectors
  181. '    END IF
  182. '    IF my > 0 THEN
  183. '        IF toggleinvertmouse = -1 THEN
  184. '            GOSUB rotate.vhat.plus: GOSUB normalize.screen.vectors
  185. '        ELSE
  186. '            GOSUB rotate.vhat.minus: GOSUB normalize.screen.vectors
  187. '        END IF
  188. '    END IF
  189. '    IF my < 0 THEN
  190. '        IF toggleinvertmouse = -1 THEN
  191. '            GOSUB rotate.vhat.minus: GOSUB normalize.screen.vectors
  192. '        ELSE
  193. '            GOSUB rotate.vhat.plus: GOSUB normalize.screen.vectors
  194. '        END IF
  195. '    END IF
  196. '    mx = 0
  197. '    my = 0
  198. 'LOOP
  200.  
  201. keyprocess:
  202. IF _KEYDOWN(119) = -1 OR _KEYDOWN(18432) = -1 THEN GOSUB strafe.camera.nhat.minus ' w or uparrow
  203. IF _KEYDOWN(115) = -1 OR _KEYDOWN(20480) = -1 THEN GOSUB strafe.camera.nhat.plus ' s or downarrow
  204. IF _KEYDOWN(97) = -1 THEN GOSUB strafe.camera.uhat.minus ' a
  205. IF _KEYDOWN(100) = -1 THEN GOSUB strafe.camera.uhat.plus ' d
  206. IF _KEYDOWN(56) = -1 THEN GOSUB rotate.vhat.plus: GOSUB normalize.screen.vectors ' 8
  207. IF _KEYDOWN(50) = -1 THEN GOSUB rotate.vhat.minus: GOSUB normalize.screen.vectors ' 2
  208. IF _KEYDOWN(19200) = -1 OR _KEYDOWN(52) = -1 THEN GOSUB rotate.uhat.minus: GOSUB normalize.screen.vectors ' 4
  209. IF _KEYDOWN(19712) = -1 OR _KEYDOWN(54) = -1 THEN GOSUB rotate.uhat.plus: GOSUB normalize.screen.vectors ' 6
  210. IF _KEYDOWN(55) = -1 THEN GOSUB rotate.clockwise ' 7
  211. IF _KEYDOWN(57) = -1 THEN GOSUB rotate.counterclockwise ' 9
  212. IF _KEYDOWN(49) = -1 THEN GOSUB rotate.uhat.minus: GOSUB normalize.screen.vectors: GOSUB rotate.clockwise ' 1
  213. IF _KEYDOWN(51) = -1 THEN GOSUB rotate.uhat.plus: GOSUB normalize.screen.vectors: GOSUB rotate.counterclockwise ' 3
  214. IF _KEYDOWN(113) = -1 THEN GOSUB strafe.camera.vhat.minus ' q
  215. IF _KEYDOWN(101) = -1 THEN GOSUB strafe.camera.vhat.plus ' e
  216.  
  218. IF (key$ <> "") THEN
  219.     SELECT CASE key$
  220.         CASE "x"
  221.             uhat(1) = 0: uhat(2) = 1: uhat(3) = 0
  222.             vhat(1) = 0: vhat(2) = 0: vhat(3) = 1
  223.         CASE "X"
  224.             uhat(1) = 0: uhat(2) = -1: uhat(3) = 0
  225.             vhat(1) = 0: vhat(2) = 0: vhat(3) = 1
  226.         CASE "y"
  227.             uhat(1) = -1: uhat(2) = 0: uhat(3) = 0
  228.             vhat(1) = 0: vhat(2) = 0: vhat(3) = 1
  229.         CASE "Y"
  230.             uhat(1) = 1: uhat(2) = 0: uhat(3) = 0
  231.             vhat(1) = 0: vhat(2) = 0: vhat(3) = 1
  232.         CASE "z"
  233.             uhat(1) = 1: uhat(2) = 0: uhat(3) = 0
  234.             vhat(1) = 0: vhat(2) = 1: vhat(3) = 0
  235.             GOSUB normalize.screen.vectors
  236.         CASE "Z"
  237.             uhat(1) = 0: uhat(2) = 1: uhat(3) = 0
  238.             vhat(1) = 1: vhat(2) = 0: vhat(3) = 0
  239.         CASE "]"
  240.             farplane(4) = farplane(4) - 1
  241.         CASE "["
  242.             farplane(4) = farplane(4) + 1
  243.         CASE " "
  244.             togglehud = -togglehud
  245.         CASE "t"
  246.             toggletimeanimate = -toggletimeanimate
  247.         CASE "i"
  248.             toggleinvertmouse = -toggleinvertmouse
  249.         CASE "v"
  250.             OPEN "snapshot-camera.txt" FOR OUTPUT AS #1
  251.             PRINT #1, camx, camy, camz
  252.             PRINT #1, uhat(1), uhat(2), uhat(3)
  253.             PRINT #1, vhat(1), vhat(2), vhat(3)
  254.             CLOSE #1
  255.         CASE CHR$(27)
  256.             SYSTEM
  257.         CASE "n"
  258.             VectorGroup(closestgroup).COMFixed = 1
  259.             FOR vectorindex = VectorGroup(closestgroup).FirstVector TO VectorGroup(closestgroup).LastVector
  260.                 vec3Dvel(vectorindex, 1) = (RND - .5) * 200
  261.                 vec3Dvel(vectorindex, 2) = (RND - .5) * 200
  262.                 vec3Dvel(vectorindex, 3) = (RND - .5) * 200
  263.             NEXT
  264.         CASE "k"
  265.             p = VectorGroup(closestgroup).Pointer
  266.             l = VectorGroup(closestgroup).Lagger
  267.             VectorGroup(l).Pointer = p
  268.             IF (p <> -999) THEN
  269.                 VectorGroup(p).Lagger = l
  270.             END IF
  271.         CASE "b"
  272.             tilesize = 5
  273.             ' Determine last object id.
  274.             p = 1
  275.             DO
  276.                 k = VectorGroup(p).Identity
  277.                 p = VectorGroup(k).Pointer
  278.                 IF (p = -999) THEN EXIT DO
  279.             LOOP
  280.             lastobjectid = k
  281.             vectorindex = VectorGroup(lastobjectid).LastVector
  282.             ' Create new group.
  283.             groupidticker = groupidticker + 1
  284.             vecgroupid = groupidticker
  285.             VectorGroup(vecgroupid).Identity = vecgroupid
  286.             VectorGroup(vecgroupid).Pointer = -999
  287.             VectorGroup(vecgroupid).Lagger = lastobjectid
  288.             VectorGroup(vecgroupid).GroupName = "Block"
  289.             VectorGroup(vecgroupid).Visible = 0
  290.             VectorGroup(vecgroupid).COMFixed = 1
  291.             VectorGroup(vecgroupid).DIMx = tilesize / 2
  292.             VectorGroup(vecgroupid).DIMy = tilesize / 2
  293.             VectorGroup(vecgroupid).DIMz = tilesize / 2
  294.             VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  295.             FOR r = 1 TO 400
  296.                 vectorindex = vectorindex + 1
  297.                 vec3Dpos(vectorindex, 1) = camx + -20 * nhat(1) + (RND - .5) * tilesize
  298.                 vec3Dpos(vectorindex, 2) = camy + -20 * nhat(2) + (RND - .5) * tilesize
  299.                 vec3Dpos(vectorindex, 3) = camz + -20 * nhat(3) + (RND - .5) * tilesize
  300.                 vec3Dvis(vectorindex) = 0
  301.                 IF RND > .5 THEN
  302.                     vec3Dcolor(vectorindex) = Lime
  303.                 ELSE
  304.                     vec3Dcolor(vectorindex) = Purple
  305.                 END IF
  306.                 GOSUB integratecom
  307.             NEXT
  308.             VectorGroup(vecgroupid).LastVector = vectorindex
  309.             GOSUB calculatecom
  310.             VectorGroup(lastobjectid).Pointer = vecgroupid
  311.     END SELECT
  314.  
  315. convert:
  316. ' Convert graphics from uv-cartesian coordinates to monitor coordinates.
  317. x0 = x: y0 = y
  318. x = x0 + midscreenx
  319. y = -y0 + midscreeny
  321.  
  322. rotate.uhat.plus:
  323. uhat(1) = uhat(1) + nhat(1) * rotspeed
  324. uhat(2) = uhat(2) + nhat(2) * rotspeed
  325. uhat(3) = uhat(3) + nhat(3) * rotspeed
  327.  
  328. rotate.uhat.minus:
  329. uhat(1) = uhat(1) - nhat(1) * rotspeed
  330. uhat(2) = uhat(2) - nhat(2) * rotspeed
  331. uhat(3) = uhat(3) - nhat(3) * rotspeed
  333.  
  334. rotate.vhat.plus:
  335. vhat(1) = vhat(1) + nhat(1) * rotspeed
  336. vhat(2) = vhat(2) + nhat(2) * rotspeed
  337. vhat(3) = vhat(3) + nhat(3) * rotspeed
  339.  
  340. rotate.vhat.minus:
  341. vhat(1) = vhat(1) - nhat(1) * rotspeed
  342. vhat(2) = vhat(2) - nhat(2) * rotspeed
  343. vhat(3) = vhat(3) - nhat(3) * rotspeed
  345.  
  346. rotate.counterclockwise:
  347. v1 = vhat(1)
  348. v2 = vhat(2)
  349. v3 = vhat(3)
  350. vhat(1) = vhat(1) + uhat(1) * rotspeed
  351. vhat(2) = vhat(2) + uhat(2) * rotspeed
  352. vhat(3) = vhat(3) + uhat(3) * rotspeed
  353. uhat(1) = uhat(1) - v1 * rotspeed
  354. uhat(2) = uhat(2) - v2 * rotspeed
  355. uhat(3) = uhat(3) - v3 * rotspeed
  357.  
  358. rotate.clockwise:
  359. v1 = vhat(1)
  360. v2 = vhat(2)
  361. v3 = vhat(3)
  362. vhat(1) = vhat(1) - uhat(1) * rotspeed
  363. vhat(2) = vhat(2) - uhat(2) * rotspeed
  364. vhat(3) = vhat(3) - uhat(3) * rotspeed
  365. uhat(1) = uhat(1) + v1 * rotspeed
  366. uhat(2) = uhat(2) + v2 * rotspeed
  367. uhat(3) = uhat(3) + v3 * rotspeed
  369.  
  370. strafe.camera.uhat.plus:
  371. camx = camx + uhat(1) * linspeed
  372. camy = camy + uhat(2) * linspeed
  373. camz = camz + uhat(3) * linspeed
  375.  
  376. strafe.camera.uhat.minus:
  377. camx = camx - uhat(1) * linspeed
  378. camy = camy - uhat(2) * linspeed
  379. camz = camz - uhat(3) * linspeed
  381.  
  382. strafe.camera.vhat.plus:
  383. camx = camx + vhat(1) * linspeed
  384. camy = camy + vhat(2) * linspeed
  385. camz = camz + vhat(3) * linspeed
  387.  
  388. strafe.camera.vhat.minus:
  389. camx = camx - vhat(1) * linspeed
  390. camy = camy - vhat(2) * linspeed
  391. camz = camz - vhat(3) * linspeed
  393.  
  394. strafe.camera.nhat.plus:
  395. camx = camx + nhat(1) * linspeed
  396. camy = camy + nhat(2) * linspeed
  397. camz = camz + nhat(3) * linspeed
  399.  
  400. strafe.camera.nhat.minus:
  401. camx = camx - nhat(1) * linspeed
  402. camy = camy - nhat(2) * linspeed
  403. camz = camz - nhat(3) * linspeed
  405.  
  406. normalize.screen.vectors:
  407. uhatmag = SQR(uhat(1) * uhat(1) + uhat(2) * uhat(2) + uhat(3) * uhat(3))
  408. uhat(1) = uhat(1) / uhatmag: uhat(2) = uhat(2) / uhatmag: uhat(3) = uhat(3) / uhatmag
  409. vhatmag = SQR(vhat(1) * vhat(1) + vhat(2) * vhat(2) + vhat(3) * vhat(3))
  410. vhat(1) = vhat(1) / vhatmag: vhat(2) = vhat(2) / vhatmag: vhat(3) = vhat(3) / vhatmag
  411. uhatdotvhat = uhat(1) * vhat(1) + uhat(2) * vhat(2) + uhat(3) * vhat(3)
  412. ' The normal vector points toward the eye.
  413. nhat(1) = uhat(2) * vhat(3) - uhat(3) * vhat(2)
  414. nhat(2) = uhat(3) * vhat(1) - uhat(1) * vhat(3)
  415. nhat(3) = uhat(1) * vhat(2) - uhat(2) * vhat(1)
  416. nhatmag = SQR(nhat(1) * nhat(1) + nhat(2) * nhat(2) + nhat(3) * nhat(3))
  417. nhat(1) = nhat(1) / nhatmag: nhat(2) = nhat(2) / nhatmag: nhat(3) = nhat(3) / nhatmag
  419.  
  420. calculate.clippingplanes:
  421. ' Calculate normal vectors to all clipping planes.
  422. h2 = screenheight / 2
  423. w2 = screenwidth / 2
  424. nearplane(1) = -nhat(1)
  425. nearplane(2) = -nhat(2)
  426. nearplane(3) = -nhat(3)
  427. farplane(1) = nhat(1)
  428. farplane(2) = nhat(2)
  429. farplane(3) = nhat(3)
  430. rightplane(1) = h2 * fovd * uhat(1) - h2 * w2 * nhat(1)
  431. rightplane(2) = h2 * fovd * uhat(2) - h2 * w2 * nhat(2)
  432. rightplane(3) = h2 * fovd * uhat(3) - h2 * w2 * nhat(3)
  433. mag = SQR(rightplane(1) * rightplane(1) + rightplane(2) * rightplane(2) + rightplane(3) * rightplane(3))
  434. rightplane(1) = rightplane(1) / mag
  435. rightplane(2) = rightplane(2) / mag
  436. rightplane(3) = rightplane(3) / mag
  437. leftplane(1) = -h2 * fovd * uhat(1) - h2 * w2 * nhat(1)
  438. leftplane(2) = -h2 * fovd * uhat(2) - h2 * w2 * nhat(2)
  439. leftplane(3) = -h2 * fovd * uhat(3) - h2 * w2 * nhat(3)
  440. mag = SQR(leftplane(1) * leftplane(1) + leftplane(2) * leftplane(2) + leftplane(3) * leftplane(3))
  441. leftplane(1) = leftplane(1) / mag
  442. leftplane(2) = leftplane(2) / mag
  443. leftplane(3) = leftplane(3) / mag
  444. topplane(1) = w2 * fovd * vhat(1) - h2 * w2 * nhat(1)
  445. topplane(2) = w2 * fovd * vhat(2) - h2 * w2 * nhat(2)
  446. topplane(3) = w2 * fovd * vhat(3) - h2 * w2 * nhat(3)
  447. mag = SQR(topplane(1) * topplane(1) + topplane(2) * topplane(2) + topplane(3) * topplane(3))
  448. topplane(1) = topplane(1) / mag
  449. topplane(2) = topplane(2) / mag
  450. topplane(3) = topplane(3) / mag
  451. bottomplane(1) = -w2 * fovd * vhat(1) - h2 * w2 * nhat(1)
  452. bottomplane(2) = -w2 * fovd * vhat(2) - h2 * w2 * nhat(2)
  453. bottomplane(3) = -w2 * fovd * vhat(3) - h2 * w2 * nhat(3)
  454. mag = SQR(bottomplane(1) * bottomplane(1) + bottomplane(2) * bottomplane(2) + bottomplane(3) * bottomplane(3))
  455. bottomplane(1) = bottomplane(1) / mag
  456. bottomplane(2) = bottomplane(2) / mag
  457. bottomplane(3) = bottomplane(3) / mag
  459.  
  460. compute.visible.groups:
  461. closestdist2 = 10000000
  462. closestgroup = 1
  463. fp42 = farplane(4) * farplane(4)
  464.  
  465. k = 1
  466. k = VectorGroup(k).Identity
  467. DO ' iterates k
  468.  
  469.     VectorGroup(k).Visible = 0
  470.  
  471.     dx = VectorGroup(k).COMx - camx
  472.     dy = VectorGroup(k).COMy - camy
  473.     dz = VectorGroup(k).COMz - camz
  474.  
  475.     dist2 = dx * dx + dy * dy + dz * dz
  476.  
  477.     IF dist2 < fp42 THEN
  478.  
  479.         groupinview = 1
  480.         IF dx * nearplane(1) + dy * nearplane(2) + dz * nearplane(3) - nearplane(4) < 0 THEN groupinview = 0
  481.         'IF dx * farplane(1) + dy * farplane(2) + dz * farplane(3) - farplane(4) < 0 THEN groupinview = 0
  482.         IF dx * rightplane(1) + dy * rightplane(2) + dz * rightplane(3) - rightplane(4) < 0 THEN groupinview = 0
  483.         IF dx * leftplane(1) + dy * leftplane(2) + dz * leftplane(3) - leftplane(4) < 0 THEN groupinview = 0
  484.         IF dx * topplane(1) + dy * topplane(2) + dz * topplane(3) - topplane(4) < 0 THEN groupinview = 0
  485.         IF dx * bottomplane(1) + dy * bottomplane(2) + dz * bottomplane(3) - bottomplane(4) < 0 THEN groupinview = 0
  486.         IF groupinview = 1 THEN
  487.  
  488.             IF (dist2 < closestdist2) THEN
  489.                 closestdist2 = dist2
  490.                 closestgroup = k
  491.             END IF
  492.  
  493.             VectorGroup(k).Visible = 1
  494.  
  495.             IF (toggletimeanimate = 1) THEN
  496.                 vecgroupid = k
  497.                 GOSUB timeanimate
  498.             END IF
  499.  
  500.             FOR i = VectorGroup(k).FirstVector TO VectorGroup(k).LastVector
  501.                 GOSUB clip.project.vectors
  502.             NEXT
  503.  
  504.         ELSE
  505.             ' Force animation regardless of clipping.
  506.             IF (VectorGroup(k).ForceAnimate = 1) THEN
  507.                 vecgroupid = k
  508.                 GOSUB timeanimate
  509.             END IF
  510.         END IF
  511.     ELSE
  512.         ' Force animation regardless of distance from camera.
  513.         IF (VectorGroup(k).ForceAnimate = 1) THEN
  514.             vecgroupid = k
  515.             GOSUB timeanimate
  516.         END IF
  517.     END IF
  518.     k = VectorGroup(k).Pointer
  519.     IF k = -999 THEN EXIT DO
  520.     k = VectorGroup(k).Identity
  523.  
  524. clip.project.vectors: ' requires i
  525. vec(i, 1) = vec3Dpos(i, 1) - camx
  526. vec(i, 2) = vec3Dpos(i, 2) - camy
  527. vec(i, 3) = vec3Dpos(i, 3) - camz
  528. fogswitch = -1
  529. vec3Dvis(i) = 0
  530. vectorinview = 1
  531. ' Perform view plane clipping.
  532. IF vec(i, 1) * nearplane(1) + vec(i, 2) * nearplane(2) + vec(i, 3) * nearplane(3) - nearplane(4) < 0 THEN vectorinview = 0
  533. IF vec(i, 1) * farplane(1) + vec(i, 2) * farplane(2) + vec(i, 3) * farplane(3) - farplane(4) < 0 THEN vectorinview = 0
  534. IF vec(i, 1) * farplane(1) + vec(i, 2) * farplane(2) + vec(i, 3) * farplane(3) - farplane(4) * .85 < 0 THEN fogswitch = 1
  535. IF vec(i, 1) * rightplane(1) + vec(i, 2) * rightplane(2) + vec(i, 3) * rightplane(3) - rightplane(4) < 0 THEN vectorinview = 0
  536. IF vec(i, 1) * leftplane(1) + vec(i, 2) * leftplane(2) + vec(i, 3) * leftplane(3) - leftplane(4) < 0 THEN vectorinview = 0
  537. IF vec(i, 1) * topplane(1) + vec(i, 2) * topplane(2) + vec(i, 3) * topplane(3) - topplane(4) < 0 THEN vectorinview = 0
  538. IF vec(i, 1) * bottomplane(1) + vec(i, 2) * bottomplane(2) + vec(i, 3) * bottomplane(3) - bottomplane(4) < 0 THEN vectorinview = 0
  539. IF vectorinview = 1 THEN
  540.     vec3Dvis(i) = 1
  541.     ' Project vectors onto the screen plane.
  542.     vec3Ddotnhat = vec(i, 1) * nhat(1) + vec(i, 2) * nhat(2) + vec(i, 3) * nhat(3)
  543.     vec2D(i, 1) = (vec(i, 1) * uhat(1) + vec(i, 2) * uhat(2) + vec(i, 3) * uhat(3)) * fovd / vec3Ddotnhat
  544.     vec2D(i, 2) = (vec(i, 1) * vhat(1) + vec(i, 2) * vhat(2) + vec(i, 3) * vhat(3)) * fovd / vec3Ddotnhat
  545.     IF fogswitch = 1 THEN vec2Dcolor(i) = Gray ELSE vec2Dcolor(i) = vec3Dcolor(i)
  548.  
  549. timeanimate: ' requires vecgroupid
  550. dt = timestep
  551.  
  552. xcom = VectorGroup(vecgroupid).COMx
  553. ycom = VectorGroup(vecgroupid).COMy
  554. zcom = VectorGroup(vecgroupid).COMz
  555. xrot = VectorGroup(vecgroupid).ROTx
  556. yrot = VectorGroup(vecgroupid).ROTy
  557. zrot = VectorGroup(vecgroupid).ROTz
  558. xrev = VectorGroup(vecgroupid).ROTx
  559. yrev = VectorGroup(vecgroupid).ROTy
  560. zrev = VectorGroup(vecgroupid).ROTz
  561. xdim = VectorGroup(vecgroupid).DIMx
  562. ydim = VectorGroup(vecgroupid).DIMy
  563. zdim = VectorGroup(vecgroupid).DIMz
  564.  
  565. IF (VectorGroup(vecgroupid).COMFixed = 0) THEN GOSUB resetcom
  566.  
  567. FOR vectorindex = VectorGroup(vecgroupid).FirstVector TO VectorGroup(vecgroupid).LastVector
  568.  
  569.     ' Linear velocity update
  570.     ax = vec3Dacc(vectorindex, 1)
  571.     ay = vec3Dacc(vectorindex, 2)
  572.     az = vec3Dacc(vectorindex, 3)
  573.     IF (ax <> 0) THEN vec3Dvel(vectorindex, 1) = vec3Dvel(vectorindex, 1) + ax * dt
  574.     IF (ay <> 0) THEN vec3Dvel(vectorindex, 2) = vec3Dvel(vectorindex, 2) + ay * dt
  575.     IF (az <> 0) THEN vec3Dvel(vectorindex, 3) = vec3Dvel(vectorindex, 3) + az * dt
  576.  
  577.     ' Linear position update with periodic boundaries inside group dimension
  578.     vx = vec3Dvel(vectorindex, 1)
  579.     vy = vec3Dvel(vectorindex, 2)
  580.     vz = vec3Dvel(vectorindex, 3)
  581.     IF (vx <> 0) THEN
  582.         px = vec3Dpos(vectorindex, 1) + vx * dt
  583.         IF ABS(px - xcom) > xdim THEN
  584.             IF (px > xcom) THEN
  585.                 px = xcom - xdim
  586.             ELSE
  587.                 px = xcom + xdim
  588.             END IF
  589.         END IF
  590.         vec3Dpos(vectorindex, 1) = px
  591.     END IF
  592.     IF (vy <> 0) THEN
  593.         py = vec3Dpos(vectorindex, 2) + vy * dt
  594.         IF ABS(py - ycom) > ydim THEN
  595.             IF (py > ycom) THEN
  596.                 py = ycom - ydim
  597.             ELSE
  598.                 py = ycom + ydim
  599.             END IF
  600.         END IF
  601.         vec3Dpos(vectorindex, 2) = py
  602.     END IF
  603.     IF (vz <> 0) THEN
  604.         pz = vec3Dpos(vectorindex, 3) + vz * dt
  605.         IF ABS(pz - zcom) > zdim THEN
  606.             IF (pz > zcom) THEN
  607.                 pz = zcom - zdim
  608.             ELSE
  609.                 pz = zcom + zdim
  610.             END IF
  611.         END IF
  612.         vec3Dpos(vectorindex, 3) = pz
  613.     END IF
  614.  
  615.     ' Rotation update
  616.     IF (xrot <> 0) THEN
  617.         anv = vec3Danv(vectorindex, 1)
  618.         yy = vec3Dpos(vectorindex, 2) - yrot
  619.         zz = vec3Dpos(vectorindex, 3) - zrot
  620.         y = yy * COS(timestep * anv) - zz * SIN(timestep * anv)
  621.         z = yy * SIN(timestep * anv) + zz * COS(timestep * anv)
  622.         vec3Dpos(vectorindex, 2) = y + yrot
  623.         vec3Dpos(vectorindex, 3) = z + zrot
  624.     END IF
  625.     IF (yrot <> 0) THEN
  626.         anv = vec3Danv(vectorindex, 2)
  627.         xx = vec3Dpos(vectorindex, 1) - xrot
  628.         zz = vec3Dpos(vectorindex, 3) - zrot
  629.         x = xx * COS(timestep * anv) + zz * SIN(timestep * anv)
  630.         z = -xx * SIN(timestep * anv) + zz * COS(timestep * anv)
  631.         vec3Dpos(vectorindex, 1) = x + xrot
  632.         vec3Dpos(vectorindex, 3) = z + zrot
  633.     END IF
  634.     IF (zrot <> 0) THEN
  635.         anv = vec3Danv(vectorindex, 3)
  636.         xx = vec3Dpos(vectorindex, 1) - xrot
  637.         yy = vec3Dpos(vectorindex, 2) - yrot
  638.         x = xx * COS(timestep * anv) - yy * SIN(timestep * anv)
  639.         y = xx * SIN(timestep * anv) + yy * COS(timestep * anv)
  640.         vec3Dpos(vectorindex, 1) = x + xrot
  641.         vec3Dpos(vectorindex, 2) = y + yrot
  642.     END IF
  643.  
  644.     ' Revolution update
  645.     IF (xrev <> 0) THEN
  646.         anv = xrev
  647.         yy = vec3Dpos(vectorindex, 2) - ycom
  648.         zz = vec3Dpos(vectorindex, 3) - zcom
  649.         y = yy * COS(timestep * anv) - zz * SIN(timestep * anv)
  650.         z = yy * SIN(timestep * anv) + zz * COS(timestep * anv)
  651.         vec3Dpos(vectorindex, 2) = y + ycom
  652.         vec3Dpos(vectorindex, 3) = z + zcom
  653.     END IF
  654.     IF (yrev <> 0) THEN
  655.         anv = yrev
  656.         xx = vec3Dpos(vectorindex, 1) - xcom
  657.         zz = vec3Dpos(vectorindex, 3) - zcom
  658.         x = xx * COS(timestep * anv) + zz * SIN(timestep * anv)
  659.         z = -xx * SIN(timestep * anv) + zz * COS(timestep * anv)
  660.         vec3Dpos(vectorindex, 1) = x + xcom
  661.         vec3Dpos(vectorindex, 3) = z + zcom
  662.     END IF
  663.     IF (zrev <> 0) THEN
  664.         anv = zrev
  665.         xx = vec3Dpos(vectorindex, 1) - xcom
  666.         yy = vec3Dpos(vectorindex, 2) - ycom
  667.         x = xx * COS(timestep * anv) - yy * SIN(timestep * anv)
  668.         y = xx * SIN(timestep * anv) + yy * COS(timestep * anv)
  669.         vec3Dpos(vectorindex, 1) = x + xcom
  670.         vec3Dpos(vectorindex, 2) = y + ycom
  671.     END IF
  672.  
  673.     IF (VectorGroup(vecgroupid).COMFixed = 0) THEN GOSUB integratecom
  675. IF (VectorGroup(vecgroupid).COMFixed = 0) THEN GOSUB calculatecom
  677.  
  678. integratecom: ' requires vecgroupid
  679. VectorGroup(vecgroupid).COMx = vec3Dpos(vectorindex, 1) + VectorGroup(vecgroupid).COMx
  680. VectorGroup(vecgroupid).COMy = vec3Dpos(vectorindex, 2) + VectorGroup(vecgroupid).COMy
  681. VectorGroup(vecgroupid).COMz = vec3Dpos(vectorindex, 3) + VectorGroup(vecgroupid).COMz
  683.  
  684. calculatecom: ' requires vecgroupid
  685. f = 1 + VectorGroup(vecgroupid).LastVector - VectorGroup(vecgroupid).FirstVector
  686. VectorGroup(vecgroupid).COMx = VectorGroup(vecgroupid).COMx / f
  687. VectorGroup(vecgroupid).COMy = VectorGroup(vecgroupid).COMy / f
  688. VectorGroup(vecgroupid).COMz = VectorGroup(vecgroupid).COMz / f
  690.  
  691. resetcom: ' requires vecgroupid
  692. VectorGroup(vecgroupid).COMx = 0
  693. VectorGroup(vecgroupid).COMy = 0
  694. VectorGroup(vecgroupid).COMz = 0
  696.  
  697. plot.visible.vectors:
  698. GOSUB plot.vectors
  699. GOSUB plot.hud
  701.  
  702. plot.vectors:
  704. numgroupvisible = 0
  705. numvectorvisible = 0
  706. k = 1
  707. k = VectorGroup(k).Identity
  709.     IF (VectorGroup(k).Visible = 1) THEN
  710.         numgroupvisible = numgroupvisible + 1
  711.         FOR i = VectorGroup(k).FirstVector TO VectorGroup(k).LastVector - 1
  712.             IF (vec3Dvis(i) = 1) THEN
  713.                 numvectorvisible = numvectorvisible + 1
  714.                 IF k = closestgroup THEN col = Yellow ELSE col = vec2Dcolor(i)
  715.                 x = vec2D(i, 1): y = vec2D(i, 2): GOSUB convert: x1 = x: y1 = y
  716.                 x = vec2D(i + 1, 1): y = vec2D(i + 1, 2): GOSUB convert: x2 = x: y2 = y
  717.                 IF ((x2 - x1) * (x2 - x1) + (y2 - y1) * (y2 - y1)) < 225 THEN
  718.                     LINE (x1, y1)-(x2, y2), col
  719.                 ELSE
  720.                     CIRCLE (x1, y1), 1, col
  721.                 END IF
  722.             END IF
  723.         NEXT
  724.     END IF
  725.     k = VectorGroup(k).Pointer
  726.     IF k = -999 THEN EXIT DO
  727.     k = VectorGroup(k).Identity
  730.  
  731. plot.hud:
  732. ' Redraw compass.
  733. x = 30 * (xhat(1) * uhat(1) + xhat(2) * uhat(2) + xhat(3) * uhat(3)): y = 30 * (xhat(1) * vhat(1) + xhat(2) * vhat(2) + xhat(3) * vhat(3)): GOSUB convert
  734. LINE (midscreenx, midscreeny)-(x, y), Red
  735. x = 30 * (yhat(1) * uhat(1) + yhat(2) * uhat(2) + yhat(3) * uhat(3)): y = 30 * (yhat(1) * vhat(1) + yhat(2) * vhat(2) + yhat(3) * vhat(3)): GOSUB convert
  736. LINE (midscreenx, midscreeny)-(x, y), Green
  737. x = 30 * (zhat(1) * uhat(1) + zhat(2) * uhat(2) + zhat(3) * uhat(3)): y = 30 * (zhat(1) * vhat(1) + zhat(2) * vhat(2) + zhat(3) * vhat(3)): GOSUB convert
  738. LINE (midscreenx, midscreeny)-(x, y), Blue
  739. IF togglehud = 1 THEN
  740.     COLOR LimeGreen
  741.     LOCATE 2, 2: PRINT "- View Info -"
  742.     COLOR DarkKhaki
  743.     LOCATE 3, 2: PRINT "FPS:"; fpsreport
  744.     LOCATE 4, 2: PRINT "Vectors:"; numvectorvisible
  745.     LOCATE 5, 2: PRINT "Groups:"; numgroupvisible
  746.     LOCATE 6, 2: PRINT "Depth:"; -farplane(4)
  747.     LOCATE 7, 2: PRINT "Adjust via [ ]"
  748.     COLOR LimeGreen
  749.     LOCATE 9, 2: PRINT "- Camera -"
  750.     COLOR DarkKhaki
  751.     LOCATE 10, 2: PRINT INT(camx); INT(camy); INT(camz)
  752.     COLOR LimeGreen
  753.     LOCATE 12, 2: PRINT "- Closest: -"
  754.     COLOR DarkKhaki
  755.     LOCATE 13, 2: PRINT LTRIM$(RTRIM$(VectorGroup(closestgroup).GroupName))
  756.     COLOR LimeGreen
  757.     a$ = "MOVE - ALIGN": LOCATE 2, screenwidth / 8 - LEN(a$): PRINT a$
  758.     COLOR DarkKhaki
  759.     a$ = "q w e - x y z": LOCATE 3, screenwidth / 8 - LEN(a$): PRINT a$
  760.     a$ = "a s d - X Y Z": LOCATE 4, screenwidth / 8 - LEN(a$): PRINT a$
  761.     a$ = "i = invert ms": LOCATE 5, screenwidth / 8 - LEN(a$): PRINT a$
  762.     COLOR LimeGreen
  763.     a$ = "- ROTATE -": LOCATE 7, screenwidth / 8 - LEN(a$): PRINT a$
  764.     COLOR DarkKhaki
  765.     a$ = "7 8 9 Mouse": LOCATE 8, screenwidth / 8 - LEN(a$): PRINT a$
  766.     a$ = "4 5 6   +  ": LOCATE 9, screenwidth / 8 - LEN(a$): PRINT a$
  767.     a$ = "1 2 3 Wheel": LOCATE 10, screenwidth / 8 - LEN(a$): PRINT a$
  768.     COLOR LimeGreen
  769.     a$ = "- CONTROL -": LOCATE 12, screenwidth / 8 - LEN(a$): PRINT a$
  770.     COLOR DarkKhaki
  771.     a$ = "t = Stop time": LOCATE 13, screenwidth / 8 - LEN(a$): PRINT a$
  772.     a$ = "b = Create": LOCATE 14, screenwidth / 8 - LEN(a$): PRINT a$
  773.     a$ = "n = Destroy": LOCATE 15, screenwidth / 8 - LEN(a$): PRINT a$
  774.     a$ = "k = Delete": LOCATE 16, screenwidth / 8 - LEN(a$): PRINT a$
  775.     COLOR LimeGreen
  776.     a$ = "SPACE = Hide Info": LOCATE (screenheight / 16) - 3, (screenwidth / 8) / 2 - LEN(a$) / 2: PRINT a$
  778.     COLOR LimeGreen
  779.     a$ = "You See: " + LTRIM$(RTRIM$(VectorGroup(closestgroup).GroupName)): LOCATE (screenheight / 16) - 3, (screenwidth / 8) / 2 - LEN(a$) / 2: PRINT a$
  782.  
  783. 'groupidfromname: ' requires n$, returns k
  784. 'k = 1
  785. 'k = VectorGroup(k).Identity
  786. 'DO ' iterates k
  787. '    IF n$ = LTRIM$(RTRIM$(VectorGroup(k).GroupName)) THEN EXIT DO
  788. '    k = VectorGroup(k).Pointer
  789. '    IF k = -999 THEN EXIT DO
  790. '    k = VectorGroup(k).Identity
  791. 'LOOP
  792. 'RETURN
  793.  
  794. ' Data.
  795.  
  796. initialize.objects:
  797. vectorindex = 0
  798. groupidticker = 0
  799. gridsize = 550
  800. tilesize = 15
  801.  
  802. '__AAA
  803. groupidticker = groupidticker + 1
  804. vecgroupid = groupidticker
  805. VectorGroup(vecgroupid).Identity = vecgroupid
  806. VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  807. VectorGroup(vecgroupid).Lagger = vecgroupid - 1 ' Fancy way to say 0.
  808. VectorGroup(vecgroupid).GroupName = "__AAA"
  809. VectorGroup(vecgroupid).Visible = 0
  810. VectorGroup(vecgroupid).COMFixed = 1
  811. VectorGroup(vecgroupid).DIMx = 5
  812. VectorGroup(vecgroupid).DIMy = 5
  813. VectorGroup(vecgroupid).DIMz = 5
  814. VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  815. FOR r = 1 TO 1
  816.     vectorindex = vectorindex + 1
  817.     vec3Dpos(vectorindex, 1) = 0
  818.     vec3Dpos(vectorindex, 2) = 0
  819.     vec3Dpos(vectorindex, 3) = -1000
  820.     vec3Dcolor(vectorindex) = White
  821.     GOSUB integratecom
  823. VectorGroup(vecgroupid).LastVector = vectorindex
  824. GOSUB calculatecom
  825.  
  826. 'Dirt
  827. h = 5
  828. FOR w = 1 TO 5
  829.     FOR u = -gridsize TO gridsize STEP tilesize
  830.         FOR v = -gridsize TO gridsize STEP tilesize
  831.             groupidticker = groupidticker + 1
  832.             vecgroupid = groupidticker
  833.             VectorGroup(vecgroupid).Identity = vecgroupid
  834.             VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  835.             VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  836.             VectorGroup(vecgroupid).GroupName = "Dirt"
  837.             VectorGroup(vecgroupid).Visible = 0
  838.             VectorGroup(vecgroupid).COMFixed = 1
  839.             VectorGroup(vecgroupid).DIMx = 35
  840.             VectorGroup(vecgroupid).DIMy = 35
  841.             VectorGroup(vecgroupid).DIMz = 35
  842.             VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  843.             FOR i = u TO u + tilesize STEP h
  844.                 FOR j = v TO v + tilesize STEP h
  845.                     IF RND > 1 - w / 5 THEN
  846.                         vectorindex = vectorindex + 1
  847.                         vec3Dpos(vectorindex, 1) = i + RND * h - RND * h
  848.                         vec3Dpos(vectorindex, 2) = j + RND * h - RND * h
  849.                         vec3Dpos(vectorindex, 3) = -(w - 1) * 70 - RND * 70
  850.                         vec3Dvis(vectorindex) = 0
  851.                         IF RND > .5 THEN
  852.                             vec3Dcolor(vectorindex) = DarkGoldenRod
  853.                         ELSE
  854.                             IF RND > .5 THEN
  855.                                 vec3Dcolor(vectorindex) = SaddleBrown
  856.                             ELSE
  857.                                 vec3Dcolor(vectorindex) = Sienna
  858.                             END IF
  859.                         END IF
  860.                         GOSUB integratecom
  861.                     END IF
  862.                 NEXT
  863.             NEXT
  864.             VectorGroup(vecgroupid).LastVector = vectorindex
  865.             GOSUB calculatecom
  866.         NEXT
  867.     NEXT
  869.  
  870. 'Grass and Puddles
  871. h = 2
  872. FOR u = -gridsize TO gridsize STEP tilesize
  873.     FOR v = -gridsize TO gridsize STEP tilesize
  874.         groupidticker = groupidticker + 1
  875.         vecgroupid = groupidticker
  876.         VectorGroup(vecgroupid).Identity = vecgroupid
  877.         VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  878.         VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  879.         VectorGroup(vecgroupid).GroupName = "Grass and Puddles"
  880.         VectorGroup(vecgroupid).Visible = 0
  881.         VectorGroup(vecgroupid).COMFixed = 1
  882.         VectorGroup(vecgroupid).DIMx = tilesize / 2
  883.         VectorGroup(vecgroupid).DIMy = tilesize / 2
  884.         VectorGroup(vecgroupid).DIMz = 3
  885.         VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  886.         FOR i = u TO u + tilesize STEP h
  887.             FOR j = v TO v + tilesize STEP h
  888.                 vectorindex = vectorindex + 1
  889.                 vec3Dpos(vectorindex, 1) = i + RND * h - RND * h
  890.                 vec3Dpos(vectorindex, 2) = j + RND * h - RND * h
  891.                 vec3Dpos(vectorindex, 3) = .5 + 1 * COS((i - 15) * .08) - 1 * COS((j - 6) * .12)
  892.                 vec3Dvis(vectorindex) = 0
  893.                 IF vec3Dpos(vectorindex, 3) > 0 THEN
  894.                     IF RND > .5 THEN
  895.                         vec3Dcolor(vectorindex) = Green
  896.                     ELSE
  897.                         vec3Dcolor(vectorindex) = ForestGreen
  898.                     END IF
  899.                 ELSE
  900.                     vec3Dvel(vectorindex, 1) = (RND - .5) * 20
  901.                     vec3Dvel(vectorindex, 2) = (RND - .5) * 20
  902.                     vec3Dvel(vectorindex, 3) = (RND - .5) * 20
  903.                     IF RND > .2 THEN
  904.                         vec3Dcolor(vectorindex) = LightSeaGreen
  905.                     ELSE
  906.                         vec3Dcolor(vectorindex) = Blue
  907.                     END IF
  908.                 END IF
  909.                 GOSUB integratecom
  910.             NEXT
  911.         NEXT
  912.         VectorGroup(vecgroupid).LastVector = vectorindex
  913.         GOSUB calculatecom
  914.     NEXT
  916.  
  917. 'Grave
  918. thickness = 2.5
  919. span = 20
  920. height = 30
  921. crux = 22
  922. FOR xloc = -90 TO -290 STEP -60
  923.     FOR yloc = 0 TO 180 STEP 45
  924.         FOR k = 0 TO height
  925.             groupidticker = groupidticker + 1
  926.             vecgroupid = groupidticker
  927.             VectorGroup(vecgroupid).Identity = vecgroupid
  928.             VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  929.             VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  930.             VectorGroup(vecgroupid).GroupName = "Grave"
  931.             VectorGroup(vecgroupid).Visible = 0
  932.             VectorGroup(vecgroupid).COMFixed = 1
  933.             VectorGroup(vecgroupid).DIMx = thickness
  934.             VectorGroup(vecgroupid).DIMy = thickness
  935.             VectorGroup(vecgroupid).DIMz = thickness
  936.             VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  937.             FOR i = -thickness TO thickness STEP thickness / 2
  938.                 FOR j = -thickness TO thickness STEP thickness / 2
  939.                     vectorindex = vectorindex + 1
  940.                     vec3Dpos(vectorindex, 1) = xloc + i + (RND - .5) * 2
  941.                     vec3Dpos(vectorindex, 2) = yloc + j + (RND - .5) * 2
  942.                     vec3Dpos(vectorindex, 3) = k + (RND - .5) * 2
  943.                     vec3Dvis(vectorindex) = 0
  944.                     IF RND > .5 THEN
  945.                         vec3Dcolor(vectorindex) = SlateGray
  946.                     ELSE
  947.                         vec3Dcolor(vectorindex) = DarkGray
  948.                     END IF
  949.                     GOSUB integratecom
  950.                 NEXT
  951.             NEXT
  952.             VectorGroup(vecgroupid).LastVector = vectorindex
  953.             GOSUB calculatecom
  954.         NEXT
  955.         FOR j = -span / 2 TO -thickness
  956.             groupidticker = groupidticker + 1
  957.             vecgroupid = groupidticker
  958.             VectorGroup(vecgroupid).Identity = vecgroupid
  959.             VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  960.             VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  961.             VectorGroup(vecgroupid).GroupName = "Grave"
  962.             VectorGroup(vecgroupid).Visible = 0
  963.             VectorGroup(vecgroupid).COMFixed = 1
  964.             VectorGroup(vecgroupid).DIMx = thickness
  965.             VectorGroup(vecgroupid).DIMy = thickness
  966.             VectorGroup(vecgroupid).DIMz = thickness
  967.             VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  968.             FOR k = -thickness TO thickness STEP thickness / 2
  969.                 FOR i = -thickness TO thickness STEP thickness / 2
  970.                     vectorindex = vectorindex + 1
  971.                     vec3Dpos(vectorindex, 1) = xloc + i + (RND - .5) * 2
  972.                     vec3Dpos(vectorindex, 2) = yloc + j + (RND - .5) * 2
  973.                     vec3Dpos(vectorindex, 3) = crux + k + (RND - .5) * 2
  974.                     vec3Dvis(vectorindex) = 0
  975.                     IF RND > .5 THEN
  976.                         vec3Dcolor(vectorindex) = SlateGray
  977.                     ELSE
  978.                         vec3Dcolor(vectorindex) = DarkGray
  979.                     END IF
  980.                     GOSUB integratecom
  981.                 NEXT
  982.             NEXT
  983.             VectorGroup(vecgroupid).LastVector = vectorindex
  984.             GOSUB calculatecom
  985.         NEXT
  986.         FOR j = thickness TO span / 2
  987.             groupidticker = groupidticker + 1
  988.             vecgroupid = groupidticker
  989.             VectorGroup(vecgroupid).Identity = vecgroupid
  990.             VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  991.             VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  992.             VectorGroup(vecgroupid).GroupName = "Grave"
  993.             VectorGroup(vecgroupid).Visible = 0
  994.             VectorGroup(vecgroupid).COMFixed = 1
  995.             VectorGroup(vecgroupid).DIMx = thickness
  996.             VectorGroup(vecgroupid).DIMy = thickness
  997.             VectorGroup(vecgroupid).DIMz = thickness
  998.             VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  999.             FOR k = -thickness TO thickness STEP thickness / 2
  1000.                 FOR i = -thickness TO thickness STEP thickness / 2
  1001.                     vectorindex = vectorindex + 1
  1002.                     vec3Dpos(vectorindex, 1) = xloc + i + (RND - .5) * 2
  1003.                     vec3Dpos(vectorindex, 2) = yloc + j + (RND - .5) * 2
  1004.                     vec3Dpos(vectorindex, 3) = crux + k + (RND - .5) * 2
  1005.                     vec3Dvis(vectorindex) = 0
  1006.                     IF RND > .5 THEN
  1007.                         vec3Dcolor(vectorindex) = SlateGray
  1008.                     ELSE
  1009.                         vec3Dcolor(vectorindex) = DarkGray
  1010.                     END IF
  1011.                     GOSUB integratecom
  1012.                 NEXT
  1013.             NEXT
  1014.             VectorGroup(vecgroupid).LastVector = vectorindex
  1015.             GOSUB calculatecom
  1016.         NEXT
  1017.     NEXT
  1019.  
  1020. 'Heaven's Bottom Layer
  1021. h = 2
  1022. FOR u = -gridsize TO gridsize STEP tilesize
  1023.     FOR v = -gridsize TO gridsize STEP tilesize
  1024.         groupidticker = groupidticker + 1
  1025.         vecgroupid = groupidticker
  1026.         VectorGroup(vecgroupid).Identity = vecgroupid
  1027.         VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1028.         VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1029.         VectorGroup(vecgroupid).GroupName = "Heaven's Bottom Layer"
  1030.         VectorGroup(vecgroupid).Visible = 0
  1031.         VectorGroup(vecgroupid).COMFixed = 1
  1032.         VectorGroup(vecgroupid).DIMx = tilesize / 2
  1033.         VectorGroup(vecgroupid).DIMy = tilesize / 2
  1034.         VectorGroup(vecgroupid).DIMz = 3
  1035.         VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1036.         FOR i = u TO u + tilesize STEP h
  1037.             FOR j = v TO v + tilesize STEP h
  1038.                 vectorindex = vectorindex + 1
  1039.                 vec3Dpos(vectorindex, 1) = i + RND * h - RND * h
  1040.                 vec3Dpos(vectorindex, 2) = j + RND * h - RND * h
  1041.                 vec3Dpos(vectorindex, 3) = 420 - RND
  1042.                 vec3Dvis(vectorindex) = 0
  1043.                 IF RND > .5 THEN
  1044.                     vec3Dcolor(vectorindex) = BlueViolet
  1045.                 ELSE
  1046.                     vec3Dcolor(vectorindex) = Cyan
  1047.                 END IF
  1048.                 GOSUB integratecom
  1049.             NEXT
  1050.         NEXT
  1051.         VectorGroup(vecgroupid).LastVector = vectorindex
  1052.         GOSUB calculatecom
  1053.     NEXT
  1055.  
  1056. 'Hell Spawn
  1057. FOR u = -gridsize TO gridsize STEP tilesize
  1058.     FOR v = -gridsize TO gridsize STEP tilesize
  1059.         groupidticker = groupidticker + 1
  1060.         vecgroupid = groupidticker
  1061.         VectorGroup(vecgroupid).Identity = vecgroupid
  1062.         VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1063.         VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1064.         VectorGroup(vecgroupid).GroupName = "Hell Spawn"
  1065.         VectorGroup(vecgroupid).Visible = 0
  1066.         VectorGroup(vecgroupid).COMFixed = 1
  1067.         VectorGroup(vecgroupid).DIMx = tilesize / 2
  1068.         VectorGroup(vecgroupid).DIMy = tilesize / 2
  1069.         VectorGroup(vecgroupid).DIMz = 35
  1070.         VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1071.         FOR i = u TO u + tilesize STEP tilesize / 5
  1072.             FOR j = v TO v + tilesize STEP tilesize / 5
  1073.                 vectorindex = vectorindex + 1
  1074.                 vec3Dpos(vectorindex, 1) = i + (RND - .5) * tilesize
  1075.                 vec3Dpos(vectorindex, 2) = j + (RND - .5) * tilesize
  1076.                 vec3Dpos(vectorindex, 3) = -350 - RND * 70
  1077.                 vec3Dvel(vectorindex, 1) = 0
  1078.                 vec3Dvel(vectorindex, 2) = 0
  1079.                 vec3Dvel(vectorindex, 3) = 400 * RND
  1080.                 vec3Dvis(vectorindex) = 0
  1081.                 IF RND > .2 THEN
  1082.                     vec3Dcolor(vectorindex) = Red
  1083.                 ELSE
  1084.                     vec3Dcolor(vectorindex) = DarkGoldenRod
  1085.                 END IF
  1086.                 GOSUB integratecom
  1087.             NEXT
  1088.         NEXT
  1089.         VectorGroup(vecgroupid).LastVector = vectorindex
  1090.         GOSUB calculatecom
  1091.         VectorGroup(vecgroupid).COMz = -350 - 35
  1092.     NEXT
  1094.  
  1095. 'Icewall East
  1096. h = 2
  1097. FOR u = -gridsize TO gridsize STEP tilesize
  1098.     FOR v = 0 TO 70 STEP tilesize
  1099.         groupidticker = groupidticker + 1
  1100.         vecgroupid = groupidticker
  1101.         VectorGroup(vecgroupid).Identity = vecgroupid
  1102.         VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1103.         VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1104.         VectorGroup(vecgroupid).GroupName = "Icewall East"
  1105.         VectorGroup(vecgroupid).Visible = 0
  1106.         VectorGroup(vecgroupid).COMFixed = 1
  1107.         VectorGroup(vecgroupid).DIMx = tilesize / 2
  1108.         VectorGroup(vecgroupid).DIMy = tilesize / 2
  1109.         VectorGroup(vecgroupid).DIMz = tilesize / 2
  1110.         VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1111.         FOR i = u TO u + tilesize STEP h
  1112.             FOR j = v TO v + tilesize STEP h
  1113.                 vectorindex = vectorindex + 1
  1114.                 vec3Dpos(vectorindex, 1) = gridsize + tilesize / 2
  1115.                 vec3Dpos(vectorindex, 2) = i + RND * h - RND * h
  1116.                 vec3Dpos(vectorindex, 3) = j + RND * h - RND * h
  1117.                 vec3Dvis(vectorindex) = 0
  1118.                 IF RND > .5 THEN
  1119.                     vec3Dcolor(vectorindex) = White
  1120.                 ELSE
  1121.                     vec3Dcolor(vectorindex) = Ivory
  1122.                 END IF
  1123.                 GOSUB integratecom
  1124.             NEXT
  1125.         NEXT
  1126.         VectorGroup(vecgroupid).LastVector = vectorindex
  1127.         GOSUB calculatecom
  1128.     NEXT
  1130.  
  1131. 'Icewall South
  1132. h = 2
  1133. FOR u = -gridsize TO gridsize STEP tilesize
  1134.     FOR v = 0 TO 70 STEP tilesize
  1135.         groupidticker = groupidticker + 1
  1136.         vecgroupid = groupidticker
  1137.         VectorGroup(vecgroupid).Identity = vecgroupid
  1138.         VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1139.         VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1140.         VectorGroup(vecgroupid).GroupName = "Icewall South"
  1141.         VectorGroup(vecgroupid).Visible = 0
  1142.         VectorGroup(vecgroupid).COMFixed = 1
  1143.         VectorGroup(vecgroupid).DIMx = tilesize / 2
  1144.         VectorGroup(vecgroupid).DIMy = tilesize / 2
  1145.         VectorGroup(vecgroupid).DIMz = tilesize / 2
  1146.         VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1147.         FOR i = u TO u + tilesize STEP h
  1148.             FOR j = v TO v + tilesize STEP h
  1149.                 vectorindex = vectorindex + 1
  1150.                 vec3Dpos(vectorindex, 1) = -gridsize
  1151.                 vec3Dpos(vectorindex, 2) = i + RND * h - RND * h
  1152.                 vec3Dpos(vectorindex, 3) = j + RND * h - RND * h
  1153.                 vec3Dvis(vectorindex) = 0
  1154.                 IF RND > .5 THEN
  1155.                     vec3Dcolor(vectorindex) = White
  1156.                 ELSE
  1157.                     vec3Dcolor(vectorindex) = Ivory
  1158.                 END IF
  1159.                 GOSUB integratecom
  1160.             NEXT
  1161.         NEXT
  1162.         VectorGroup(vecgroupid).LastVector = vectorindex
  1163.         GOSUB calculatecom
  1164.     NEXT
  1166.  
  1167. 'Icewall North
  1168. h = 2
  1169. FOR u = -gridsize TO gridsize STEP tilesize
  1170.     FOR v = 0 TO 70 STEP tilesize
  1171.         groupidticker = groupidticker + 1
  1172.         vecgroupid = groupidticker
  1173.         VectorGroup(vecgroupid).Identity = vecgroupid
  1174.         VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1175.         VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1176.         VectorGroup(vecgroupid).GroupName = "Icewall North"
  1177.         VectorGroup(vecgroupid).Visible = 0
  1178.         VectorGroup(vecgroupid).COMFixed = 1
  1179.         VectorGroup(vecgroupid).DIMx = tilesize / 2
  1180.         VectorGroup(vecgroupid).DIMy = tilesize / 2
  1181.         VectorGroup(vecgroupid).DIMz = tilesize / 2
  1182.         VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1183.         FOR i = u TO u + tilesize STEP h
  1184.             FOR j = v TO v + tilesize STEP h
  1185.                 vectorindex = vectorindex + 1
  1186.                 vec3Dpos(vectorindex, 1) = i + RND * h - RND * h
  1187.                 vec3Dpos(vectorindex, 2) = gridsize + tilesize / 2
  1188.                 vec3Dpos(vectorindex, 3) = j + RND * h - RND * h
  1189.                 vec3Dvis(vectorindex) = 0
  1190.                 IF RND > .5 THEN
  1191.                     vec3Dcolor(vectorindex) = White
  1192.                 ELSE
  1193.                     vec3Dcolor(vectorindex) = Ivory
  1194.                 END IF
  1195.                 GOSUB integratecom
  1196.             NEXT
  1197.         NEXT
  1198.         VectorGroup(vecgroupid).LastVector = vectorindex
  1199.         GOSUB calculatecom
  1200.     NEXT
  1202.  
  1203. 'Icewall West
  1204. h = 2
  1205. FOR u = -gridsize TO gridsize STEP tilesize
  1206.     FOR v = 0 TO 70 STEP tilesize
  1207.         groupidticker = groupidticker + 1
  1208.         vecgroupid = groupidticker
  1209.         VectorGroup(vecgroupid).Identity = vecgroupid
  1210.         VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1211.         VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1212.         VectorGroup(vecgroupid).GroupName = "Icewall West"
  1213.         VectorGroup(vecgroupid).Visible = 0
  1214.         VectorGroup(vecgroupid).COMFixed = 1
  1215.         VectorGroup(vecgroupid).DIMx = tilesize / 2
  1216.         VectorGroup(vecgroupid).DIMy = tilesize / 2
  1217.         VectorGroup(vecgroupid).DIMz = tilesize / 2
  1218.         VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1219.         FOR i = u TO u + tilesize STEP h
  1220.             FOR j = v TO v + tilesize STEP h
  1221.                 vectorindex = vectorindex + 1
  1222.                 vec3Dpos(vectorindex, 1) = i + RND * h - RND * h
  1223.                 vec3Dpos(vectorindex, 2) = -gridsize
  1224.                 vec3Dpos(vectorindex, 3) = j + RND * h - RND * h
  1225.                 vec3Dvis(vectorindex) = 0
  1226.                 IF RND > .5 THEN
  1227.                     vec3Dcolor(vectorindex) = White
  1228.                 ELSE
  1229.                     vec3Dcolor(vectorindex) = Ivory
  1230.                 END IF
  1231.                 GOSUB integratecom
  1232.             NEXT
  1233.         NEXT
  1234.         VectorGroup(vecgroupid).LastVector = vectorindex
  1235.         GOSUB calculatecom
  1236.     NEXT
  1238.  
  1239. 'Lake of Fire
  1240. h = 2
  1241. FOR u = -gridsize TO gridsize STEP tilesize
  1242.     FOR v = -gridsize TO gridsize STEP tilesize
  1243.         groupidticker = groupidticker + 1
  1244.         vecgroupid = groupidticker
  1245.         VectorGroup(vecgroupid).Identity = vecgroupid
  1246.         VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1247.         VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1248.         VectorGroup(vecgroupid).GroupName = "Lake of Fire"
  1249.         VectorGroup(vecgroupid).Visible = 0
  1250.         VectorGroup(vecgroupid).COMFixed = 1
  1251.         VectorGroup(vecgroupid).DIMx = tilesize / 2
  1252.         VectorGroup(vecgroupid).DIMy = tilesize / 2
  1253.         VectorGroup(vecgroupid).DIMz = tilesize / 2
  1254.         VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1255.         FOR i = u TO u + tilesize STEP h
  1256.             FOR j = v TO v + tilesize STEP h
  1257.                 vectorindex = vectorindex + 1
  1258.                 vec3Dpos(vectorindex, 1) = i + RND * h - RND * h
  1259.                 vec3Dpos(vectorindex, 2) = j + RND * h - RND * h
  1260.                 vec3Dpos(vectorindex, 3) = -350 - 70 - RND
  1261.                 vec3Dvis(vectorindex) = 0
  1262.                 IF RND > .2 THEN
  1263.                     vec3Dcolor(vectorindex) = Red
  1264.                 ELSE
  1265.                     vec3Dcolor(vectorindex) = Indigo
  1266.                 END IF
  1267.                 GOSUB integratecom
  1268.             NEXT
  1269.         NEXT
  1270.         VectorGroup(vecgroupid).LastVector = vectorindex
  1271.         GOSUB calculatecom
  1272.     NEXT
  1274.  
  1275. 'Megalith
  1276. ctrx = -90
  1277. ctry = -320
  1278. ctrz = 4
  1279. w = 8
  1280. h = 256
  1281. dens = 100
  1282. FOR k = 1 TO h STEP w
  1283.     FOR i = -h / 20 + k / 20 TO h / 20 - k / 20 STEP w
  1284.         FOR j = -h / 20 + k / 20 TO h / 20 - k / 20 STEP w
  1285.             groupidticker = groupidticker + 1
  1286.             vecgroupid = groupidticker
  1287.             VectorGroup(vecgroupid).Identity = vecgroupid
  1288.             VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1289.             VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1290.             VectorGroup(vecgroupid).GroupName = "Megalith"
  1291.             VectorGroup(vecgroupid).Visible = 0
  1292.             VectorGroup(vecgroupid).COMFixed = 1
  1293.             VectorGroup(vecgroupid).DIMx = w / 2
  1294.             VectorGroup(vecgroupid).DIMy = w / 2
  1295.             VectorGroup(vecgroupid).DIMz = w / 2
  1296.             VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1297.             FOR q = 1 TO dens
  1298.                 vectorindex = vectorindex + 1
  1299.                 vec3Dpos(vectorindex, 1) = ctrx + i + (RND - .5) * w
  1300.                 vec3Dpos(vectorindex, 2) = ctry + j + (RND - .5) * w
  1301.                 vec3Dpos(vectorindex, 3) = ctrz + k + (RND - .5) * w
  1302.                 vec3Dvis(vectorindex) = 0
  1303.                 IF RND > .5 THEN
  1304.                     vec3Dcolor(vectorindex) = Cyan
  1305.                 ELSE
  1306.                     vec3Dcolor(vectorindex) = Teal
  1307.                 END IF
  1308.                 GOSUB integratecom
  1309.             NEXT
  1310.             VectorGroup(vecgroupid).LastVector = vectorindex
  1311.             GOSUB calculatecom
  1312.         NEXT
  1313.     NEXT
  1315.  
  1316. 'Pyramid
  1317. ctrx = -90
  1318. ctry = -120
  1319. ctrz = 4
  1320. w = 8
  1321. h = 56
  1322. dens = 50
  1323. FOR k = 1 TO h STEP w
  1324.     FOR i = -h / 2 + k / 2 TO h / 2 - k / 2 STEP w
  1325.         FOR j = -h / 2 + k / 2 TO h / 2 - k / 2 STEP w
  1326.             groupidticker = groupidticker + 1
  1327.             vecgroupid = groupidticker
  1328.             VectorGroup(vecgroupid).Identity = vecgroupid
  1329.             VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1330.             VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1331.             VectorGroup(vecgroupid).GroupName = "Pyramid"
  1332.             VectorGroup(vecgroupid).Visible = 0
  1333.             VectorGroup(vecgroupid).COMFixed = 1
  1334.             VectorGroup(vecgroupid).DIMx = tilesize / 2
  1335.             VectorGroup(vecgroupid).DIMy = tilesize / 2
  1336.             VectorGroup(vecgroupid).DIMz = tilesize / 2
  1337.             VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1338.             FOR q = 1 TO dens
  1339.                 vectorindex = vectorindex + 1
  1340.                 vec3Dpos(vectorindex, 1) = ctrx + i + (RND - .5) * w
  1341.                 vec3Dpos(vectorindex, 2) = ctry + j + (RND - .5) * w
  1342.                 vec3Dpos(vectorindex, 3) = ctrz + k + (RND - .5) * w
  1343.                 vec3Dvis(vectorindex) = 0
  1344.                 IF RND > .5 THEN
  1345.                     vec3Dcolor(vectorindex) = DarkGoldenRod
  1346.                 ELSE
  1347.                     vec3Dcolor(vectorindex) = GoldenRod
  1348.                 END IF
  1349.                 GOSUB integratecom
  1350.             NEXT
  1351.             VectorGroup(vecgroupid).LastVector = vectorindex
  1352.             GOSUB calculatecom
  1353.         NEXT
  1354.     NEXT
  1356.  
  1357. 'Rain
  1358. FOR u = -gridsize TO gridsize STEP tilesize
  1359.     FOR v = -gridsize TO gridsize STEP tilesize
  1360.         groupidticker = groupidticker + 1
  1361.         vecgroupid = groupidticker
  1362.         VectorGroup(vecgroupid).Identity = vecgroupid
  1363.         VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1364.         VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1365.         VectorGroup(vecgroupid).GroupName = "Rain"
  1366.         VectorGroup(vecgroupid).Visible = 0
  1367.         VectorGroup(vecgroupid).COMFixed = 1
  1368.         VectorGroup(vecgroupid).DIMx = tilesize / 2
  1369.         VectorGroup(vecgroupid).DIMy = tilesize / 2
  1370.         VectorGroup(vecgroupid).DIMz = 35
  1371.         VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1372.         FOR i = u TO u + tilesize STEP tilesize '/ 3
  1373.             FOR j = v TO v + tilesize STEP tilesize '/ 3
  1374.                 vectorindex = vectorindex + 1
  1375.                 vec3Dpos(vectorindex, 1) = i + (RND - .5) * tilesize
  1376.                 vec3Dpos(vectorindex, 2) = j + (RND - .5) * tilesize
  1377.                 vec3Dpos(vectorindex, 3) = RND * 70
  1378.                 vec3Dvel(vectorindex, 1) = 0
  1379.                 vec3Dvel(vectorindex, 2) = 0
  1380.                 vec3Dvel(vectorindex, 3) = -400 * RND
  1381.                 vec3Dvis(vectorindex) = 0
  1382.                 IF RND > 5 THEN
  1383.                     vec3Dcolor(vectorindex) = Aquamarine
  1384.                 ELSE
  1385.                     vec3Dcolor(vectorindex) = DodgerBlue
  1386.                 END IF
  1387.                 GOSUB integratecom
  1388.             NEXT
  1389.         NEXT
  1390.         VectorGroup(vecgroupid).LastVector = vectorindex
  1391.         GOSUB calculatecom
  1392.         VectorGroup(vecgroupid).COMz = 35
  1393.     NEXT
  1395.  
  1396. 'Sky
  1397. h = 2
  1398. FOR u = -gridsize TO gridsize STEP tilesize
  1399.     FOR v = -gridsize TO gridsize STEP tilesize
  1400.         groupidticker = groupidticker + 1
  1401.         vecgroupid = groupidticker
  1402.         VectorGroup(vecgroupid).Identity = vecgroupid
  1403.         VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1404.         VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1405.         VectorGroup(vecgroupid).GroupName = "Sky"
  1406.         VectorGroup(vecgroupid).Visible = 0
  1407.         VectorGroup(vecgroupid).COMFixed = 1
  1408.         VectorGroup(vecgroupid).DIMx = tilesize / 2
  1409.         VectorGroup(vecgroupid).DIMy = tilesize / 2
  1410.         VectorGroup(vecgroupid).DIMz = 3
  1411.         VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1412.         FOR i = u TO u + tilesize STEP h
  1413.             FOR j = v TO v + tilesize STEP h
  1414.                 vectorindex = vectorindex + 1
  1415.                 vec3Dpos(vectorindex, 1) = i + (RND - RND) * h
  1416.                 vec3Dpos(vectorindex, 2) = j + (RND - RND) * h
  1417.                 vec3Dpos(vectorindex, 3) = 70 + (RND - RND) * h
  1418.                 vec3Dvel(vectorindex, 1) = (RND - RND) * 2
  1419.                 vec3Dvel(vectorindex, 2) = (RND - RND) * 2
  1420.                 vec3Dvel(vectorindex, 3) = (RND - RND) * 2
  1421.                 vec3Danv(vectorindex, 1) = 0
  1422.                 vec3Danv(vectorindex, 2) = 0
  1423.                 vec3Danv(vectorindex, 3) = 0
  1424.                 vec3Dvis(vectorindex) = 0
  1425.                 IF RND > .5 THEN
  1426.                     vec3Dcolor(vectorindex) = Snow
  1427.                 ELSE
  1428.                     vec3Dcolor(vectorindex) = RoyalBlue
  1429.                 END IF
  1430.                 GOSUB integratecom
  1431.             NEXT
  1432.         NEXT
  1433.         VectorGroup(vecgroupid).LastVector = vectorindex
  1434.         GOSUB calculatecom
  1435.     NEXT
  1437.  
  1438. 'Stars
  1439. h = 5
  1440. FOR w = 1 TO 5
  1441.     FOR u = -gridsize TO gridsize STEP tilesize
  1442.         FOR v = -gridsize TO gridsize STEP tilesize
  1443.             groupidticker = groupidticker + 1
  1444.             vecgroupid = groupidticker
  1445.             VectorGroup(vecgroupid).Identity = vecgroupid
  1446.             VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1447.             VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1448.             VectorGroup(vecgroupid).GroupName = "Stars"
  1449.             VectorGroup(vecgroupid).Visible = 0
  1450.             VectorGroup(vecgroupid).COMFixed = 1
  1451.             VectorGroup(vecgroupid).DIMx = tilesize / 2
  1452.             VectorGroup(vecgroupid).DIMy = tilesize / 2
  1453.             VectorGroup(vecgroupid).DIMz = tilesize / 2
  1454.             VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1455.             FOR i = u TO u + tilesize STEP h
  1456.                 FOR j = v TO v + tilesize STEP h
  1457.                     IF RND > 1 - w / 5 THEN
  1458.                         vectorindex = vectorindex + 1
  1459.                         vec3Dpos(vectorindex, 1) = i + RND * h - RND * h
  1460.                         vec3Dpos(vectorindex, 2) = j + RND * h - RND * h
  1461.                         vec3Dpos(vectorindex, 3) = w * 70 + RND * 70
  1462.                         vec3Dvis(vectorindex) = 0
  1463.                         IF RND > .5 THEN
  1464.                             vec3Dcolor(vectorindex) = GhostWhite
  1465.                         ELSE
  1466.                             IF RND > .5 THEN
  1467.                                 vec3Dcolor(vectorindex) = White
  1468.                             ELSE
  1469.                                 vec3Dcolor(vectorindex) = DarkGray
  1470.                             END IF
  1471.                         END IF
  1472.                         GOSUB integratecom
  1473.                     END IF
  1474.                 NEXT
  1475.             NEXT
  1476.             VectorGroup(vecgroupid).LastVector = vectorindex
  1477.             GOSUB calculatecom
  1478.         NEXT
  1479.     NEXT
  1481.  
  1482. 'Sun
  1483. radius = 10
  1484. dx = .0628
  1485. dy = .0628
  1486. xl = 0: xr = 2 * pi
  1487. yl = 0: yr = pi
  1488. xrange = 1 + INT((-xl + xr) / dx)
  1489. yrange = 1 + INT((-yl + yr) / dy)
  1490. FOR i = 1 TO xrange STEP 10
  1491.     FOR j = 1 TO yrange STEP 10
  1492.         groupidticker = groupidticker + 1
  1493.         vecgroupid = groupidticker
  1494.         VectorGroup(vecgroupid).Identity = vecgroupid
  1495.         VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1496.         VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1497.         VectorGroup(vecgroupid).GroupName = "Sun"
  1498.         VectorGroup(vecgroupid).Visible = 0
  1499.         VectorGroup(vecgroupid).COMFixed = 1
  1500.         VectorGroup(vecgroupid).DIMx = radius
  1501.         VectorGroup(vecgroupid).DIMy = radius
  1502.         VectorGroup(vecgroupid).DIMz = radius
  1503.         VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1504.         FOR u = i TO i + 10 STEP 1
  1505.             FOR v = j TO j + 10 STEP 1
  1506.                 vectorindex = vectorindex + 1
  1507.                 theta = u * dx - dx
  1508.                 phi = v * dy - dy
  1509.                 vec3Dpos(vectorindex, 1) = radius * SIN(phi) * COS(theta)
  1510.                 vec3Dpos(vectorindex, 2) = radius * SIN(phi) * SIN(theta)
  1511.                 vec3Dpos(vectorindex, 3) = 90 + radius * COS(phi)
  1512.                 vec3Dvis(vectorindex) = 0
  1513.                 IF RND > .5 THEN
  1514.                     vec3Dcolor(vectorindex) = Sunglow
  1515.                 ELSE
  1516.                     vec3Dcolor(vectorindex) = SunsetOrange
  1517.                 END IF
  1518.                 GOSUB integratecom
  1519.             NEXT
  1520.         NEXT
  1521.         GOSUB integratecom
  1522.         VectorGroup(vecgroupid).LastVector = vectorindex
  1523.         GOSUB calculatecom
  1524.     NEXT
  1526.  
  1527. 'Moon
  1528. radius = 4
  1529. au = 60
  1530. dx = (2 * pi / radius) * .05
  1531. dy = (2 * pi / radius) * .05
  1532. xl = 0: xr = 2 * pi
  1533. yl = 0: yr = pi
  1534. xrange = 1 + INT((-xl + xr) / dx)
  1535. yrange = 1 + INT((-yl + yr) / dy)
  1536. groupidticker = groupidticker + 1
  1537. vecgroupid = groupidticker
  1538. VectorGroup(vecgroupid).Identity = vecgroupid
  1539. VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1540. VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1541. VectorGroup(vecgroupid).GroupName = "Moon"
  1542. VectorGroup(vecgroupid).Visible = 0
  1543. VectorGroup(vecgroupid).ForceAnimate = 1
  1544. VectorGroup(vecgroupid).COMFixed = 0
  1545. VectorGroup(vecgroupid).ROTx = 0
  1546. VectorGroup(vecgroupid).ROTy = 0
  1547. VectorGroup(vecgroupid).ROTz = 90
  1548. VectorGroup(vecgroupid).REVx = 1.5
  1549. VectorGroup(vecgroupid).REVy = 0
  1550. VectorGroup(vecgroupid).REVz = 0
  1551. VectorGroup(vecgroupid).DIMx = 2 * radius + 1
  1552. VectorGroup(vecgroupid).DIMy = 2 * radius + 1
  1553. VectorGroup(vecgroupid).DIMz = 2 * radius + 1
  1554. VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1555. FOR i = 1 TO xrange
  1556.     FOR j = 1 TO yrange
  1557.         vectorindex = vectorindex + 1
  1558.         theta = i * dx - dx
  1559.         phi = j * dy - dy
  1560.         vec3Dpos(vectorindex, 1) = au + radius * SIN(phi) * COS(theta)
  1561.         vec3Dpos(vectorindex, 2) = radius * SIN(phi) * SIN(theta)
  1562.         vec3Dpos(vectorindex, 3) = 90 + radius * COS(phi)
  1563.         vec3Danv(vectorindex, 1) = 0
  1564.         vec3Danv(vectorindex, 2) = 0
  1565.         vec3Danv(vectorindex, 3) = 1.5
  1566.         vec3Dvis(vectorindex) = 0
  1567.         IF RND > .5 THEN
  1568.             vec3Dcolor(vectorindex) = Gray
  1569.         ELSE
  1570.             vec3Dcolor(vectorindex) = PaleGoldenRod
  1571.         END IF
  1572.         GOSUB integratecom
  1573.         VectorGroup(vecgroupid).LastVector = vectorindex
  1574.         GOSUB calculatecom
  1575.     NEXT
  1577.  
  1578. 'Waves or Particles? (1)
  1579. FOR i = 1 TO 5 STEP 1
  1580.     FOR k = 1 TO 5 STEP 1
  1581.         groupidticker = groupidticker + 1
  1582.         vecgroupid = groupidticker
  1583.         VectorGroup(vecgroupid).Identity = vecgroupid
  1584.         VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1585.         VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1586.         VectorGroup(vecgroupid).GroupName = "Waves or Particles?"
  1587.         VectorGroup(vecgroupid).Visible = 0
  1588.         VectorGroup(vecgroupid).COMFixed = 1
  1589.         VectorGroup(vecgroupid).DIMx = 4
  1590.         VectorGroup(vecgroupid).DIMy = 1
  1591.         VectorGroup(vecgroupid).DIMz = 4
  1592.         VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1593.         FOR u = i TO i + 1 STEP .05
  1594.             FOR v = k TO k + 1 STEP .05
  1595.                 vectorindex = vectorindex + 1
  1596.                 vec3Dpos(vectorindex, 1) = 70 + 7 * u
  1597.                 vec3Dpos(vectorindex, 2) = 80 + 1 * COS((u ^ 2 - v ^ 2))
  1598.                 vec3Dpos(vectorindex, 3) = 10 + 7 * v
  1599.                 vec3Dvis(vectorindex) = 0
  1600.                 IF vec3Dpos(vectorindex, 2) < 80 THEN
  1601.                     vec3Dcolor(vectorindex) = DarkBlue
  1602.                 ELSE
  1603.                     vec3Dcolor(vectorindex) = DeepPink
  1604.                 END IF
  1605.                 GOSUB integratecom
  1606.             NEXT
  1607.         NEXT
  1608.         VectorGroup(vecgroupid).LastVector = vectorindex
  1609.         GOSUB calculatecom
  1610.     NEXT
  1612.  
  1613. 'Waves or Particles? (2)
  1614. FOR i = 1 TO 5 STEP 1
  1615.     FOR k = 1 TO 5 STEP 1
  1616.         groupidticker = groupidticker + 1
  1617.         vecgroupid = groupidticker
  1618.         VectorGroup(vecgroupid).Identity = vecgroupid
  1619.         VectorGroup(vecgroupid).Pointer = vecgroupid + 1
  1620.         VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1621.         VectorGroup(vecgroupid).GroupName = "Particles or Waves?"
  1622.         VectorGroup(vecgroupid).Visible = 0
  1623.         VectorGroup(vecgroupid).COMFixed = 1
  1624.         VectorGroup(vecgroupid).DIMx = 4
  1625.         VectorGroup(vecgroupid).DIMy = 1
  1626.         VectorGroup(vecgroupid).DIMz = 4
  1627.         VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1628.         FOR u = i TO i + 1 STEP .05
  1629.             FOR v = k TO k + 1 STEP .05
  1630.                 vectorindex = vectorindex + 1
  1631.                 vec3Dpos(vectorindex, 1) = -7 * u
  1632.                 vec3Dpos(vectorindex, 2) = 80 + 1 * COS(2 * ((u - 7) ^ 2 - (v - 5) ^ 2))
  1633.                 vec3Dpos(vectorindex, 3) = 10 + 7 * v
  1634.                 vec3Dvis(vectorindex) = 0
  1635.                 IF vec3Dpos(vectorindex, 2) < 80 THEN
  1636.                     vec3Dcolor(vectorindex) = Magenta
  1637.                 ELSE
  1638.                     vec3Dcolor(vectorindex) = Chocolate
  1639.                 END IF
  1640.                 GOSUB integratecom
  1641.             NEXT
  1642.         NEXT
  1643.         VectorGroup(vecgroupid).LastVector = vectorindex
  1644.         GOSUB calculatecom
  1645.     NEXT
  1647.  
  1648. '__ZZZ
  1649. groupidticker = groupidticker + 1
  1650. vecgroupid = groupidticker
  1651. VectorGroup(vecgroupid).Identity = vecgroupid
  1652. VectorGroup(vecgroupid).Pointer = -999
  1653. VectorGroup(vecgroupid).Lagger = vecgroupid - 1
  1654. VectorGroup(vecgroupid).GroupName = "__ZZZ"
  1655. VectorGroup(vecgroupid).COMFixed = 1
  1656. VectorGroup(vecgroupid).FirstVector = vectorindex + 1
  1657. FOR r = 1 TO 1
  1658.     vectorindex = vectorindex + 1
  1659.     vec3Dpos(vectorindex, 1) = 0
  1660.     vec3Dpos(vectorindex, 2) = 0
  1661.     vec3Dpos(vectorindex, 3) = -1000
  1662.     vec3Dcolor(vectorindex) = White
  1663.     GOSUB integratecom
  1665. VectorGroup(vecgroupid).LastVector = vectorindex
  1666. GOSUB calculatecom
  1668.  

screenshot.png
* Sanctum.bas (Filesize: 62.07 KB, Downloads: 450)
« Last Edit: August 25, 2019, 12:10:02 am by The Librarian »
Logged
Pages: [1]   Go Up
« previous next »