Flight simulator

[MasterGy]

Hello !
I'm thinking of a game. I’ve tried a shooter-like game before, now I want something that shows the speed. I was thinking of a fighter simulator. (I really want a car game, but to this day I can't approach the physics of the car)
Because the plane is fast and a trajectory flies over quickly, I generated the seemingly random ground from a function.
I searched the net for a solution and found a very interesting solution to the real-time endless trajectory: Perlin noise. I used this.
Check it out, it’s interesting, and it can be used for a lot of things.
https://www.arendpeter.com/Perlin_Noise.html
We're just racing for now. By turning the whole thing left and right, we turn and “turn” with the A-D keys. If you become an enemy, it may be interesting to shoot him.
The source code is concise, but whoever cares will disassemble it.

Code: QB64: [Select]

1. 'mikrolepes
2.  
3. Const pip180 = 3.141592 / 180
4. Const perlin_pers = 7 'amplitudo
5. Const perlin_level = 5 'level of resolution wave
6. Const perlin_multiplier = 0.006
7. Const perlin_smooth_noise_interpolation_type = 0 '0-none '1-standard
8. Const perlin_flat_soil_limit = .8
9. monx = _DesktopWidth * .7: mony = monx / 16 * 9 'screen size
10.  
11. map_shadow = 18
12. position_speed = .7
13. map_resolution = 70
14.  
15. map_dimxy = 1200
16. map_dimz = 250
17.  
18. map_zoom_xy = 6
19. map_zoom_distance = 15
20. mouse_sens_xy = .01
21. mouse_sens_z = .01
22.  
23. cam(4) = pip180 * 270
24.  
25. me(0) = 2000: me(1) = 2000
26.  
27. 'create pseudo random buffer
28. Dim Shared noise_rand_c: noise_rand_c = 10000: Dim Shared noise_rand(noise_rand_c - 1): For t = 0 To noise_rand_c - 1: noise_rand(t) = Rnd(1): Next t
29.  
30. 'create texture
31. Dim color_temp As _Integer64: Dim Shared mapz_min, mapz_max
32. text_size = 8: text_size_marg = 8: text_deep = 40: text_height_scale = 15: Dim texture(text_height_scale - 1, text_deep - 1)
33. For t = 0 To text_deep - 1: For t2 = 0 To text_height_scale - 1
34.     dark = 1 - (1 / text_deep * t)
35.     color_temp = _RGB32((255 * dark), (100 * dark), ((255 / text_height_scale * t2)) * dark)
36.     temp = _NewImage(text_size, text_size, 32): _Dest temp: Cls: marg = text_size * text_size_marg / 100
37. Line (marg, marg)-(text_size - marg, text_size - marg), color_temp, BF: texture(t2, t) = _CopyImage(temp, 33): _FreeImage temp: Next t2, t
38.  
39. 'find min/max Z
40. mapz_min = 9999999: mapz_max = -mapz_min: For t = 0 To 1999: c = Perlin2D(1000 * Rnd(1), 1000 * Rnd(1))
41.     If mapz_min > c Then mapz_min = c
42.     If mapz_max < c Then mapz_max = c
43. Next t
44.  
45.  
46. 'fill map-buffer
47. Dim map(map_resolution - 1, map_resolution - 1, 5): perl_setx = Int(me(0)): perl_sety = Int(me(1))
48. For map_x = 0 To map_resolution - 1: For map_y = 0 To map_resolution - 1: map_z = Perlin2D(perl_setx + map_x, perl_sety + map_y)
49. map((map_x + perl_setx) Mod map_resolution, (map_y + perl_sety) Mod map_resolution, 0) = 1 / (mapz_max - mapz_min) * (map_z - mapz_min): Next map_y, map_x
50.  
51.  
52. mon = _NewImage(monx, mony, 32): Screen mon: _FullScreen: _MouseHide
53.  
54. 'MAP array
55. '0-perlin Z-data
56. '1-maptriangle calculate X
57. '2-maptriangle calculate Y
58. '3-maptriangle calculate Z
59. '4-distance from me (color)
60.  
61.  
62.  
63. Do: _Limit 40
64.  
65.     mousex = 0: mousey = 0: While _MouseInput: mousex = mousex + _MouseMovementX: mousey = mousey + _MouseMovementY: Wend
66.     'cam(3) = cam(3) + mousex * mouse_sens_xy
67.     cam(4) = cam(4) + mousey * mouse_sens_z
68.     '    cam(3) = mouse_sens * .1
69.     rotx = rotx + mousex / 100
70.     max_incli = 60: If Abs(rotx) > (max_incli * pip180) Then rotx = max_incli * pip180 * Sgn(rotx)
71.     cam(3) = cam(3) + rotx / 40
72.  
73.     'control position
74.     kw = _KeyDown(119)
75.     kw = -1
76.     ks = _KeyDown(115): ka = _KeyDown(97): kd = _KeyDown(100): new_direction = (Abs(ka Or kd Or kw) Or -Abs(ks)) * position_speed
77.  
78.     deg_XY = -90 * Abs(ka) + 90 * Abs(kd)
79.  
80.     szog_xy = cam(3) + deg_XY * pip180
81.  
82.     '   szog_z = cam(4)
83.  
84.     me(0) = me(0) - Sin(szog_xy) * Cos(szog_z) * new_direction
85.     me(1) = me(1) - Cos(szog_xy) * Cos(szog_z) * new_direction
86.     '    me(2) = me(2) + SIN(szog_z) * new_direction '* ABS(ka = 0 AND kd = 0)
87.  
88.     '    cam(2) = me(2)
89.     '    LOCATE 1, 1: PRINT cam(2)
90.  
91.  
92.  
93.     'replace the missing row in the buffer line Y
94.     Do While Int(Int(me(0))) <> perl_setx: dir = -Sgn(Int(me(0)) - perl_setx): perl_setx = perl_setx - dir
95.         For map_y = 0 To map_resolution - 1: temp = perl_setx + Abs(dir = -1) * (map_resolution - 1): map_z = Perlin2D(temp, perl_sety + map_y)
96.     map(temp Mod map_resolution, (map_y + perl_sety) Mod map_resolution, 0) = 1 / (mapz_max - mapz_min) * (map_z - mapz_min): Next map_y: Loop
97.  
98.     'replace the missing row in the buffer line X
99.     Do While Int(Int(me(1))) <> perl_sety: dir = -Sgn(Int(me(1)) - perl_sety): perl_sety = perl_sety - dir
100.         For map_x = 0 To map_resolution - 1: temp = perl_sety + Abs(dir = -1) * (map_resolution - 1): map_z = Perlin2D(perl_setx + map_x, temp)
101.     map((map_x + perl_setx) Mod map_resolution, temp Mod map_resolution, 0) = 1 / (mapz_max - mapz_min) * (map_z - mapz_min): Next map_x: Loop
102.  
103.  
104.  
105.     cosrotz = Cos(cam(3)): sinrotz = Sin(cam(3)): cosrotx = Cos(cam(4)): sinrotx = Sin(cam(4))
106.  
107.     'calculating position and textures
108.     '    e1 = perl_setx - INT(perl_setx)
109.  
110.     '   e2 = perl_sety - INT(perl_sety)
111.     For map_x = 0 To map_resolution - 1: px = -map_dimxy / 2 + map_dimxy / map_resolution * (map_x - e1) - cam(0)
112.         For map_y = 0 To map_resolution - 1: py = -map_dimxy / 2 + map_dimxy / map_resolution * (map_y - e2) - cam(1)
113.             map_z = map((map_x + perl_setx) Mod map_resolution, (map_y + perl_sety) Mod map_resolution, 0)
114.             If map_z > perlin_flat_soil_limit Then map_z = perlin_flat_soil_limit
115.             pz2 = map_dimz * map_z - cam(2)
116.             px3 = px * cosrotz - py * sinrotz: py2 = px * sinrotz + py * cosrotz: py3 = py2 * cosrotx - pz2 * sinrotx: pz3 = py2 * sinrotx + pz2 * cosrotx
117.             px4 = (px3 * Cos(rotx)) - (py3 * Sin(rotx)): py4 = (px3 * Sin(rotx)) + (py3 * Cos(rotx))
118.             map(map_x, map_y, 1) = -px4 * map_zoom_xy: map(map_x, map_y, 2) = -py4 * map_zoom_xy: map(map_x, map_y, 3) = -pz3 * map_zoom_distance
119.             map(map_x, map_y, 4) = Int(Sqr(px3 * px3 + py3 * py3 + pz3 * pz3) / map_shadow)
120.             If map(map_x, map_y, 4) > text_deep - 1 Then map(map_x, map_y, 4) = text_deep - 1
121.             map(map_x, map_y, 5) = Int(text_height_scale * map_z)
122.             If map(map_x, map_y, 5) > text_height_scale - 1 Then map(map_x, map_y, 5) = text_height_scale - 1
123.             If map(map_x, map_y, 5) < 0 Then map(map_x, map_y, 5) = 0
124.     Next map_y, map_x
125.  
126.     'do maptriangle from squares !
127.     For map_x = 0 To map_resolution - 2: For map_y = 0 To map_resolution - 2
128.         m0x = map(map_x, map_y, 1): m0y = map(map_x, map_y, 2): m0z = map(map_x, map_y, 3)
129.         m1x = map(map_x + 1, map_y, 1): m1y = map(map_x + 1, map_y, 2): m1z = map(map_x + 1, map_y, 3)
130.         m2x = map(map_x, map_y + 1, 1): m2y = map(map_x, map_y + 1, 2): m2z = map(map_x, map_y + 1, 3)
131.         m3x = map(map_x + 1, map_y + 1, 1): m3y = map(map_x + 1, map_y + 1, 2): m3z = map(map_x + 1, map_y + 1, 3)
132.         atexture = texture(map(map_x, map_y, 5), map(map_x, map_y, 4))
133.         _MapTriangle (0, 0)-(text_size - 1, 0)-(0, text_size - 1), atexture To(m0x, m0y, m0z)-(m1x, m1y, m1z)-(m2x, m2y, m2z)
134.         _MapTriangle (0, 0)-(text_size - 1, 0)-(0, text_size - 1), atexture To(m3x, m3y, m3z)-(m1x, m1y, m1z)-(m2x, m2y, m2z)
135.     Next map_y, map_x
136.  
137.     _Display: Cls
138.  
139. Loop
140.  
141.  
142.  
143.  
144.  
145. Function noise (x, y) 'reading pseudo random buffer
146.     x2 = Int(x): y2 = Int(y): a = Int(Abs(x2) + Abs(y2)) Mod noise_rand_c: b = Abs(x2) Mod noise_rand_c: c = noise_rand(a) + noise_rand(b): noise = c - Int(c)
147. End Function
148.  
149. Function SmoothNoise (x, y)
150.     Select Case perlin_smooth_noise_interpolation_type
151.         Case 0: SmoothNoise = noise(x, y)
152.         Case 1: corners = (noise(x - 1, y - 1) + noise(x + 1, y - 1) + noise(x - 1, y + 1) + noise(x + 1, y + 1)) / 16
153.             sides = (noise(x - 1, y) + noise(x + 1, y) + noise(x, y - 1) + noise(x, y + 1)) / 8: center = noise(x, y) / 4: SmoothNoise = corners + sides + center
154.     End Select
155. End Function
156.  
157. Function Perlin2D (x, y): total = 0: For t = 0 To perlin_level - 1: frequency = 2 ^ t * perlin_multiplier: amplitude = perlin_pers ^ t
158.     total = total + InterpolatedNoise(x * frequency, y * frequency) * amplitude: Next t: Perlin2D = total
159. End Function
160.  
161. Function InterpolatedNoise (x, y): integer_X = Int(x): fractional_X = x - integer_X: integer_y = Int(y): fractional_Y = y - integer_y
162.     v1 = SmoothNoise(integer_X, integer_y): v2 = SmoothNoise(integer_X + 1, integer_y)
163.     v3 = SmoothNoise(integer_X, integer_y + 1): v4 = SmoothNoise(integer_X + 1, integer_y + 1)
164.     i1 = Interpolate(v1, v2, fractional_X): i2 = Interpolate(v3, v4, fractional_X): InterpolatedNoise = Interpolate(i1, i2, fractional_Y)
165. End Function
166.  
167. Function Interpolate (a, b, x): Interpolate = a * (1 - x) + b * x: End Function
168.  
169.  

[bplus]

More genius work from MasterGy! Nice :)

[johnno56]

Very nicely done! So quick and smooth. Reminds me of the old DOS game 'Magic Carpet'... Looking forward to the next iteration...

[Qwerkey]

More genius work from MasterGy! Nice :)

@bplus is not wrong.  Very impressive (up to the standard of @Ashish with _GL) terrain.  Don't fancy trying to land a plane on that terrain, though!

If you just let this code run, it eventually crashes (aeronautical pun!).

[bplus]

How long for you before it crashes? @Qwerkey

Heck of a line (101):

Code: QB64: [Select]

1. map((map_x + perl_setx) Mod map_resolution, temp Mod map_resolution, 0) = 1 / (mapz_max - mapz_min) * (map_z - mapz_min): Next map_x: Loop

Map is 3 dim's

Code: QB64: [Select]

1. Dim map(map_resolution - 1, map_resolution - 1, 5)
2.  

I see only 1 comma inside a () inside line 101???

[MasterGy]

Hello !
Thanks for trying it! I didn’t know the “Magic Carpet” game, but I watched it and really like the scenery! I imagined a similar graphic, I wish it would look like mine.
Yes, you are right that real physics, real fighter control, proper runway would be hard to put together from that. I don't strive for 100% reality, but rather I just want to show the feeling of speed, flight, speed, that would give the gaming experience.
I always develop in QB64 1.3 because it creates an exe quickly. When do you experience a crash? I haven't had such a problem yet.
Bplus, the MAP array records a lot. Calculating an elevation point with “perlin” takes a long time, so the MAP array is a buffer. As we move on the track, it always counts only the next X or Y row, then adds it to the MAP array and shifts the entire array by 1.
The position of the points, the distance of the point from us, the height value, and more can be used to determine objects, such as trees, columns, flowers, based on different conditions. To be or not to be in a particular place.
Thank you Qwerky for making a mistake. In vain do I think, I do not understand why the address of the array may have fallen out of range. I put INT in it, I hope it filters out the trouble.
I refined it a bit more. I want it to look like “Magic Carpet”.

Code: QB64: [Select]

1. 'mikrolepes
2.  
3. CONST pip180 = 3.141592 / 180
4. CONST perlin_pers = 7 'amplitudo
5. CONST perlin_level = 5 'level of resolution wave
6. CONST perlin_multiplier = 0.006
7. CONST perlin_smooth_noise_interpolation_type = 0 '0-none '1-standard
8. CONST perlin_flat_soil_limit = .8
9. DIM SHARED rotating(2), cosrotz, sinrotz, cosrotx, sinrotx, map_zoom_xy, map_zoom_distance, rotx
10.  
11. monx = _DESKTOPWIDTH * .7: mony = monx / 16 * 9 'screen size
12.  
13. map_shadow = 15
14. position_speed = .6
15. map_resolution = 60
16.  
17. map_dimxy = 1200
18. map_dimz = 250
19.  
20. map_zoom_xy = 6
21. map_zoom_distance = 15
22. mouse_sens_xy = .01
23. mouse_sens_z = .01
24.  
25. cam(4) = pip180 * 270
26.  
27. me(0) = 2000: me(1) = 2000
28. DIM clmn(1999, 5): ration_column = 0.001
29.  
30.  
31. 'create pseudo random buffer
32. DIM SHARED noise_rand_c: noise_rand_c = 10000: DIM SHARED noise_rand(noise_rand_c - 1): FOR t = 0 TO noise_rand_c - 1: noise_rand(t) = RND(1): NEXT t
33.  
34.  
35.  
36.  
37.  
38.  
39. 'create texture
40. DIM color_temp AS _INTEGER64: DIM SHARED mapz_min, mapz_max
41. text_size = 8: text_size_marg = 8: text_deep = 40: text_height_scale = 15: DIM texture(text_height_scale - 1, text_deep - 1)
42. FOR t = 0 TO text_deep - 1: FOR t2 = 0 TO text_height_scale - 1: dark = 1 - (1 / (text_deep - 1) * t)
43.         color_temp = _RGB32((255 * dark), (100 * dark), ((255 / text_height_scale * t2)) * dark)
44.         temp = _NEWIMAGE(text_size, text_size, 32): _DEST temp: CLS: marg = text_size * text_size_marg / 100
45. LINE (marg, marg)-(text_size - marg, text_size - marg), color_temp, BF: texture(t2, t) = _COPYIMAGE(temp, 33): _FREEIMAGE temp: NEXT t2, t
46.  
47.  
48.  
49. 'find min/max Z
50. mapz_min = 9999999: mapz_max = -mapz_min: FOR t = 0 TO 1999: c = Perlin2D(1000 * RND(1), 1000 * RND(1))
51.     IF mapz_min > c THEN mapz_min = c
52.     IF mapz_max < c THEN mapz_max = c
53. NEXT t
54.  
55.  
56. 'fill map-buffer
57. DIM map(map_resolution - 1, map_resolution - 1, 9): perl_setx = INT(me(0)): perl_sety = INT(me(1))
58. FOR map_x = 0 TO map_resolution - 1: FOR map_y = 0 TO map_resolution - 1: map_z = Perlin2D(perl_setx + map_x, perl_sety + map_y)
59.         map((map_x + perl_setx) MOD map_resolution, (map_y + perl_sety) MOD map_resolution, 0) = 1 / (mapz_max - mapz_min) * (map_z - mapz_min)
60. map((map_x + perl_setx) MOD map_resolution, (map_y + perl_sety) MOD map_resolution, 6) = noise(map_x, map_y) < ration_column: NEXT map_y, map_x
61.  
62. mon = _NEWIMAGE(monx, mony, 32): SCREEN mon: _FULLSCREEN: _MOUSEHIDE
63.  
64.  
65. 'MAP array
66. '0-perlin Z-data
67. '1-maptriangle calculate X
68. '2-maptriangle calculate Y
69. '3-maptriangle calculate Z
70. '4-distance from me (color)
71. '5-texture height scale
72. '6-is there a column ?
73.  
74. DO: _LIMIT 40
75.  
76.     mousex = 0: mousey = 0: WHILE _MOUSEINPUT: mousex = mousex + _MOUSEMOVEMENTX: mousey = mousey + _MOUSEMOVEMENTY: WEND
77.     'cam(3) = cam(3) + mousex * mouse_sens_xy
78.     cam(4) = cam(4) + mousey * mouse_sens_z
79.     '    cam(3) = mouse_sens * .1
80.     rotx = rotx + mousex / 150
81.     max_incli = 60: IF ABS(rotx) > (max_incli * pip180) THEN rotx = max_incli * pip180 * SGN(rotx)
82.     cam(3) = cam(3) + rotx / 30
83.  
84.     'control position
85.     kw = _KEYDOWN(119)
86.     '    kw = -1
87.     ks = _KEYDOWN(115): ka = _KEYDOWN(97): kd = _KEYDOWN(100): new_direction = (ABS(ka OR kd OR kw) OR -ABS(ks)) * position_speed
88.  
89.     deg_XY = -90 * ABS(ka) + 90 * ABS(kd)
90.  
91.     szog_xy = cam(3) + deg_XY * pip180
92.  
93.     '   szog_z = cam(4)
94.  
95.     me(0) = me(0) - SIN(szog_xy) * COS(szog_z) * new_direction
96.     me(1) = me(1) - COS(szog_xy) * COS(szog_z) * new_direction
97.     '    me(2) = me(2) + SIN(szog_z) * new_direction '* ABS(ka = 0 AND kd = 0)
98.  
99.     '    cam(2) = me(2)
100.     '    LOCATE 1, 1: PRINT cam(2)
101.  
102.  
103.  
104.     'replace the missing row in the buffer line Y
105.     DO WHILE INT(INT(me(0))) <> perl_setx: dir = -SGN(INT(me(0)) - perl_setx): perl_setx = perl_setx - dir
106.         FOR map_y = 0 TO map_resolution - 1: temp = perl_setx + ABS(dir = -1) * (map_resolution - 1): map_z = Perlin2D(temp, perl_sety + map_y)
107.             map(temp MOD map_resolution, (map_y + perl_sety) MOD map_resolution, 0) = 1 / (mapz_max - mapz_min) * (map_z - mapz_min)
108.     map(INT(temp MOD map_resolution), INT((map_y + perl_sety) MOD map_resolution), 6) = noise(temp, perl_sety + map_y) < ration_column: NEXT map_y: LOOP
109.  
110.     'replace the missing row in the buffer line X
111.     DO WHILE INT(INT(me(1))) <> perl_sety: dir = -SGN(INT(me(1)) - perl_sety): perl_sety = perl_sety - dir
112.         FOR map_x = 0 TO map_resolution - 1: temp = perl_sety + ABS(dir = -1) * (map_resolution - 1): map_z = Perlin2D(perl_setx + map_x, temp)
113.             map(INT((map_x + perl_setx) MOD map_resolution), INT(temp MOD map_resolution), 0) = 1 / (mapz_max - mapz_min) * (map_z - mapz_min)
114.     map((map_x + perl_setx) MOD map_resolution, temp MOD map_resolution, 6) = noise(perl_setx + map_x, temp) < ration_column: NEXT map_x: LOOP
115.  
116.     cosrotz = COS(cam(3)): sinrotz = SIN(cam(3)): cosrotx = COS(cam(4)): sinrotx = SIN(cam(4)) 'to rotating angles
117.  
118.     'calculating position and textures
119.     clmn(0, 0) = 0 'reset column counter
120.     FOR map_x = 0 TO map_resolution - 1: px = -map_dimxy / 2 + map_dimxy / map_resolution * (map_x - (me(0) - INT(me(0)))) - cam(0)
121.         FOR map_y = 0 TO map_resolution - 1: py = -map_dimxy / 2 + map_dimxy / map_resolution * (map_y - (me(1) - INT(me(1)))) - cam(1)
122.             read_mapx = INT((map_x + perl_setx) MOD map_resolution): read_mapy = INT((map_y + perl_sety) MOD map_resolution)
123.             map_z = map(read_mapx, read_mapy, 0): IF map_z > perlin_flat_soil_limit THEN map_z = perlin_flat_soil_limit
124.             pz2 = map_dimz * map_z - cam(2): rotate px, py, pz2
125.             map(map_x, map_y, 1) = rotating(0): map(map_x, map_y, 2) = rotating(1): map(map_x, map_y, 3) = rotating(2)
126.             map(map_x, map_y, 4) = INT(SQR(px * px + py * py) / map_shadow)
127.             IF map(map_x, map_y, 4) > text_deep - 1 THEN map(map_x, map_y, 4) = text_deep - 1
128.             map(map_x, map_y, 5) = INT(text_height_scale * map_z)
129.             IF map(map_x, map_y, 5) > text_height_scale - 1 THEN map(map_x, map_y, 5) = text_height_scale - 1
130.             IF map(map_x, map_y, 5) < 0 THEN map(map_x, map_y, 5) = 0
131.             IF map(read_mapx, read_mapy, 6) THEN clmn(clmn(0, 0) + 1, 2) = map(map_x, map_y, 4): clmn(clmn(0, 0) + 1, 0) = px: clmn(clmn(0, 0) + 1, 1) = py: clmn(0, 0) = clmn(0, 0) + 1
132.  
133.     NEXT map_y, map_x
134.  
135.     'do maptriangle from squares !
136.     FOR map_x = 0 TO map_resolution - 2: FOR map_y = 0 TO map_resolution - 2
137.             m0x = map(map_x, map_y, 1): m0y = map(map_x, map_y, 2): m0z = map(map_x, map_y, 3)
138.             m1x = map(map_x + 1, map_y, 1): m1y = map(map_x + 1, map_y, 2): m1z = map(map_x + 1, map_y, 3)
139.             m2x = map(map_x, map_y + 1, 1): m2y = map(map_x, map_y + 1, 2): m2z = map(map_x, map_y + 1, 3)
140.             m3x = map(map_x + 1, map_y + 1, 1): m3y = map(map_x + 1, map_y + 1, 2): m3z = map(map_x + 1, map_y + 1, 3)
141.             atexture = texture(map(map_x, map_y, 5), map(map_x, map_y, 4))
142.             _MAPTRIANGLE (0, 0)-(text_size - 1, 0)-(0, text_size - 1), atexture TO(m0x, m0y, m0z)-(m1x, m1y, m1z)-(m2x, m2y, m2z)
143.             _MAPTRIANGLE (0, 0)-(text_size - 1, 0)-(0, text_size - 1), atexture TO(m3x, m3y, m3z)-(m1x, m1y, m1z)-(m2x, m2y, m2z)
144.     NEXT map_y, map_x
145.  
146.     'draw columns
147.     ac_ang = 5: ac_rad = 20: column_height = 800: column_rot = column_rot + 1
148.     FOR ac = 0 TO clmn(0, 0) - 1: FOR t = 0 TO ac_ang - 1: temp = (360 / ac_ang * t + column_rot) * pip180
149.             px = clmn(ac + 1, 0) + SIN(temp) * ac_rad: py = clmn(ac + 1, 1) + COS(temp) * ac_rad
150.             pz2 = 400: rotate px, py, pz2: ac(t, 0, 0) = rotating(0): ac(t, 1, 0) = rotating(1): ac(t, 2, 0) = rotating(2)
151.         pz2 = -300: rotate px, py, pz2: ac(t, 0, 1) = rotating(0): ac(t, 1, 1) = rotating(1): ac(t, 2, 1) = rotating(2): NEXT t
152.  
153.         FOR t = 0 TO ac_ang - 1: t2 = (t + 1) MOD ac_ang
154.             m0x = ac(t, 0, 0): m0y = ac(t, 1, 0): m0z = ac(t, 2, 0): m1x = ac(t2, 0, 0): m1y = ac(t2, 1, 0): m1z = ac(t2, 2, 0)
155.             m2x = ac(t, 0, 1): m2y = ac(t, 1, 1): m2z = ac(t, 2, 1): m3x = ac(t2, 0, 1): m3y = ac(t2, 1, 1): m3z = ac(t2, 2, 1)
156.             atexture = texture(0, clmn(ac + 1, 2))
157.             _MAPTRIANGLE (0, 0)-(text_size - 1, 0)-(0, text_size - 1), atexture TO(m0x, m0y, m0z)-(m1x, m1y, m1z)-(m2x, m2y, m2z)
158.             _MAPTRIANGLE (0, 0)-(text_size - 1, 0)-(0, text_size - 1), atexture TO(m3x, m3y, m3z)-(m1x, m1y, m1z)-(m2x, m2y, m2z)
159.     NEXT t, ac
160.  
161.  
162.     _DISPLAY: CLS , _RGB32(15, 0, 0)
163.  
164.  
165. LOOP
166.  
167.  
168.  
169.  
170.  
171. FUNCTION noise (x, y) 'reading pseudo random buffer
172.     x2 = INT(x): y2 = INT(y): a = INT(ABS(x2) + ABS(y2)) MOD noise_rand_c: b = ABS(x2) MOD noise_rand_c: c = noise_rand(a) + noise_rand(b): noise = c - INT(c)
173. END FUNCTION
174.  
175. FUNCTION SmoothNoise (x, y)
176.     SELECT CASE perlin_smooth_noise_interpolation_type
177.         CASE 0: SmoothNoise = noise(x, y)
178.         CASE 1: corners = (noise(x - 1, y - 1) + noise(x + 1, y - 1) + noise(x - 1, y + 1) + noise(x + 1, y + 1)) / 16
179.             sides = (noise(x - 1, y) + noise(x + 1, y) + noise(x, y - 1) + noise(x, y + 1)) / 8: center = noise(x, y) / 4: SmoothNoise = corners + sides + center
180.     END SELECT
181. END FUNCTION
182.  
183. FUNCTION Perlin2D (x, y): total = 0: FOR t = 0 TO perlin_level - 1: frequency = 2 ^ t * perlin_multiplier: amplitude = perlin_pers ^ t
184.     total = total + InterpolatedNoise(x * frequency, y * frequency) * amplitude: NEXT t: Perlin2D = total
185. END FUNCTION
186.  
187. FUNCTION InterpolatedNoise (x, y): integer_X = INT(x): fractional_X = x - integer_X: integer_y = INT(y): fractional_Y = y - integer_y
188.     v1 = SmoothNoise(integer_X, integer_y): v2 = SmoothNoise(integer_X + 1, integer_y)
189.     v3 = SmoothNoise(integer_X, integer_y + 1): v4 = SmoothNoise(integer_X + 1, integer_y + 1)
190.     i1 = Interpolate(v1, v2, fractional_X): i2 = Interpolate(v3, v4, fractional_X): InterpolatedNoise = Interpolate(i1, i2, fractional_Y)
191. END FUNCTION
192.  
193. FUNCTION Interpolate (a, b, x): Interpolate = a * (1 - x) + b * x: END FUNCTION
194.  
195. SUB rotate (px, py, pz2)
196.     px3 = px * cosrotz - py * sinrotz: py2 = px * sinrotz + py * cosrotz: py3 = py2 * cosrotx - pz2 * sinrotx: pz3 = py2 * sinrotx + pz2 * cosrotx
197.     px4 = (px3 * COS(rotx)) - (py3 * SIN(rotx)): py4 = (px3 * SIN(rotx)) + (py3 * COS(rotx))
198.     rotating(0) = -px4 * map_zoom_xy: rotating(1) = -py4 * map_zoom_xy: rotating(2) = -pz3 * map_zoom_distance
199. END SUB
200.  
201.  

[bplus]

OK I ran this for 1 hr 20 mins while lunched, read paper did puzzles and walked dog and no problems when I came back to it (2nd version).

[johnno56]

First. Nice improvements. Manual control. Sweet. The rotating "tree trunks" were a surprise... Cool.

This may be a silly question (bplus... don't even think about it... lol) but is it just me that has difficulty exiting the program? I can alt+tab to select QB64 IDE and terminate it, but the demo is still running... I should have read the listing more carefully and inserted an "escape clause"... Moo ha ha...

Still a very cool demo...

[Dav]

What a great demo.  Pretty amazing work!  Small code but does so much.  Well done.

- Dav

[_vince]

nice work, looks great!  I am interested in this and have been modifying your program, also have some game ideas

[MasterGy]

Thanks !
Thanks Dav for trying it and the others too!
I’m very curious about your ideas, Vince, and I’m curious about what you modified on it. Where can I see it?
I changed a lot of it. I worked a lot of hours to make it enjoyable to manage. The hills just throw it and slow it down, and it bounces off the columns. When steering, you can feel that it is not just the direction that matters, but the thrust behind us. I wonder what you say about management, I had a lot of work in it. It feels a bit like drifting. You go in one direction at high speed and suddenly I put the force in another direction, but the inertia and thrust play enjoyably. In flight, W is the gas, S is the brake, and all directions can be achieved with the sensual movement of the mouse. If the speed is 0, it will land on the ground by pressing S continuously.
Press F to toggle between flight and walk mode.
Currently, the goal is to get to the “target”.
There’s a lot of computation and I’ve managed to simplify a few things, so far, my otherwise slow machine is comfortable with 40FPS.
I hope that next time a heavenly war can be fought. I thought a lot about how to entrust the control of an enemy machine to the program.

download (2Mb)
https://drive.google.com/file/d/1DWwX1ayPSQLx8rjgWZXhJvx_uCLN9auk/view?usp=sharing

[Ashish]

@MasterGy
Great Job!! Looks similar to a "textured" version of this - https://www.qb64.org/forum/index.php?topic=1139.msg103744#msg103744 which has heavily used  Perlin noise.
If you are interested, you can read this article - https://www.redblobgames.com/maps/terrain-from-noise/

[MasterGy]

Ashish, I'm glad you brought it to my attention! I used to watch Openworld and I really liked the graphics. Then I saw that there was opengl in it and I was scared and I didn’t even look further. :)
Now I understand that you also made the map with perlin noise, and even the textures, only with a different approach! I really like you to assign textures to height values. This is how I do the placement of the columns. All I miss is that if you’re already making the map with Perlin noise anyway, why haven’t you made it to infinite expanse. Originally, I was looking for this solution because I wanted an infinite map that would be generated in real time.
What you wrote as a link demonstrates this Perlin noise very well! I like the way the website shows how the map changes when we change a value.
And your opengl solution is a special congratulations, I still don't understand the whole thing :(

[Dav]

This is looking great, @MasterGy!  Kind of reminds me of one of my favorite DOS games, Terminal Velocity.  I'll be following this game development!

- Dav

[MasterGy]

I promised that I would show the sequel when there was a heavenly battle. This, in turn, required a solution that would take the plane to its destination as efficiently as possible. I have to show it because it’s very spectacular! Everything is the same as before. Press F to toggle between walking, flying, and autopilot. The autopilot goes to the target in a superhuman way, elegantly bypassing the columns at maximum speed.

Code: QB64: [Select]

1. CONST pip180 = 3.141592 / 180
2. CONST perlin_pers = 7 'amplitudo
3. CONST perlin_level = 5 'level of resolution wave
4. 'CONST perlin_multiplier = 0.006
5. CONST perlin_multiplier = 0.008
6. CONST perlin_smooth_noise_interpolation_type = 0 '0-none '1-standard
7. CONST perlin_flat_soil_limit_low = .8 'high value is low area (negate)
8. CONST perlin_flat_soil_limit_high = .1
9. CONST flight_array_size = 20: DIM flight_save(flight_array_size - 1)
10. DIM SHARED rotating(2), cosrotz, sinrotz, cosrotx, sinrotx, map_zoom_xy, map_zoom_distance, see_point, flight(19, flight_array_size), me(19), rotx_bevel
11. DIM SHARED mouse_sens_xy, mouse_sens_z, position_speed, clmn_rad
12. DIM SHARED map_x, map_y, perl_setx, perl_sety, map_resolution, map_z, map_dimz, map_dimxy, flight_high_max
13. monx = 1366: mony = monx / 16 * 9 'screen size
14.  
15. RANDOMIZE TIMER
16. win_marg = monx / 25
17.  
18. flight_high_max = -70
19.  
20. text_deep = 30 'make deep shadow pictures
21.  
22.  
23. map_center = 100000 'start of perlin-map
24. text_field$ = "need\field2.bmp"
25. text_field$ = "need\field1.bmp" 'field map texture
26. text_multi = 40 'zoom to texture
27.  
28. position_speed = .9
29. map_resolution = 60 'see map resolution
30. DIM SHARED map_resper2: map_resper2 = INT(map_resolution / 2)
31.  
32.  
33. DIM SHARED max_incli: max_incli = 65 'flight max rotation
34.  
35. map_dimxy = 1200 * .8
36. map_dimz = 250 * .8
37.  
38. map_zoom_xy = 6 * .2
39. map_zoom_distance = 15 * .2
40.  
41.  
42. mouse_sens_xy = 1.6
43. mouse_sens_z = 1
44.  
45. me(6) = pip180 * 270
46.  
47. me(0) = map_center: me(1) = map_center
48.  
49. control_type$(0) = "walking": control_type$(1) = "flying": control_type$(2) = "flying/autopilot"
50. DIM color_temp(9) AS _INTEGER64
51. targ_text = _LOADIMAGE("need\target.jpg", 33)
52.  
53.  
54.  
55.  
56.  
57. DIM clmn(1999, 5): clmn_rad = 30: ration_column = 0.0003: ration_column = 0.002
58.  
59.  
60. 'install radar
61. rad_ts = 600: rad_size = monx / 5.5: r1 = rad_ts / 2 * .96: r2 = rad_ts / 2 * .80
62. temp = _NEWIMAGE(rad_ts, rad_ts, 32): _DEST temp: CLS , 0: color_temp(0) = _RGB32(200, 200, 200)
63. CIRCLE (rad_ts / 2, rad_ts / 2), r1, color_temp(0): CIRCLE (rad_ts / 2, rad_ts / 2), r2, color_temp(0)
64. PAINT (rad_ts / 2 + (r1 + r2) / 2, rad_ts / 2), color_temp(0), color_temp(0): PAINT (rad_ts / 2, rad_ts / 2), _RGBA(100, 100, 255, 150), color_temp(0)
65. rad_text = _COPYIMAGE(temp, 33): _FREEIMAGE temp
66.  
67. 'radar object
68. rado_ts = 150: temp = _NEWIMAGE(rado_ts, rado_ts, 32): _DEST temp: CLS , 0: color_temp(0) = _RGB32(180, 0, 0)
69. CIRCLE (rado_ts / 2, rado_ts / 2), rado_ts / 2 * .98, color_temp(0)
70. PAINT (rado_ts / 2, rado_ts / 2), color_temp(0), color_temp(0): rado_text = _COPYIMAGE(temp, 33): _FREEIMAGE temp
71.  
72. 'make sky
73. da = 8: db = 8: sky_points = da * db: DIM sky_points(sky_points - 1, 9), sq(sky_points - 1, 7): sky_image = _LOADIMAGE("need\sky.jpg", 33)
74. FOR da2 = 0 TO da - 1: dega = 360 / (da - 1) * da2 * pip180: FOR db2 = 0 TO db - 1: degb = 180 / (db - 1) * db2 * pip180: ss = 1400
75. ap = da2 * db + db2: sky_points(ap, 0) = SIN(degb) * COS(dega) * ss: sky_points(ap, 1) = SIN(degb) * SIN(dega) * ss: sky_points(ap, 2) = COS(degb) * ss: NEXT db2, da2
76. FOR da2 = 0 TO da - 2: FOR db2 = 0 TO db - 2: sqa = da2 * db + db2: sq(sqa, 0) = sqa: sq(sqa, 1) = sq(sqa, 0) + 1: sq(sqa, 2) = sq(sqa, 0) + db: sq(sqa, 3) = sq(sqa, 2) + 1
77.         sq(sqa, 4) = _WIDTH(sky_image) - (_WIDTH(sky_image) / (da - 1) * da2) - 1: sq(sqa, 5) = _WIDTH(sky_image) - (_WIDTH(sky_image) / (da - 1) * (da2 + 1)) - 1
78. sq(sqa, 6) = INT(_HEIGHT(sky_image) / (db - 1) * db2): sq(sqa, 7) = INT(_HEIGHT(sky_image) / (db - 1) * (db2 + 1)): NEXT db2, da2
79.  
80. 'create pseudo random buffer
81. DIM SHARED noise_rand_c: noise_rand_c = 10000: DIM SHARED noise_rand(noise_rand_c - 1): FOR t = 0 TO noise_rand_c - 1: noise_rand(t) = RND(1): NEXT t
82.  
83. 'create texture
84. DIM temp_color(9) AS _INTEGER64: DIM SHARED mapz_min, mapz_max: DIM texture(text_deep - 1, 5): temp = _LOADIMAGE(text_field$, 32)
85. FOR t = 0 TO text_deep - 1: dark = (1 - (1 / (text_deep - 1) * t)): transp = dark * 2000: IF transp > 255 THEN transp = 255
86.     temp2 = _NEWIMAGE(_WIDTH(temp), _HEIGHT(temp), 32): _SOURCE temp: _DEST temp2: CLS 0, _RGBA32(0, 0, 0, transp)
87.     FOR tx = 0 TO _WIDTH(temp) - 1: FOR ty = 0 TO _HEIGHT(temp) - 1: temp_color(0) = POINT(tx, ty)
88.             c1 = _RED32(temp_color(0)) * dark: c2 = _GREEN32(temp_color(0)) * dark: c3 = _BLUE32(temp_color(0)) * dark
89.     PSET (tx, ty), _RGBA32(c1, c2, c3, transp): NEXT ty, tx: texture(t, 0) = _COPYIMAGE(temp2, 33): _FREEIMAGE temp2
90. NEXT t
91.  
92. 'find min/max Z
93. DIM SHARED mapz_multiplier: mapz_min = 9999999: mapz_max = -mapz_min: FOR t = 0 TO 2999: c = Perlin2D_original(5000 * RND(1), 5000 * RND(1))
94.     mapz_min = c * ABS(mapz_min > c) + mapz_min * (ABS(mapz_min > c) XOR 1): mapz_max = c * ABS(mapz_max < c) + mapz_max * (ABS(mapz_max < c) XOR 1)
95. NEXT t: mapz_multiplier = 1 / (mapz_max - mapz_min)
96.  
97.  
98. 'fill map-buffer
99. DIM SHARED map(map_resolution - 1, map_resolution - 1, 19): perl_setx = INT(me(0)): perl_sety = INT(me(1))
100. FOR map_x = 0 TO map_resolution - 1: FOR map_y = 0 TO map_resolution - 1
101.         map((map_x + perl_setx) MOD map_resolution, (map_y + perl_sety) MOD map_resolution, 0) = perlin2d(perl_setx + map_x, perl_sety + map_y)
102. map((map_x + perl_setx) MOD map_resolution, (map_y + perl_sety) MOD map_resolution, 6) = noise(map_x, map_y) < ration_column: NEXT map_y, map_x
103.  
104. 'calculating shadows
105. FOR map_x = 0 TO map_resolution - 1: FOR map_y = 0 TO map_resolution - 1
106.         dis = INT((text_deep - 1) / map_resper2 * SQR((map_x - map_resper2) ^ 2 + (map_y - map_resper2) ^ 2)): IF dis < 0 OR dis > text_deep - 1 THEN _CONTINUE
107. map(map_x, map_y, 9) = 1: map(map_x, map_y, 8) = dis: NEXT map_y, map_x
108.  
109.  
110.  
111.  
112. 'MAP array
113. '0-perlin Z-data
114. '1-maptriangle calculate X
115. '2-maptriangle calculate Y
116. '3-maptriangle calculate Z
117. '4-distance from me (color)
118. '5-texture height scale
119. '6-is there a column ?
120. '7-is the point visible?
121. '8-shadow-table
122. '9-not used area signal
123.  
124. 'ME array
125. '0-me X location
126. '1-me Y location
127. '2-me Z location
128. '3-vector_x
129. '4-vector_y
130. '5-me XY angle CAM3
131. '6-me XY CAM4
132. '7-me kanyarodas merteke
133.  
134. 'FLIGHT array
135. '0-active
136. '1-rotx - kanyarodas merteke
137. '2-rotx_bevel
138. '3-location X
139. '4-location Y
140. '5-location Z
141. '6-vector X
142. '7-vector Y
143. '8-W gas
144. '9-S brake
145. '10-mouseX
146. '11-mouseZ
147. '12 cam4  /me6
148. '13 cam3 /me5
149. '15 actual speed
150. '16 column (XY) crash /last step
151. '17 land (Z) crash /last step
152.  
153.  
154. auto_precalc = 28
155. auto_mouse_scale = 140
156. auto_try_resolution = 22
157. DIM flight_auto(auto_try_resolution - 1, flight_array_size - 1)
158. DIM auto_efficiency(auto_try_resolution - 1, 5) '0- 1 if crash          '1- target-flight distance
159.  
160.  
161.  
162.  
163. GOSUB new_target
164. mon = _NEWIMAGE(monx, mony, 32): SCREEN mon: _FULLSCREEN _SQUAREPIXELS , _SMOOTH: _MOUSEHIDE: _DISPLAYORDER _HARDWARE , _SOFTWARE
165. control_type = 0
166.  
167.  
168.  
169.  
170. _PRINTSTRING (0, 0), "moving:mouse+WASD jump:space walking/fly: F"
171. DO: _LIMIT 40
172.     IF _KEYDOWN(27) THEN SYSTEM
173.     kf = _KEYDOWN(102): k_space = _KEYDOWN(32)
174.     kw = _KEYDOWN(119): ks = _KEYDOWN(115): mousex = 0: mousey = 0: WHILE _MOUSEINPUT: mousex = mousex + _MOUSEMOVEMENTX: mousey = mousey + _MOUSEMOVEMENTY: WEND
175.  
176.     IF ut_kf = 0 AND kf THEN control_type = (control_type + 1) MOD 3
177.     ut_kf = kf
178.  
179.     walk_speed = .18
180.  
181.  
182.     _PRINTSTRING (500, 0), "control type : " + control_type$(control_type) + "           "
183.  
184.     SELECT CASE control_type '(0-walking 1-flying)
185.         CASE 0
186.             me(3) = 0: me(4) = 0: me(7) = 0: flight(0, 1) = 0
187.             me(5) = me(5) + mousex * mouse_sens_xy / 200: me(6) = me(6) + mousey * mouse_sens_xy / 200
188.             kw = _KEYDOWN(119): ks = _KEYDOWN(115): ka = _KEYDOWN(97): kd = _KEYDOWN(100)
189.             go = ABS(ka OR kd OR kw OR ks): direction = (-90 * ABS(ka) + 90 * ABS(kd) + 180 * ABS(ks)) * go * pip180
190.             go_x = -(SIN(direction + me(5)) * walk_speed) * go: go_y = -(COS(direction + me(5)) * walk_speed) * go
191.             IF crash_to_column(me(0) + go_x, me(1) + go_y) THEN go_x = 0: go_y = 0
192.  
193.             me(0) = me(0) + go_x: me(1) = me(1) + go_y
194.  
195.             IF jump_cf = 0 THEN
196.                 IF (map_deep(me(0), me(1)) - .18) * map_dimz > me(2) THEN
197.                     walk_c_down = walk_c_down + .2: me(2) = me(2) + walk_c_down: free_jump = 0
198.                 ELSE walk_c_down = 0: me(2) = (map_deep(me(0), me(1)) - .18) * map_dimz: free_jump = 1
199.                 END IF
200.             END IF
201.  
202.             IF k_space AND free_jump THEN jump_cf = 15
203.             rotx_bevel = rotx_bevel * .95: me(2) = me(2) - jump_cf / 4: jump_cf = jump_cf - 1: IF jump_cf < 0 THEN jump_cf = 0
204.  
205.  
206.         CASE 1
207.             'contact between me and flight0
208.             flight(0, 2) = rotx_bevel: flight(0, 8) = ABS(kw): flight(0, 9) = ABS(ks): flight(0, 10) = mousex * mouse_sens_xy: flight(0, 11) = mousey * mouse_sens_z
209.             flight(0, 3) = me(0): flight(0, 4) = me(1): flight(0, 5) = me(2): flight(0, 6) = me(3): flight(0, 7) = me(4): flight(0, 12) = me(6): flight(0, 13) = me(5)
210.             flight_control 0, 0
211.             me(0) = flight(0, 3): me(1) = flight(0, 4): me(2) = flight(0, 5): me(3) = flight(0, 6): me(4) = flight(0, 7): me(5) = flight(0, 13): me(6) = flight(0, 12)
212.             rotx_bevel = flight(0, 2)
213.         CASE 2 'autopilot
214.  
215.             flight(0, 2) = rotx_bevel: flight(0, 8) = ABS(kw): flight(0, 9) = ABS(ks): flight(0, 10) = mousex * mouse_sens_xy: flight(0, 11) = mousey * mouse_sens_z
216.             flight(0, 3) = me(0): flight(0, 4) = me(1): flight(0, 5) = me(2): flight(0, 6) = me(3): flight(0, 7) = me(4): flight(0, 12) = me(6): flight(0, 13) = me(5)
217.  
218.             FOR t = 0 TO flight_array_size - 1: flight_save(t) = flight(0, t): NEXT t
219.  
220.             'calculation of options
221.             FOR a_try = 0 TO auto_try_resolution - 1: auto_efficiency(a_try, 0) = 0
222.                 FOR t = 0 TO flight_array_size - 1: flight(0, t) = flight_save(t): NEXT t
223.                 flight(0, 8) = 1: flight(0, 10) = -auto_mouse_scale + (auto_mouse_scale * 2) / (auto_try_resolution - 1) * a_try
224.                 FOR a_pre = 1 TO auto_precalc: flight_control 0, 1: IF flight(0, 16) THEN EXIT FOR
225.                 NEXT a_pre: IF flight(0, 16) THEN auto_efficiency(a_try, 0) = 1: _CONTINUE
226.             auto_efficiency(a_try, 1) = ABS((targ(0) - flight(0, 3)) ^ 2 + (targ(1) - flight(0, 4)) ^ 2): NEXT a_try
227.  
228.             'choosing the most effective option
229.             adis = 999999: auto_opt = -1: FOR t = 0 TO auto_try_resolution - 1
230.                 IF auto_efficiency(t, 0) = 0 AND auto_efficiency(t, 1) < adis THEN adis = auto_efficiency(t, 1): auto_opt = t
231.             NEXT t: IF auto_opt = -1 THEN auto_opt = INT(auto_try_resolution * RND(1)): auto_rand = auto_rand + 1
232.  
233.             FOR t = 0 TO flight_array_size - 1: flight(0, t) = flight_save(t): NEXT t
234.             flight(0, 8) = 1: flight(0, 10) = -auto_mouse_scale + (auto_mouse_scale * 2) / (auto_try_resolution - 1) * auto_opt: flight_control 0, 0
235.             me(0) = flight(0, 3): me(1) = flight(0, 4): me(2) = flight(0, 5): me(3) = flight(0, 6): me(4) = flight(0, 7): me(5) = flight(0, 13): me(6) = flight(0, 12)
236.             rotx_bevel = flight(0, 2)
237.  
238.     END SELECT
239.  
240.  
241.  
242.     'replace the missing row in the buffer line Y
243.     DO WHILE INT(INT(me(0))) <> perl_setx: dir = -SGN(INT(me(0)) - perl_setx): perl_setx = perl_setx - dir
244.         temp = INT(perl_setx + ABS(dir = -1) * (map_resolution - 1)): temp_m = INT(temp MOD map_resolution)
245.         FOR map_y = 0 TO map_resolution - 1: map(temp_m, INT(map_y + perl_sety) MOD map_resolution, 0) = perlin2d(temp, perl_sety + map_y)
246.     map(temp_m, INT((map_y + perl_sety) MOD map_resolution), 6) = noise(temp, perl_sety + map_y) < ration_column: NEXT map_y: LOOP
247.  
248.     'replace the missing row in the buffer line X
249.     DO WHILE INT(INT(me(1))) <> perl_sety: dir = -SGN(INT(me(1)) - perl_sety): perl_sety = perl_sety - dir
250.         temp = perl_sety + ABS(dir = -1) * (map_resolution - 1): temp_m = INT(temp MOD map_resolution)
251.         FOR map_x = 0 TO map_resolution - 1: map(INT((map_x + perl_setx) MOD map_resolution), temp_m, 0) = perlin2d(perl_setx + map_x, temp)
252.     map(INT(map_x + perl_setx) MOD map_resolution, temp_m, 6) = noise(perl_setx + map_x, temp) < ration_column: NEXT map_x: LOOP
253.  
254.     cosrotz = COS(me(5)): sinrotz = SIN(me(5)): cosrotx = COS(me(6)): sinrotx = SIN(me(6)) 'to rotating angles
255.  
256.  
257.     'draw sky
258.     FOR actual_point = 0 TO sky_points - 1: rotate sky_points(actual_point, 0), sky_points(actual_point, 1), sky_points(actual_point, 2), 0
259.     sky_points(actual_point, 4) = rotating(0): sky_points(actual_point, 5) = rotating(1): sky_points(actual_point, 6) = rotating(2): NEXT actual_point
260.     FOR asq = 0 TO sky_points - 1
261.         wx0 = sky_points(sq(asq, 0), 4): wy0 = sky_points(sq(asq, 0), 5): wz0 = sky_points(sq(asq, 0), 6)
262.         wx1 = sky_points(sq(asq, 1), 4): wy1 = sky_points(sq(asq, 1), 5): wz1 = sky_points(sq(asq, 1), 6)
263.         wx2 = sky_points(sq(asq, 2), 4): wy2 = sky_points(sq(asq, 2), 5): wz2 = sky_points(sq(asq, 2), 6)
264.         wx3 = sky_points(sq(asq, 3), 4): wy3 = sky_points(sq(asq, 3), 5): wz3 = sky_points(sq(asq, 3), 6)
265.         sy0 = sq(asq, 6): sx0 = sq(asq, 4): sy1 = sq(asq, 7): sx1 = sq(asq, 4): sy2 = sq(asq, 6): sx2 = sq(asq, 5): sy3 = sq(asq, 7): sx3 = sq(asq, 5)
266.         _MAPTRIANGLE (sx0, sy0)-(sx1, sy1)-(sx2, sy2), sky_image TO(wx0, wy0, wz0)-(wx1, wy1, wz1)-(wx2, wy2, wz2), , _SMOOTH
267.         _MAPTRIANGLE (sx3, sy3)-(sx1, sy1)-(sx2, sy2), sky_image TO(wx3, wy3, wz3)-(wx1, wy1, wz1)-(wx2, wy2, wz2), , _SMOOTH
268.     NEXT asq
269.  
270.     'calculating position and textures
271.     clmn(0, 0) = 0 'reset column counter
272.     FOR map_x = 0 TO map_resolution - 1: px = -map_dimxy / 2 + map_dimxy / map_resolution * (map_x - (me(0) - INT(me(0))))
273.         FOR map_y = 0 TO map_resolution - 1: IF map(map_x, map_y, 9) = 0 THEN _CONTINUE
274.             py = -map_dimxy / 2 + map_dimxy / map_resolution * (map_y - (me(1) - INT(me(1))))
275.             read_mapx = INT((map_x + perl_setx) MOD map_resolution): read_mapy = INT((map_y + perl_sety) MOD map_resolution): map_z = map(read_mapx, read_mapy, 0)
276.             map(map_x, map_y, 7) = 0: pz2 = map_dimz * map_z - me(2): rotate px, py, pz2, 1: IF see_point = 0 THEN _CONTINUE
277.             map(map_x, map_y, 1) = rotating(0): map(map_x, map_y, 2) = rotating(1): map(map_x, map_y, 3) = rotating(2): map(map_x, map_y, 7) = 1
278.             map(map_x, map_y, 5) = INT(text_height_scale * map_z): IF map(map_x, map_y, 5) > text_height_scale - 1 THEN map(map_x, map_y, 5) = text_height_scale - 1
279.             IF map(map_x, map_y, 5) < 0 THEN map(map_x, map_y, 5) = 0
280.             IF map(read_mapx, read_mapy, 6) THEN clmn(clmn(0, 0) + 1, 2) = map(map_x, map_y, 8): clmn(clmn(0, 0) + 1, 0) = px: clmn(clmn(0, 0) + 1, 1) = py: clmn(0, 0) = clmn(0, 0) + 1
281.     NEXT map_y, map_x
282.  
283.  
284.     'do maptriangle from squares !
285.     FOR map_x = 0 TO map_resolution - 2: FOR map_y = 0 TO map_resolution - 2
286.             IF (map(map_x, map_y, 7) AND map(map_x + 1, map_y, 7) AND map(map_x, map_y + 1, 7) AND map(map_x + 1, map_y + 1, 7)) = 0 THEN _CONTINUE
287.             m0x = map(map_x, map_y, 1): m0y = map(map_x, map_y, 2): m0z = map(map_x, map_y, 3)
288.             m1x = map(map_x + 1, map_y, 1): m1y = map(map_x + 1, map_y, 2): m1z = map(map_x + 1, map_y, 3)
289.             m2x = map(map_x, map_y + 1, 1): m2y = map(map_x, map_y + 1, 2): m2z = map(map_x, map_y + 1, 3)
290.             m3x = map(map_x + 1, map_y + 1, 1): m3y = map(map_x + 1, map_y + 1, 2): m3z = map(map_x + 1, map_y + 1, 3)
291.             atexture = texture(map(map_x, map_y, 8), 0)
292.  
293.             sx1 = (perl_setx + map_x) * text_multi: sy1 = (perl_sety + map_y) * text_multi: sx2 = sx1 + text_multi: sy2 = sy1 + text_multi
294.  
295.             _MAPTRIANGLE (sx1, sy1)-(sx2, sy1)-(sx1, sy2), atexture TO(m0x, m0y, m0z)-(m1x, m1y, m1z)-(m2x, m2y, m2z)
296.             _MAPTRIANGLE (sx2, sy2)-(sx2, sy1)-(sx1, sy2), atexture TO(m3x, m3y, m3z)-(m1x, m1y, m1z)-(m2x, m2y, m2z)
297.     NEXT map_y, map_x
298.  
299.     'draw columns
300.     ac_ang = 11: column_height = 800: column_rot = column_rot + 1: FOR ac = 0 TO clmn(0, 0) - 1: FOR t = 0 TO ac_ang - 1
301.             temp = (360 / ac_ang * t + column_rot) * pip180: px = clmn(ac + 1, 0) + SIN(temp) * clmn_rad: py = clmn(ac + 1, 1) + COS(temp) * clmn_rad
302.             pz2 = 300 - me(2): rotate px, py, pz2, 0: ac(t, 0, 0) = rotating(0): ac(t, 1, 0) = rotating(1): ac(t, 2, 0) = rotating(2)
303.         pz2 = -150 - me(2): rotate px, py, pz2, 0: ac(t, 0, 1) = rotating(0): ac(t, 1, 1) = rotating(1): ac(t, 2, 1) = rotating(2): NEXT t
304.  
305.         FOR t = 0 TO ac_ang - 1: t2 = (t + 1) MOD ac_ang
306.             m0x = ac(t, 0, 0): m0y = ac(t, 1, 0): m0z = ac(t, 2, 0): m1x = ac(t2, 0, 0): m1y = ac(t2, 1, 0): m1z = ac(t2, 2, 0)
307.             m2x = ac(t, 0, 1): m2y = ac(t, 1, 1): m2z = ac(t, 2, 1): m3x = ac(t2, 0, 1): m3y = ac(t2, 1, 1): m3z = ac(t2, 2, 1)
308.             atexture = texture(clmn(ac + 1, 2), 0): dd = 4: pm = 40 * dd: pm2 = 310 * dd: sx1 = pm * t: sy1 = pm2: sx2 = sx1 + pm: sy2 = sy1 + pm2
309.             _MAPTRIANGLE (sx1, sy1)-(sx2, sy1)-(sx1, sy2), atexture TO(m0x, m0y, m0z)-(m1x, m1y, m1z)-(m2x, m2y, m2z)
310.             _MAPTRIANGLE (sx2, sy2)-(sx2, sy1)-(sx1, sy2), atexture TO(m3x, m3y, m3z)-(m1x, m1y, m1z)-(m2x, m2y, m2z)
311.  
312.     NEXT t, ac
313.  
314.  
315.  
316.     'target draw
317.     dis = SQR((targ(0) - me(0)) ^ 2 + (targ(1) - me(1)) ^ 2): _PRINTSTRING (0, 20), "target distance:" + STR$(dis)
318.  
319.     IF dis < map_resolution / 2 * 1.5 THEN
320.         ee = map_dimxy / map_resolution: targ_size = 50: rotate (targ(0) - me(0)) * ee, (targ(1) - me(1)) * ee, (targ(2) - me(2)) * 1, 0
321.         m0x = rotating(0) - targ_size: m0y = rotating(1) - targ_size: m0z = rotating(2): m1x = rotating(0) + targ_size: m1y = rotating(1) - targ_size: m1z = rotating(2)
322.         m2x = rotating(0) - targ_size: m2y = rotating(1) + targ_size: m2z = rotating(2): m3x = rotating(0) + targ_size: m3y = rotating(1) + targ_size: m3z = rotating(2)
323.         _MAPTRIANGLE (0, 500 - 1)-(500 - 1, 500 - 1)-(0, 0), targ_text TO(m0x, m0y, m0z)-(m1x, m1y, m1z)-(m2x, m2y, m2z)
324.         _MAPTRIANGLE (500 - 1, 0)-(500 - 1, 500 - 1)-(0, 0), targ_text TO(m3x, m3y, m3z)-(m1x, m1y, m1z)-(m2x, m2y, m2z)
325.     END IF
326.  
327.     IF new_target_signal AND INT(RND(1) * 2) THEN _PRINTSTRING (600, 20), "NEW TARGET !!!!!"
328.     new_target_signal = new_target_signal - 1: IF new_target_signal < 0 THEN new_target_signal = 0
329.  
330.  
331.  
332.     IF dis < 4 OR bill$ = "l" THEN GOSUB new_target
333.  
334.  
335.     'draw radar
336.     m0x = win_marg: m1x = m0x + rad_size: m0y = mony - rad_size - win_marg: m1y = m0y + rad_size
337.     _MAPTRIANGLE (0, 0)-(rad_ts - 1, 0)-(0, rad_ts - 1), rad_text TO(m0x, m0y)-(m1x, m0y)-(m0x, m1y), , _SMOOTH
338.     _MAPTRIANGLE (rad_ts - 1, rad_ts - 1)-(rad_ts - 1, 0)-(0, rad_ts - 1), rad_text TO(m1x, m1y)-(m1x, m0y)-(m0x, m1y), , _SMOOTH
339.  
340.     'draw objects in radar
341.     rad_obj_multi = .8: rad_obj_limit = rad_size / 2 * .9: angle = me(5) + (degree(targ(1) - me(1), targ(0) - me(0)) - 90) * pip180
342.     dis = SQR((targ(0) - me(0)) ^ 2 + (targ(1) - me(1)) ^ 2) * rad_obj_multi: IF dis > rad_obj_limit THEN dis = rad_obj_limit
343.     x = SIN(angle) * dis + win_marg + rad_size / 2: y = mony - win_marg - rad_size / 2 + COS(angle) * dis: ros = 10
344.     _MAPTRIANGLE (0, 0)-(rado_ts - 1, 0)-(0, rado_ts - 1), rado_text TO(x - ros, y - ros)-(x + ros, y - ros)-(x - ros, y + ros), , _SMOOTH
345.     _MAPTRIANGLE (rado_ts - 1, rado_ts - 1)-(rado_ts - 1, 0)-(0, rado_ts - 1), rado_text TO(x + ros, y + ros)-(x + ros, y - ros)-(x - ros, y + ros), , _SMOOTH
346.  
347.  
348.  
349.  
350.  
351.     _DISPLAY: _PRINTSTRING (0, 20), SPACE$(160)
352.  
353.  
354. LOOP
355.  
356. new_target: l = 200: targ(0) = map_center + l / 2 + l * RND(1): targ(1) = map_center + l / 2 + l * RND(1): targ(2) = (map_deep(targ(0), targ(1)) - .1 - targ_size / 2)
357. new_target_signal = 160: RETURN
358.  
359.  
360.  
361.  
362. FUNCTION noise (x, y) 'reading pseudo random buffer
363.     x2 = INT(x): y2 = INT(y): a = INT(ABS(x2) + ABS(y2)) MOD noise_rand_c: b = ABS(x2) MOD noise_rand_c: c = noise_rand(a) + noise_rand(b): noise = c - INT(c)
364. END FUNCTION
365.  
366. FUNCTION SmoothNoise (x, y)
367.     SELECT CASE perlin_smooth_noise_interpolation_type
368.         CASE 0: SmoothNoise = noise(x, y)
369.         CASE 1: corners = (noise(x - 1, y - 1) + noise(x + 1, y - 1) + noise(x - 1, y + 1) + noise(x + 1, y + 1)) / 16
370.             sides = (noise(x - 1, y) + noise(x + 1, y) + noise(x, y - 1) + noise(x, y + 1)) / 8: center = noise(x, y) / 4: SmoothNoise = corners + sides + center
371. END SELECT: END FUNCTION
372.  
373. FUNCTION Perlin2D_original (x, y): total = 0: FOR t = 0 TO perlin_level - 1: frequency = 2 ^ t * perlin_multiplier: amplitude = perlin_pers ^ t
374. total = total + InterpolatedNoise(x * frequency, y * frequency) * amplitude: NEXT t: Perlin2D_original = total: END FUNCTION
375.  
376. FUNCTION perlin2d (x, y): perlin2dx = mapz_multiplier * (Perlin2D_original(x, y) - mapz_min)
377.     IF perlin2dx > perlin_flat_soil_limit_low THEN perlin2dx = perlin_flat_soil_limit_low
378.     IF perlin2dx < perlin_flat_soil_limit_high THEN perlin2dx = perlin_flat_soil_limit_high
379.     perlin2d = perlin2dx
380. END FUNCTION
381.  
382. FUNCTION InterpolatedNoise (x, y): integer_X = INT(x): fractional_X = x - integer_X: integer_y = INT(y): fractional_Y = y - integer_y
383.     v1 = SmoothNoise(integer_X, integer_y): v2 = SmoothNoise(integer_X + 1, integer_y)
384.     v3 = SmoothNoise(integer_X, integer_y + 1): v4 = SmoothNoise(integer_X + 1, integer_y + 1)
385. i1 = Interpolate(v1, v2, fractional_X): i2 = Interpolate(v3, v4, fractional_X): InterpolatedNoise = Interpolate(i1, i2, fractional_Y): END FUNCTION
386.  
387. FUNCTION Interpolate (a, b, x): Interpolate = a * (1 - x) + b * x: END FUNCTION
388.  
389. SUB rotate (px, py, pz2, see_analysis)
390.     px3 = px * cosrotz - py * sinrotz: py2 = px * sinrotz + py * cosrotz: py3 = py2 * cosrotx - pz2 * sinrotx: pz3 = -(py2 * sinrotx + pz2 * cosrotx)
391.     see_point = pz3 < 200: IF see_point = 0 AND see_analysis THEN EXIT SUB
392.     rotating(0) = -(px3 * COS(rotx_bevel) - py3 * SIN(rotx_bevel)) * map_zoom_xy: rotating(1) = -(px3 * SIN(rotx_bevel) + py3 * COS(rotx_bevel)) * map_zoom_xy
393. rotating(2) = pz3 * map_zoom_distance: END SUB
394.  
395.  
396. SUB flight_control (af, test)
397.  
398.     rotx = flight(af, 1): me0 = flight(af, 3): me1 = flight(af, 4): me2 = flight(af, 5): me3 = flight(af, 6): me4 = flight(af, 7): kw = ABS(flight(af, 8))
399.     ks = ABS(flight(af, 9)): mousey = flight(af, 11): mousex = flight(af, 10): me6 = flight(af, 12): me5 = flight(af, 13): rotx_bevel3 = flight(af, 2)
400.  
401.     rotx = (rotx + mousex * .005) * .95: IF ABS(rotx) > (max_incli * pip180) THEN rotx = max_incli * pip180 * SGN(rotx)
402.  
403.     actual_speed = 1 / position_speed * SQR(me3 * me3 + me4 * me4): rotx_bevel2 = rotx * actual_speed * kw
404.     rotx_bevel3 = rotx_bevel3 + (rotx_bevel2 - rotx_bevel3) * .05: me5 = me5 + rotx / 20: me6 = me6 + mousey / 100: lim_fy = 65: temp = me6 / pip180
405.     IF temp > 270 + lim_fy THEN temp = 270 + lim_fy
406.     IF temp < 270 - lim_fy THEN temp = 270 - lim_fy
407.     me6 = temp * pip180: me2 = me2 + COS(-me6) * actual_speed * 5 * kw + 1 * ABS(ks AND actual_speed < .05) 'calculating Z
408.     actual_deep = (map_deep(me0, me1) - .1) * map_dimz
409.     IF me2 > actual_deep - .1 THEN me2 = actual_deep - .1 - actual_speed / 4: me3 = me3 * .9: me4 = me4 * .9: flight(af, 17) = 1 ELSE flight(af, 17) = 0
410.     IF me2 < flight_high_max THEN me2 = flight_high_max
411.     position_accel = .02: ww = .9999: szog_xy = me5 + rotx 'calculating XY
412.  
413.     new_vec_x = -SIN(szog_xy) * position_accel * kw: new_vec_y = -COS(szog_xy) * position_accel * kw: vec_x = (me3 * ww + new_vec_x): vec_y = (me4 * ww + new_vec_y)
414.     vec_sum = SQR(vec_x * vec_x + vec_y * vec_y): IF vec_sum > position_speed THEN vec_sum = position_speed / vec_sum ELSE vec_sum = 1
415.     me3 = vec_x * vec_sum: me4 = vec_y * vec_sum: me0 = me0 + me3: me1 = me1 + me4: IF ks THEN deacc = .96 ELSE deacc = .99: IF kw THEN deacc = 1
416.     me3 = me3 * deacc: me4 = me4 * deacc: flight(af, 16) = 0
417.     IF crash_to_column(me0, me1) THEN
418.         flight(af, 16) = 1
419.         IF test = 0 THEN me3 = -me3 * .7: me4 = -me4 * .7: DO: m0 = m0 + me3: me1 = me1 + me4: LOOP WHILE crash_to_column(me0, me1)
420.     END IF
421.     flight(af, 1) = rotx: flight(af, 3) = me0: flight(af, 4) = me1: flight(af, 5) = me2: flight(af, 6) = me3: flight(af, 7) = me4: flight(af, 2) = rotx_bevel3
422.     flight(af, 12) = me6: flight(af, 13) = me5: flight(af, 15) = actual_speed
423.  
424. END SUB
425.  
426. FUNCTION map_deep (x, y): read_mapx1 = INT(x + map_resper2) MOD map_resolution: read_mapx2 = INT(x + map_resper2 + 1) MOD map_resolution
427.     read_mapy1 = INT(y + map_resper2) MOD map_resolution: read_mapy2 = INT(y + map_resper2 + 1) MOD map_resolution
428.     mapx1 = Interpolate(map(read_mapx1, read_mapy1, 0), map(read_mapx2, read_mapy1, 0), me(0) - INT(me(0)))
429.     mapx2 = Interpolate(map(read_mapx1, read_mapy2, 0), map(read_mapx2, read_mapy2, 0), me(0) - INT(me(0)))
430. map_deep = Interpolate(mapx1, mapx2, me(1) - INT(me(1))): END FUNCTION
431.  
432.  
433. FUNCTION crash_to_column (x, y): pixs_column = INT(map_dimxy / map_resolution / clmn_rad) + 2
434.     FOR tx = -pixs_column TO pixs_column: FOR ty = -pixs_column TO pixs_column
435.             read_mapx = INT(x + map_resper2 + tx) MOD map_resolution: read_mapy = INT(y + map_resper2 + ty) MOD map_resolution
436.             IF map(read_mapx, read_mapy, 6) THEN crash_to_column = 1: EXIT FUNCTION
437. NEXT ty, tx: END FUNCTION
438.  
439. FUNCTION degree (a, b): degreex = ATN((a + .00001) / (b + .00001)) / pip180: degreex = degreex - 180 * ABS(0 > b): degreex = degreex - 360 * (degreex < 0)
440. degree = degreex: END FUNCTION
441.  
442.