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