Heightmaps: Fix uninitialised values in mgv5/mgv6. findGroundLevel: Return -MAP_GENER...
[oweals/minetest.git] / src / mapgen_v5.cpp
1 /*
2 Minetest
3 Copyright (C) 2010-2013 kwolekr, Ryan Kwolek <kwolekr@minetest.net>
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 "content_sao.h"
28 #include "nodedef.h"
29 #include "voxelalgorithms.h"
30 #include "profiler.h"
31 #include "settings.h" // For g_settings
32 #include "main.h" // For g_profiler
33 #include "emerge.h"
34 #include "dungeongen.h"
35 #include "cavegen.h"
36 #include "treegen.h"
37 #include "mg_biome.h"
38 #include "mg_ore.h"
39 #include "mg_decoration.h"
40 #include "mapgen_v5.h"
41 #include "util/directiontables.h"
42
43
44 FlagDesc flagdesc_mapgen_v5[] = {
45         {NULL,         0}
46 };
47
48
49 MapgenV5::MapgenV5(int mapgenid, MapgenParams *params, EmergeManager *emerge)
50         : Mapgen(mapgenid, params, emerge)
51 {
52         this->m_emerge = emerge;
53         this->bmgr     = emerge->biomemgr;
54
55         // amount of elements to skip for the next index
56         // for noise/height/biome maps (not vmanip)
57         this->ystride = csize.X;
58         this->zstride = csize.X * (csize.Y + 2);
59
60         this->biomemap  = new u8[csize.X * csize.Z];
61         this->heightmap = new s16[csize.X * csize.Z];
62
63         MapgenV5Params *sp = (MapgenV5Params *)params->sparams;
64         this->spflags      = sp->spflags;
65
66         // Terrain noise
67         noise_filler_depth = new Noise(&sp->np_filler_depth, seed, csize.X, csize.Z);
68         noise_factor       = new Noise(&sp->np_factor,       seed, csize.X, csize.Z);
69         noise_height       = new Noise(&sp->np_height,       seed, csize.X, csize.Z);
70
71         // 3D terrain noise
72         noise_cave1        = new Noise(&sp->np_cave1,   seed, csize.X, csize.Y + 2, csize.Z);
73         noise_cave2        = new Noise(&sp->np_cave2,   seed, csize.X, csize.Y + 2, csize.Z);
74         noise_ground       = new Noise(&sp->np_ground,  seed, csize.X, csize.Y + 2, csize.Z);
75
76         // Biome noise
77         noise_heat         = new Noise(&params->np_biome_heat,     seed, csize.X, csize.Z);
78         noise_humidity     = new Noise(&params->np_biome_humidity, seed, csize.X, csize.Z);
79
80         //// Resolve nodes to be used
81         INodeDefManager *ndef = emerge->ndef;
82
83         c_stone           = ndef->getId("mapgen_stone");
84         c_dirt            = ndef->getId("mapgen_dirt");
85         c_dirt_with_grass = ndef->getId("mapgen_dirt_with_grass");
86         c_sand            = ndef->getId("mapgen_sand");
87         c_water_source    = ndef->getId("mapgen_water_source");
88         c_lava_source     = ndef->getId("mapgen_lava_source");
89         c_gravel          = ndef->getId("mapgen_gravel");
90         c_cobble          = ndef->getId("mapgen_cobble");
91         c_ice             = ndef->getId("default:ice");
92         c_mossycobble     = ndef->getId("mapgen_mossycobble");
93         c_sandbrick       = ndef->getId("mapgen_sandstonebrick");
94         c_stair_cobble    = ndef->getId("mapgen_stair_cobble");
95         c_stair_sandstone = ndef->getId("mapgen_stair_sandstone");
96         if (c_ice == CONTENT_IGNORE)
97                 c_ice = CONTENT_AIR;
98         if (c_mossycobble == CONTENT_IGNORE)
99                 c_mossycobble = c_cobble;
100         if (c_sandbrick == CONTENT_IGNORE)
101                 c_sandbrick = c_desert_stone;
102         if (c_stair_cobble == CONTENT_IGNORE)
103                 c_stair_cobble = c_cobble;
104         if (c_stair_sandstone == CONTENT_IGNORE)
105                 c_stair_sandstone = c_sandbrick;
106 }
107
108
109 MapgenV5::~MapgenV5()
110 {
111         delete noise_filler_depth;
112         delete noise_factor;
113         delete noise_height;
114         delete noise_cave1;
115         delete noise_cave2;
116         delete noise_ground;
117
118         delete noise_heat;
119         delete noise_humidity;
120
121         delete[] heightmap;
122         delete[] biomemap;
123 }
124
125
126 MapgenV5Params::MapgenV5Params()
127 {
128         spflags = 0;
129
130         np_filler_depth = NoiseParams(0, 1,  v3f(150, 150, 150), 261,    4, 0.7,  2.0);
131         np_factor       = NoiseParams(0, 1,  v3f(250, 250, 250), 920381, 3, 0.45, 2.0);
132         np_height       = NoiseParams(0, 10, v3f(250, 250, 250), 84174,  4, 0.5,  2.0);
133         np_cave1        = NoiseParams(0, 12, v3f(50,  50,  50),  52534,  4, 0.5,  2.0);
134         np_cave2        = NoiseParams(0, 12, v3f(50,  50,  50),  10325,  4, 0.5,  2.0);
135         np_ground       = NoiseParams(0, 40, v3f(80,  80,  80),  983240, 4, 0.55, 2.0, NOISE_FLAG_EASED);
136 }
137
138
139 //#define CAVE_NOISE_SCALE 12.0
140 //#define CAVE_NOISE_THRESHOLD (1.5/CAVE_NOISE_SCALE) = 0.125
141
142
143 void MapgenV5Params::readParams(const Settings *settings)
144 {
145         settings->getFlagStrNoEx("mgv5_spflags", spflags, flagdesc_mapgen_v5);
146
147         settings->getNoiseParams("mgv5_np_filler_depth", np_filler_depth);
148         settings->getNoiseParams("mgv5_np_factor",       np_factor);
149         settings->getNoiseParams("mgv5_np_height",       np_height);
150         settings->getNoiseParams("mgv5_np_cave1",        np_cave1);
151         settings->getNoiseParams("mgv5_np_cave2",        np_cave2);
152         settings->getNoiseParams("mgv5_np_ground",       np_ground);
153 }
154
155
156 void MapgenV5Params::writeParams(Settings *settings) const
157 {
158         settings->setFlagStr("mgv5_spflags", spflags, flagdesc_mapgen_v5, (u32)-1);
159
160         settings->setNoiseParams("mgv5_np_filler_depth", np_filler_depth);
161         settings->setNoiseParams("mgv5_np_factor",       np_factor);
162         settings->setNoiseParams("mgv5_np_height",       np_height);
163         settings->setNoiseParams("mgv5_np_cave1",        np_cave1);
164         settings->setNoiseParams("mgv5_np_cave2",        np_cave2);
165         settings->setNoiseParams("mgv5_np_ground",       np_ground);
166 }
167
168
169 int MapgenV5::getGroundLevelAtPoint(v2s16 p)
170 {
171         //TimeTaker t("getGroundLevelAtPoint", NULL, PRECISION_MICRO);
172
173         float f = 0.55 + NoisePerlin2D(&noise_factor->np, p.X, p.Y, seed);
174         if (f < 0.01)
175                 f = 0.01;
176         else if (f >= 1.0)
177                 f *= 1.6;
178         float h = water_level + NoisePerlin2D(&noise_height->np, p.X, p.Y, seed);
179
180         s16 search_top = water_level + 15;
181         s16 search_base = water_level;
182
183         s16 level = -31000;
184         for (s16 y = search_top; y >= search_base; y--) {
185                 float n_ground = NoisePerlin3D(&noise_ground->np, p.X, y, p.Y, seed);
186                 if (n_ground * f > y - h) {
187                         if (y >= search_top - 7)
188                                 break;
189                         else
190                                 level = y;
191                                 break;
192                 }
193         }
194
195         //printf("getGroundLevelAtPoint: %dus\n", t.stop());
196         return level;
197 }
198
199
200 void MapgenV5::makeChunk(BlockMakeData *data)
201 {
202         // Pre-conditions
203         assert(data->vmanip);
204         assert(data->nodedef);
205         assert(data->blockpos_requested.X >= data->blockpos_min.X &&
206                 data->blockpos_requested.Y >= data->blockpos_min.Y &&
207                 data->blockpos_requested.Z >= data->blockpos_min.Z);
208         assert(data->blockpos_requested.X <= data->blockpos_max.X &&
209                 data->blockpos_requested.Y <= data->blockpos_max.Y &&
210                 data->blockpos_requested.Z <= data->blockpos_max.Z);
211
212         generating = true;
213         vm   = data->vmanip;
214         ndef = data->nodedef;
215         //TimeTaker t("makeChunk");
216
217         v3s16 blockpos_min = data->blockpos_min;
218         v3s16 blockpos_max = data->blockpos_max;
219         node_min = blockpos_min * MAP_BLOCKSIZE;
220         node_max = (blockpos_max + v3s16(1, 1, 1)) * MAP_BLOCKSIZE - v3s16(1, 1, 1);
221         full_node_min = (blockpos_min - 1) * MAP_BLOCKSIZE;
222         full_node_max = (blockpos_max + 2) * MAP_BLOCKSIZE - v3s16(1, 1, 1);
223
224         // Create a block-specific seed
225         blockseed = getBlockSeed2(full_node_min, seed);
226
227         // Make some noise
228         calculateNoise();
229
230         // Generate base terrain
231         s16 stone_surface_max_y = generateBaseTerrain();
232
233         // Create heightmap
234         updateHeightmap(node_min, node_max);
235
236         // Create biomemap at heightmap surface
237         bmgr->calcBiomes(csize.X, csize.Z, noise_heat->result,
238                 noise_humidity->result, heightmap, biomemap);
239
240         // Actually place the biome-specific nodes
241         generateBiomes(noise_heat->result, noise_humidity->result);
242
243         // Generate caves
244         if ((flags & MG_CAVES) && (stone_surface_max_y >= node_min.Y))
245                 generateCaves(stone_surface_max_y);
246
247         // Generate dungeons and desert temples
248         if ((flags & MG_DUNGEONS) && (stone_surface_max_y >= node_min.Y)) {
249                 DungeonGen dgen(this, NULL);
250                 dgen.generate(blockseed, full_node_min, full_node_max);
251         }
252
253         // Generate the registered decorations
254         m_emerge->decomgr->placeAllDecos(this, blockseed, node_min, node_max);
255
256         // Generate the registered ores
257         m_emerge->oremgr->placeAllOres(this, blockseed, node_min, node_max);
258
259         // Sprinkle some dust on top after everything else was generated
260         dustTopNodes();
261
262         //printf("makeChunk: %dms\n", t.stop());
263
264         // Add top and bottom side of water to transforming_liquid queue
265         updateLiquid(&data->transforming_liquid, full_node_min, full_node_max);
266
267         // Calculate lighting
268         if (flags & MG_LIGHT) {
269                 calcLighting(node_min - v3s16(0, 1, 0), node_max + v3s16(0, 1, 0),
270                         full_node_min, full_node_max);
271         }
272
273         this->generating = false;
274 }
275
276
277 void MapgenV5::calculateNoise()
278 {
279         //TimeTaker t("calculateNoise", NULL, PRECISION_MICRO);
280         int x = node_min.X;
281         int y = node_min.Y - 1;
282         int z = node_min.Z;
283
284         noise_factor->perlinMap2D(x, z);
285         noise_height->perlinMap2D(x, z);
286         noise_ground->perlinMap3D(x, y, z);
287
288         if (flags & MG_CAVES) {
289                 noise_cave1->perlinMap3D(x, y, z);
290                 noise_cave2->perlinMap3D(x, y, z);
291         }
292
293         if (node_max.Y >= water_level) {
294                 noise_filler_depth->perlinMap2D(x, z);
295                 noise_heat->perlinMap2D(x, z);
296                 noise_humidity->perlinMap2D(x, z);
297         }
298
299         //printf("calculateNoise: %dus\n", t.stop());
300 }
301
302
303 //bool is_cave(u32 index) {
304 //      double d1 = contour(noise_cave1->result[index]);
305 //      double d2 = contour(noise_cave2->result[index]);
306 //      return d1*d2 > CAVE_NOISE_THRESHOLD;
307 //}
308
309 //bool val_is_ground(v3s16 p, u32 index, u32 index2d) {
310 //      double f = 0.55 + noise_factor->result[index2d];
311 //      if(f < 0.01)
312 //              f = 0.01;
313 //      else if(f >= 1.0)
314 //              f *= 1.6;
315 //      double h = WATER_LEVEL + 10 * noise_height->result[index2d];
316 //      return (noise_ground->result[index] * f > (double)p.Y - h);
317 //}
318
319
320 int MapgenV5::generateBaseTerrain()
321 {
322         u32 index = 0;
323         u32 index2d = 0;
324         int stone_surface_max_y = -MAP_GENERATION_LIMIT;
325
326         for (s16 z=node_min.Z; z<=node_max.Z; z++) {
327                 for (s16 y=node_min.Y - 1; y<=node_max.Y + 1; y++) {
328                         u32 i = vm->m_area.index(node_min.X, y, z);
329                         for (s16 x=node_min.X; x<=node_max.X; x++, i++, index++, index2d++) {
330                                 if (vm->m_data[i].getContent() != CONTENT_IGNORE)
331                                         continue;
332
333                                 float f = 0.55 + noise_factor->result[index2d];
334                                 if (f < 0.01)
335                                         f = 0.01;
336                                 else if (f >= 1.0)
337                                         f *= 1.6;
338                                 float h = noise_height->result[index2d];
339
340                                 if (noise_ground->result[index] * f < y - h) {
341                                         if (y <= water_level)
342                                                 vm->m_data[i] = MapNode(c_water_source);
343                                         else
344                                                 vm->m_data[i] = MapNode(CONTENT_AIR);
345                                 } else {
346                                         vm->m_data[i] = MapNode(c_stone);
347                                         if (y > stone_surface_max_y)
348                                                 stone_surface_max_y = y;
349                                 }
350                         }
351                         index2d = index2d - ystride;
352                 }
353                 index2d = index2d + ystride;
354         }
355
356         return stone_surface_max_y;
357 }
358
359
360 void MapgenV5::generateBiomes(float *heat_map, float *humidity_map)
361 {
362         if (node_max.Y < water_level)
363                 return;
364
365         MapNode n_air(CONTENT_AIR);
366         MapNode n_stone(c_stone);
367         MapNode n_water(c_water_source);
368
369         v3s16 em = vm->m_area.getExtent();
370         u32 index = 0;
371
372         for (s16 z = node_min.Z; z <= node_max.Z; z++)
373         for (s16 x = node_min.X; x <= node_max.X; x++, index++) {
374                 Biome *biome = NULL;
375                 s16 dfiller = 0;
376                 s16 y0_top = 0;
377                 s16 y0_filler = 0;
378                 s16 depth_water_top = 0;
379
380                 s16 nplaced = 0;
381                 u32 i = vm->m_area.index(x, node_max.Y, z);
382
383                 content_t c_above = vm->m_data[i + em.X].getContent();
384                 bool have_air = c_above == CONTENT_AIR;
385
386                 for (s16 y = node_max.Y; y >= node_min.Y; y--) {
387                         content_t c = vm->m_data[i].getContent();
388
389                         if (c != CONTENT_IGNORE && c != CONTENT_AIR && (y == node_max.Y || have_air)) {
390                                 biome           = bmgr->getBiome(heat_map[index], humidity_map[index], y);
391                                 dfiller         = biome->depth_filler + noise_filler_depth->result[index];
392                                 y0_top          = biome->depth_top;
393                                 y0_filler       = biome->depth_top + dfiller;
394                                 depth_water_top = biome->depth_water_top;
395                         }
396
397                         if (c == c_stone && have_air) {
398                                 content_t c_below = vm->m_data[i - em.X].getContent();
399
400                                 if (c_below != CONTENT_AIR) {
401                                         if (nplaced < y0_top) {
402                                                 vm->m_data[i] = MapNode(biome->c_top);
403                                                 nplaced++;
404                                         } else if (nplaced < y0_filler && nplaced >= y0_top) {
405                                                 vm->m_data[i] = MapNode(biome->c_filler);
406                                                 nplaced++;
407                                         } else if (c == c_stone) {
408                                                 have_air = false;
409                                                 nplaced  = 0;
410                                                 vm->m_data[i] = MapNode(biome->c_stone);
411                                         } else {
412                                                 have_air = false;
413                                                 nplaced  = 0;
414                                         }
415                                 } else if (c == c_stone) {
416                                         have_air = false;
417                                         nplaced = 0;
418                                         vm->m_data[i] = MapNode(biome->c_stone);
419                                 }
420                         } else if (c == c_stone) {
421                                 have_air = false;
422                                 nplaced = 0;
423                                 vm->m_data[i] = MapNode(biome->c_stone);
424                         } else if (c == c_water_source) {
425                                 have_air = true;
426                                 nplaced = 0;
427                                 if(y > water_level - depth_water_top)
428                                         vm->m_data[i] = MapNode(biome->c_water_top);
429                                 else
430                                         vm->m_data[i] = MapNode(biome->c_water);
431                         } else if (c == CONTENT_AIR) {
432                                 have_air = true;
433                                 nplaced = 0;
434                         }
435
436                         vm->m_area.add_y(em, i, -1);
437                 }
438         }
439 }
440
441
442 void MapgenV5::generateCaves(int max_stone_y)
443 {
444         u32 index = 0;
445         u32 index2d = 0;
446
447         for (s16 z=node_min.Z; z<=node_max.Z; z++) {
448                 for (s16 y=node_min.Y - 1; y<=node_max.Y + 1; y++) {
449                         u32 i = vm->m_area.index(node_min.X, y, z);
450                         for (s16 x=node_min.X; x<=node_max.X; x++, i++, index++, index2d++) {
451                                 Biome *biome = (Biome *)bmgr->get(biomemap[index2d]);
452                                 content_t c = vm->m_data[i].getContent();
453                                 if (c == CONTENT_AIR
454                                                 || (y <= water_level
455                                                 && c != biome->c_stone
456                                                 && c != c_stone))
457                                         continue;
458
459                                 float d1 = contour(noise_cave1->result[index]);
460                                 float d2 = contour(noise_cave2->result[index]);
461                                 if (d1*d2 > 0.125)
462                                         vm->m_data[i] = MapNode(CONTENT_AIR);
463                         }
464                         index2d = index2d - ystride;
465                 }
466                 index2d = index2d + ystride;
467         }
468
469         if (node_max.Y > -256)
470                 return;
471
472         PseudoRandom ps(blockseed + 21343);
473         u32 bruises_count = (ps.range(1, 5) == 1) ? ps.range(1, 2) : 0;
474         for (u32 i = 0; i < bruises_count; i++) {
475                 CaveV5 cave(this, &ps);
476                 cave.makeCave(node_min, node_max, max_stone_y);
477         }
478 }
479
480
481 void MapgenV5::dustTopNodes()
482 {
483         if (node_max.Y < water_level)
484                 return;
485
486         v3s16 em = vm->m_area.getExtent();
487         u32 index = 0;
488
489         for (s16 z = node_min.Z; z <= node_max.Z; z++)
490         for (s16 x = node_min.X; x <= node_max.X; x++, index++) {
491                 Biome *biome = (Biome *)bmgr->get(biomemap[index]);
492
493                 if (biome->c_dust == CONTENT_IGNORE)
494                         continue;
495
496                 s16 y = node_max.Y + 1;
497                 u32 vi = vm->m_area.index(x, y, z);
498                 for (; y >= node_min.Y; y--) {
499                         if (vm->m_data[vi].getContent() != CONTENT_AIR)
500                                 break;
501
502                         vm->m_area.add_y(em, vi, -1);
503                 }
504
505                 content_t c = vm->m_data[vi].getContent();
506                 if (!ndef->get(c).buildable_to && c != CONTENT_IGNORE
507                                 && c != biome->c_dust) {
508                         if (y == node_max.Y + 1)
509                                 continue;
510
511                         vm->m_area.add_y(em, vi, 1);
512                         vm->m_data[vi] = MapNode(biome->c_dust);
513                 }
514         }
515 }
516