Cavegen: Place correct biome surface in tunnel entrances
[oweals/minetest.git] / src / cavegen.cpp
1 /*
2 Minetest
3 Copyright (C) 2010-2013 celeron55, Perttu Ahola <celeron55@gmail.com>
4
5 This program is free software; you can redistribute it and/or modify
6 it under the terms of the GNU Lesser General Public License as published by
7 the Free Software Foundation; either version 2.1 of the License, or
8 (at your option) any later version.
9
10 This program is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
13 GNU Lesser General Public License for more details.
14
15 You should have received a copy of the GNU Lesser General Public License along
16 with this program; if not, write to the Free Software Foundation, Inc.,
17 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
18 */
19
20 #include "util/numeric.h"
21 #include "map.h"
22 #include "mapgen.h"
23 #include "mapgen_v5.h"
24 #include "mapgen_v6.h"
25 #include "mapgen_v7.h"
26 #include "mg_biome.h"
27 #include "cavegen.h"
28
29 static NoiseParams nparams_caveliquids(0, 1, v3f(150.0, 150.0, 150.0), 776, 3, 0.6, 2.0);
30
31
32 ////
33 //// CavesNoiseIntersection
34 ////
35
36 CavesNoiseIntersection::CavesNoiseIntersection(
37         INodeDefManager *nodedef, BiomeManager *biomemgr, v3s16 chunksize,
38         NoiseParams *np_cave1, NoiseParams *np_cave2, s32 seed, float cave_width)
39 {
40         assert(nodedef);
41         assert(biomemgr);
42
43         m_ndef = nodedef;
44         m_bmgr = biomemgr;
45
46         m_csize = chunksize;
47         m_cave_width = cave_width;
48
49         m_ystride    = m_csize.X;
50         m_zstride_1d = m_csize.X * (m_csize.Y + 1);
51
52         // Noises are created using 1-down overgeneration
53         // A Nx-by-1-by-Nz-sized plane is at the bottom of the desired for
54         // re-carving the solid overtop placed for blocking sunlight
55         noise_cave1 = new Noise(np_cave1, seed, m_csize.X, m_csize.Y + 1, m_csize.Z);
56         noise_cave2 = new Noise(np_cave2, seed, m_csize.X, m_csize.Y + 1, m_csize.Z);
57 }
58
59
60 CavesNoiseIntersection::~CavesNoiseIntersection()
61 {
62         delete noise_cave1;
63         delete noise_cave2;
64 }
65
66
67 void CavesNoiseIntersection::generateCaves(MMVManip *vm,
68         v3s16 nmin, v3s16 nmax, u8 *biomemap)
69 {
70         assert(vm);
71         assert(biomemap);
72
73         noise_cave1->perlinMap3D(nmin.X, nmin.Y - 1, nmin.Z);
74         noise_cave2->perlinMap3D(nmin.X, nmin.Y - 1, nmin.Z);
75
76         v3s16 em = vm->m_area.getExtent();
77         u32 index2d = 0;  // Biomemap index
78
79         for (s16 z = nmin.Z; z <= nmax.Z; z++)
80         for (s16 x = nmin.X; x <= nmax.X; x++, index2d++) {
81                 bool column_is_open = false;  // Is column open to overground
82                 bool is_under_river = false;  // Is column under river water
83                 bool is_under_tunnel = false;  // Is tunnel or is under tunnel
84                 // Indexes at column top
85                 u32 vi = vm->m_area.index(x, nmax.Y, z);
86                 u32 index3d = (z - nmin.Z) * m_zstride_1d + m_csize.Y * m_ystride +
87                         (x - nmin.X);  // 3D noise index
88                 // Biome of column
89                 Biome *biome = (Biome *)m_bmgr->getRaw(biomemap[index2d]);
90                 u16 depth_top = biome->depth_top;
91                 u16 base_filler = depth_top + biome->depth_filler;
92                 u16 depth_riverbed = biome->depth_riverbed;
93                 u16 nplaced = 0;
94                 // Don't excavate the overgenerated stone at nmax.Y + 1,
95                 // this creates a 'roof' over the tunnel, preventing light in
96                 // tunnels at mapchunk borders when generating mapchunks upwards.
97                 // This 'roof' is removed when the mapchunk above is generated.
98                 for (s16 y = nmax.Y; y >= nmin.Y - 1; y--,
99                                 index3d -= m_ystride,
100                                 vm->m_area.add_y(em, vi, -1)) {
101
102                         content_t c = vm->m_data[vi].getContent();
103                         if (c == CONTENT_AIR || c == biome->c_water_top ||
104                                         c == biome->c_water) {
105                                 column_is_open = true;
106                                 continue;
107                         } else if (c == biome->c_river_water) {
108                                 column_is_open = true;
109                                 is_under_river = true;
110                                 continue;
111                         }
112                         // Ground
113                         float d1 = contour(noise_cave1->result[index3d]);
114                         float d2 = contour(noise_cave2->result[index3d]);
115
116                         if (d1 * d2 > m_cave_width && m_ndef->get(c).is_ground_content) {
117                                 // In tunnel and ground content, excavate
118                                 vm->m_data[vi] = MapNode(CONTENT_AIR);
119                                 is_under_tunnel = true;
120                         } else if (column_is_open && is_under_tunnel &&
121                                         (c == biome->c_stone || c == biome->c_filler)) {
122                                 // Tunnel entrance floor, place biome surface nodes
123                                 if (is_under_river) {
124                                         if (nplaced < depth_riverbed) {
125                                                 vm->m_data[vi] = MapNode(biome->c_riverbed);
126                                                 nplaced++;
127                                         } else {
128                                                 // Disable top/filler placement
129                                                 column_is_open = false;
130                                                 is_under_river = false;
131                                                 is_under_tunnel = false;
132                                         }
133                                 } else if (nplaced < depth_top) {
134                                         vm->m_data[vi] = MapNode(biome->c_top);
135                                         nplaced++;
136                                 } else if (nplaced < base_filler) {
137                                         vm->m_data[vi] = MapNode(biome->c_filler);
138                                         nplaced++;
139                                 } else {
140                                         // Disable top/filler placement
141                                         column_is_open = false;
142                                         is_under_tunnel = false;
143                                 }
144                         } else {
145                                 // Not tunnel or tunnel entrance floor
146                                 column_is_open = false;
147                         }
148                 }
149         }
150 }
151
152
153 ////
154 //// CavesRandomWalk
155 ////
156
157 CavesRandomWalk::CavesRandomWalk(
158         INodeDefManager *ndef,
159         GenerateNotifier *gennotify,
160         s32 seed,
161         int water_level,
162         content_t water_source,
163         content_t lava_source)
164 {
165         assert(ndef);
166
167         this->ndef           = ndef;
168         this->gennotify      = gennotify;
169         this->seed           = seed;
170         this->water_level    = water_level;
171         this->np_caveliquids = &nparams_caveliquids;
172         this->lava_depth     = DEFAULT_LAVA_DEPTH;
173
174         c_water_source = water_source;
175         if (c_water_source == CONTENT_IGNORE)
176                 c_water_source = ndef->getId("mapgen_water_source");
177         if (c_water_source == CONTENT_IGNORE)
178                 c_water_source = CONTENT_AIR;
179
180         c_lava_source = lava_source;
181         if (c_lava_source == CONTENT_IGNORE)
182                 c_lava_source = ndef->getId("mapgen_lava_source");
183         if (c_lava_source == CONTENT_IGNORE)
184                 c_lava_source = CONTENT_AIR;
185 }
186
187
188 void CavesRandomWalk::makeCave(MMVManip *vm, v3s16 nmin, v3s16 nmax,
189         PseudoRandom *ps, bool is_large_cave, int max_stone_height, s16 *heightmap)
190 {
191         assert(vm);
192         assert(ps);
193
194         this->vm         = vm;
195         this->ps         = ps;
196         this->node_min   = nmin;
197         this->node_max   = nmax;
198         this->heightmap  = heightmap;
199         this->large_cave = is_large_cave;
200
201         this->ystride = nmax.X - nmin.X + 1;
202
203         // Set initial parameters from randomness
204         int dswitchint = ps->range(1, 14);
205         flooded = ps->range(1, 2) == 2;
206
207         if (large_cave) {
208                 part_max_length_rs = ps->range(2, 4);
209                 tunnel_routepoints = ps->range(5, ps->range(15, 30));
210                 min_tunnel_diameter = 5;
211                 max_tunnel_diameter = ps->range(7, ps->range(8, 24));
212         } else {
213                 part_max_length_rs = ps->range(2, 9);
214                 tunnel_routepoints = ps->range(10, ps->range(15, 30));
215                 min_tunnel_diameter = 2;
216                 max_tunnel_diameter = ps->range(2, 6);
217         }
218
219         large_cave_is_flat = (ps->range(0, 1) == 0);
220
221         main_direction = v3f(0, 0, 0);
222
223         // Allowed route area size in nodes
224         ar = node_max - node_min + v3s16(1, 1, 1);
225         // Area starting point in nodes
226         of = node_min;
227
228         // Allow a bit more
229         //(this should be more than the maximum radius of the tunnel)
230         const s16 insure = 10;
231         s16 more = MYMAX(MAP_BLOCKSIZE - max_tunnel_diameter / 2 - insure, 1);
232         ar += v3s16(1, 0, 1) * more * 2;
233         of -= v3s16(1, 0, 1) * more;
234
235         route_y_min = 0;
236         // Allow half a diameter + 7 over stone surface
237         route_y_max = -of.Y + max_stone_y + max_tunnel_diameter / 2 + 7;
238
239         // Limit maximum to area
240         route_y_max = rangelim(route_y_max, 0, ar.Y - 1);
241
242         if (large_cave) {
243                 s16 minpos = 0;
244                 if (node_min.Y < water_level && node_max.Y > water_level) {
245                         minpos = water_level - max_tunnel_diameter / 3 - of.Y;
246                         route_y_max = water_level + max_tunnel_diameter / 3 - of.Y;
247                 }
248                 route_y_min = ps->range(minpos, minpos + max_tunnel_diameter);
249                 route_y_min = rangelim(route_y_min, 0, route_y_max);
250         }
251
252         s16 route_start_y_min = route_y_min;
253         s16 route_start_y_max = route_y_max;
254
255         route_start_y_min = rangelim(route_start_y_min, 0, ar.Y - 1);
256         route_start_y_max = rangelim(route_start_y_max, route_start_y_min, ar.Y - 1);
257
258         // Randomize starting position
259         orp.Z = (float)(ps->next() % ar.Z) + 0.5f;
260         orp.Y = (float)(ps->range(route_start_y_min, route_start_y_max)) + 0.5f;
261         orp.X = (float)(ps->next() % ar.X) + 0.5f;
262
263         // Add generation notify begin event
264         if (gennotify) {
265                 v3s16 abs_pos(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
266                 GenNotifyType notifytype = large_cave ?
267                         GENNOTIFY_LARGECAVE_BEGIN : GENNOTIFY_CAVE_BEGIN;
268                 gennotify->addEvent(notifytype, abs_pos);
269         }
270
271         // Generate some tunnel starting from orp
272         for (u16 j = 0; j < tunnel_routepoints; j++)
273                 makeTunnel(j % dswitchint == 0);
274
275         // Add generation notify end event
276         if (gennotify) {
277                 v3s16 abs_pos(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
278                 GenNotifyType notifytype = large_cave ?
279                         GENNOTIFY_LARGECAVE_END : GENNOTIFY_CAVE_END;
280                 gennotify->addEvent(notifytype, abs_pos);
281         }
282 }
283
284
285 void CavesRandomWalk::makeTunnel(bool dirswitch)
286 {
287         if (dirswitch && !large_cave) {
288                 main_direction.Z = ((float)(ps->next() % 20) - (float)10) / 10;
289                 main_direction.Y = ((float)(ps->next() % 20) - (float)10) / 30;
290                 main_direction.X = ((float)(ps->next() % 20) - (float)10) / 10;
291
292                 main_direction *= (float)ps->range(0, 10) / 10;
293         }
294
295         // Randomize size
296         s16 min_d = min_tunnel_diameter;
297         s16 max_d = max_tunnel_diameter;
298         rs = ps->range(min_d, max_d);
299         s16 rs_part_max_length_rs = rs * part_max_length_rs;
300
301         v3s16 maxlen;
302         if (large_cave) {
303                 maxlen = v3s16(
304                         rs_part_max_length_rs,
305                         rs_part_max_length_rs / 2,
306                         rs_part_max_length_rs
307                 );
308         } else {
309                 maxlen = v3s16(
310                         rs_part_max_length_rs,
311                         ps->range(1, rs_part_max_length_rs),
312                         rs_part_max_length_rs
313                 );
314         }
315
316         v3f vec;
317         // Jump downward sometimes
318         if (!large_cave && ps->range(0, 12) == 0) {
319                 vec.Z = (float)(ps->next() % (maxlen.Z * 1)) - (float)maxlen.Z / 2;
320                 vec.Y = (float)(ps->next() % (maxlen.Y * 2)) - (float)maxlen.Y;
321                 vec.X = (float)(ps->next() % (maxlen.X * 1)) - (float)maxlen.X / 2;
322         } else {
323                 vec.Z = (float)(ps->next() % (maxlen.Z * 1)) - (float)maxlen.Z / 2;
324                 vec.Y = (float)(ps->next() % (maxlen.Y * 1)) - (float)maxlen.Y / 2;
325                 vec.X = (float)(ps->next() % (maxlen.X * 1)) - (float)maxlen.X / 2;
326         }
327
328         // Do not make caves that are above ground.
329         // It is only necessary to check the startpoint and endpoint.
330         v3s16 p1 = v3s16(orp.X, orp.Y, orp.Z) + of + rs / 2;
331         v3s16 p2 = v3s16(vec.X, vec.Y, vec.Z) + p1;
332         if (isPosAboveSurface(p1) || isPosAboveSurface(p2))
333                 return;
334
335         vec += main_direction;
336
337         v3f rp = orp + vec;
338         if (rp.X < 0)
339                 rp.X = 0;
340         else if (rp.X >= ar.X)
341                 rp.X = ar.X - 1;
342
343         if (rp.Y < route_y_min)
344                 rp.Y = route_y_min;
345         else if (rp.Y >= route_y_max)
346                 rp.Y = route_y_max - 1;
347
348         if (rp.Z < 0)
349                 rp.Z = 0;
350         else if (rp.Z >= ar.Z)
351                 rp.Z = ar.Z - 1;
352
353         vec = rp - orp;
354
355         float veclen = vec.getLength();
356         if (veclen < 0.05f)
357                 veclen = 1.0f;
358
359         // Every second section is rough
360         bool randomize_xz = (ps->range(1, 2) == 1);
361
362         // Carve routes
363         for (float f = 0.f; f < 1.0f; f += 1.0f / veclen)
364                 carveRoute(vec, f, randomize_xz);
365
366         orp = rp;
367 }
368
369
370 void CavesRandomWalk::carveRoute(v3f vec, float f, bool randomize_xz)
371 {
372         MapNode airnode(CONTENT_AIR);
373         MapNode waternode(c_water_source);
374         MapNode lavanode(c_lava_source);
375
376         v3s16 startp(orp.X, orp.Y, orp.Z);
377         startp += of;
378
379         float nval = NoisePerlin3D(np_caveliquids, startp.X,
380                 startp.Y, startp.Z, seed);
381         MapNode liquidnode = (nval < 0.40f && node_max.Y < lava_depth) ?
382                 lavanode : waternode;
383
384         v3f fp = orp + vec * f;
385         fp.X += 0.1f * ps->range(-10, 10);
386         fp.Z += 0.1f * ps->range(-10, 10);
387         v3s16 cp(fp.X, fp.Y, fp.Z);
388
389         s16 d0 = -rs / 2;
390         s16 d1 = d0 + rs;
391         if (randomize_xz) {
392                 d0 += ps->range(-1, 1);
393                 d1 += ps->range(-1, 1);
394         }
395
396         bool flat_cave_floor = !large_cave && ps->range(0, 2) == 2;
397
398         for (s16 z0 = d0; z0 <= d1; z0++) {
399                 s16 si = rs / 2 - MYMAX(0, abs(z0) - rs / 7 - 1);
400                 for (s16 x0 = -si - ps->range(0,1); x0 <= si - 1 + ps->range(0,1); x0++) {
401                         s16 maxabsxz = MYMAX(abs(x0), abs(z0));
402
403                         s16 si2 = rs / 2 - MYMAX(0, maxabsxz - rs / 7 - 1);
404
405                         for (s16 y0 = -si2; y0 <= si2; y0++) {
406                                 // Make better floors in small caves
407                                 if (flat_cave_floor && y0 <= -rs / 2 && rs <= 7)
408                                         continue;
409
410                                 if (large_cave_is_flat) {
411                                         // Make large caves not so tall
412                                         if (rs > 7 && abs(y0) >= rs / 3)
413                                                 continue;
414                                 }
415
416                                 v3s16 p(cp.X + x0, cp.Y + y0, cp.Z + z0);
417                                 p += of;
418
419                                 if (vm->m_area.contains(p) == false)
420                                         continue;
421
422                                 u32 i = vm->m_area.index(p);
423                                 content_t c = vm->m_data[i].getContent();
424                                 if (!ndef->get(c).is_ground_content)
425                                         continue;
426
427                                 if (large_cave) {
428                                         int full_ymin = node_min.Y - MAP_BLOCKSIZE;
429                                         int full_ymax = node_max.Y + MAP_BLOCKSIZE;
430
431                                         if (flooded && full_ymin < water_level && full_ymax > water_level)
432                                                 vm->m_data[i] = (p.Y <= water_level) ? waternode : airnode;
433                                         else if (flooded && full_ymax < water_level)
434                                                 vm->m_data[i] = (p.Y < startp.Y - 4) ? liquidnode : airnode;
435                                         else
436                                                 vm->m_data[i] = airnode;
437                                 } else {
438                                         if (c == CONTENT_IGNORE)
439                                                 continue;
440
441                                         vm->m_data[i] = airnode;
442                                         vm->m_flags[i] |= VMANIP_FLAG_CAVE;
443                                 }
444                         }
445                 }
446         }
447 }
448
449
450 inline bool CavesRandomWalk::isPosAboveSurface(v3s16 p)
451 {
452         if (heightmap != NULL &&
453                         p.Z >= node_min.Z && p.Z <= node_max.Z &&
454                         p.X >= node_min.X && p.X <= node_max.X) {
455                 u32 index = (p.Z - node_min.Z) * ystride + (p.X - node_min.X);
456                 if (heightmap[index] < p.Y)
457                         return true;
458         } else if (p.Y > water_level) {
459                 return true;
460         }
461
462         return false;
463 }
464
465
466 ////
467 //// CavesV6
468 ////
469
470 CavesV6::CavesV6(INodeDefManager *ndef, GenerateNotifier *gennotify,
471         int water_level, content_t water_source, content_t lava_source)
472 {
473         assert(ndef);
474
475         this->ndef        = ndef;
476         this->gennotify   = gennotify;
477         this->water_level = water_level;
478
479         c_water_source = water_source;
480         if (c_water_source == CONTENT_IGNORE)
481                 c_water_source = ndef->getId("mapgen_water_source");
482         if (c_water_source == CONTENT_IGNORE)
483                 c_water_source = CONTENT_AIR;
484
485         c_lava_source = lava_source;
486         if (c_lava_source == CONTENT_IGNORE)
487                 c_lava_source = ndef->getId("mapgen_lava_source");
488         if (c_lava_source == CONTENT_IGNORE)
489                 c_lava_source = CONTENT_AIR;
490 }
491
492
493 void CavesV6::makeCave(MMVManip *vm, v3s16 nmin, v3s16 nmax,
494         PseudoRandom *ps, PseudoRandom *ps2,
495         bool is_large_cave, int max_stone_height, s16 *heightmap)
496 {
497         assert(vm);
498         assert(ps);
499         assert(ps2);
500
501         this->vm         = vm;
502         this->ps         = ps;
503         this->ps2        = ps2;
504         this->node_min   = nmin;
505         this->node_max   = nmax;
506         this->heightmap  = heightmap;
507         this->large_cave = is_large_cave;
508
509         this->ystride = nmax.X - nmin.X + 1;
510
511         // Set initial parameters from randomness
512         min_tunnel_diameter = 2;
513         max_tunnel_diameter = ps->range(2, 6);
514         int dswitchint      = ps->range(1, 14);
515         if (large_cave) {
516                 part_max_length_rs  = ps->range(2, 4);
517                 tunnel_routepoints  = ps->range(5, ps->range(15, 30));
518                 min_tunnel_diameter = 5;
519                 max_tunnel_diameter = ps->range(7, ps->range(8, 24));
520         } else {
521                 part_max_length_rs = ps->range(2, 9);
522                 tunnel_routepoints = ps->range(10, ps->range(15, 30));
523         }
524         large_cave_is_flat = (ps->range(0, 1) == 0);
525
526         main_direction = v3f(0, 0, 0);
527
528         // Allowed route area size in nodes
529         ar = node_max - node_min + v3s16(1, 1, 1);
530         // Area starting point in nodes
531         of = node_min;
532
533         // Allow a bit more
534         //(this should be more than the maximum radius of the tunnel)
535         const s16 max_spread_amount = MAP_BLOCKSIZE;
536         const s16 insure = 10;
537         s16 more = MYMAX(max_spread_amount - max_tunnel_diameter / 2 - insure, 1);
538         ar += v3s16(1, 0, 1) * more * 2;
539         of -= v3s16(1, 0, 1) * more;
540
541         route_y_min = 0;
542         // Allow half a diameter + 7 over stone surface
543         route_y_max = -of.Y + max_stone_height + max_tunnel_diameter / 2 + 7;
544
545         // Limit maximum to area
546         route_y_max = rangelim(route_y_max, 0, ar.Y - 1);
547
548         if (large_cave) {
549                 s16 minpos = 0;
550                 if (node_min.Y < water_level && node_max.Y > water_level) {
551                         minpos = water_level - max_tunnel_diameter / 3 - of.Y;
552                         route_y_max = water_level + max_tunnel_diameter / 3 - of.Y;
553                 }
554                 route_y_min = ps->range(minpos, minpos + max_tunnel_diameter);
555                 route_y_min = rangelim(route_y_min, 0, route_y_max);
556         }
557
558         s16 route_start_y_min = route_y_min;
559         s16 route_start_y_max = route_y_max;
560
561         route_start_y_min = rangelim(route_start_y_min, 0, ar.Y - 1);
562         route_start_y_max = rangelim(route_start_y_max, route_start_y_min, ar.Y - 1);
563
564         // Randomize starting position
565         orp.Z = (float)(ps->next() % ar.Z) + 0.5f;
566         orp.Y = (float)(ps->range(route_start_y_min, route_start_y_max)) + 0.5f;
567         orp.X = (float)(ps->next() % ar.X) + 0.5f;
568
569         // Add generation notify begin event
570         if (gennotify != NULL) {
571                 v3s16 abs_pos(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
572                 GenNotifyType notifytype = large_cave ?
573                         GENNOTIFY_LARGECAVE_BEGIN : GENNOTIFY_CAVE_BEGIN;
574                 gennotify->addEvent(notifytype, abs_pos);
575         }
576
577         // Generate some tunnel starting from orp
578         for (u16 j = 0; j < tunnel_routepoints; j++)
579                 makeTunnel(j % dswitchint == 0);
580
581         // Add generation notify end event
582         if (gennotify != NULL) {
583                 v3s16 abs_pos(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
584                 GenNotifyType notifytype = large_cave ?
585                         GENNOTIFY_LARGECAVE_END : GENNOTIFY_CAVE_END;
586                 gennotify->addEvent(notifytype, abs_pos);
587         }
588 }
589
590
591 void CavesV6::makeTunnel(bool dirswitch)
592 {
593         if (dirswitch && !large_cave) {
594                 main_direction.Z = ((float)(ps->next() % 20) - (float)10) / 10;
595                 main_direction.Y = ((float)(ps->next() % 20) - (float)10) / 30;
596                 main_direction.X = ((float)(ps->next() % 20) - (float)10) / 10;
597
598                 main_direction *= (float)ps->range(0, 10) / 10;
599         }
600
601         // Randomize size
602         s16 min_d = min_tunnel_diameter;
603         s16 max_d = max_tunnel_diameter;
604         rs = ps->range(min_d, max_d);
605         s16 rs_part_max_length_rs = rs * part_max_length_rs;
606
607         v3s16 maxlen;
608         if (large_cave) {
609                 maxlen = v3s16(
610                         rs_part_max_length_rs,
611                         rs_part_max_length_rs / 2,
612                         rs_part_max_length_rs
613                 );
614         } else {
615                 maxlen = v3s16(
616                         rs_part_max_length_rs,
617                         ps->range(1, rs_part_max_length_rs),
618                         rs_part_max_length_rs
619                 );
620         }
621
622         v3f vec;
623         vec.Z = (float)(ps->next() % maxlen.Z) - (float)maxlen.Z / 2;
624         vec.Y = (float)(ps->next() % maxlen.Y) - (float)maxlen.Y / 2;
625         vec.X = (float)(ps->next() % maxlen.X) - (float)maxlen.X / 2;
626
627         // Jump downward sometimes
628         if (!large_cave && ps->range(0, 12) == 0) {
629                 vec.Z = (float)(ps->next() % maxlen.Z) - (float)maxlen.Z / 2;
630                 vec.Y = (float)(ps->next() % (maxlen.Y * 2)) - (float)maxlen.Y;
631                 vec.X = (float)(ps->next() % maxlen.X) - (float)maxlen.X / 2;
632         }
633
634         // Do not make caves that are entirely above ground, to fix shadow bugs
635         // caused by overgenerated large caves.
636         // It is only necessary to check the startpoint and endpoint.
637         v3s16 p1 = v3s16(orp.X, orp.Y, orp.Z) + of + rs / 2;
638         v3s16 p2 = v3s16(vec.X, vec.Y, vec.Z) + p1;
639
640         // If startpoint and endpoint are above ground, disable placement of nodes
641         // in carveRoute while still running all PseudoRandom calls to ensure caves
642         // are consistent with existing worlds.
643         bool tunnel_above_ground =
644                 p1.Y > getSurfaceFromHeightmap(p1) &&
645                 p2.Y > getSurfaceFromHeightmap(p2);
646
647         vec += main_direction;
648
649         v3f rp = orp + vec;
650         if (rp.X < 0)
651                 rp.X = 0;
652         else if (rp.X >= ar.X)
653                 rp.X = ar.X - 1;
654
655         if (rp.Y < route_y_min)
656                 rp.Y = route_y_min;
657         else if (rp.Y >= route_y_max)
658                 rp.Y = route_y_max - 1;
659
660         if (rp.Z < 0)
661                 rp.Z = 0;
662         else if (rp.Z >= ar.Z)
663                 rp.Z = ar.Z - 1;
664
665         vec = rp - orp;
666
667         float veclen = vec.getLength();
668         // As odd as it sounds, veclen is *exactly* 0.0 sometimes, causing a FPE
669         if (veclen < 0.05f)
670                 veclen = 1.0f;
671
672         // Every second section is rough
673         bool randomize_xz = (ps2->range(1, 2) == 1);
674
675         // Carve routes
676         for (float f = 0.f; f < 1.0f; f += 1.0f / veclen)
677                 carveRoute(vec, f, randomize_xz, tunnel_above_ground);
678
679         orp = rp;
680 }
681
682
683 void CavesV6::carveRoute(v3f vec, float f, bool randomize_xz,
684         bool tunnel_above_ground)
685 {
686         MapNode airnode(CONTENT_AIR);
687         MapNode waternode(c_water_source);
688         MapNode lavanode(c_lava_source);
689
690         v3s16 startp(orp.X, orp.Y, orp.Z);
691         startp += of;
692
693         v3f fp = orp + vec * f;
694         fp.X += 0.1f * ps->range(-10, 10);
695         fp.Z += 0.1f * ps->range(-10, 10);
696         v3s16 cp(fp.X, fp.Y, fp.Z);
697
698         s16 d0 = -rs / 2;
699         s16 d1 = d0 + rs;
700         if (randomize_xz) {
701                 d0 += ps->range(-1, 1);
702                 d1 += ps->range(-1, 1);
703         }
704
705         for (s16 z0 = d0; z0 <= d1; z0++) {
706                 s16 si = rs / 2 - MYMAX(0, abs(z0) - rs / 7 - 1);
707                 for (s16 x0 = -si - ps->range(0,1); x0 <= si - 1 + ps->range(0,1); x0++) {
708                         if (tunnel_above_ground)
709                                 continue;
710
711                         s16 maxabsxz = MYMAX(abs(x0), abs(z0));
712                         s16 si2 = rs / 2 - MYMAX(0, maxabsxz - rs / 7 - 1);
713                         for (s16 y0 = -si2; y0 <= si2; y0++) {
714                                 if (large_cave_is_flat) {
715                                         // Make large caves not so tall
716                                         if (rs > 7 && abs(y0) >= rs / 3)
717                                                 continue;
718                                 }
719
720                                 v3s16 p(cp.X + x0, cp.Y + y0, cp.Z + z0);
721                                 p += of;
722
723                                 if (vm->m_area.contains(p) == false)
724                                         continue;
725
726                                 u32 i = vm->m_area.index(p);
727                                 content_t c = vm->m_data[i].getContent();
728                                 if (!ndef->get(c).is_ground_content)
729                                         continue;
730
731                                 if (large_cave) {
732                                         int full_ymin = node_min.Y - MAP_BLOCKSIZE;
733                                         int full_ymax = node_max.Y + MAP_BLOCKSIZE;
734
735                                         if (full_ymin < water_level && full_ymax > water_level) {
736                                                 vm->m_data[i] = (p.Y <= water_level) ? waternode : airnode;
737                                         } else if (full_ymax < water_level) {
738                                                 vm->m_data[i] = (p.Y < startp.Y - 2) ? lavanode : airnode;
739                                         } else {
740                                                 vm->m_data[i] = airnode;
741                                         }
742                                 } else {
743                                         if (c == CONTENT_IGNORE || c == CONTENT_AIR)
744                                                 continue;
745
746                                         vm->m_data[i] = airnode;
747                                         vm->m_flags[i] |= VMANIP_FLAG_CAVE;
748                                 }
749                         }
750                 }
751         }
752 }
753
754
755 inline s16 CavesV6::getSurfaceFromHeightmap(v3s16 p)
756 {
757         if (heightmap != NULL &&
758                         p.Z >= node_min.Z && p.Z <= node_max.Z &&
759                         p.X >= node_min.X && p.X <= node_max.X) {
760                 u32 index = (p.Z - node_min.Z) * ystride + (p.X - node_min.X);
761                 return heightmap[index];
762         } else {
763                 return water_level;
764         }
765 }