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 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         // Carve routes
 197         for (float f = 0; f < 1.0; f += 1.0 / veclen)
 198                 carveRoute(vec, f, randomize_xz);
 199
 200         orp = rp;
 201 }
 202
 203
 204 void CaveV5::carveRoute(v3f vec, float f, bool randomize_xz) {
 205         MapNode airnode(CONTENT_AIR);
 206         MapNode waternode(c_water_source);
 207         MapNode lavanode(c_lava_source);
 208
 209         v3s16 startp(orp.X, orp.Y, orp.Z);
 210         startp += of;
 211
 212         float nval = NoisePerlin3D(np_caveliquids, startp.X,
 213                 startp.Y, startp.Z, mg->seed);
 214         MapNode liquidnode = nval < 0.40 ? lavanode : waternode;
 215
 216         v3f fp = orp + vec * f;
 217         fp.X += 0.1 * ps->range(-10, 10);
 218         fp.Z += 0.1 * ps->range(-10, 10);
 219         v3s16 cp(fp.X, fp.Y, fp.Z);
 220
 221         s16 d0 = -rs/2;
 222         s16 d1 = d0 + rs;
 223         if (randomize_xz) {
 224                 d0 += ps->range(-1, 1);
 225                 d1 += ps->range(-1, 1);
 226         }
 227
 228         for (s16 z0 = d0; z0 <= d1; z0++) {
 229                 s16 si = rs / 2 - MYMAX(0, abs(z0) - rs / 7 - 1);
 230                 for (s16 x0 = -si - ps->range(0,1); x0 <= si - 1 + ps->range(0,1); x0++) {
 231                         s16 maxabsxz = MYMAX(abs(x0), abs(z0));
 232
 233                         s16 si2 = rs / 2 - MYMAX(0, maxabsxz - rs / 7 - 1);
 234
 235                         for (s16 y0 = -si2; y0 <= si2; y0++) {
 236                                 if (large_cave_is_flat) {
 237                                         // Make large caves not so tall
 238                                         if (rs > 7 && abs(y0) >= rs / 3)
 239                                                 continue;
 240                                 }
 241
 242                                 v3s16 p(cp.X + x0, cp.Y + y0, cp.Z + z0);
 243                                 p += of;
 244
 245                                 if (vm->m_area.contains(p) == false)
 246                                         continue;
 247
 248                                 u32 i = vm->m_area.index(p);
 249                                 content_t c = vm->m_data[i].getContent();
 250                                 if (!ndef->get(c).is_ground_content)
 251                                         continue;
 252
 253                                 int full_ymin = node_min.Y - MAP_BLOCKSIZE;
 254                                 int full_ymax = node_max.Y + MAP_BLOCKSIZE;
 255
 256                                 if (flooded && full_ymin < water_level &&
 257                                                 full_ymax > water_level)
 258                                         vm->m_data[i] = (p.Y <= water_level) ?
 259                                                 waternode : airnode;
 260                                 else if (flooded && full_ymax < water_level)
 261                                         vm->m_data[i] = (p.Y < startp.Y - 4) ?
 262                                                 liquidnode : airnode;
 263                                 else
 264                                         vm->m_data[i] = airnode;
 265                         }
 266                 }
 267         }
 268 }
 269
 270
 271 ///////////////////////////////////////// Caves V6
 272
 273
 274 CaveV6::CaveV6(MapgenV6 *mg, PseudoRandom *ps, PseudoRandom *ps2, bool is_large_cave) {
 275         this->mg             = mg;
 276         this->vm             = mg->vm;
 277         this->ndef           = mg->ndef;
 278         this->water_level    = mg->water_level;
 279         this->large_cave     = is_large_cave;
 280         this->ps             = ps;
 281         this->ps2            = ps2;
 282         this->c_water_source = mg->c_water_source;
 283         this->c_lava_source  = mg->c_lava_source;
 284
 285         min_tunnel_diameter = 2;
 286         max_tunnel_diameter = ps->range(2, 6);
 287         dswitchint          = ps->range(1, 14);
 288         flooded             = true;
 289
 290         if (large_cave) {
 291                 part_max_length_rs  = ps->range(2,4);
 292                 tunnel_routepoints  = ps->range(5, ps->range(15,30));
 293                 min_tunnel_diameter = 5;
 294                 max_tunnel_diameter = ps->range(7, ps->range(8,24));
 295         } else {
 296                 part_max_length_rs = ps->range(2,9);
 297                 tunnel_routepoints = ps->range(10, ps->range(15,30));
 298         }
 299
 300         large_cave_is_flat = (ps->range(0,1) == 0);
 301 }
 302
 303
 304 void CaveV6::makeCave(v3s16 nmin, v3s16 nmax, int max_stone_height) {
 305         node_min = nmin;
 306         node_max = nmax;
 307         max_stone_y = max_stone_height;
 308         main_direction = v3f(0, 0, 0);
 309
 310         // Allowed route area size in nodes
 311         ar = node_max - node_min + v3s16(1, 1, 1);
 312         // Area starting point in nodes
 313         of = node_min;
 314
 315         // Allow a bit more
 316         //(this should be more than the maximum radius of the tunnel)
 317         const s16 max_spread_amount = MAP_BLOCKSIZE;
 318         s16 insure = 10;
 319         s16 more = MYMAX(max_spread_amount - max_tunnel_diameter / 2 - insure, 1);
 320         ar += v3s16(1,0,1) * more * 2;
 321         of -= v3s16(1,0,1) * more;
 322
 323         route_y_min = 0;
 324         // Allow half a diameter + 7 over stone surface
 325         route_y_max = -of.Y + max_stone_y + max_tunnel_diameter / 2 + 7;
 326
 327         // Limit maximum to area
 328         route_y_max = rangelim(route_y_max, 0, ar.Y - 1);
 329
 330         if (large_cave) {
 331                 s16 min = 0;
 332                 if (node_min.Y < water_level && node_max.Y > water_level) {
 333                         min = water_level - max_tunnel_diameter/3 - of.Y;
 334                         route_y_max = water_level + max_tunnel_diameter/3 - of.Y;
 335                 }
 336                 route_y_min = ps->range(min, min + max_tunnel_diameter);
 337                 route_y_min = rangelim(route_y_min, 0, route_y_max);
 338         }
 339
 340         s16 route_start_y_min = route_y_min;
 341         s16 route_start_y_max = route_y_max;
 342
 343         route_start_y_min = rangelim(route_start_y_min, 0, ar.Y-1);
 344         route_start_y_max = rangelim(route_start_y_max, route_start_y_min, ar.Y-1);
 345
 346         // Randomize starting position
 347         orp = v3f(
 348                 (float)(ps->next() % ar.X) + 0.5,
 349                 (float)(ps->range(route_start_y_min, route_start_y_max)) + 0.5,
 350                 (float)(ps->next() % ar.Z) + 0.5
 351         );
 352
 353         // Add generation notify begin event
 354         v3s16 abs_pos(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 355         GenNotifyType notifytype = large_cave ?
 356                 GENNOTIFY_LARGECAVE_BEGIN : GENNOTIFY_CAVE_BEGIN;
 357         mg->gennotify.addEvent(notifytype, abs_pos);
 358
 359         // Generate some tunnel starting from orp
 360         for (u16 j = 0; j < tunnel_routepoints; j++)
 361                 makeTunnel(j % dswitchint == 0);
 362
 363         // Add generation notify end event
 364         abs_pos = v3s16(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 365         notifytype = large_cave ?
 366                 GENNOTIFY_LARGECAVE_END : GENNOTIFY_CAVE_END;
 367         mg->gennotify.addEvent(notifytype, abs_pos);
 368 }
 369
 370
 371 void CaveV6::makeTunnel(bool dirswitch) {
 372         if (dirswitch && !large_cave) {
 373                 main_direction = v3f(
 374                         ((float)(ps->next() % 20) - (float)10) / 10,
 375                         ((float)(ps->next() % 20) - (float)10) / 30,
 376                         ((float)(ps->next() % 20) - (float)10) / 10
 377                 );
 378                 main_direction *= (float)ps->range(0, 10) / 10;
 379         }
 380
 381         // Randomize size
 382         s16 min_d = min_tunnel_diameter;
 383         s16 max_d = max_tunnel_diameter;
 384         rs = ps->range(min_d, max_d);
 385         s16 rs_part_max_length_rs = rs * part_max_length_rs;
 386
 387         v3s16 maxlen;
 388         if (large_cave) {
 389                 maxlen = v3s16(
 390                         rs_part_max_length_rs,
 391                         rs_part_max_length_rs / 2,
 392                         rs_part_max_length_rs
 393                 );
 394         } else {
 395                 maxlen = v3s16(
 396                         rs_part_max_length_rs,
 397                         ps->range(1, rs_part_max_length_rs),
 398                         rs_part_max_length_rs
 399                 );
 400         }
 401
 402         v3f vec(
 403                 (float)(ps->next() % maxlen.X) - (float)maxlen.X / 2,
 404                 (float)(ps->next() % maxlen.Y) - (float)maxlen.Y / 2,
 405                 (float)(ps->next() % maxlen.Z) - (float)maxlen.Z / 2
 406         );
 407
 408         // Jump downward sometimes
 409         if (!large_cave && ps->range(0, 12) == 0) {
 410                 vec = v3f(
 411                         (float)(ps->next() % maxlen.X) - (float)maxlen.X / 2,
 412                         (float)(ps->next() % (maxlen.Y * 2)) - (float)maxlen.Y,
 413                         (float)(ps->next() % maxlen.Z) - (float)maxlen.Z / 2
 414                 );
 415         }
 416
 417         // Do not make caves that are entirely above ground.
 418         // It is only necessary to check the startpoint and endpoint.
 419         v3s16 orpi(orp.X, orp.Y, orp.Z);
 420         v3s16 veci(vec.X, vec.Y, vec.Z);
 421         s16 h1;
 422         s16 h2;
 423
 424         v3s16 p1 = orpi + veci + of + rs / 2;
 425         if (p1.Z >= node_min.Z && p1.Z <= node_max.Z &&
 426                         p1.X >= node_min.X && p1.X <= node_max.X) {
 427                 u32 index1 = (p1.Z - node_min.Z) * mg->ystride +
 428                         (p1.X - node_min.X);
 429                 h1 = mg->heightmap[index1];
 430         } else {
 431                 h1 = water_level; // If not in heightmap
 432         }
 433
 434         v3s16 p2 = orpi + of + rs / 2;
 435         if (p2.Z >= node_min.Z && p2.Z <= node_max.Z &&
 436                         p2.X >= node_min.X && p2.X <= node_max.X) {
 437                 u32 index2 = (p2.Z - node_min.Z) * mg->ystride +
 438                         (p2.X - node_min.X);
 439                 h2 = mg->heightmap[index2];
 440         } else {
 441                 h2 = water_level;
 442         }
 443
 444         if (p1.Y > h1 && p2.Y > h2) // If startpoint and endpoint are above ground
 445                 return;
 446
 447         vec += main_direction;
 448
 449         v3f rp = orp + vec;
 450         if (rp.X < 0)
 451                 rp.X = 0;
 452         else if (rp.X >= ar.X)
 453                 rp.X = ar.X - 1;
 454
 455         if (rp.Y < route_y_min)
 456                 rp.Y = route_y_min;
 457         else if (rp.Y >= route_y_max)
 458                 rp.Y = route_y_max - 1;
 459
 460         if (rp.Z < 0)
 461                 rp.Z = 0;
 462         else if (rp.Z >= ar.Z)
 463                 rp.Z = ar.Z - 1;
 464
 465         vec = rp - orp;
 466
 467         float veclen = vec.getLength();
 468         // As odd as it sounds, veclen is *exactly* 0.0 sometimes, causing a FPE
 469         if (veclen < 0.05)
 470                 veclen = 1.0;
 471
 472         // Every second section is rough
 473         bool randomize_xz = (ps2->range(1, 2) == 1);
 474
 475         // Carve routes
 476         for (float f = 0; f < 1.0; f += 1.0 / veclen)
 477                 carveRoute(vec, f, randomize_xz);
 478
 479         orp = rp;
 480 }
 481
 482
 483 void CaveV6::carveRoute(v3f vec, float f, bool randomize_xz) {
 484         MapNode airnode(CONTENT_AIR);
 485         MapNode waternode(c_water_source);
 486         MapNode lavanode(c_lava_source);
 487
 488         v3s16 startp(orp.X, orp.Y, orp.Z);
 489         startp += of;
 490
 491         v3f fp = orp + vec * f;
 492         fp.X += 0.1 * ps->range(-10, 10);
 493         fp.Z += 0.1 * ps->range(-10, 10);
 494         v3s16 cp(fp.X, fp.Y, fp.Z);
 495
 496         s16 d0 = -rs/2;
 497         s16 d1 = d0 + rs;
 498         if (randomize_xz) {
 499                 d0 += ps->range(-1, 1);
 500                 d1 += ps->range(-1, 1);
 501         }
 502
 503         for (s16 z0 = d0; z0 <= d1; z0++) {
 504                 s16 si = rs / 2 - MYMAX(0, abs(z0) - rs / 7 - 1);
 505                 for (s16 x0 = -si - ps->range(0,1); x0 <= si - 1 + ps->range(0,1); x0++) {
 506                         s16 maxabsxz = MYMAX(abs(x0), abs(z0));
 507                         s16 si2 = rs / 2 - MYMAX(0, maxabsxz - rs / 7 - 1);
 508                         for (s16 y0 = -si2; y0 <= si2; y0++) {
 509                                 if (large_cave_is_flat) {
 510                                         // Make large caves not so tall
 511                                         if (rs > 7 && abs(y0) >= rs / 3)
 512                                                 continue;
 513                                 }
 514
 515                                 v3s16 p(cp.X + x0, cp.Y + y0, cp.Z + z0);
 516                                 p += of;
 517
 518                                 if (vm->m_area.contains(p) == false)
 519                                         continue;
 520
 521                                 u32 i = vm->m_area.index(p);
 522                                 content_t c = vm->m_data[i].getContent();
 523                                 if (!ndef->get(c).is_ground_content)
 524                                         continue;
 525
 526                                 if (large_cave) {
 527                                         int full_ymin = node_min.Y - MAP_BLOCKSIZE;
 528                                         int full_ymax = node_max.Y + MAP_BLOCKSIZE;
 529
 530                                         if (flooded && full_ymin < water_level &&
 531                                                         full_ymax > water_level) {
 532                                                 vm->m_data[i] = (p.Y <= water_level) ?
 533                                                         waternode : airnode;
 534                                         } else if (flooded && full_ymax < water_level) {
 535                                                 vm->m_data[i] = (p.Y < startp.Y - 2) ?
 536                                                         lavanode : airnode;
 537                                         } else {
 538                                                 vm->m_data[i] = airnode;
 539                                         }
 540                                 } else {
 541                                         if (c == CONTENT_IGNORE || c == CONTENT_AIR)
 542                                                 continue;
 543
 544                                         vm->m_data[i] = airnode;
 545                                         vm->m_flags[i] |= VMANIP_FLAG_CAVE;
 546                                 }
 547                         }
 548                 }
 549         }
 550 }
 551
 552
 553 ///////////////////////////////////////// Caves V7
 554
 555
 556 CaveV7::CaveV7(MapgenV7 *mg, PseudoRandom *ps) {
 557         this->mg             = mg;
 558         this->vm             = mg->vm;
 559         this->ndef           = mg->ndef;
 560         this->water_level    = mg->water_level;
 561         this->ps             = ps;
 562         this->c_water_source = mg->c_water_source;
 563         this->c_lava_source  = mg->c_lava_source;
 564         this->c_ice          = mg->c_ice;
 565         this->np_caveliquids = &nparams_caveliquids;
 566
 567         dswitchint = ps->range(1, 14);
 568         flooded    = ps->range(1, 2) == 2;
 569
 570         part_max_length_rs  = ps->range(2, 4);
 571         tunnel_routepoints  = ps->range(5, ps->range(15, 30));
 572         min_tunnel_diameter = 5;
 573         max_tunnel_diameter = ps->range(7, ps->range(8, 24));
 574
 575         large_cave_is_flat = (ps->range(0, 1) == 0);
 576 }
 577
 578
 579 void CaveV7::makeCave(v3s16 nmin, v3s16 nmax, int max_stone_height) {
 580         node_min = nmin;
 581         node_max = nmax;
 582         max_stone_y = max_stone_height;
 583         main_direction = v3f(0, 0, 0);
 584
 585         // Allowed route area size in nodes
 586         ar = node_max - node_min + v3s16(1, 1, 1);
 587         // Area starting point in nodes
 588         of = node_min;
 589
 590         // Allow a bit more
 591         //(this should be more than the maximum radius of the tunnel)
 592         s16 insure = 10;
 593         s16 more = MYMAX(MAP_BLOCKSIZE - max_tunnel_diameter / 2 - insure, 1);
 594         ar += v3s16(1,0,1) * more * 2;
 595         of -= v3s16(1,0,1) * more;
 596
 597         route_y_min = 0;
 598         // Allow half a diameter + 7 over stone surface
 599         route_y_max = -of.Y + max_stone_y + max_tunnel_diameter / 2 + 7;
 600
 601         // Limit maximum to area
 602         route_y_max = rangelim(route_y_max, 0, ar.Y - 1);
 603
 604         s16 min = 0;
 605         if (node_min.Y < water_level && node_max.Y > water_level) {
 606                 min = water_level - max_tunnel_diameter/3 - of.Y;
 607                 route_y_max = water_level + max_tunnel_diameter/3 - of.Y;
 608         }
 609         route_y_min = ps->range(min, min + max_tunnel_diameter);
 610         route_y_min = rangelim(route_y_min, 0, route_y_max);
 611
 612         s16 route_start_y_min = route_y_min;
 613         s16 route_start_y_max = route_y_max;
 614
 615         route_start_y_min = rangelim(route_start_y_min, 0, ar.Y - 1);
 616         route_start_y_max = rangelim(route_start_y_max, route_start_y_min, ar.Y - 1);
 617
 618         // Randomize starting position
 619         orp = v3f(
 620                 (float)(ps->next() % ar.X) + 0.5,
 621                 (float)(ps->range(route_start_y_min, route_start_y_max)) + 0.5,
 622                 (float)(ps->next() % ar.Z) + 0.5
 623         );
 624
 625         // Add generation notify begin event
 626         v3s16 abs_pos(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 627         GenNotifyType notifytype = GENNOTIFY_LARGECAVE_BEGIN;
 628         mg->gennotify.addEvent(notifytype, abs_pos);
 629
 630         // Generate some tunnel starting from orp
 631         for (u16 j = 0; j < tunnel_routepoints; j++)
 632                 makeTunnel(j % dswitchint == 0);
 633
 634         // Add generation notify end event
 635         abs_pos = v3s16(of.X + orp.X, of.Y + orp.Y, of.Z + orp.Z);
 636         notifytype = GENNOTIFY_LARGECAVE_END;
 637         mg->gennotify.addEvent(notifytype, abs_pos);
 638 }
 639
 640
 641 void CaveV7::makeTunnel(bool dirswitch) {
 642
 643         // Randomize size
 644         s16 min_d = min_tunnel_diameter;
 645         s16 max_d = max_tunnel_diameter;
 646         rs = ps->range(min_d, max_d);
 647         s16 rs_part_max_length_rs = rs * part_max_length_rs;
 648
 649         v3s16 maxlen;
 650         maxlen = v3s16(
 651                 rs_part_max_length_rs,
 652                 rs_part_max_length_rs / 2,
 653                 rs_part_max_length_rs
 654         );
 655
 656         v3f vec;
 657         // Jump downward sometimes
 658         vec = v3f(
 659                 (float)(ps->next() % maxlen.X) - (float)maxlen.X / 2,
 660                 (float)(ps->next() % maxlen.Y) - (float)maxlen.Y / 2,
 661                 (float)(ps->next() % maxlen.Z) - (float)maxlen.Z / 2
 662         );
 663
 664         // Do not make caves that are above ground.
 665         // It is only necessary to check the startpoint and endpoint.
 666         v3s16 orpi(orp.X, orp.Y, orp.Z);
 667         v3s16 veci(vec.X, vec.Y, vec.Z);
 668         v3s16 p;
 669
 670         p = orpi + veci + of + rs / 2;
 671         if (p.Z >= node_min.Z && p.Z <= node_max.Z &&
 672                         p.X >= node_min.X && p.X <= node_max.X) {
 673                 u32 index = (p.Z - node_min.Z) * mg->ystride + (p.X - node_min.X);
 674                 s16 h = mg->ridge_heightmap[index];
 675                 if (h < p.Y)
 676                         return;
 677         } else if (p.Y > water_level) {
 678                 return; // If it's not in our heightmap, use a simple heuristic
 679         }
 680
 681         p = orpi + of + rs / 2;
 682         if (p.Z >= node_min.Z && p.Z <= node_max.Z &&
 683                         p.X >= node_min.X && p.X <= node_max.X) {
 684                 u32 index = (p.Z - node_min.Z) * mg->ystride + (p.X - node_min.X);
 685                 s16 h = mg->ridge_heightmap[index];
 686                 if (h < p.Y)
 687                         return;
 688         } else if (p.Y > water_level) {
 689                 return;
 690         }
 691
 692         vec += main_direction;
 693
 694         v3f rp = orp + vec;
 695         if (rp.X < 0)
 696                 rp.X = 0;
 697         else if (rp.X >= ar.X)
 698                 rp.X = ar.X - 1;
 699
 700         if (rp.Y < route_y_min)
 701                 rp.Y = route_y_min;
 702         else if (rp.Y >= route_y_max)
 703                 rp.Y = route_y_max - 1;
 704
 705         if (rp.Z < 0)
 706                 rp.Z = 0;
 707         else if (rp.Z >= ar.Z)
 708                 rp.Z = ar.Z - 1;
 709
 710         vec = rp - orp;
 711
 712         float veclen = vec.getLength();
 713         if (veclen < 0.05)
 714                 veclen = 1.0;
 715
 716         // Every second section is rough
 717         bool randomize_xz = (ps->range(1, 2) == 1);
 718
 719         // Carve routes
 720         for (float f = 0; f < 1.0; f += 1.0 / veclen)
 721                 carveRoute(vec, f, randomize_xz);
 722
 723         orp = rp;
 724 }
 725
 726
 727 void CaveV7::carveRoute(v3f vec, float f, bool randomize_xz) {
 728         MapNode airnode(CONTENT_AIR);
 729         MapNode waternode(c_water_source);
 730         MapNode lavanode(c_lava_source);
 731
 732         v3s16 startp(orp.X, orp.Y, orp.Z);
 733         startp += of;
 734
 735         float nval = NoisePerlin3D(np_caveliquids, startp.X,
 736                 startp.Y, startp.Z, mg->seed);
 737         MapNode liquidnode = (nval < 0.40 && node_max.Y < MGV7_LAVA_DEPTH) ?
 738                 lavanode : waternode;
 739
 740         v3f fp = orp + vec * f;
 741         fp.X += 0.1 * ps->range(-10, 10);
 742         fp.Z += 0.1 * ps->range(-10, 10);
 743         v3s16 cp(fp.X, fp.Y, fp.Z);
 744
 745         s16 d0 = -rs/2;
 746         s16 d1 = d0 + rs;
 747         if (randomize_xz) {
 748                 d0 += ps->range(-1, 1);
 749                 d1 += ps->range(-1, 1);
 750         }
 751
 752         for (s16 z0 = d0; z0 <= d1; z0++) {
 753                 s16 si = rs / 2 - MYMAX(0, abs(z0) - rs / 7 - 1);
 754                 for (s16 x0 = -si - ps->range(0,1); x0 <= si - 1 + ps->range(0,1); x0++) {
 755                         s16 maxabsxz = MYMAX(abs(x0), abs(z0));
 756
 757                         s16 si2 = rs / 2 - MYMAX(0, maxabsxz - rs / 7 - 1);
 758
 759                         for (s16 y0 = -si2; y0 <= si2; y0++) {
 760                                 if (large_cave_is_flat) {
 761                                         // Make large caves not so tall
 762                                         if (rs > 7 && abs(y0) >= rs / 3)
 763                                                 continue;
 764                                 }
 765
 766                                 v3s16 p(cp.X + x0, cp.Y + y0, cp.Z + z0);
 767                                 p += of;
 768
 769                                 if (vm->m_area.contains(p) == false)
 770                                         continue;
 771
 772                                 u32 i = vm->m_area.index(p);
 773                                 content_t c = vm->m_data[i].getContent();
 774                                 if (!ndef->get(c).is_ground_content)
 775                                         continue;
 776
 777                                 int full_ymin = node_min.Y - MAP_BLOCKSIZE;
 778                                 int full_ymax = node_max.Y + MAP_BLOCKSIZE;
 779
 780                                 if (flooded && full_ymin < water_level &&
 781                                                 full_ymax > water_level)
 782                                         vm->m_data[i] = (p.Y <= water_level) ?
 783                                                 waternode : airnode;
 784                                 else if (flooded && full_ymax < water_level)
 785                                         vm->m_data[i] = (p.Y < startp.Y - 4) ?
 786                                                 liquidnode : airnode;
 787                                 else
 788                                         vm->m_data[i] = airnode;
 789                         }
 790                 }
 791         }
 792 }
 793