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