src/mapgen_v6.cpp

   1 /*
   2 Minetest
   3 Copyright (C) 2010-2015 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
  21 #include "mapgen.h"
  22 #include "voxel.h"
  23 #include "noise.h"
  24 #include "mapblock.h"
  25 #include "mapnode.h"
  26 #include "map.h"
  27 //#include "serverobject.h"
  28 #include "content_sao.h"
  29 #include "nodedef.h"
  30 #include "voxelalgorithms.h"
  31 //#include "profiler.h" // For TimeTaker
  32 #include "settings.h" // For g_settings
  33 #include "emerge.h"
  34 #include "dungeongen.h"
  35 #include "cavegen.h"
  36 #include "treegen.h"
  37 #include "mg_ore.h"
  38 #include "mg_decoration.h"
  39 #include "mapgen_v6.h"
  40
  41
  42 FlagDesc flagdesc_mapgen_v6[] = {
  43         {"jungles",    MGV6_JUNGLES},
  44         {"biomeblend", MGV6_BIOMEBLEND},
  45         {"mudflow",    MGV6_MUDFLOW},
  46         {"snowbiomes", MGV6_SNOWBIOMES},
  47         {"flat",       MGV6_FLAT},
  48         {"trees",      MGV6_TREES},
  49         {NULL,         0}
  50 };
  51
  52
  53 /////////////////////////////////////////////////////////////////////////////
  54
  55
  56 MapgenV6::MapgenV6(int mapgenid, MapgenV6Params *params, EmergeManager *emerge)
  57         : Mapgen(mapgenid, params, emerge)
  58 {
  59         this->m_emerge = emerge;
  60         this->ystride = csize.X; //////fix this
  61
  62         this->heightmap = new s16[csize.X * csize.Z];
  63
  64         this->spflags     = params->spflags;
  65         this->freq_desert = params->freq_desert;
  66         this->freq_beach  = params->freq_beach;
  67
  68         np_cave        = &params->np_cave;
  69         np_humidity    = &params->np_humidity;
  70         np_trees       = &params->np_trees;
  71         np_apple_trees = &params->np_apple_trees;
  72
  73         //// Create noise objects
  74         noise_terrain_base   = new Noise(&params->np_terrain_base,   seed, csize.X, csize.Y);
  75         noise_terrain_higher = new Noise(&params->np_terrain_higher, seed, csize.X, csize.Y);
  76         noise_steepness      = new Noise(&params->np_steepness,      seed, csize.X, csize.Y);
  77         noise_height_select  = new Noise(&params->np_height_select,  seed, csize.X, csize.Y);
  78         noise_mud            = new Noise(&params->np_mud,            seed, csize.X, csize.Y);
  79         noise_beach          = new Noise(&params->np_beach,          seed, csize.X, csize.Y);
  80         noise_biome          = new Noise(&params->np_biome,          seed,
  81                         csize.X + 2 * MAP_BLOCKSIZE, csize.Y + 2 * MAP_BLOCKSIZE);
  82         noise_humidity       = new Noise(&params->np_humidity,       seed,
  83                         csize.X + 2 * MAP_BLOCKSIZE, csize.Y + 2 * MAP_BLOCKSIZE);
  84
  85         //// Resolve nodes to be used
  86         INodeDefManager *ndef = emerge->ndef;
  87
  88         c_stone           = ndef->getId("mapgen_stone");
  89         c_dirt            = ndef->getId("mapgen_dirt");
  90         c_dirt_with_grass = ndef->getId("mapgen_dirt_with_grass");
  91         c_sand            = ndef->getId("mapgen_sand");
  92         c_water_source    = ndef->getId("mapgen_water_source");
  93         c_lava_source     = ndef->getId("mapgen_lava_source");
  94         c_gravel          = ndef->getId("mapgen_gravel");
  95         c_desert_stone    = ndef->getId("mapgen_desert_stone");
  96         c_desert_sand     = ndef->getId("mapgen_desert_sand");
  97         c_dirt_with_snow  = ndef->getId("mapgen_dirt_with_snow");
  98         c_snow            = ndef->getId("mapgen_snow");
  99         c_snowblock       = ndef->getId("mapgen_snowblock");
 100         c_ice             = ndef->getId("mapgen_ice");
 101
 102         if (c_gravel == CONTENT_IGNORE)
 103                 c_gravel = c_stone;
 104         if (c_desert_stone == CONTENT_IGNORE)
 105                 c_desert_stone = c_stone;
 106         if (c_desert_sand == CONTENT_IGNORE)
 107                 c_desert_sand = c_sand;
 108         if (c_dirt_with_snow == CONTENT_IGNORE)
 109                 c_dirt_with_snow = c_dirt_with_grass;
 110         if (c_snow == CONTENT_IGNORE)
 111                 c_snow = CONTENT_AIR;
 112         if (c_snowblock == CONTENT_IGNORE)
 113                 c_snowblock = c_dirt_with_grass;
 114         if (c_ice == CONTENT_IGNORE)
 115                 c_ice = c_water_source;
 116
 117         c_cobble             = ndef->getId("mapgen_cobble");
 118         c_mossycobble        = ndef->getId("mapgen_mossycobble");
 119         c_stair_cobble       = ndef->getId("mapgen_stair_cobble");
 120         c_stair_desert_stone = ndef->getId("mapgen_stair_desert_stone");
 121
 122         if (c_mossycobble == CONTENT_IGNORE)
 123                 c_mossycobble = c_cobble;
 124         if (c_stair_cobble == CONTENT_IGNORE)
 125                 c_stair_cobble = c_cobble;
 126         if (c_stair_desert_stone == CONTENT_IGNORE)
 127                 c_stair_desert_stone = c_desert_stone;
 128 }
 129
 130
 131 MapgenV6::~MapgenV6()
 132 {
 133         delete noise_terrain_base;
 134         delete noise_terrain_higher;
 135         delete noise_steepness;
 136         delete noise_height_select;
 137         delete noise_mud;
 138         delete noise_beach;
 139         delete noise_biome;
 140         delete noise_humidity;
 141
 142         delete[] heightmap;
 143 }
 144
 145
 146 MapgenV6Params::MapgenV6Params()
 147 {
 148         spflags = MGV6_JUNGLES | MGV6_SNOWBIOMES | MGV6_TREES |
 149                 MGV6_BIOMEBLEND | MGV6_MUDFLOW;
 150
 151         freq_desert = 0.45;
 152         freq_beach  = 0.15;
 153
 154         np_terrain_base   = NoiseParams(-4,   20.0, v3f(250.0, 250.0, 250.0), 82341,  5, 0.6,  2.0);
 155         np_terrain_higher = NoiseParams(20,   16.0, v3f(500.0, 500.0, 500.0), 85039,  5, 0.6,  2.0);
 156         np_steepness      = NoiseParams(0.85, 0.5,  v3f(125.0, 125.0, 125.0), -932,   5, 0.7,  2.0);
 157         np_height_select  = NoiseParams(0,    1.0,  v3f(250.0, 250.0, 250.0), 4213,   5, 0.69, 2.0);
 158         np_mud            = NoiseParams(4,    2.0,  v3f(200.0, 200.0, 200.0), 91013,  3, 0.55, 2.0);
 159         np_beach          = NoiseParams(0,    1.0,  v3f(250.0, 250.0, 250.0), 59420,  3, 0.50, 2.0);
 160         np_biome          = NoiseParams(0,    1.0,  v3f(500.0, 500.0, 500.0), 9130,   3, 0.50, 2.0);
 161         np_cave           = NoiseParams(6,    6.0,  v3f(250.0, 250.0, 250.0), 34329,  3, 0.50, 2.0);
 162         np_humidity       = NoiseParams(0.5,  0.5,  v3f(500.0, 500.0, 500.0), 72384,  3, 0.50, 2.0);
 163         np_trees          = NoiseParams(0,    1.0,  v3f(125.0, 125.0, 125.0), 2,      4, 0.66, 2.0);
 164         np_apple_trees    = NoiseParams(0,    1.0,  v3f(100.0, 100.0, 100.0), 342902, 3, 0.45, 2.0);
 165 }
 166
 167
 168 void MapgenV6Params::readParams(const Settings *settings)
 169 {
 170         settings->getFlagStrNoEx("mgv6_spflags", spflags, flagdesc_mapgen_v6);
 171         settings->getFloatNoEx("mgv6_freq_desert", freq_desert);
 172         settings->getFloatNoEx("mgv6_freq_beach",  freq_beach);
 173
 174         settings->getNoiseParams("mgv6_np_terrain_base",   np_terrain_base);
 175         settings->getNoiseParams("mgv6_np_terrain_higher", np_terrain_higher);
 176         settings->getNoiseParams("mgv6_np_steepness",      np_steepness);
 177         settings->getNoiseParams("mgv6_np_height_select",  np_height_select);
 178         settings->getNoiseParams("mgv6_np_mud",            np_mud);
 179         settings->getNoiseParams("mgv6_np_beach",          np_beach);
 180         settings->getNoiseParams("mgv6_np_biome",          np_biome);
 181         settings->getNoiseParams("mgv6_np_cave",           np_cave);
 182         settings->getNoiseParams("mgv6_np_humidity",       np_humidity);
 183         settings->getNoiseParams("mgv6_np_trees",          np_trees);
 184         settings->getNoiseParams("mgv6_np_apple_trees",    np_apple_trees);
 185 }
 186
 187
 188 void MapgenV6Params::writeParams(Settings *settings) const
 189 {
 190         settings->setFlagStr("mgv6_spflags", spflags, flagdesc_mapgen_v6, U32_MAX);
 191         settings->setFloat("mgv6_freq_desert", freq_desert);
 192         settings->setFloat("mgv6_freq_beach",  freq_beach);
 193
 194         settings->setNoiseParams("mgv6_np_terrain_base",   np_terrain_base);
 195         settings->setNoiseParams("mgv6_np_terrain_higher", np_terrain_higher);
 196         settings->setNoiseParams("mgv6_np_steepness",      np_steepness);
 197         settings->setNoiseParams("mgv6_np_height_select",  np_height_select);
 198         settings->setNoiseParams("mgv6_np_mud",            np_mud);
 199         settings->setNoiseParams("mgv6_np_beach",          np_beach);
 200         settings->setNoiseParams("mgv6_np_biome",          np_biome);
 201         settings->setNoiseParams("mgv6_np_cave",           np_cave);
 202         settings->setNoiseParams("mgv6_np_humidity",       np_humidity);
 203         settings->setNoiseParams("mgv6_np_trees",          np_trees);
 204         settings->setNoiseParams("mgv6_np_apple_trees",    np_apple_trees);
 205 }
 206
 207
 208 //////////////////////// Some helper functions for the map generator
 209
 210 // Returns Y one under area minimum if not found
 211 s16 MapgenV6::find_stone_level(v2s16 p2d)
 212 {
 213         v3s16 em = vm->m_area.getExtent();
 214         s16 y_nodes_max = vm->m_area.MaxEdge.Y;
 215         s16 y_nodes_min = vm->m_area.MinEdge.Y;
 216         u32 i = vm->m_area.index(p2d.X, y_nodes_max, p2d.Y);
 217         s16 y;
 218
 219         for (y = y_nodes_max; y >= y_nodes_min; y--) {
 220                 content_t c = vm->m_data[i].getContent();
 221                 if (c != CONTENT_IGNORE && (c == c_stone || c == c_desert_stone))
 222                         break;
 223
 224                 vm->m_area.add_y(em, i, -1);
 225         }
 226         return (y >= y_nodes_min) ? y : y_nodes_min - 1;
 227 }
 228
 229
 230 // Required by mapgen.h
 231 bool MapgenV6::block_is_underground(u64 seed, v3s16 blockpos)
 232 {
 233         /*s16 minimum_groundlevel = (s16)get_sector_minimum_ground_level(
 234                         seed, v2s16(blockpos.X, blockpos.Z));*/
 235         // Nah, this is just a heuristic, just return something
 236         s16 minimum_groundlevel = water_level;
 237
 238         if(blockpos.Y * MAP_BLOCKSIZE + MAP_BLOCKSIZE <= minimum_groundlevel)
 239                 return true;
 240         else
 241                 return false;
 242 }
 243
 244
 245 //////////////////////// Base terrain height functions
 246
 247 float MapgenV6::baseTerrainLevel(float terrain_base, float terrain_higher,
 248         float steepness, float height_select)
 249 {
 250         float base   = 1 + terrain_base;
 251         float higher = 1 + terrain_higher;
 252
 253         // Limit higher ground level to at least base
 254         if(higher < base)
 255                 higher = base;
 256
 257         // Steepness factor of cliffs
 258         float b = steepness;
 259         b = rangelim(b, 0.0, 1000.0);
 260         b = 5 * b * b * b * b * b * b * b;
 261         b = rangelim(b, 0.5, 1000.0);
 262
 263         // Values 1.5...100 give quite horrible looking slopes
 264         if (b > 1.5 && b < 100.0)
 265                 b = (b < 10.0) ? 1.5 : 100.0;
 266
 267         float a_off = -0.20; // Offset to more low
 268         float a = 0.5 + b * (a_off + height_select);
 269         a = rangelim(a, 0.0, 1.0); // Limit
 270
 271         return base * (1.0 - a) + higher * a;
 272 }
 273
 274
 275 float MapgenV6::baseTerrainLevelFromNoise(v2s16 p)
 276 {
 277         if (spflags & MGV6_FLAT)
 278                 return water_level;
 279
 280         float terrain_base   = NoisePerlin2D_PO(&noise_terrain_base->np,
 281                                                         p.X, 0.5, p.Y, 0.5, seed);
 282         float terrain_higher = NoisePerlin2D_PO(&noise_terrain_higher->np,
 283                                                         p.X, 0.5, p.Y, 0.5, seed);
 284         float steepness      = NoisePerlin2D_PO(&noise_steepness->np,
 285                                                         p.X, 0.5, p.Y, 0.5, seed);
 286         float height_select  = NoisePerlin2D_PO(&noise_height_select->np,
 287                                                         p.X, 0.5, p.Y, 0.5, seed);
 288
 289         return baseTerrainLevel(terrain_base, terrain_higher,
 290                                                         steepness, height_select);
 291 }
 292
 293
 294 float MapgenV6::baseTerrainLevelFromMap(v2s16 p)
 295 {
 296         int index = (p.Y - node_min.Z) * ystride + (p.X - node_min.X);
 297         return baseTerrainLevelFromMap(index);
 298 }
 299
 300
 301 float MapgenV6::baseTerrainLevelFromMap(int index)
 302 {
 303         if (spflags & MGV6_FLAT)
 304                 return water_level;
 305
 306         float terrain_base   = noise_terrain_base->result[index];
 307         float terrain_higher = noise_terrain_higher->result[index];
 308         float steepness      = noise_steepness->result[index];
 309         float height_select  = noise_height_select->result[index];
 310
 311         return baseTerrainLevel(terrain_base, terrain_higher,
 312                                                         steepness, height_select);
 313 }
 314
 315
 316 s16 MapgenV6::find_ground_level_from_noise(u64 seed, v2s16 p2d, s16 precision)
 317 {
 318         return baseTerrainLevelFromNoise(p2d) + MGV6_AVERAGE_MUD_AMOUNT;
 319 }
 320
 321
 322 int MapgenV6::getGroundLevelAtPoint(v2s16 p)
 323 {
 324         return baseTerrainLevelFromNoise(p) + MGV6_AVERAGE_MUD_AMOUNT;
 325 }
 326
 327
 328 int MapgenV6::getSpawnLevelAtPoint(v2s16 p)
 329 {
 330         s16 level_at_point = baseTerrainLevelFromNoise(p) + MGV6_AVERAGE_MUD_AMOUNT;
 331         if (level_at_point <= water_level ||
 332                         level_at_point > water_level + 16)
 333                 return MAX_MAP_GENERATION_LIMIT;  // Unsuitable spawn point
 334         else
 335                 return level_at_point;
 336 }
 337
 338
 339 //////////////////////// Noise functions
 340
 341 float MapgenV6::getMudAmount(v2s16 p)
 342 {
 343         int index = (p.Y - node_min.Z) * ystride + (p.X - node_min.X);
 344         return getMudAmount(index);
 345 }
 346
 347
 348 bool MapgenV6::getHaveBeach(v2s16 p)
 349 {
 350         int index = (p.Y - node_min.Z) * ystride + (p.X - node_min.X);
 351         return getHaveBeach(index);
 352 }
 353
 354
 355 BiomeV6Type MapgenV6::getBiome(v2s16 p)
 356 {
 357         int index = (p.Y - full_node_min.Z) * (ystride + 2 * MAP_BLOCKSIZE)
 358                         + (p.X - full_node_min.X);
 359         return getBiome(index, p);
 360 }
 361
 362
 363 float MapgenV6::getHumidity(v2s16 p)
 364 {
 365         /*double noise = noise2d_perlin(
 366                 0.5+(float)p.X/500, 0.5+(float)p.Y/500,
 367                 seed+72384, 4, 0.66);
 368         noise = (noise + 1.0)/2.0;*/
 369
 370         int index = (p.Y - full_node_min.Z) * (ystride + 2 * MAP_BLOCKSIZE)
 371                         + (p.X - full_node_min.X);
 372         float noise = noise_humidity->result[index];
 373
 374         if (noise < 0.0)
 375                 noise = 0.0;
 376         if (noise > 1.0)
 377                 noise = 1.0;
 378         return noise;
 379 }
 380
 381
 382 float MapgenV6::getTreeAmount(v2s16 p)
 383 {
 384         /*double noise = noise2d_perlin(
 385                         0.5+(float)p.X/125, 0.5+(float)p.Y/125,
 386                         seed+2, 4, 0.66);*/
 387
 388         float noise = NoisePerlin2D(np_trees, p.X, p.Y, seed);
 389         float zeroval = -0.39;
 390         if (noise < zeroval)
 391                 return 0;
 392         else
 393                 return 0.04 * (noise - zeroval) / (1.0 - zeroval);
 394 }
 395
 396
 397 bool MapgenV6::getHaveAppleTree(v2s16 p)
 398 {
 399         /*is_apple_tree = noise2d_perlin(
 400                 0.5+(float)p.X/100, 0.5+(float)p.Z/100,
 401                 data->seed+342902, 3, 0.45) > 0.2;*/
 402
 403         float noise = NoisePerlin2D(np_apple_trees, p.X, p.Y, seed);
 404
 405         return noise > 0.2;
 406 }
 407
 408
 409 float MapgenV6::getMudAmount(int index)
 410 {
 411         if (spflags & MGV6_FLAT)
 412                 return MGV6_AVERAGE_MUD_AMOUNT;
 413
 414         /*return ((float)AVERAGE_MUD_AMOUNT + 2.0 * noise2d_perlin(
 415                         0.5+(float)p.X/200, 0.5+(float)p.Y/200,
 416                         seed+91013, 3, 0.55));*/
 417
 418         return noise_mud->result[index];
 419 }
 420
 421
 422 bool MapgenV6::getHaveBeach(int index)
 423 {
 424         // Determine whether to have sand here
 425         /*double sandnoise = noise2d_perlin(
 426                         0.2+(float)p2d.X/250, 0.7+(float)p2d.Y/250,
 427                         seed+59420, 3, 0.50);*/
 428
 429         float sandnoise = noise_beach->result[index];
 430         return (sandnoise > freq_beach);
 431 }
 432
 433
 434 BiomeV6Type MapgenV6::getBiome(int index, v2s16 p)
 435 {
 436         // Just do something very simple as for now
 437         /*double d = noise2d_perlin(
 438                         0.6+(float)p2d.X/250, 0.2+(float)p2d.Y/250,
 439                         seed+9130, 3, 0.50);*/
 440
 441         float d = noise_biome->result[index];
 442         float h = noise_humidity->result[index];
 443
 444         if (spflags & MGV6_SNOWBIOMES) {
 445                 float blend = (spflags & MGV6_BIOMEBLEND) ? noise2d(p.X, p.Y, seed) / 40 : 0;
 446
 447                 if (d > MGV6_FREQ_HOT + blend) {
 448                         if (h > MGV6_FREQ_JUNGLE + blend)
 449                                 return BT_JUNGLE;
 450                         else
 451                                 return BT_DESERT;
 452                 } else if (d < MGV6_FREQ_SNOW + blend) {
 453                         if (h > MGV6_FREQ_TAIGA + blend)
 454                                 return BT_TAIGA;
 455                         else
 456                                 return BT_TUNDRA;
 457                 } else {
 458                         return BT_NORMAL;
 459                 }
 460         } else {
 461                 if (d > freq_desert)
 462                         return BT_DESERT;
 463
 464                 if ((spflags & MGV6_BIOMEBLEND) && (d > freq_desert - 0.10) &&
 465                                 ((noise2d(p.X, p.Y, seed) + 1.0) > (freq_desert - d) * 20.0))
 466                         return BT_DESERT;
 467
 468                 if ((spflags & MGV6_JUNGLES) && h > 0.75)
 469                         return BT_JUNGLE;
 470                 else
 471                         return BT_NORMAL;
 472         }
 473 }
 474
 475
 476 u32 MapgenV6::get_blockseed(u64 seed, v3s16 p)
 477 {
 478         s32 x = p.X, y = p.Y, z = p.Z;
 479         return (u32)(seed % 0x100000000ULL) + z * 38134234 + y * 42123 + x * 23;
 480 }
 481
 482
 483 //////////////////////// Map generator
 484
 485 void MapgenV6::makeChunk(BlockMakeData *data)
 486 {
 487         // Pre-conditions
 488         assert(data->vmanip);
 489         assert(data->nodedef);
 490         assert(data->blockpos_requested.X >= data->blockpos_min.X &&
 491                 data->blockpos_requested.Y >= data->blockpos_min.Y &&
 492                 data->blockpos_requested.Z >= data->blockpos_min.Z);
 493         assert(data->blockpos_requested.X <= data->blockpos_max.X &&
 494                 data->blockpos_requested.Y <= data->blockpos_max.Y &&
 495                 data->blockpos_requested.Z <= data->blockpos_max.Z);
 496
 497         this->generating = true;
 498         this->vm   = data->vmanip;
 499         this->ndef = data->nodedef;
 500
 501         // Hack: use minimum block coords for old code that assumes a single block
 502         v3s16 blockpos_min = data->blockpos_min;
 503         v3s16 blockpos_max = data->blockpos_max;
 504
 505         // Area of central chunk
 506         node_min = blockpos_min * MAP_BLOCKSIZE;
 507         node_max = (blockpos_max + v3s16(1, 1, 1)) * MAP_BLOCKSIZE - v3s16(1, 1, 1);
 508
 509         // Full allocated area
 510         full_node_min = (blockpos_min - 1) * MAP_BLOCKSIZE;
 511         full_node_max = (blockpos_max + 2) * MAP_BLOCKSIZE - v3s16(1, 1, 1);
 512
 513         central_area_size = node_max - node_min + v3s16(1, 1, 1);
 514         assert(central_area_size.X == central_area_size.Z);
 515
 516         int volume_blocks = (blockpos_max.X - blockpos_min.X + 1)
 517                                           * (blockpos_max.Y - blockpos_min.Y + 1)
 518                                           * (blockpos_max.Z - blockpos_max.Z + 1);
 519
 520         volume_nodes = volume_blocks *
 521                 MAP_BLOCKSIZE * MAP_BLOCKSIZE * MAP_BLOCKSIZE;
 522
 523         // Create a block-specific seed
 524         blockseed = get_blockseed(data->seed, full_node_min);
 525
 526         // Make some noise
 527         calculateNoise();
 528
 529         // Maximum height of the stone surface and obstacles.
 530         // This is used to guide the cave generation
 531         s16 stone_surface_max_y;
 532
 533         // Generate general ground level to full area
 534         stone_surface_max_y = generateGround();
 535
 536         // Create initial heightmap to limit caves
 537         updateHeightmap(node_min, node_max);
 538
 539         const s16 max_spread_amount = MAP_BLOCKSIZE;
 540         // Limit dirt flow area by 1 because mud is flown into neighbors.
 541         s16 mudflow_minpos = -max_spread_amount + 1;
 542         s16 mudflow_maxpos = central_area_size.X + max_spread_amount - 2;
 543
 544         // Loop this part, it will make stuff look older and newer nicely
 545         const u32 age_loops = 2;
 546         for (u32 i_age = 0; i_age < age_loops; i_age++) { // Aging loop
 547                 // Make caves (this code is relatively horrible)
 548                 if (flags & MG_CAVES)
 549                         generateCaves(stone_surface_max_y);
 550
 551                 // Add mud to the central chunk
 552                 addMud();
 553
 554                 // Flow mud away from steep edges
 555                 if (spflags & MGV6_MUDFLOW)
 556                         flowMud(mudflow_minpos, mudflow_maxpos);
 557
 558         }
 559
 560         // Update heightmap after mudflow
 561         updateHeightmap(node_min, node_max);
 562
 563         // Add dungeons
 564         if ((flags & MG_DUNGEONS) && (stone_surface_max_y >= node_min.Y)) {
 565                 DungeonParams dp;
 566
 567                 dp.seed = seed;
 568                 dp.c_water       = c_water_source;
 569                 dp.c_river_water = c_water_source;
 570                 dp.rooms_min     = 2;
 571                 dp.rooms_max     = 16;
 572                 dp.y_min         = -MAX_MAP_GENERATION_LIMIT;
 573                 dp.y_max         = MAX_MAP_GENERATION_LIMIT;
 574                 dp.np_density    = NoiseParams(0.9, 0.5, v3f(500.0, 500.0, 500.0), 0, 2, 0.8, 2.0);
 575                 dp.np_alt_wall   = NoiseParams(-0.4, 1.0, v3f(40.0, 40.0, 40.0), 32474, 6, 1.1, 2.0);
 576
 577                 if (getBiome(0, v2s16(node_min.X, node_min.Z)) == BT_DESERT) {
 578                         dp.c_wall     = c_desert_stone;
 579                         dp.c_alt_wall = CONTENT_IGNORE;
 580                         dp.c_stair    = c_stair_desert_stone;
 581
 582                         dp.diagonal_dirs = true;
 583                         dp.holesize      = v3s16(2, 3, 2);
 584                         dp.roomsize      = v3s16(2, 5, 2);
 585                         dp.notifytype    = GENNOTIFY_TEMPLE;
 586                 } else {
 587                         dp.c_wall     = c_cobble;
 588                         dp.c_alt_wall = c_mossycobble;
 589                         dp.c_stair    = c_stair_cobble;
 590
 591                         dp.diagonal_dirs = false;
 592                         dp.holesize      = v3s16(1, 2, 1);
 593                         dp.roomsize      = v3s16(0, 0, 0);
 594                         dp.notifytype    = GENNOTIFY_DUNGEON;
 595                 }
 596
 597                 DungeonGen dgen(ndef, &gennotify, &dp);
 598                 dgen.generate(vm, blockseed, full_node_min, full_node_max);
 599         }
 600
 601         // Add top and bottom side of water to transforming_liquid queue
 602         updateLiquid(&data->transforming_liquid, full_node_min, full_node_max);
 603
 604         // Add surface nodes
 605         growGrass();
 606
 607         // Generate some trees, and add grass, if a jungle
 608         if (spflags & MGV6_TREES)
 609                 placeTreesAndJungleGrass();
 610
 611         // Generate the registered decorations
 612         if (flags & MG_DECORATIONS)
 613                 m_emerge->decomgr->placeAllDecos(this, blockseed, node_min, node_max);
 614
 615         // Generate the registered ores
 616         m_emerge->oremgr->placeAllOres(this, blockseed, node_min, node_max);
 617
 618         // Calculate lighting
 619         if (flags & MG_LIGHT)
 620                 calcLighting(node_min - v3s16(1, 1, 1) * MAP_BLOCKSIZE,
 621                         node_max + v3s16(1, 0, 1) * MAP_BLOCKSIZE,
 622                         full_node_min, full_node_max);
 623
 624         this->generating = false;
 625 }
 626
 627
 628 void MapgenV6::calculateNoise()
 629 {
 630         int x = node_min.X;
 631         int z = node_min.Z;
 632         int fx = full_node_min.X;
 633         int fz = full_node_min.Z;
 634
 635         if (!(spflags & MGV6_FLAT)) {
 636                 noise_terrain_base->perlinMap2D_PO(x, 0.5, z, 0.5);
 637                 noise_terrain_higher->perlinMap2D_PO(x, 0.5, z, 0.5);
 638                 noise_steepness->perlinMap2D_PO(x, 0.5, z, 0.5);
 639                 noise_height_select->perlinMap2D_PO(x, 0.5, z, 0.5);
 640                 noise_mud->perlinMap2D_PO(x, 0.5, z, 0.5);
 641         }
 642
 643         noise_beach->perlinMap2D_PO(x, 0.2, z, 0.7);
 644
 645         noise_biome->perlinMap2D_PO(fx, 0.6, fz, 0.2);
 646         noise_humidity->perlinMap2D_PO(fx, 0.0, fz, 0.0);
 647         // Humidity map does not need range limiting 0 to 1,
 648         // only humidity at point does
 649 }
 650
 651
 652 int MapgenV6::generateGround()
 653 {
 654         //TimeTaker timer1("Generating ground level");
 655         MapNode n_air(CONTENT_AIR), n_water_source(c_water_source);
 656         MapNode n_stone(c_stone), n_desert_stone(c_desert_stone);
 657         MapNode n_ice(c_ice);
 658         int stone_surface_max_y = -MAX_MAP_GENERATION_LIMIT;
 659
 660         u32 index = 0;
 661         for (s16 z = node_min.Z; z <= node_max.Z; z++)
 662         for (s16 x = node_min.X; x <= node_max.X; x++, index++) {
 663                 // Surface height
 664                 s16 surface_y = (s16)baseTerrainLevelFromMap(index);
 665
 666                 // Log it
 667                 if (surface_y > stone_surface_max_y)
 668                         stone_surface_max_y = surface_y;
 669
 670                 BiomeV6Type bt = getBiome(v2s16(x, z));
 671
 672                 // Fill ground with stone
 673                 v3s16 em = vm->m_area.getExtent();
 674                 u32 i = vm->m_area.index(x, node_min.Y, z);
 675                 for (s16 y = node_min.Y; y <= node_max.Y; y++) {
 676                         if (vm->m_data[i].getContent() == CONTENT_IGNORE) {
 677                                 if (y <= surface_y) {
 678                                         vm->m_data[i] = (y >= MGV6_DESERT_STONE_BASE
 679                                                         && bt == BT_DESERT) ?
 680                                                 n_desert_stone : n_stone;
 681                                 } else if (y <= water_level) {
 682                                         vm->m_data[i] = (y >= MGV6_ICE_BASE
 683                                                         && bt == BT_TUNDRA) ?
 684                                                 n_ice : n_water_source;
 685                                 } else {
 686                                         vm->m_data[i] = n_air;
 687                                 }
 688                         }
 689                         vm->m_area.add_y(em, i, 1);
 690                 }
 691         }
 692
 693         return stone_surface_max_y;
 694 }
 695
 696
 697 void MapgenV6::addMud()
 698 {
 699         // 15ms @cs=8
 700         //TimeTaker timer1("add mud");
 701         MapNode n_dirt(c_dirt), n_gravel(c_gravel);
 702         MapNode n_sand(c_sand), n_desert_sand(c_desert_sand);
 703         MapNode addnode;
 704
 705         u32 index = 0;
 706         for (s16 z = node_min.Z; z <= node_max.Z; z++)
 707         for (s16 x = node_min.X; x <= node_max.X; x++, index++) {
 708                 // Randomize mud amount
 709                 s16 mud_add_amount = getMudAmount(index) / 2.0 + 0.5;
 710
 711                 // Find ground level
 712                 s16 surface_y = find_stone_level(v2s16(x, z)); /////////////////optimize this!
 713
 714                 // Handle area not found
 715                 if (surface_y == vm->m_area.MinEdge.Y - 1)
 716                         continue;
 717
 718                 BiomeV6Type bt = getBiome(v2s16(x, z));
 719                 addnode = (bt == BT_DESERT) ? n_desert_sand : n_dirt;
 720
 721                 if (bt == BT_DESERT && surface_y + mud_add_amount <= water_level + 1) {
 722                         addnode = n_sand;
 723                 } else if (mud_add_amount <= 0) {
 724                         mud_add_amount = 1 - mud_add_amount;
 725                         addnode = n_gravel;
 726                 } else if (bt != BT_DESERT && getHaveBeach(index) &&
 727                                 surface_y + mud_add_amount <= water_level + 2) {
 728                         addnode = n_sand;
 729                 }
 730
 731                 if ((bt == BT_DESERT || bt == BT_TUNDRA) && surface_y > 20)
 732                         mud_add_amount = MYMAX(0, mud_add_amount - (surface_y - 20) / 5);
 733
 734                 /* If topmost node is grass, change it to mud.  It might be if it was
 735                 // flown to there from a neighboring chunk and then converted.
 736                 u32 i = vm->m_area.index(x, surface_y, z);
 737                 if (vm->m_data[i].getContent() == c_dirt_with_grass)
 738                         vm->m_data[i] = n_dirt;*/
 739
 740                 // Add mud on ground
 741                 s16 mudcount = 0;
 742                 v3s16 em = vm->m_area.getExtent();
 743                 s16 y_start = surface_y + 1;
 744                 u32 i = vm->m_area.index(x, y_start, z);
 745                 for (s16 y = y_start; y <= node_max.Y; y++) {
 746                         if (mudcount >= mud_add_amount)
 747                                 break;
 748
 749                         vm->m_data[i] = addnode;
 750                         mudcount++;
 751
 752                         vm->m_area.add_y(em, i, 1);
 753                 }
 754         }
 755 }
 756
 757
 758 void MapgenV6::flowMud(s16 &mudflow_minpos, s16 &mudflow_maxpos)
 759 {
 760         // 340ms @cs=8
 761         //TimeTaker timer1("flow mud");
 762
 763         // Iterate a few times
 764         for (s16 k = 0; k < 3; k++) {
 765                 for (s16 z = mudflow_minpos; z <= mudflow_maxpos; z++)
 766                 for (s16 x = mudflow_minpos; x <= mudflow_maxpos; x++) {
 767                         // Invert coordinates every 2nd iteration
 768                         if (k % 2 == 0) {
 769                                 x = mudflow_maxpos - (x - mudflow_minpos);
 770                                 z = mudflow_maxpos - (z - mudflow_minpos);
 771                         }
 772
 773                         // Node position in 2d
 774                         v2s16 p2d = v2s16(node_min.X, node_min.Z) + v2s16(x, z);
 775
 776                         v3s16 em = vm->m_area.getExtent();
 777                         u32 i = vm->m_area.index(p2d.X, node_max.Y, p2d.Y);
 778                         s16 y = node_max.Y;
 779
 780                         while (y >= node_min.Y) {
 781
 782                         for (;; y--) {
 783                                 MapNode *n = NULL;
 784                                 // Find mud
 785                                 for (; y >= node_min.Y; y--) {
 786                                         n = &vm->m_data[i];
 787                                         if (n->getContent() == c_dirt ||
 788                                                         n->getContent() == c_dirt_with_grass ||
 789                                                         n->getContent() == c_gravel)
 790                                                 break;
 791
 792                                         vm->m_area.add_y(em, i, -1);
 793                                 }
 794
 795                                 // Stop if out of area
 796                                 //if(vmanip.m_area.contains(i) == false)
 797                                 if (y < node_min.Y)
 798                                         break;
 799
 800                                 if (n->getContent() == c_dirt ||
 801                                                 n->getContent() == c_dirt_with_grass) {
 802                                         // Make it exactly mud
 803                                         n->setContent(c_dirt);
 804
 805                                         // Don't flow it if the stuff under it is not mud
 806                                         {
 807                                                 u32 i2 = i;
 808                                                 vm->m_area.add_y(em, i2, -1);
 809                                                 // Cancel if out of area
 810                                                 if (vm->m_area.contains(i2) == false)
 811                                                         continue;
 812                                                 MapNode *n2 = &vm->m_data[i2];
 813                                                 if (n2->getContent() != c_dirt &&
 814                                                                 n2->getContent() != c_dirt_with_grass)
 815                                                         continue;
 816                                         }
 817                                 }
 818
 819                                 v3s16 dirs4[4] = {
 820                                         v3s16(0, 0, 1), // back
 821                                         v3s16(1, 0, 0), // right
 822                                         v3s16(0, 0, -1), // front
 823                                         v3s16(-1, 0, 0), // left
 824                                 };
 825
 826                                 // Check that upper is air or doesn't exist.
 827                                 // Cancel dropping if upper keeps it in place
 828                                 u32 i3 = i;
 829                                 vm->m_area.add_y(em, i3, 1);
 830                                 if (vm->m_area.contains(i3) == true &&
 831                                                 ndef->get(vm->m_data[i3]).walkable)
 832                                         continue;
 833
 834                                 // Drop mud on side
 835                                 for (u32 di = 0; di < 4; di++) {
 836                                         v3s16 dirp = dirs4[di];
 837                                         u32 i2 = i;
 838                                         // Move to side
 839                                         vm->m_area.add_p(em, i2, dirp);
 840                                         // Fail if out of area
 841                                         if (vm->m_area.contains(i2) == false)
 842                                                 continue;
 843                                         // Check that side is air
 844                                         MapNode *n2 = &vm->m_data[i2];
 845                                         if (ndef->get(*n2).walkable)
 846                                                 continue;
 847                                         // Check that under side is air
 848                                         vm->m_area.add_y(em, i2, -1);
 849                                         if (vm->m_area.contains(i2) == false)
 850                                                 continue;
 851                                         n2 = &vm->m_data[i2];
 852                                         if (ndef->get(*n2).walkable)
 853                                                 continue;
 854                                         // Loop further down until not air
 855                                         bool dropped_to_unknown = false;
 856                                         do {
 857                                                 vm->m_area.add_y(em, i2, -1);
 858                                                 n2 = &vm->m_data[i2];
 859                                                 // if out of known area
 860                                                 if (vm->m_area.contains(i2) == false ||
 861                                                                 n2->getContent() == CONTENT_IGNORE) {
 862                                                         dropped_to_unknown = true;
 863                                                         break;
 864                                                 }
 865                                         } while (ndef->get(*n2).walkable == false);
 866                                         // Loop one up so that we're in air
 867                                         vm->m_area.add_y(em, i2, 1);
 868                                         n2 = &vm->m_data[i2];
 869
 870                                         bool old_is_water = (n->getContent() == c_water_source);
 871                                         // Move mud to new place
 872                                         if (!dropped_to_unknown) {
 873                                                 *n2 = *n;
 874                                                 // Set old place to be air (or water)
 875                                                 if (old_is_water)
 876                                                         *n = MapNode(c_water_source);
 877                                                 else
 878                                                         *n = MapNode(CONTENT_AIR);
 879                                         }
 880
 881                                         // Done
 882                                         break;
 883                                 }
 884                         }
 885                         }
 886                 }
 887         }
 888 }
 889
 890
 891 void MapgenV6::placeTreesAndJungleGrass()
 892 {
 893         //TimeTaker t("placeTrees");
 894         if (node_max.Y < water_level)
 895                 return;
 896
 897         PseudoRandom grassrandom(blockseed + 53);
 898         content_t c_junglegrass = ndef->getId("mapgen_junglegrass");
 899         // if we don't have junglegrass, don't place cignore... that's bad
 900         if (c_junglegrass == CONTENT_IGNORE)
 901                 c_junglegrass = CONTENT_AIR;
 902         MapNode n_junglegrass(c_junglegrass);
 903         v3s16 em = vm->m_area.getExtent();
 904
 905         // Divide area into parts
 906         s16 div = 8;
 907         s16 sidelen = central_area_size.X / div;
 908         double area = sidelen * sidelen;
 909
 910         // N.B.  We must add jungle grass first, since tree leaves will
 911         // obstruct the ground, giving us a false ground level
 912         for (s16 z0 = 0; z0 < div; z0++)
 913         for (s16 x0 = 0; x0 < div; x0++) {
 914                 // Center position of part of division
 915                 v2s16 p2d_center(
 916                         node_min.X + sidelen / 2 + sidelen * x0,
 917                         node_min.Z + sidelen / 2 + sidelen * z0
 918                 );
 919                 // Minimum edge of part of division
 920                 v2s16 p2d_min(
 921                         node_min.X + sidelen * x0,
 922                         node_min.Z + sidelen * z0
 923                 );
 924                 // Maximum edge of part of division
 925                 v2s16 p2d_max(
 926                         node_min.X + sidelen + sidelen * x0 - 1,
 927                         node_min.Z + sidelen + sidelen * z0 - 1
 928                 );
 929
 930                 // Get biome at center position of part of division
 931                 BiomeV6Type bt = getBiome(p2d_center);
 932
 933                 // Amount of trees
 934                 u32 tree_count;
 935                 if (bt == BT_JUNGLE || bt == BT_TAIGA || bt == BT_NORMAL) {
 936                         tree_count = area * getTreeAmount(p2d_center);
 937                         if (bt == BT_JUNGLE)
 938                                 tree_count *= 4;
 939                 } else {
 940                         tree_count = 0;
 941                 }
 942
 943                 // Add jungle grass
 944                 if (bt == BT_JUNGLE) {
 945                         float humidity = getHumidity(p2d_center);
 946                         u32 grass_count = 5 * humidity * tree_count;
 947                         for (u32 i = 0; i < grass_count; i++) {
 948                                 s16 x = grassrandom.range(p2d_min.X, p2d_max.X);
 949                                 s16 z = grassrandom.range(p2d_min.Y, p2d_max.Y);
 950                                 int mapindex = central_area_size.X * (z - node_min.Z)
 951                                                                 + (x - node_min.X);
 952                                 s16 y = heightmap[mapindex];
 953                                 if (y < water_level)
 954                                         continue;
 955
 956                                 u32 vi = vm->m_area.index(x, y, z);
 957                                 // place on dirt_with_grass, since we know it is exposed to sunlight
 958                                 if (vm->m_data[vi].getContent() == c_dirt_with_grass) {
 959                                         vm->m_area.add_y(em, vi, 1);
 960                                         vm->m_data[vi] = n_junglegrass;
 961                                 }
 962                         }
 963                 }
 964
 965                 // Put trees in random places on part of division
 966                 for (u32 i = 0; i < tree_count; i++) {
 967                         s16 x = myrand_range(p2d_min.X, p2d_max.X);
 968                         s16 z = myrand_range(p2d_min.Y, p2d_max.Y);
 969                         int mapindex = central_area_size.X * (z - node_min.Z)
 970                                                         + (x - node_min.X);
 971                         s16 y = heightmap[mapindex];
 972                         // Don't make a tree under water level
 973                         // Don't make a tree so high that it doesn't fit
 974                         if (y < water_level || y > node_max.Y - 6)
 975                                 continue;
 976
 977                         v3s16 p(x, y, z);
 978                         // Trees grow only on mud and grass and snowblock
 979                         {
 980                                 u32 i = vm->m_area.index(p);
 981                                 content_t c = vm->m_data[i].getContent();
 982                                 if (c != c_dirt &&
 983                                                 c != c_dirt_with_grass &&
 984                                                 c != c_dirt_with_snow &&
 985                                                 c != c_snowblock)
 986                                         continue;
 987                         }
 988                         p.Y++;
 989
 990                         // Make a tree
 991                         if (bt == BT_JUNGLE) {
 992                                 treegen::make_jungletree(*vm, p, ndef, myrand());
 993                         } else if (bt == BT_TAIGA) {
 994                                 treegen::make_pine_tree(*vm, p - v3s16(0, 1, 0), ndef, myrand());
 995                         } else if (bt == BT_NORMAL) {
 996                                 bool is_apple_tree = (myrand_range(0, 3) == 0) &&
 997                                                         getHaveAppleTree(v2s16(x, z));
 998                                 treegen::make_tree(*vm, p, is_apple_tree, ndef, myrand());
 999                         }
1000                 }
1001         }
1002         //printf("placeTreesAndJungleGrass: %dms\n", t.stop());
1003 }
1004
1005
1006 void MapgenV6::growGrass() // Add surface nodes
1007 {
1008         MapNode n_dirt_with_grass(c_dirt_with_grass);
1009         MapNode n_dirt_with_snow(c_dirt_with_snow);
1010         MapNode n_snowblock(c_snowblock);
1011         MapNode n_snow(c_snow);
1012         v3s16 em = vm->m_area.getExtent();
1013
1014         u32 index = 0;
1015         for (s16 z = full_node_min.Z; z <= full_node_max.Z; z++)
1016         for (s16 x = full_node_min.X; x <= full_node_max.X; x++, index++) {
1017                 // Find the lowest surface to which enough light ends up to make
1018                 // grass grow.  Basically just wait until not air and not leaves.
1019                 s16 surface_y = 0;
1020                 {
1021                         u32 i = vm->m_area.index(x, node_max.Y, z);
1022                         s16 y;
1023                         // Go to ground level
1024                         for (y = node_max.Y; y >= full_node_min.Y; y--) {
1025                                 MapNode &n = vm->m_data[i];
1026                                 if (ndef->get(n).param_type != CPT_LIGHT ||
1027                                                 ndef->get(n).liquid_type != LIQUID_NONE ||
1028                                                 n.getContent() == c_ice)
1029                                         break;
1030                                 vm->m_area.add_y(em, i, -1);
1031                         }
1032                         surface_y = (y >= full_node_min.Y) ? y : full_node_min.Y;
1033                 }
1034
1035                 BiomeV6Type bt = getBiome(index, v2s16(x, z));
1036                 u32 i = vm->m_area.index(x, surface_y, z);
1037                 content_t c = vm->m_data[i].getContent();
1038                 if (surface_y >= water_level - 20) {
1039                         if (bt == BT_TAIGA && c == c_dirt) {
1040                                 vm->m_data[i] = n_snowblock;
1041                                 vm->m_area.add_y(em, i, -1);
1042                                 vm->m_data[i] = n_dirt_with_snow;
1043                         } else if (bt == BT_TUNDRA) {
1044                                 if (c == c_dirt) {
1045                                         vm->m_data[i] = n_dirt_with_snow;
1046                                 } else if (c == c_stone && surface_y < node_max.Y) {
1047                                         vm->m_area.add_y(em, i, 1);
1048                                         vm->m_data[i] = n_snow;
1049                                 }
1050                         } else if (c == c_dirt) {
1051                                 vm->m_data[i] = n_dirt_with_grass;
1052                         }
1053                 }
1054         }
1055 }
1056
1057
1058 void MapgenV6::generateCaves(int max_stone_y)
1059 {
1060         float cave_amount = NoisePerlin2D(np_cave, node_min.X, node_min.Y, seed);
1061         int volume_nodes = (node_max.X - node_min.X + 1) *
1062                                            (node_max.Y - node_min.Y + 1) * MAP_BLOCKSIZE;
1063         cave_amount = MYMAX(0.0, cave_amount);
1064         u32 caves_count = cave_amount * volume_nodes / 50000;
1065         u32 bruises_count = 1;
1066         PseudoRandom ps(blockseed + 21343);
1067         PseudoRandom ps2(blockseed + 1032);
1068
1069         if (ps.range(1, 6) == 1)
1070                 bruises_count = ps.range(0, ps.range(0, 2));
1071
1072         if (getBiome(v2s16(node_min.X, node_min.Z)) == BT_DESERT) {
1073                 caves_count   /= 3;
1074                 bruises_count /= 3;
1075         }
1076
1077         for (u32 i = 0; i < caves_count + bruises_count; i++) {
1078                 CavesV6 cave(ndef, &gennotify, water_level, c_water_source, c_lava_source);
1079
1080                 bool large_cave = (i >= caves_count);
1081                 cave.makeCave(vm, node_min, node_max, &ps, &ps2,
1082                         large_cave, max_stone_y, heightmap);
1083         }
1084 }