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