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 ///////////////////////////////////////// Caves V5
  32
  33
  34 CaveV5::CaveV5(MapgenV5 *mg, PseudoRandom *ps) {
  35         this->mg   = mg;
  36         this->vm   = mg->vm;
  37         this->ndef = mg->ndef;
  38         this->water_level = mg->water_level;
  39         this->ps = ps;
  40         this->c_water_source = mg->c_water_source;
  41         this->c_lava_source  = mg->c_lava_source;
  42         this->c_ice          = mg->c_ice;
  43         this->np_caveliquids = &nparams_caveliquids;
  44
  45         dswitchint = ps->range(1, 14);
  46         flooded    = ps->range(1, 2) == 2;
  47
  48         part_max_length_rs = ps->range(2, 4);
  49         tunnel_routepoints = ps->range(5, ps->range(15, 30));
  50         min_tunnel_diameter = 5;
  51         max_tunnel_diameter = ps->range(7, ps->range(8, 24));
  52
  53         large_cave_is_flat = (ps->range(0, 1) == 0);
  54 }
  55
  56
  57 void CaveV5::makeCave(v3s16 nmin, v3s16 nmax, int max_stone_height) {
  58         node_min = nmin;
  59         node_max = nmax;
  60         main_direction = v3f(0, 0, 0);
  61
  62         // Allowed route area size in nodes
  63         ar = node_max - node_min + v3s16(1, 1, 1);
  64         // Area starting point in nodes
  65         of = node_min;
  66
  67         // Allow a bit more
  68         //(this should be more than the maximum radius of the tunnel)
  69         s16 insure = 10;
  70         s16 more = MYMAX(MAP_BLOCKSIZE - max_tunnel_diameter / 2 - insure, 1);
  71         ar += v3s16(1,0,1) * more * 2;
  72         of -= v3s16(1,0,1) * more;
  73
  74         route_y_min = 0;
  75         // Allow half a diameter + 7 over stone surface
  76         route_y_max = -of.Y + max_stone_y + max_tunnel_diameter / 2 + 7;
  77
  78         // Limit maximum to area
  79         route_y_max = rangelim(route_y_max, 0, ar.Y - 1);
  80
  81         s16 min = 0;
  82                 if (node_min.Y < water_level && node_max.Y > water_level) {
  83                         min = water_level - max_tunnel_diameter/3 - of.Y;
  84                         route_y_max = water_level + max_tunnel_diameter/3 - of.Y;
  85                 }
  86         route_y_min = ps->range(min, min + max_tunnel_diameter);
  87         route_y_min = rangelim(route_y_min, 0, route_y_max);
  88
  89         s16 route_start_y_min = route_y_min;
  90         s16 route_start_y_max = route_y_max;
  91
  92         route_start_y_min = rangelim(route_start_y_min, 0, ar.Y - 1);
  93         route_start_y_max = rangelim(route_start_y_max, route_start_y_min, ar.Y - 1);
  94
  95         // Randomize starting position
  96         orp = v3f(
  97                 (float)(ps->next() % ar.X) + 0.5,
  98                 (float)(ps->range(route_start_y_min, route_start_y_max)) + 0.5,
  99                 (float)(ps->next() % ar.Z) + 0.5
 100         );
 101
 102         // Add generation notify begin event
 103         v3s16 abs_pos(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 104         GenNotifyType notifytype = GENNOTIFY_LARGECAVE_BEGIN;
 105         mg->gennotify.addEvent(notifytype, abs_pos);
 106
 107         // Generate some tunnel starting from orp
 108         for (u16 j = 0; j < tunnel_routepoints; j++)
 109                 makeTunnel(j % dswitchint == 0);
 110
 111         // Add generation notify end event
 112         abs_pos = v3s16(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 113         notifytype = GENNOTIFY_LARGECAVE_END;
 114         mg->gennotify.addEvent(notifytype, abs_pos);
 115 }
 116
 117
 118 void CaveV5::makeTunnel(bool dirswitch) {
 119
 120         // Randomize size
 121         s16 min_d = min_tunnel_diameter;
 122         s16 max_d = max_tunnel_diameter;
 123         rs = ps->range(min_d, max_d);
 124         s16 rs_part_max_length_rs = rs * part_max_length_rs;
 125
 126         v3s16 maxlen;
 127         maxlen = v3s16(
 128                 rs_part_max_length_rs,
 129                 rs_part_max_length_rs / 2,
 130                 rs_part_max_length_rs
 131         );
 132
 133         v3f vec;
 134         // Jump downward sometimes
 135         vec = v3f(
 136                 (float)(ps->next() % maxlen.X) - (float)maxlen.X / 2,
 137                 (float)(ps->next() % maxlen.Y) - (float)maxlen.Y / 2,
 138                 (float)(ps->next() % maxlen.Z) - (float)maxlen.Z / 2
 139         );
 140
 141         // Do not make large caves that are above ground.
 142         // It is only necessary to check the startpoint and endpoint.
 143         v3s16 orpi(orp.X, orp.Y, orp.Z);
 144         v3s16 veci(vec.X, vec.Y, vec.Z);
 145         v3s16 p;
 146
 147         p = orpi + veci + of + rs / 2;
 148         if (p.Z >= node_min.Z && p.Z <= node_max.Z &&
 149                         p.X >= node_min.X && p.X <= node_max.X) {
 150                 u32 index = (p.Z - node_min.Z) * mg->ystride + (p.X - node_min.X);
 151                 s16 h = mg->heightmap[index];
 152                 if (h < p.Y)
 153                         return;
 154         } else if (p.Y > water_level) {
 155                 return; // If it's not in our heightmap, use a simple heuristic
 156         }
 157
 158         p = orpi + of + rs / 2;
 159         if (p.Z >= node_min.Z && p.Z <= node_max.Z &&
 160                 p.X >= node_min.X && p.X <= node_max.X) {
 161                 u32 index = (p.Z - node_min.Z) * mg->ystride + (p.X - node_min.X);
 162                 s16 h = mg->heightmap[index];
 163                 if (h < p.Y)
 164                         return;
 165         } else if (p.Y > water_level) {
 166                 return;
 167         }
 168
 169         vec += main_direction;
 170
 171         v3f rp = orp + vec;
 172         if (rp.X < 0)
 173                 rp.X = 0;
 174         else if (rp.X >= ar.X)
 175                 rp.X = ar.X - 1;
 176
 177         if (rp.Y < route_y_min)
 178                 rp.Y = route_y_min;
 179         else if (rp.Y >= route_y_max)
 180                 rp.Y = route_y_max - 1;
 181
 182         if (rp.Z < 0)
 183                 rp.Z = 0;
 184         else if (rp.Z >= ar.Z)
 185                 rp.Z = ar.Z - 1;
 186
 187         vec = rp - orp;
 188
 189         float veclen = vec.getLength();
 190         if (veclen < 0.05)
 191                 veclen = 1.0;
 192
 193         // Every second section is rough
 194         bool randomize_xz = (ps->range(1, 2) == 1);
 195
 196         // Make a ravine every once in a while if it's long enough
 197         //float xylen = vec.X * vec.X + vec.Z * vec.Z;
 198         //disable ravines for now
 199         bool is_ravine = false; //(xylen > 500.0) && !large_cave && (ps->range(1, 8) == 1);
 200
 201         // Carve routes
 202         for (float f = 0; f < 1.0; f += 1.0 / veclen)
 203                 carveRoute(vec, f, randomize_xz, is_ravine);
 204
 205         orp = rp;
 206 }
 207
 208
 209 void CaveV5::carveRoute(v3f vec, float f, bool randomize_xz, bool is_ravine) {
 210         MapNode airnode(CONTENT_AIR);
 211         MapNode waternode(c_water_source);
 212         MapNode lavanode(c_lava_source);
 213
 214         v3s16 startp(orp.X, orp.Y, orp.Z);
 215         startp += of;
 216
 217         float nval = NoisePerlin3D(np_caveliquids, startp.X,
 218                                                         startp.Y, startp.Z, mg->seed);
 219         MapNode liquidnode = nval < 0.40 ? lavanode : waternode;
 220
 221         v3f fp = orp + vec * f;
 222         fp.X += 0.1 * ps->range(-10, 10);
 223         fp.Z += 0.1 * ps->range(-10, 10);
 224         v3s16 cp(fp.X, fp.Y, fp.Z);
 225
 226         s16 d0 = -rs/2;
 227         s16 d1 = d0 + rs;
 228         if (randomize_xz) {
 229                 d0 += ps->range(-1, 1);
 230                 d1 += ps->range(-1, 1);
 231         }
 232
 233         bool should_make_cave_hole = ps->range(1, 10) == 1;
 234
 235         for (s16 z0 = d0; z0 <= d1; z0++) {
 236                 s16 si = rs / 2 - MYMAX(0, abs(z0) - rs / 7 - 1);
 237                 for (s16 x0 = -si - ps->range(0,1); x0 <= si - 1 + ps->range(0,1); x0++) {
 238                         s16 maxabsxz = MYMAX(abs(x0), abs(z0));
 239
 240                         s16 si2 = is_ravine ? MYMIN(ps->range(25, 26), ar.Y) :
 241                                                                  rs / 2 - MYMAX(0, maxabsxz - rs / 7 - 1);
 242
 243                         for (s16 y0 = -si2; y0 <= si2; y0++) {
 244                                 if (large_cave_is_flat) {
 245                                         // Make large caves not so tall
 246                                         if (rs > 7 && abs(y0) >= rs / 3)
 247                                                 continue;
 248                                 }
 249
 250                                 v3s16 p(cp.X + x0, cp.Y + y0, cp.Z + z0);
 251                                 p += of;
 252
 253                                 if (!is_ravine && mg->heightmap && should_make_cave_hole &&
 254                                         p.X <= node_max.X && p.Z <= node_max.Z) {
 255                                         int maplen = node_max.X - node_min.X + 1;
 256                                         int idx = (p.Z - node_min.Z) * maplen + (p.X - node_min.X);
 257                                         if (p.Y >= mg->heightmap[idx] - 2)
 258                                                 continue;
 259                                 }
 260
 261                                 if (vm->m_area.contains(p) == false)
 262                                         continue;
 263
 264                                 u32 i = vm->m_area.index(p);
 265
 266                                 // Don't replace air, water, lava, or ice
 267                                 content_t c = vm->m_data[i].getContent();
 268                                 if (!ndef->get(c).is_ground_content || c == CONTENT_AIR ||
 269                                         c == c_water_source || c == c_lava_source || c == c_ice)
 270                                         continue;
 271
 272                                 int full_ymin = node_min.Y - MAP_BLOCKSIZE;
 273                                 int full_ymax = node_max.Y + MAP_BLOCKSIZE;
 274
 275                                 if (flooded && full_ymin < water_level && full_ymax > water_level)
 276                                         vm->m_data[i] = (p.Y <= water_level) ? waternode : airnode;
 277                                 else if (flooded && full_ymax < water_level)
 278                                         vm->m_data[i] = (p.Y < startp.Y - 4) ? liquidnode : airnode;
 279                                 else
 280                                         vm->m_data[i] = airnode;
 281                         }
 282                 }
 283         }
 284 }
 285
 286
 287 ///////////////////////////////////////// Caves V6
 288
 289
 290 CaveV6::CaveV6(MapgenV6 *mg, PseudoRandom *ps, PseudoRandom *ps2, bool is_large_cave) {
 291         this->mg   = mg;
 292         this->vm   = mg->vm;
 293         this->ndef = mg->ndef;
 294         this->water_level = mg->water_level;
 295         this->large_cave = is_large_cave;
 296         this->ps  = ps;
 297         this->ps2 = ps2;
 298         this->c_water_source = mg->c_water_source;
 299         this->c_lava_source  = mg->c_lava_source;
 300
 301         min_tunnel_diameter = 2;
 302         max_tunnel_diameter = ps->range(2, 6);
 303         dswitchint = ps->range(1, 14);
 304         flooded = true;
 305
 306         if (large_cave) {
 307                 part_max_length_rs = ps->range(2,4);
 308                 tunnel_routepoints = ps->range(5, ps->range(15,30));
 309                 min_tunnel_diameter = 5;
 310                 max_tunnel_diameter = ps->range(7, ps->range(8,24));
 311         } else {
 312                 part_max_length_rs = ps->range(2,9);
 313                 tunnel_routepoints = ps->range(10, ps->range(15,30));
 314         }
 315
 316         large_cave_is_flat = (ps->range(0,1) == 0);
 317 }
 318
 319
 320 void CaveV6::makeCave(v3s16 nmin, v3s16 nmax, int max_stone_height) {
 321         node_min = nmin;
 322         node_max = nmax;
 323         max_stone_y = max_stone_height;
 324         main_direction = v3f(0, 0, 0);
 325
 326         // Allowed route area size in nodes
 327         ar = node_max - node_min + v3s16(1, 1, 1);
 328         // Area starting point in nodes
 329         of = node_min;
 330
 331         // Allow a bit more
 332         //(this should be more than the maximum radius of the tunnel)
 333         const s16 max_spread_amount = MAP_BLOCKSIZE;
 334         s16 insure = 10;
 335         s16 more = MYMAX(max_spread_amount - max_tunnel_diameter / 2 - insure, 1);
 336         ar += v3s16(1,0,1) * more * 2;
 337         of -= v3s16(1,0,1) * more;
 338
 339         route_y_min = 0;
 340         // Allow half a diameter + 7 over stone surface
 341         route_y_max = -of.Y + max_stone_y + max_tunnel_diameter / 2 + 7;
 342
 343         // Limit maximum to area
 344         route_y_max = rangelim(route_y_max, 0, ar.Y - 1);
 345
 346         if (large_cave) {
 347                 s16 min = 0;
 348                 if (node_min.Y < water_level && node_max.Y > water_level) {
 349                         min = water_level - max_tunnel_diameter/3 - of.Y;
 350                         route_y_max = water_level + max_tunnel_diameter/3 - of.Y;
 351                 }
 352                 route_y_min = ps->range(min, min + max_tunnel_diameter);
 353                 route_y_min = rangelim(route_y_min, 0, route_y_max);
 354         }
 355
 356         s16 route_start_y_min = route_y_min;
 357         s16 route_start_y_max = route_y_max;
 358
 359         route_start_y_min = rangelim(route_start_y_min, 0, ar.Y-1);
 360         route_start_y_max = rangelim(route_start_y_max, route_start_y_min, ar.Y-1);
 361
 362         // Randomize starting position
 363         orp = v3f(
 364                 (float)(ps->next() % ar.X) + 0.5,
 365                 (float)(ps->range(route_start_y_min, route_start_y_max)) + 0.5,
 366                 (float)(ps->next() % ar.Z) + 0.5
 367         );
 368
 369         // Add generation notify begin event
 370         v3s16 abs_pos(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 371         GenNotifyType notifytype = large_cave ?
 372                 GENNOTIFY_LARGECAVE_BEGIN : GENNOTIFY_CAVE_BEGIN;
 373         mg->gennotify.addEvent(notifytype, abs_pos);
 374
 375         // Generate some tunnel starting from orp
 376         for (u16 j = 0; j < tunnel_routepoints; j++)
 377                 makeTunnel(j % dswitchint == 0);
 378
 379         // Add generation notify end event
 380         abs_pos = v3s16(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 381         notifytype = large_cave ?
 382                 GENNOTIFY_LARGECAVE_END : GENNOTIFY_CAVE_END;
 383         mg->gennotify.addEvent(notifytype, abs_pos);
 384 }
 385
 386
 387 void CaveV6::makeTunnel(bool dirswitch) {
 388         if (dirswitch && !large_cave) {
 389                 main_direction = v3f(
 390                         ((float)(ps->next() % 20) - (float)10) / 10,
 391                         ((float)(ps->next() % 20) - (float)10) / 30,
 392                         ((float)(ps->next() % 20) - (float)10) / 10
 393                 );
 394                 main_direction *= (float)ps->range(0, 10) / 10;
 395         }
 396
 397         // Randomize size
 398         s16 min_d = min_tunnel_diameter;
 399         s16 max_d = max_tunnel_diameter;
 400         rs = ps->range(min_d, max_d);
 401         s16 rs_part_max_length_rs = rs * part_max_length_rs;
 402
 403         v3s16 maxlen;
 404         if (large_cave) {
 405                 maxlen = v3s16(
 406                         rs_part_max_length_rs,
 407                         rs_part_max_length_rs / 2,
 408                         rs_part_max_length_rs
 409                 );
 410         } else {
 411                 maxlen = v3s16(
 412                         rs_part_max_length_rs,
 413                         ps->range(1, rs_part_max_length_rs),
 414                         rs_part_max_length_rs
 415                 );
 416         }
 417
 418         v3f vec(
 419                 (float)(ps->next() % maxlen.X) - (float)maxlen.X / 2,
 420                 (float)(ps->next() % maxlen.Y) - (float)maxlen.Y / 2,
 421                 (float)(ps->next() % maxlen.Z) - (float)maxlen.Z / 2
 422         );
 423
 424         // Jump downward sometimes
 425         if (!large_cave && ps->range(0, 12) == 0) {
 426                 vec = v3f(
 427                         (float)(ps->next() % maxlen.X) - (float)maxlen.X / 2,
 428                         (float)(ps->next() % (maxlen.Y * 2)) - (float)maxlen.Y,
 429                         (float)(ps->next() % maxlen.Z) - (float)maxlen.Z / 2
 430                 );
 431         }
 432
 433         // Do not make large caves that are entirely above ground.
 434         // It is only necessary to check the startpoint and endpoint.
 435         if (large_cave) {
 436                 v3s16 orpi(orp.X, orp.Y, orp.Z);
 437                 v3s16 veci(vec.X, vec.Y, vec.Z);
 438                 s16 h1;
 439                 s16 h2;
 440
 441                 v3s16 p1 = orpi + veci + of + rs / 2;
 442                 if (p1.Z >= node_min.Z && p1.Z <= node_max.Z &&
 443                                 p1.X >= node_min.X && p1.X <= node_max.X) {
 444                         u32 index1 = (p1.Z - node_min.Z) * mg->ystride + (p1.X - node_min.X);
 445                         h1 = mg->heightmap[index1];
 446                 } else {
 447                         h1 = water_level; // If not in heightmap
 448                 }
 449
 450                 v3s16 p2 = orpi + of + rs / 2;
 451                 if (p2.Z >= node_min.Z && p2.Z <= node_max.Z &&
 452                                 p2.X >= node_min.X && p2.X <= node_max.X) {
 453                         u32 index2 = (p2.Z - node_min.Z) * mg->ystride + (p2.X - node_min.X);
 454                         h2 = mg->heightmap[index2];
 455                 } else {
 456                         h2 = water_level;
 457                 }
 458
 459                 if (p1.Y > h1 && p2.Y > h2) // If startpoint and endpoint are above ground
 460                         return;
 461         }
 462
 463         vec += main_direction;
 464
 465         v3f rp = orp + vec;
 466         if (rp.X < 0)
 467                 rp.X = 0;
 468         else if (rp.X >= ar.X)
 469                 rp.X = ar.X - 1;
 470
 471         if (rp.Y < route_y_min)
 472                 rp.Y = route_y_min;
 473         else if (rp.Y >= route_y_max)
 474                 rp.Y = route_y_max - 1;
 475
 476         if (rp.Z < 0)
 477                 rp.Z = 0;
 478         else if (rp.Z >= ar.Z)
 479                 rp.Z = ar.Z - 1;
 480
 481         vec = rp - orp;
 482
 483         float veclen = vec.getLength();
 484         // As odd as it sounds, veclen is *exactly* 0.0 sometimes, causing a FPE
 485         if (veclen < 0.05)
 486                 veclen = 1.0;
 487
 488         // Every second section is rough
 489         bool randomize_xz = (ps2->range(1, 2) == 1);
 490
 491         // Carve routes
 492         for (float f = 0; f < 1.0; f += 1.0 / veclen)
 493                 carveRoute(vec, f, randomize_xz);
 494
 495         orp = rp;
 496 }
 497
 498
 499 void CaveV6::carveRoute(v3f vec, float f, bool randomize_xz) {
 500         MapNode airnode(CONTENT_AIR);
 501         MapNode waternode(c_water_source);
 502         MapNode lavanode(c_lava_source);
 503
 504         v3s16 startp(orp.X, orp.Y, orp.Z);
 505         startp += of;
 506
 507         v3f fp = orp + vec * f;
 508         fp.X += 0.1 * ps->range(-10, 10);
 509         fp.Z += 0.1 * ps->range(-10, 10);
 510         v3s16 cp(fp.X, fp.Y, fp.Z);
 511
 512         s16 d0 = -rs/2;
 513         s16 d1 = d0 + rs;
 514         if (randomize_xz) {
 515                 d0 += ps->range(-1, 1);
 516                 d1 += ps->range(-1, 1);
 517         }
 518
 519         for (s16 z0 = d0; z0 <= d1; z0++) {
 520                 s16 si = rs / 2 - MYMAX(0, abs(z0) - rs / 7 - 1);
 521                 for (s16 x0 = -si - ps->range(0,1); x0 <= si - 1 + ps->range(0,1); x0++) {
 522                         s16 maxabsxz = MYMAX(abs(x0), abs(z0));
 523                         s16 si2 = rs / 2 - MYMAX(0, maxabsxz - rs / 7 - 1);
 524                         for (s16 y0 = -si2; y0 <= si2; y0++) {
 525                                 if (large_cave_is_flat) {
 526                                         // Make large caves not so tall
 527                                         if (rs > 7 && abs(y0) >= rs / 3)
 528                                                 continue;
 529                                 }
 530
 531                                 v3s16 p(cp.X + x0, cp.Y + y0, cp.Z + z0);
 532                                 p += of;
 533
 534                                 if (vm->m_area.contains(p) == false)
 535                                         continue;
 536
 537                                 u32 i = vm->m_area.index(p);
 538                                 content_t c = vm->m_data[i].getContent();
 539                                 if (!ndef->get(c).is_ground_content)
 540                                         continue;
 541
 542                                 if (large_cave) {
 543                                         int full_ymin = node_min.Y - MAP_BLOCKSIZE;
 544                                         int full_ymax = node_max.Y + MAP_BLOCKSIZE;
 545
 546                                         if (flooded && full_ymin < water_level && full_ymax > water_level) {
 547                                                 vm->m_data[i] = (p.Y <= water_level) ? waternode : airnode;
 548                                         } else if (flooded && full_ymax < water_level) {
 549                                                 vm->m_data[i] = (p.Y < startp.Y - 2) ? lavanode : airnode;
 550                                         } else {
 551                                                 vm->m_data[i] = airnode;
 552                                         }
 553                                 } else {
 554                                         // Don't replace air or water or lava or ignore
 555                                         if (c == CONTENT_IGNORE || c == CONTENT_AIR ||
 556                                                 c == c_water_source || c == c_lava_source)
 557                                                 continue;
 558
 559                                         vm->m_data[i] = airnode;
 560                                         vm->m_flags[i] |= VMANIP_FLAG_CAVE;
 561                                 }
 562                         }
 563                 }
 564         }
 565 }
 566
 567
 568 ///////////////////////////////////////// Caves V7
 569
 570
 571 CaveV7::CaveV7(MapgenV7 *mg, PseudoRandom *ps) {
 572         this->mg   = mg;
 573         this->vm   = mg->vm;
 574         this->ndef = mg->ndef;
 575         this->water_level = mg->water_level;
 576         this->ps = ps;
 577         this->c_water_source = mg->c_water_source;
 578         this->c_lava_source  = mg->c_lava_source;
 579         this->c_ice          = mg->c_ice;
 580         this->np_caveliquids = &nparams_caveliquids;
 581
 582         dswitchint = ps->range(1, 14);
 583         flooded    = ps->range(1, 2) == 2;
 584
 585         part_max_length_rs = ps->range(2, 4);
 586         tunnel_routepoints = ps->range(5, ps->range(15, 30));
 587         min_tunnel_diameter = 5;
 588         max_tunnel_diameter = ps->range(7, ps->range(8, 24));
 589
 590         large_cave_is_flat = (ps->range(0, 1) == 0);
 591 }
 592
 593
 594 void CaveV7::makeCave(v3s16 nmin, v3s16 nmax, int max_stone_height) {
 595         node_min = nmin;
 596         node_max = nmax;
 597         max_stone_y = max_stone_height;
 598         main_direction = v3f(0, 0, 0);
 599
 600         // Allowed route area size in nodes
 601         ar = node_max - node_min + v3s16(1, 1, 1);
 602         // Area starting point in nodes
 603         of = node_min;
 604
 605         // Allow a bit more
 606         //(this should be more than the maximum radius of the tunnel)
 607         s16 insure = 10;
 608         s16 more = MYMAX(MAP_BLOCKSIZE - max_tunnel_diameter / 2 - insure, 1);
 609         ar += v3s16(1,0,1) * more * 2;
 610         of -= v3s16(1,0,1) * more;
 611
 612         route_y_min = 0;
 613         // Allow half a diameter + 7 over stone surface
 614         route_y_max = -of.Y + max_stone_y + max_tunnel_diameter / 2 + 7;
 615
 616         // Limit maximum to area
 617         route_y_max = rangelim(route_y_max, 0, ar.Y - 1);
 618
 619         s16 min = 0;
 620         if (node_min.Y < water_level && node_max.Y > water_level) {
 621                 min = water_level - max_tunnel_diameter/3 - of.Y;
 622                 route_y_max = water_level + max_tunnel_diameter/3 - of.Y;
 623         }
 624         route_y_min = ps->range(min, min + max_tunnel_diameter);
 625         route_y_min = rangelim(route_y_min, 0, route_y_max);
 626
 627         s16 route_start_y_min = route_y_min;
 628         s16 route_start_y_max = route_y_max;
 629
 630         route_start_y_min = rangelim(route_start_y_min, 0, ar.Y - 1);
 631         route_start_y_max = rangelim(route_start_y_max, route_start_y_min, ar.Y - 1);
 632
 633         // Randomize starting position
 634         orp = v3f(
 635                 (float)(ps->next() % ar.X) + 0.5,
 636                 (float)(ps->range(route_start_y_min, route_start_y_max)) + 0.5,
 637                 (float)(ps->next() % ar.Z) + 0.5
 638         );
 639
 640         // Add generation notify begin event
 641         v3s16 abs_pos(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 642         GenNotifyType notifytype = GENNOTIFY_LARGECAVE_BEGIN;
 643         mg->gennotify.addEvent(notifytype, abs_pos);
 644
 645         // Generate some tunnel starting from orp
 646         for (u16 j = 0; j < tunnel_routepoints; j++)
 647                 makeTunnel(j % dswitchint == 0);
 648
 649         // Add generation notify end event
 650         abs_pos = v3s16(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 651         notifytype = GENNOTIFY_LARGECAVE_END;
 652         mg->gennotify.addEvent(notifytype, abs_pos);
 653 }
 654
 655
 656 void CaveV7::makeTunnel(bool dirswitch) {
 657
 658         // Randomize size
 659         s16 min_d = min_tunnel_diameter;
 660         s16 max_d = max_tunnel_diameter;
 661         rs = ps->range(min_d, max_d);
 662         s16 rs_part_max_length_rs = rs * part_max_length_rs;
 663
 664         v3s16 maxlen;
 665         maxlen = v3s16(
 666                 rs_part_max_length_rs,
 667                 rs_part_max_length_rs / 2,
 668                 rs_part_max_length_rs
 669         );
 670
 671         v3f vec;
 672         // Jump downward sometimes
 673         vec = v3f(
 674                 (float)(ps->next() % maxlen.X) - (float)maxlen.X / 2,
 675                 (float)(ps->next() % maxlen.Y) - (float)maxlen.Y / 2,
 676                 (float)(ps->next() % maxlen.Z) - (float)maxlen.Z / 2
 677         );
 678
 679         // Do not make large caves that are above ground.
 680         // It is only necessary to check the startpoint and endpoint.
 681         v3s16 orpi(orp.X, orp.Y, orp.Z);
 682         v3s16 veci(vec.X, vec.Y, vec.Z);
 683         v3s16 p;
 684
 685         p = orpi + veci + of + rs / 2;
 686         if (p.Z >= node_min.Z && p.Z <= node_max.Z &&
 687                         p.X >= node_min.X && p.X <= node_max.X) {
 688                 u32 index = (p.Z - node_min.Z) * mg->ystride + (p.X - node_min.X);
 689                 s16 h = mg->ridge_heightmap[index];
 690                 if (h < p.Y)
 691                         return;
 692         } else if (p.Y > water_level) {
 693                 return; // If it's not in our heightmap, use a simple heuristic
 694         }
 695
 696         p = orpi + of + rs / 2;
 697         if (p.Z >= node_min.Z && p.Z <= node_max.Z &&
 698                         p.X >= node_min.X && p.X <= node_max.X) {
 699                 u32 index = (p.Z - node_min.Z) * mg->ystride + (p.X - node_min.X);
 700                 s16 h = mg->ridge_heightmap[index];
 701                 if (h < p.Y)
 702                         return;
 703         } else if (p.Y > water_level) {
 704                 return;
 705         }
 706
 707         vec += main_direction;
 708
 709         v3f rp = orp + vec;
 710         if (rp.X < 0)
 711                 rp.X = 0;
 712         else if (rp.X >= ar.X)
 713                 rp.X = ar.X - 1;
 714
 715         if (rp.Y < route_y_min)
 716                 rp.Y = route_y_min;
 717         else if (rp.Y >= route_y_max)
 718                 rp.Y = route_y_max - 1;
 719
 720         if (rp.Z < 0)
 721                 rp.Z = 0;
 722         else if (rp.Z >= ar.Z)
 723                 rp.Z = ar.Z - 1;
 724
 725         vec = rp - orp;
 726
 727         float veclen = vec.getLength();
 728         if (veclen < 0.05)
 729                 veclen = 1.0;
 730
 731         // Every second section is rough
 732         bool randomize_xz = (ps->range(1, 2) == 1);
 733
 734         // Make a ravine every once in a while if it's long enough
 735         //float xylen = vec.X * vec.X + vec.Z * vec.Z;
 736         //disable ravines for now
 737         bool is_ravine = false; //(xylen > 500.0) && !large_cave && (ps->range(1, 8) == 1);
 738
 739         // Carve routes
 740         for (float f = 0; f < 1.0; f += 1.0 / veclen)
 741                 carveRoute(vec, f, randomize_xz, is_ravine);
 742
 743         orp = rp;
 744 }
 745
 746
 747 void CaveV7::carveRoute(v3f vec, float f, bool randomize_xz, bool is_ravine) {
 748         MapNode airnode(CONTENT_AIR);
 749         MapNode waternode(c_water_source);
 750         MapNode lavanode(c_lava_source);
 751
 752         v3s16 startp(orp.X, orp.Y, orp.Z);
 753         startp += of;
 754
 755         float nval = NoisePerlin3D(np_caveliquids, startp.X,
 756                                                         startp.Y, startp.Z, mg->seed);
 757         MapNode liquidnode = (nval < 0.40 && node_max.Y < -256) ? lavanode : waternode;
 758
 759         v3f fp = orp + vec * f;
 760         fp.X += 0.1 * ps->range(-10, 10);
 761         fp.Z += 0.1 * ps->range(-10, 10);
 762         v3s16 cp(fp.X, fp.Y, fp.Z);
 763
 764         s16 d0 = -rs/2;
 765         s16 d1 = d0 + rs;
 766         if (randomize_xz) {
 767                 d0 += ps->range(-1, 1);
 768                 d1 += ps->range(-1, 1);
 769         }
 770
 771         bool should_make_cave_hole = ps->range(1, 10) == 1;
 772
 773         for (s16 z0 = d0; z0 <= d1; z0++) {
 774                 s16 si = rs / 2 - MYMAX(0, abs(z0) - rs / 7 - 1);
 775                 for (s16 x0 = -si - ps->range(0,1); x0 <= si - 1 + ps->range(0,1); x0++) {
 776                         s16 maxabsxz = MYMAX(abs(x0), abs(z0));
 777
 778                         s16 si2 = is_ravine ? MYMIN(ps->range(25, 26), ar.Y) :
 779                                                                  rs / 2 - MYMAX(0, maxabsxz - rs / 7 - 1);
 780
 781                         for (s16 y0 = -si2; y0 <= si2; y0++) {
 782                                 if (large_cave_is_flat) {
 783                                         // Make large caves not so tall
 784                                         if (rs > 7 && abs(y0) >= rs / 3)
 785                                                 continue;
 786                                 }
 787
 788                                 v3s16 p(cp.X + x0, cp.Y + y0, cp.Z + z0);
 789                                 p += of;
 790
 791                                 if (!is_ravine && mg->heightmap && should_make_cave_hole &&
 792                                         p.X <= node_max.X && p.Z <= node_max.Z) {
 793                                         int maplen = node_max.X - node_min.X + 1;
 794                                         int idx = (p.Z - node_min.Z) * maplen + (p.X - node_min.X);
 795                                         if (p.Y >= mg->heightmap[idx] - 2)
 796                                                 continue;
 797                                 }
 798
 799                                 if (vm->m_area.contains(p) == false)
 800                                         continue;
 801
 802                                 u32 i = vm->m_area.index(p);
 803
 804                                 // Don't replace air, water, lava, or ice
 805                                 content_t c = vm->m_data[i].getContent();
 806                                 if (!ndef->get(c).is_ground_content || c == CONTENT_AIR ||
 807                                         c == c_water_source || c == c_lava_source || c == c_ice)
 808                                         continue;
 809
 810                                 int full_ymin = node_min.Y - MAP_BLOCKSIZE;
 811                                 int full_ymax = node_max.Y + MAP_BLOCKSIZE;
 812
 813                                 if (flooded && full_ymin < water_level && full_ymax > water_level)
 814                                         vm->m_data[i] = (p.Y <= water_level) ? waternode : airnode;
 815                                 else if (flooded && full_ymax < water_level)
 816                                         vm->m_data[i] = (p.Y < startp.Y - 4) ? liquidnode : airnode;
 817                                 else
 818                                         vm->m_data[i] = airnode;
 819                         }
 820                 }
 821         }
 822 }
 823