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 "cavegen.h"
  27
  28 NoiseParams nparams_caveliquids(0, 1, v3f(150.0, 150.0, 150.0), 776, 3, 0.6, 2.0);
  29
  30 ////
  31 //// CavesRandomWalk
  32 ////
  33
  34 CavesRandomWalk::CavesRandomWalk(
  35         INodeDefManager *ndef,
  36         GenerateNotifier *gennotify,
  37         int seed,
  38         int water_level,
  39         content_t water_source,
  40         content_t lava_source)
  41 {
  42         assert(ndef);
  43
  44         this->ndef           = ndef;
  45         this->gennotify      = gennotify;
  46         this->seed           = seed;
  47         this->water_level    = water_level;
  48         this->np_caveliquids = &nparams_caveliquids;
  49         this->lava_depth     = DEFAULT_LAVA_DEPTH;
  50
  51         c_water_source = water_source;
  52         if (c_water_source == CONTENT_IGNORE)
  53                 c_water_source = ndef->getId("mapgen_water_source");
  54         if (c_water_source == CONTENT_IGNORE)
  55                 c_water_source = CONTENT_AIR;
  56
  57         c_lava_source = lava_source;
  58         if (c_lava_source == CONTENT_IGNORE)
  59                 c_lava_source = ndef->getId("mapgen_lava_source");
  60         if (c_lava_source == CONTENT_IGNORE)
  61                 c_lava_source = CONTENT_AIR;
  62 }
  63
  64
  65 void CavesRandomWalk::makeCave(MMVManip *vm, v3s16 nmin, v3s16 nmax,
  66         PseudoRandom *ps, bool is_large_cave, int max_stone_height, s16 *heightmap)
  67 {
  68         assert(vm);
  69         assert(ps);
  70
  71         this->vm         = vm;
  72         this->ps         = ps;
  73         this->node_min   = nmin;
  74         this->node_max   = nmax;
  75         this->heightmap  = heightmap;
  76         this->large_cave = is_large_cave;
  77
  78         this->ystride = nmax.X - nmin.X + 1;
  79
  80         // Set initial parameters from randomness
  81         dswitchint = ps->range(1, 14);
  82         flooded    = ps->range(1, 2) == 2;
  83
  84         if (large_cave) {
  85                 part_max_length_rs = ps->range(2, 4);
  86                 tunnel_routepoints = ps->range(5, ps->range(15, 30));
  87                 min_tunnel_diameter = 5;
  88                 max_tunnel_diameter = ps->range(7, ps->range(8, 24));
  89         } else {
  90                 part_max_length_rs = ps->range(2, 9);
  91                 tunnel_routepoints = ps->range(10, ps->range(15, 30));
  92                 min_tunnel_diameter = 2;
  93                 max_tunnel_diameter = ps->range(2, 6);
  94         }
  95
  96         large_cave_is_flat = (ps->range(0, 1) == 0);
  97
  98         main_direction = v3f(0, 0, 0);
  99
 100         // Allowed route area size in nodes
 101         ar = node_max - node_min + v3s16(1, 1, 1);
 102         // Area starting point in nodes
 103         of = node_min;
 104
 105         // Allow a bit more
 106         //(this should be more than the maximum radius of the tunnel)
 107         s16 insure = 10;
 108         s16 more = MYMAX(MAP_BLOCKSIZE - max_tunnel_diameter / 2 - insure, 1);
 109         ar += v3s16(1,0,1) * more * 2;
 110         of -= v3s16(1,0,1) * more;
 111
 112         route_y_min = 0;
 113         // Allow half a diameter + 7 over stone surface
 114         route_y_max = -of.Y + max_stone_y + max_tunnel_diameter / 2 + 7;
 115
 116         // Limit maximum to area
 117         route_y_max = rangelim(route_y_max, 0, ar.Y - 1);
 118
 119         if (large_cave) {
 120                 s16 minpos = 0;
 121                 if (node_min.Y < water_level && node_max.Y > water_level) {
 122                         minpos = water_level - max_tunnel_diameter / 3 - of.Y;
 123                         route_y_max = water_level + max_tunnel_diameter / 3 - of.Y;
 124                 }
 125                 route_y_min = ps->range(minpos, minpos + max_tunnel_diameter);
 126                 route_y_min = rangelim(route_y_min, 0, route_y_max);
 127         }
 128
 129         s16 route_start_y_min = route_y_min;
 130         s16 route_start_y_max = route_y_max;
 131
 132         route_start_y_min = rangelim(route_start_y_min, 0, ar.Y - 1);
 133         route_start_y_max = rangelim(route_start_y_max, route_start_y_min, ar.Y - 1);
 134
 135         // Randomize starting position
 136         orp.Z = (float)(ps->next() % ar.Z) + 0.5;
 137         orp.Y = (float)(ps->range(route_start_y_min, route_start_y_max)) + 0.5;
 138         orp.X = (float)(ps->next() % ar.X) + 0.5;
 139
 140         // Add generation notify begin event
 141         if (gennotify) {
 142                 v3s16 abs_pos(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 143                 GenNotifyType notifytype = large_cave ?
 144                         GENNOTIFY_LARGECAVE_BEGIN : GENNOTIFY_CAVE_BEGIN;
 145                 gennotify->addEvent(notifytype, abs_pos);
 146         }
 147
 148         // Generate some tunnel starting from orp
 149         for (u16 j = 0; j < tunnel_routepoints; j++)
 150                 makeTunnel(j % dswitchint == 0);
 151
 152         // Add generation notify end event
 153         if (gennotify) {
 154                 v3s16 abs_pos(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 155                 GenNotifyType notifytype = large_cave ?
 156                         GENNOTIFY_LARGECAVE_END : GENNOTIFY_CAVE_END;
 157                 gennotify->addEvent(notifytype, abs_pos);
 158         }
 159 }
 160
 161
 162 void CavesRandomWalk::makeTunnel(bool dirswitch)
 163 {
 164         if (dirswitch && !large_cave) {
 165                 main_direction.Z = ((float)(ps->next() % 20) - (float)10) / 10;
 166                 main_direction.Y = ((float)(ps->next() % 20) - (float)10) / 30;
 167                 main_direction.X = ((float)(ps->next() % 20) - (float)10) / 10;
 168
 169                 main_direction *= (float)ps->range(0, 10) / 10;
 170         }
 171
 172         // Randomize size
 173         s16 min_d = min_tunnel_diameter;
 174         s16 max_d = max_tunnel_diameter;
 175         rs = ps->range(min_d, max_d);
 176         s16 rs_part_max_length_rs = rs * part_max_length_rs;
 177
 178         v3s16 maxlen;
 179         if (large_cave) {
 180                 maxlen = v3s16(
 181                         rs_part_max_length_rs,
 182                         rs_part_max_length_rs / 2,
 183                         rs_part_max_length_rs
 184                 );
 185         } else {
 186                 maxlen = v3s16(
 187                         rs_part_max_length_rs,
 188                         ps->range(1, rs_part_max_length_rs),
 189                         rs_part_max_length_rs
 190                 );
 191         }
 192
 193         v3f vec;
 194         // Jump downward sometimes
 195         if (!large_cave && ps->range(0, 12) == 0) {
 196                 vec.Z = (float)(ps->next() % (maxlen.Z * 1)) - (float)maxlen.Z / 2;
 197                 vec.Y = (float)(ps->next() % (maxlen.Y * 2)) - (float)maxlen.Y;
 198                 vec.X = (float)(ps->next() % (maxlen.X * 1)) - (float)maxlen.X / 2;
 199         } else {
 200                 vec.Z = (float)(ps->next() % (maxlen.Z * 1)) - (float)maxlen.Z / 2;
 201                 vec.Y = (float)(ps->next() % (maxlen.Y * 1)) - (float)maxlen.Y / 2;
 202                 vec.X = (float)(ps->next() % (maxlen.X * 1)) - (float)maxlen.X / 2;
 203         }
 204
 205         // Do not make caves that are above ground.
 206         // It is only necessary to check the startpoint and endpoint.
 207         v3s16 p1 = v3s16(orp.X, orp.Y, orp.Z) + of + rs / 2;
 208         v3s16 p2 = v3s16(vec.X, vec.Y, vec.Z) + p1;
 209         if (isPosAboveSurface(p1) || isPosAboveSurface(p2))
 210                 return;
 211
 212         vec += main_direction;
 213
 214         v3f rp = orp + vec;
 215         if (rp.X < 0)
 216                 rp.X = 0;
 217         else if (rp.X >= ar.X)
 218                 rp.X = ar.X - 1;
 219
 220         if (rp.Y < route_y_min)
 221                 rp.Y = route_y_min;
 222         else if (rp.Y >= route_y_max)
 223                 rp.Y = route_y_max - 1;
 224
 225         if (rp.Z < 0)
 226                 rp.Z = 0;
 227         else if (rp.Z >= ar.Z)
 228                 rp.Z = ar.Z - 1;
 229
 230         vec = rp - orp;
 231
 232         float veclen = vec.getLength();
 233         if (veclen < 0.05)
 234                 veclen = 1.0;
 235
 236         // Every second section is rough
 237         bool randomize_xz = (ps->range(1, 2) == 1);
 238
 239         // Carve routes
 240         for (float f = 0; f < 1.0; f += 1.0 / veclen)
 241                 carveRoute(vec, f, randomize_xz);
 242
 243         orp = rp;
 244 }
 245
 246
 247 void CavesRandomWalk::carveRoute(v3f vec, float f, bool randomize_xz)
 248 {
 249         MapNode airnode(CONTENT_AIR);
 250         MapNode waternode(c_water_source);
 251         MapNode lavanode(c_lava_source);
 252
 253         v3s16 startp(orp.X, orp.Y, orp.Z);
 254         startp += of;
 255
 256         float nval = NoisePerlin3D(np_caveliquids, startp.X,
 257                 startp.Y, startp.Z, seed);
 258         MapNode liquidnode = (nval < 0.40 && node_max.Y < lava_depth) ?
 259                 lavanode : waternode;
 260
 261         v3f fp = orp + vec * f;
 262         fp.X += 0.1 * ps->range(-10, 10);
 263         fp.Z += 0.1 * ps->range(-10, 10);
 264         v3s16 cp(fp.X, fp.Y, fp.Z);
 265
 266         s16 d0 = -rs/2;
 267         s16 d1 = d0 + rs;
 268         if (randomize_xz) {
 269                 d0 += ps->range(-1, 1);
 270                 d1 += ps->range(-1, 1);
 271         }
 272
 273         bool flat_cave_floor = !large_cave && ps->range(0, 2) == 2;
 274
 275         for (s16 z0 = d0; z0 <= d1; z0++) {
 276                 s16 si = rs / 2 - MYMAX(0, abs(z0) - rs / 7 - 1);
 277                 for (s16 x0 = -si - ps->range(0,1); x0 <= si - 1 + ps->range(0,1); x0++) {
 278                         s16 maxabsxz = MYMAX(abs(x0), abs(z0));
 279
 280                         s16 si2 = rs / 2 - MYMAX(0, maxabsxz - rs / 7 - 1);
 281
 282                         for (s16 y0 = -si2; y0 <= si2; y0++) {
 283                                 // Make better floors in small caves
 284                                 if (flat_cave_floor && y0 <= -rs / 2 && rs <= 7)
 285                                         continue;
 286
 287                                 if (large_cave_is_flat) {
 288                                         // Make large caves not so tall
 289                                         if (rs > 7 && abs(y0) >= rs / 3)
 290                                                 continue;
 291                                 }
 292
 293                                 v3s16 p(cp.X + x0, cp.Y + y0, cp.Z + z0);
 294                                 p += of;
 295
 296                                 if (vm->m_area.contains(p) == false)
 297                                         continue;
 298
 299                                 u32 i = vm->m_area.index(p);
 300                                 content_t c = vm->m_data[i].getContent();
 301                                 if (!ndef->get(c).is_ground_content)
 302                                         continue;
 303
 304                                 if (large_cave) {
 305                                         int full_ymin = node_min.Y - MAP_BLOCKSIZE;
 306                                         int full_ymax = node_max.Y + MAP_BLOCKSIZE;
 307
 308                                         if (flooded && full_ymin < water_level && full_ymax > water_level)
 309                                                 vm->m_data[i] = (p.Y <= water_level) ? waternode : airnode;
 310                                         else if (flooded && full_ymax < water_level)
 311                                                 vm->m_data[i] = (p.Y < startp.Y - 4) ? liquidnode : airnode;
 312                                         else
 313                                                 vm->m_data[i] = airnode;
 314                                 } else {
 315                                         if (c == CONTENT_IGNORE)
 316                                                 continue;
 317
 318                                         vm->m_data[i] = airnode;
 319                                         vm->m_flags[i] |= VMANIP_FLAG_CAVE;
 320                                 }
 321                         }
 322                 }
 323         }
 324 }
 325
 326
 327 inline bool CavesRandomWalk::isPosAboveSurface(v3s16 p)
 328 {
 329         if (heightmap != NULL &&
 330                         p.Z >= node_min.Z && p.Z <= node_max.Z &&
 331                         p.X >= node_min.X && p.X <= node_max.X) {
 332                 u32 index = (p.Z - node_min.Z) * ystride + (p.X - node_min.X);
 333                 if (heightmap[index] < p.Y)
 334                         return true;
 335         } else if (p.Y > water_level) {
 336                 return true;
 337         }
 338
 339         return false;
 340 }
 341
 342
 343 ////
 344 //// CavesV6
 345 ////
 346
 347 CavesV6::CavesV6(INodeDefManager *ndef, GenerateNotifier *gennotify,
 348         int water_level, content_t water_source, content_t lava_source)
 349 {
 350         assert(ndef);
 351
 352         this->ndef        = ndef;
 353         this->gennotify   = gennotify;
 354         this->water_level = water_level;
 355
 356         c_water_source = water_source;
 357         if (c_water_source == CONTENT_IGNORE)
 358                 c_water_source = ndef->getId("mapgen_water_source");
 359         if (c_water_source == CONTENT_IGNORE)
 360                 c_water_source = CONTENT_AIR;
 361
 362         c_lava_source = lava_source;
 363         if (c_lava_source == CONTENT_IGNORE)
 364                 c_lava_source = ndef->getId("mapgen_lava_source");
 365         if (c_lava_source == CONTENT_IGNORE)
 366                 c_lava_source = CONTENT_AIR;
 367 }
 368
 369
 370 void CavesV6::makeCave(MMVManip *vm, v3s16 nmin, v3s16 nmax,
 371         PseudoRandom *ps, PseudoRandom *ps2,
 372         bool is_large_cave, int max_stone_height, s16 *heightmap)
 373 {
 374         assert(vm);
 375         assert(ps);
 376         assert(ps2);
 377
 378         this->vm         = vm;
 379         this->ps         = ps;
 380         this->ps2        = ps2;
 381         this->node_min   = nmin;
 382         this->node_max   = nmax;
 383         this->heightmap  = heightmap;
 384         this->large_cave = is_large_cave;
 385
 386         this->ystride = nmax.X - nmin.X + 1;
 387
 388         // Set initial parameters from randomness
 389         min_tunnel_diameter = 2;
 390         max_tunnel_diameter = ps->range(2, 6);
 391         dswitchint          = ps->range(1, 14);
 392         if (large_cave) {
 393                 part_max_length_rs  = ps->range(2, 4);
 394                 tunnel_routepoints  = ps->range(5, ps->range(15, 30));
 395                 min_tunnel_diameter = 5;
 396                 max_tunnel_diameter = ps->range(7, ps->range(8, 24));
 397         } else {
 398                 part_max_length_rs = ps->range(2, 9);
 399                 tunnel_routepoints = ps->range(10, ps->range(15, 30));
 400         }
 401         large_cave_is_flat = (ps->range(0, 1) == 0);
 402
 403         main_direction = v3f(0, 0, 0);
 404
 405         // Allowed route area size in nodes
 406         ar = node_max - node_min + v3s16(1, 1, 1);
 407         // Area starting point in nodes
 408         of = node_min;
 409
 410         // Allow a bit more
 411         //(this should be more than the maximum radius of the tunnel)
 412         const s16 max_spread_amount = MAP_BLOCKSIZE;
 413         s16 insure = 10;
 414         s16 more = MYMAX(max_spread_amount - max_tunnel_diameter / 2 - insure, 1);
 415         ar += v3s16(1, 0, 1) * more * 2;
 416         of -= v3s16(1, 0, 1) * more;
 417
 418         route_y_min = 0;
 419         // Allow half a diameter + 7 over stone surface
 420         route_y_max = -of.Y + max_stone_height + max_tunnel_diameter / 2 + 7;
 421
 422         // Limit maximum to area
 423         route_y_max = rangelim(route_y_max, 0, ar.Y - 1);
 424
 425         if (large_cave) {
 426                 s16 minpos = 0;
 427                 if (node_min.Y < water_level && node_max.Y > water_level) {
 428                         minpos = water_level - max_tunnel_diameter / 3 - of.Y;
 429                         route_y_max = water_level + max_tunnel_diameter / 3 - of.Y;
 430                 }
 431                 route_y_min = ps->range(minpos, minpos + max_tunnel_diameter);
 432                 route_y_min = rangelim(route_y_min, 0, route_y_max);
 433         }
 434
 435         s16 route_start_y_min = route_y_min;
 436         s16 route_start_y_max = route_y_max;
 437
 438         route_start_y_min = rangelim(route_start_y_min, 0, ar.Y - 1);
 439         route_start_y_max = rangelim(route_start_y_max, route_start_y_min, ar.Y - 1);
 440
 441         // Randomize starting position
 442         orp.Z = (float)(ps->next() % ar.Z) + 0.5;
 443         orp.Y = (float)(ps->range(route_start_y_min, route_start_y_max)) + 0.5;
 444         orp.X = (float)(ps->next() % ar.X) + 0.5;
 445
 446         // Add generation notify begin event
 447         if (gennotify != NULL) {
 448                 v3s16 abs_pos(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 449                 GenNotifyType notifytype = large_cave ?
 450                         GENNOTIFY_LARGECAVE_BEGIN : GENNOTIFY_CAVE_BEGIN;
 451                 gennotify->addEvent(notifytype, abs_pos);
 452         }
 453
 454         // Generate some tunnel starting from orp
 455         for (u16 j = 0; j < tunnel_routepoints; j++)
 456                 makeTunnel(j % dswitchint == 0);
 457
 458         // Add generation notify end event
 459         if (gennotify != NULL) {
 460                 v3s16 abs_pos(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 461                 GenNotifyType notifytype = large_cave ?
 462                         GENNOTIFY_LARGECAVE_END : GENNOTIFY_CAVE_END;
 463                 gennotify->addEvent(notifytype, abs_pos);
 464         }
 465 }
 466
 467
 468 void CavesV6::makeTunnel(bool dirswitch)
 469 {
 470         if (dirswitch && !large_cave) {
 471                 main_direction.Z = ((float)(ps->next() % 20) - (float)10) / 10;
 472                 main_direction.Y = ((float)(ps->next() % 20) - (float)10) / 30;
 473                 main_direction.X = ((float)(ps->next() % 20) - (float)10) / 10;
 474
 475                 main_direction *= (float)ps->range(0, 10) / 10;
 476         }
 477
 478         // Randomize size
 479         s16 min_d = min_tunnel_diameter;
 480         s16 max_d = max_tunnel_diameter;
 481         rs = ps->range(min_d, max_d);
 482         s16 rs_part_max_length_rs = rs * part_max_length_rs;
 483
 484         v3s16 maxlen;
 485         if (large_cave) {
 486                 maxlen = v3s16(
 487                         rs_part_max_length_rs,
 488                         rs_part_max_length_rs / 2,
 489                         rs_part_max_length_rs
 490                 );
 491         } else {
 492                 maxlen = v3s16(
 493                         rs_part_max_length_rs,
 494                         ps->range(1, rs_part_max_length_rs),
 495                         rs_part_max_length_rs
 496                 );
 497         }
 498
 499         v3f vec;
 500         vec.Z = (float)(ps->next() % maxlen.Z) - (float)maxlen.Z / 2;
 501         vec.Y = (float)(ps->next() % maxlen.Y) - (float)maxlen.Y / 2;
 502         vec.X = (float)(ps->next() % maxlen.X) - (float)maxlen.X / 2;
 503
 504         // Jump downward sometimes
 505         if (!large_cave && ps->range(0, 12) == 0) {
 506                 vec.Z = (float)(ps->next() % maxlen.Z) - (float)maxlen.Z / 2;
 507                 vec.Y = (float)(ps->next() % (maxlen.Y * 2)) - (float)maxlen.Y;
 508                 vec.X = (float)(ps->next() % maxlen.X) - (float)maxlen.X / 2;
 509         }
 510
 511         // Do not make caves that are entirely above ground, to fix shadow bugs
 512         // caused by overgenerated large caves.
 513         // It is only necessary to check the startpoint and endpoint.
 514         v3s16 p1 = v3s16(orp.X, orp.Y, orp.Z) + of + rs / 2;
 515         v3s16 p2 = v3s16(vec.X, vec.Y, vec.Z) + p1;
 516
 517         // If startpoint and endpoint are above ground, disable placement of nodes
 518         // in carveRoute while still running all PseudoRandom calls to ensure caves
 519         // are consistent with existing worlds.
 520         bool tunnel_above_ground =
 521                 p1.Y > getSurfaceFromHeightmap(p1) &&
 522                 p2.Y > getSurfaceFromHeightmap(p2);
 523
 524         vec += main_direction;
 525
 526         v3f rp = orp + vec;
 527         if (rp.X < 0)
 528                 rp.X = 0;
 529         else if (rp.X >= ar.X)
 530                 rp.X = ar.X - 1;
 531
 532         if (rp.Y < route_y_min)
 533                 rp.Y = route_y_min;
 534         else if (rp.Y >= route_y_max)
 535                 rp.Y = route_y_max - 1;
 536
 537         if (rp.Z < 0)
 538                 rp.Z = 0;
 539         else if (rp.Z >= ar.Z)
 540                 rp.Z = ar.Z - 1;
 541
 542         vec = rp - orp;
 543
 544         float veclen = vec.getLength();
 545         // As odd as it sounds, veclen is *exactly* 0.0 sometimes, causing a FPE
 546         if (veclen < 0.05)
 547                 veclen = 1.0;
 548
 549         // Every second section is rough
 550         bool randomize_xz = (ps2->range(1, 2) == 1);
 551
 552         // Carve routes
 553         for (float f = 0; f < 1.0; f += 1.0 / veclen)
 554                 carveRoute(vec, f, randomize_xz, tunnel_above_ground);
 555
 556         orp = rp;
 557 }
 558
 559
 560 void CavesV6::carveRoute(v3f vec, float f, bool randomize_xz,
 561         bool tunnel_above_ground)
 562 {
 563         MapNode airnode(CONTENT_AIR);
 564         MapNode waternode(c_water_source);
 565         MapNode lavanode(c_lava_source);
 566
 567         v3s16 startp(orp.X, orp.Y, orp.Z);
 568         startp += of;
 569
 570         v3f fp = orp + vec * f;
 571         fp.X += 0.1 * ps->range(-10, 10);
 572         fp.Z += 0.1 * ps->range(-10, 10);
 573         v3s16 cp(fp.X, fp.Y, fp.Z);
 574
 575         s16 d0 = -rs / 2;
 576         s16 d1 = d0 + rs;
 577         if (randomize_xz) {
 578                 d0 += ps->range(-1, 1);
 579                 d1 += ps->range(-1, 1);
 580         }
 581
 582         for (s16 z0 = d0; z0 <= d1; z0++) {
 583                 s16 si = rs / 2 - MYMAX(0, abs(z0) - rs / 7 - 1);
 584                 for (s16 x0 = -si - ps->range(0,1); x0 <= si - 1 + ps->range(0,1); x0++) {
 585                         if (tunnel_above_ground)
 586                                 continue;
 587
 588                         s16 maxabsxz = MYMAX(abs(x0), abs(z0));
 589                         s16 si2 = rs / 2 - MYMAX(0, maxabsxz - rs / 7 - 1);
 590                         for (s16 y0 = -si2; y0 <= si2; y0++) {
 591                                 if (large_cave_is_flat) {
 592                                         // Make large caves not so tall
 593                                         if (rs > 7 && abs(y0) >= rs / 3)
 594                                                 continue;
 595                                 }
 596
 597                                 v3s16 p(cp.X + x0, cp.Y + y0, cp.Z + z0);
 598                                 p += of;
 599
 600                                 if (vm->m_area.contains(p) == false)
 601                                         continue;
 602
 603                                 u32 i = vm->m_area.index(p);
 604                                 content_t c = vm->m_data[i].getContent();
 605                                 if (!ndef->get(c).is_ground_content)
 606                                         continue;
 607
 608                                 if (large_cave) {
 609                                         int full_ymin = node_min.Y - MAP_BLOCKSIZE;
 610                                         int full_ymax = node_max.Y + MAP_BLOCKSIZE;
 611
 612                                         if (full_ymin < water_level && full_ymax > water_level) {
 613                                                 vm->m_data[i] = (p.Y <= water_level) ? waternode : airnode;
 614                                         } else if (full_ymax < water_level) {
 615                                                 vm->m_data[i] = (p.Y < startp.Y - 2) ? lavanode : airnode;
 616                                         } else {
 617                                                 vm->m_data[i] = airnode;
 618                                         }
 619                                 } else {
 620                                         if (c == CONTENT_IGNORE || c == CONTENT_AIR)
 621                                                 continue;
 622
 623                                         vm->m_data[i] = airnode;
 624                                         vm->m_flags[i] |= VMANIP_FLAG_CAVE;
 625                                 }
 626                         }
 627                 }
 628         }
 629 }
 630
 631
 632 inline s16 CavesV6::getSurfaceFromHeightmap(v3s16 p)
 633 {
 634         if (heightmap != NULL &&
 635                         p.Z >= node_min.Z && p.Z <= node_max.Z &&
 636                         p.X >= node_min.X && p.X <= node_max.X) {
 637                 u32 index = (p.Z - node_min.Z) * ystride + (p.X - node_min.X);
 638                 return heightmap[index];
 639         } else {
 640                 return water_level;
 641         }
 642 }