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