AR Design
UBC EML collab with UBC SALA - visualizing IoT data in AR
FastSimplexNoise.cs
Go to the documentation of this file.
1 // Copyright (c) Microsoft Corporation. All rights reserved.
2 // Licensed under the MIT License. See LICENSE in the project root for license information.
3 
4 using System;
5 
9 public class FastSimplexNoise
10 {
11  private const double STRETCH_2D = -1.0 / 4.73205;
12  private const double STRETCH_3D = -1.0 / 6.0;
13  private const double STRETCH_4D = -1.0 / 7.23607;
14  private const double SQUISH_2D = 1.0 / 2.73205;
15  private const double SQUISH_3D = 1.0 / 3.0;
16  private const double SQUISH_4D = 1.0 / 7.23607;
17  private const double NORM_2D = 1.0 / 47.0;
18  private const double NORM_3D = 1.0 / 103.0;
19  private const double NORM_4D = 1.0 / 30.0;
20 
21  private const Int64 SEEDVAL_1 = 6364136223846793005L;
22  private const Int64 SEEDVAL_2 = 1442695040888963407L;
23 
24  private byte[] perm;
25  private byte[] perm2D;
26  private byte[] perm3D;
27  private byte[] perm4D;
28 
29  private static double[] gradients2D = new double[]
30  {
31  5, 2, 2, 5,
32  -5, 2, -2, 5,
33  5, -2, 2, -5,
34  -5, -2, -2, -5,
35  };
36 
37  private static double[] gradients3D =
38  {
39  -11, 4, 4, -4, 11, 4, -4, 4, 11,
40  11, 4, 4, 4, 11, 4, 4, 4, 11,
41  -11, -4, 4, -4, -11, 4, -4, -4, 11,
42  11, -4, 4, 4, -11, 4, 4, -4, 11,
43  -11, 4, -4, -4, 11, -4, -4, 4, -11,
44  11, 4, -4, 4, 11, -4, 4, 4, -11,
45  -11, -4, -4, -4, -11, -4, -4, -4, -11,
46  11, -4, -4, 4, -11, -4, 4, -4, -11,
47  };
48 
49  private static double[] gradients4D =
50  {
51  3, 1, 1, 1, 1, 3, 1, 1, 1, 1, 3, 1, 1, 1, 1, 3,
52  -3, 1, 1, 1, -1, 3, 1, 1, -1, 1, 3, 1, -1, 1, 1, 3,
53  3, -1, 1, 1, 1, -3, 1, 1, 1, -1, 3, 1, 1, -1, 1, 3,
54  -3, -1, 1, 1, -1, -3, 1, 1, -1, -1, 3, 1, -1, -1, 1, 3,
55  3, 1, -1, 1, 1, 3, -1, 1, 1, 1, -3, 1, 1, 1, -1, 3,
56  -3, 1, -1, 1, -1, 3, -1, 1, -1, 1, -3, 1, -1, 1, -1, 3,
57  3, -1, -1, 1, 1, -3, -1, 1, 1, -1, -3, 1, 1, -1, -1, 3,
58  -3, -1, -1, 1, -1, -3, -1, 1, -1, -1, -3, 1, -1, -1, -1, 3,
59  3, 1, 1, -1, 1, 3, 1, -1, 1, 1, 3, -1, 1, 1, 1, -3,
60  -3, 1, 1, -1, -1, 3, 1, -1, -1, 1, 3, -1, -1, 1, 1, -3,
61  3, -1, 1, -1, 1, -3, 1, -1, 1, -1, 3, -1, 1, -1, 1, -3,
62  -3, -1, 1, -1, -1, -3, 1, -1, -1, -1, 3, -1, -1, -1, 1, -3,
63  3, 1, -1, -1, 1, 3, -1, -1, 1, 1, -3, -1, 1, 1, -1, -3,
64  -3, 1, -1, -1, -1, 3, -1, -1, -1, 1, -3, -1, -1, 1, -1, -3,
65  3, -1, -1, -1, 1, -3, -1, -1, 1, -1, -3, -1, 1, -1, -1, -3,
66  -3, -1, -1, -1, -1, -3, -1, -1, -1, -1, -3, -1, -1, -1, -1, -3,
67  };
68 
69  private static int[] p2D = new int[] { 0, 0, 1, -1, 0, 0, -1, 1, 0, 2, 1, 1, 1, 2, 2, 0, 1, 2, 0, 2, 1, 0, 0, 0 };
70  private static int[] p3D = new int[] { 0, 0, 1, -1, 0, 0, 1, 0, -1, 0, 0, -1, 1, 0, 0, 0, 1, -1, 0, 0, -1, 0, 1, 0, 0, -1, 1, 0, 2, 1, 1, 0, 1, 1, 1, -1, 0, 2, 1, 0, 1, 1, 1, -1, 1, 0, 2, 0, 1, 1, 1, -1, 1, 1, 1, 3, 2, 1, 0, 3, 1, 2, 0, 1, 3, 2, 0, 1, 3, 1, 0, 2, 1, 3, 0, 2, 1, 3, 0, 1, 2, 1, 1, 1, 0, 0, 2, 2, 0, 0, 1, 1, 0, 1, 0, 2, 0, 2, 0, 1, 1, 0, 0, 1, 2, 0, 0, 2, 2, 0, 0, 0, 0, 1, 1, -1, 1, 2, 0, 0, 0, 0, 1, -1, 1, 1, 2, 0, 0, 0, 0, 1, 1, 1, -1, 2, 3, 1, 1, 1, 2, 0, 0, 2, 2, 3, 1, 1, 1, 2, 2, 0, 0, 2, 3, 1, 1, 1, 2, 0, 2, 0, 2, 1, 1, -1, 1, 2, 0, 0, 2, 2, 1, 1, -1, 1, 2, 2, 0, 0, 2, 1, -1, 1, 1, 2, 0, 0, 2, 2, 1, -1, 1, 1, 2, 0, 2, 0, 2, 1, 1, 1, -1, 2, 2, 0, 0, 2, 1, 1, 1, -1, 2, 0, 2, 0 };
71  private static int[] p4D = new int[] { 0, 0, 1, -1, 0, 0, 0, 1, 0, -1, 0, 0, 1, 0, 0, -1, 0, 0, -1, 1, 0, 0, 0, 0, 1, -1, 0, 0, 0, 1, 0, -1, 0, 0, -1, 0, 1, 0, 0, 0, -1, 1, 0, 0, 0, 0, 1, -1, 0, 0, -1, 0, 0, 1, 0, 0, -1, 0, 1, 0, 0, 0, -1, 1, 0, 2, 1, 1, 0, 0, 1, 1, 1, -1, 0, 1, 1, 1, 0, -1, 0, 2, 1, 0, 1, 0, 1, 1, -1, 1, 0, 1, 1, 0, 1, -1, 0, 2, 0, 1, 1, 0, 1, -1, 1, 1, 0, 1, 0, 1, 1, -1, 0, 2, 1, 0, 0, 1, 1, 1, -1, 0, 1, 1, 1, 0, -1, 1, 0, 2, 0, 1, 0, 1, 1, -1, 1, 0, 1, 1, 0, 1, -1, 1, 0, 2, 0, 0, 1, 1, 1, -1, 0, 1, 1, 1, 0, -1, 1, 1, 1, 4, 2, 1, 1, 0, 4, 1, 2, 1, 0, 4, 1, 1, 2, 0, 1, 4, 2, 1, 0, 1, 4, 1, 2, 0, 1, 4, 1, 1, 0, 2, 1, 4, 2, 0, 1, 1, 4, 1, 0, 2, 1, 4, 1, 0, 1, 2, 1, 4, 0, 2, 1, 1, 4, 0, 1, 2, 1, 4, 0, 1, 1, 2, 1, 2, 1, 1, 0, 0, 3, 2, 1, 0, 0, 3, 1, 2, 0, 0, 1, 2, 1, 0, 1, 0, 3, 2, 0, 1, 0, 3, 1, 0, 2, 0, 1, 2, 0, 1, 1, 0, 3, 0, 2, 1, 0, 3, 0, 1, 2, 0, 1, 2, 1, 0, 0, 1, 3, 2, 0, 0, 1, 3, 1, 0, 0, 2, 1, 2, 0, 1, 0, 1, 3, 0, 2, 0, 1, 3, 0, 1, 0, 2, 1, 2, 0, 0, 1, 1, 3, 0, 0, 2, 1, 3, 0, 0, 1, 2, 2, 3, 1, 1, 1, 0, 2, 1, 1, 1, -1, 2, 2, 0, 0, 0, 2, 3, 1, 1, 0, 1, 2, 1, 1, -1, 1, 2, 2, 0, 0, 0, 2, 3, 1, 0, 1, 1, 2, 1, -1, 1, 1, 2, 2, 0, 0, 0, 2, 3, 1, 1, 1, 0, 2, 1, 1, 1, -1, 2, 0, 2, 0, 0, 2, 3, 1, 1, 0, 1, 2, 1, 1, -1, 1, 2, 0, 2, 0, 0, 2, 3, 0, 1, 1, 1, 2, -1, 1, 1, 1, 2, 0, 2, 0, 0, 2, 3, 1, 1, 1, 0, 2, 1, 1, 1, -1, 2, 0, 0, 2, 0, 2, 3, 1, 0, 1, 1, 2, 1, -1, 1, 1, 2, 0, 0, 2, 0, 2, 3, 0, 1, 1, 1, 2, -1, 1, 1, 1, 2, 0, 0, 2, 0, 2, 3, 1, 1, 0, 1, 2, 1, 1, -1, 1, 2, 0, 0, 0, 2, 2, 3, 1, 0, 1, 1, 2, 1, -1, 1, 1, 2, 0, 0, 0, 2, 2, 3, 0, 1, 1, 1, 2, -1, 1, 1, 1, 2, 0, 0, 0, 2, 2, 1, 1, 1, -1, 0, 1, 1, 1, 0, -1, 0, 0, 0, 0, 0, 2, 1, 1, -1, 1, 0, 1, 1, 0, 1, -1, 0, 0, 0, 0, 0, 2, 1, -1, 1, 1, 0, 1, 0, 1, 1, -1, 0, 0, 0, 0, 0, 2, 1, 1, -1, 0, 1, 1, 1, 0, -1, 1, 0, 0, 0, 0, 0, 2, 1, -1, 1, 0, 1, 1, 0, 1, -1, 1, 0, 0, 0, 0, 0, 2, 1, -1, 0, 1, 1, 1, 0, -1, 1, 1, 0, 0, 0, 0, 0, 2, 1, 1, 1, -1, 0, 1, 1, 1, 0, -1, 2, 2, 0, 0, 0, 2, 1, 1, -1, 1, 0, 1, 1, 0, 1, -1, 2, 2, 0, 0, 0, 2, 1, 1, -1, 0, 1, 1, 1, 0, -1, 1, 2, 2, 0, 0, 0, 2, 1, 1, 1, -1, 0, 1, 1, 1, 0, -1, 2, 0, 2, 0, 0, 2, 1, -1, 1, 1, 0, 1, 0, 1, 1, -1, 2, 0, 2, 0, 0, 2, 1, -1, 1, 0, 1, 1, 0, 1, -1, 1, 2, 0, 2, 0, 0, 2, 1, 1, -1, 1, 0, 1, 1, 0, 1, -1, 2, 0, 0, 2, 0, 2, 1, -1, 1, 1, 0, 1, 0, 1, 1, -1, 2, 0, 0, 2, 0, 2, 1, -1, 0, 1, 1, 1, 0, -1, 1, 1, 2, 0, 0, 2, 0, 2, 1, 1, -1, 0, 1, 1, 1, 0, -1, 1, 2, 0, 0, 0, 2, 2, 1, -1, 1, 0, 1, 1, 0, 1, -1, 1, 2, 0, 0, 0, 2, 2, 1, -1, 0, 1, 1, 1, 0, -1, 1, 1, 2, 0, 0, 0, 2, 3, 1, 1, 0, 0, 0, 2, 2, 0, 0, 0, 2, 1, 1, 1, -1, 3, 1, 0, 1, 0, 0, 2, 0, 2, 0, 0, 2, 1, 1, 1, -1, 3, 1, 0, 0, 1, 0, 2, 0, 0, 2, 0, 2, 1, 1, 1, -1, 3, 1, 1, 0, 0, 0, 2, 2, 0, 0, 0, 2, 1, 1, -1, 1, 3, 1, 0, 1, 0, 0, 2, 0, 2, 0, 0, 2, 1, 1, -1, 1, 3, 1, 0, 0, 0, 1, 2, 0, 0, 0, 2, 2, 1, 1, -1, 1, 3, 1, 1, 0, 0, 0, 2, 2, 0, 0, 0, 2, 1, -1, 1, 1, 3, 1, 0, 0, 1, 0, 2, 0, 0, 2, 0, 2, 1, -1, 1, 1, 3, 1, 0, 0, 0, 1, 2, 0, 0, 0, 2, 2, 1, -1, 1, 1, 3, 1, 0, 1, 0, 0, 2, 0, 2, 0, 0, 2, -1, 1, 1, 1, 3, 1, 0, 0, 1, 0, 2, 0, 0, 2, 0, 2, -1, 1, 1, 1, 3, 1, 0, 0, 0, 1, 2, 0, 0, 0, 2, 2, -1, 1, 1, 1, 3, 3, 2, 1, 0, 0, 3, 1, 2, 0, 0, 4, 1, 1, 1, 1, 3, 3, 2, 0, 1, 0, 3, 1, 0, 2, 0, 4, 1, 1, 1, 1, 3, 3, 0, 2, 1, 0, 3, 0, 1, 2, 0, 4, 1, 1, 1, 1, 3, 3, 2, 0, 0, 1, 3, 1, 0, 0, 2, 4, 1, 1, 1, 1, 3, 3, 0, 2, 0, 1, 3, 0, 1, 0, 2, 4, 1, 1, 1, 1, 3, 3, 0, 0, 2, 1, 3, 0, 0, 1, 2, 4, 1, 1, 1, 1, 3, 3, 2, 1, 0, 0, 3, 1, 2, 0, 0, 2, 1, 1, 1, -1, 3, 3, 2, 0, 1, 0, 3, 1, 0, 2, 0, 2, 1, 1, 1, -1, 3, 3, 0, 2, 1, 0, 3, 0, 1, 2, 0, 2, 1, 1, 1, -1, 3, 3, 2, 1, 0, 0, 3, 1, 2, 0, 0, 2, 1, 1, -1, 1, 3, 3, 2, 0, 0, 1, 3, 1, 0, 0, 2, 2, 1, 1, -1, 1, 3, 3, 0, 2, 0, 1, 3, 0, 1, 0, 2, 2, 1, 1, -1, 1, 3, 3, 2, 0, 1, 0, 3, 1, 0, 2, 0, 2, 1, -1, 1, 1, 3, 3, 2, 0, 0, 1, 3, 1, 0, 0, 2, 2, 1, -1, 1, 1, 3, 3, 0, 0, 2, 1, 3, 0, 0, 1, 2, 2, 1, -1, 1, 1, 3, 3, 0, 2, 1, 0, 3, 0, 1, 2, 0, 2, -1, 1, 1, 1, 3, 3, 0, 2, 0, 1, 3, 0, 1, 0, 2, 2, -1, 1, 1, 1, 3, 3, 0, 0, 2, 1, 3, 0, 0, 1, 2, 2, -1, 1, 1, 1 };
72  private static int[] lookupPairs2D = new int[] { 0, 1, 1, 0, 4, 1, 17, 0, 20, 2, 21, 2, 22, 5, 23, 5, 26, 4, 39, 3, 42, 4, 43, 3 };
73  private static int[] lookupPairs3D = new int[] { 0, 2, 1, 1, 2, 2, 5, 1, 6, 0, 7, 0, 32, 2, 34, 2, 129, 1, 133, 1, 160, 5, 161, 5, 518, 0, 519, 0, 546, 4, 550, 4, 645, 3, 647, 3, 672, 5, 673, 5, 674, 4, 677, 3, 678, 4, 679, 3, 680, 13, 681, 13, 682, 12, 685, 14, 686, 12, 687, 14, 712, 20, 714, 18, 809, 21, 813, 23, 840, 20, 841, 21, 1198, 19, 1199, 22, 1226, 18, 1230, 19, 1325, 23, 1327, 22, 1352, 15, 1353, 17, 1354, 15, 1357, 17, 1358, 16, 1359, 16, 1360, 11, 1361, 10, 1362, 11, 1365, 10, 1366, 9, 1367, 9, 1392, 11, 1394, 11, 1489, 10, 1493, 10, 1520, 8, 1521, 8, 1878, 9, 1879, 9, 1906, 7, 1910, 7, 2005, 6, 2007, 6, 2032, 8, 2033, 8, 2034, 7, 2037, 6, 2038, 7, 2039, 6 };
74  private static int[] lookupPairs4D = new int[] { 0, 3, 1, 2, 2, 3, 5, 2, 6, 1, 7, 1, 8, 3, 9, 2, 10, 3, 13, 2, 16, 3, 18, 3, 22, 1, 23, 1, 24, 3, 26, 3, 33, 2, 37, 2, 38, 1, 39, 1, 41, 2, 45, 2, 54, 1, 55, 1, 56, 0, 57, 0, 58, 0, 59, 0, 60, 0, 61, 0, 62, 0, 63, 0, 256, 3, 258, 3, 264, 3, 266, 3, 272, 3, 274, 3, 280, 3, 282, 3, 2049, 2, 2053, 2, 2057, 2, 2061, 2, 2081, 2, 2085, 2, 2089, 2, 2093, 2, 2304, 9, 2305, 9, 2312, 9, 2313, 9, 16390, 1, 16391, 1, 16406, 1, 16407, 1, 16422, 1, 16423, 1, 16438, 1, 16439, 1, 16642, 8, 16646, 8, 16658, 8, 16662, 8, 18437, 6, 18439, 6, 18469, 6, 18471, 6, 18688, 9, 18689, 9, 18690, 8, 18693, 6, 18694, 8, 18695, 6, 18696, 9, 18697, 9, 18706, 8, 18710, 8, 18725, 6, 18727, 6, 131128, 0, 131129, 0, 131130, 0, 131131, 0, 131132, 0, 131133, 0, 131134, 0, 131135, 0, 131352, 7, 131354, 7, 131384, 7, 131386, 7, 133161, 5, 133165, 5, 133177, 5, 133181, 5, 133376, 9, 133377, 9, 133384, 9, 133385, 9, 133400, 7, 133402, 7, 133417, 5, 133421, 5, 133432, 7, 133433, 5, 133434, 7, 133437, 5, 147510, 4, 147511, 4, 147518, 4, 147519, 4, 147714, 8, 147718, 8, 147730, 8, 147734, 8, 147736, 7, 147738, 7, 147766, 4, 147767, 4, 147768, 7, 147770, 7, 147774, 4, 147775, 4, 149509, 6, 149511, 6, 149541, 6, 149543, 6, 149545, 5, 149549, 5, 149558, 4, 149559, 4, 149561, 5, 149565, 5, 149566, 4, 149567, 4, 149760, 9, 149761, 9, 149762, 8, 149765, 6, 149766, 8, 149767, 6, 149768, 9, 149769, 9, 149778, 8, 149782, 8, 149784, 7, 149786, 7, 149797, 6, 149799, 6, 149801, 5, 149805, 5, 149814, 4, 149815, 4, 149816, 7, 149817, 5, 149818, 7, 149821, 5, 149822, 4, 149823, 4, 149824, 37, 149825, 37, 149826, 36, 149829, 34, 149830, 36, 149831, 34, 149832, 37, 149833, 37, 149842, 36, 149846, 36, 149848, 35, 149850, 35, 149861, 34, 149863, 34, 149865, 33, 149869, 33, 149878, 32, 149879, 32, 149880, 35, 149881, 33, 149882, 35, 149885, 33, 149886, 32, 149887, 32, 150080, 49, 150082, 48, 150088, 49, 150098, 48, 150104, 47, 150106, 47, 151873, 46, 151877, 45, 151881, 46, 151909, 45, 151913, 44, 151917, 44, 152128, 49, 152129, 46, 152136, 49, 152137, 46, 166214, 43, 166215, 42, 166230, 43, 166247, 42, 166262, 41, 166263, 41, 166466, 48, 166470, 43, 166482, 48, 166486, 43, 168261, 45, 168263, 42, 168293, 45, 168295, 42, 168512, 31, 168513, 28, 168514, 31, 168517, 28, 168518, 25, 168519, 25, 280952, 40, 280953, 39, 280954, 40, 280957, 39, 280958, 38, 280959, 38, 281176, 47, 281178, 47, 281208, 40, 281210, 40, 282985, 44, 282989, 44, 283001, 39, 283005, 39, 283208, 30, 283209, 27, 283224, 30, 283241, 27, 283256, 22, 283257, 22, 297334, 41, 297335, 41, 297342, 38, 297343, 38, 297554, 29, 297558, 24, 297562, 29, 297590, 24, 297594, 21, 297598, 21, 299365, 26, 299367, 23, 299373, 26, 299383, 23, 299389, 20, 299391, 20, 299584, 31, 299585, 28, 299586, 31, 299589, 28, 299590, 25, 299591, 25, 299592, 30, 299593, 27, 299602, 29, 299606, 24, 299608, 30, 299610, 29, 299621, 26, 299623, 23, 299625, 27, 299629, 26, 299638, 24, 299639, 23, 299640, 22, 299641, 22, 299642, 21, 299645, 20, 299646, 21, 299647, 20, 299648, 61, 299649, 60, 299650, 61, 299653, 60, 299654, 59, 299655, 59, 299656, 58, 299657, 57, 299666, 55, 299670, 54, 299672, 58, 299674, 55, 299685, 52, 299687, 51, 299689, 57, 299693, 52, 299702, 54, 299703, 51, 299704, 56, 299705, 56, 299706, 53, 299709, 50, 299710, 53, 299711, 50, 299904, 61, 299906, 61, 299912, 58, 299922, 55, 299928, 58, 299930, 55, 301697, 60, 301701, 60, 301705, 57, 301733, 52, 301737, 57, 301741, 52, 301952, 79, 301953, 79, 301960, 76, 301961, 76, 316038, 59, 316039, 59, 316054, 54, 316071, 51, 316086, 54, 316087, 51, 316290, 78, 316294, 78, 316306, 73, 316310, 73, 318085, 77, 318087, 77, 318117, 70, 318119, 70, 318336, 79, 318337, 79, 318338, 78, 318341, 77, 318342, 78, 318343, 77, 430776, 56, 430777, 56, 430778, 53, 430781, 50, 430782, 53, 430783, 50, 431000, 75, 431002, 72, 431032, 75, 431034, 72, 432809, 74, 432813, 69, 432825, 74, 432829, 69, 433032, 76, 433033, 76, 433048, 75, 433065, 74, 433080, 75, 433081, 74, 447158, 71, 447159, 68, 447166, 71, 447167, 68, 447378, 73, 447382, 73, 447386, 72, 447414, 71, 447418, 72, 447422, 71, 449189, 70, 449191, 70, 449197, 69, 449207, 68, 449213, 69, 449215, 68, 449408, 67, 449409, 67, 449410, 66, 449413, 64, 449414, 66, 449415, 64, 449416, 67, 449417, 67, 449426, 66, 449430, 66, 449432, 65, 449434, 65, 449445, 64, 449447, 64, 449449, 63, 449453, 63, 449462, 62, 449463, 62, 449464, 65, 449465, 63, 449466, 65, 449469, 63, 449470, 62, 449471, 62, 449472, 19, 449473, 19, 449474, 18, 449477, 16, 449478, 18, 449479, 16, 449480, 19, 449481, 19, 449490, 18, 449494, 18, 449496, 17, 449498, 17, 449509, 16, 449511, 16, 449513, 15, 449517, 15, 449526, 14, 449527, 14, 449528, 17, 449529, 15, 449530, 17, 449533, 15, 449534, 14, 449535, 14, 449728, 19, 449729, 19, 449730, 18, 449734, 18, 449736, 19, 449737, 19, 449746, 18, 449750, 18, 449752, 17, 449754, 17, 449784, 17, 449786, 17, 451520, 19, 451521, 19, 451525, 16, 451527, 16, 451528, 19, 451529, 19, 451557, 16, 451559, 16, 451561, 15, 451565, 15, 451577, 15, 451581, 15, 451776, 19, 451777, 19, 451784, 19, 451785, 19, 465858, 18, 465861, 16, 465862, 18, 465863, 16, 465874, 18, 465878, 18, 465893, 16, 465895, 16, 465910, 14, 465911, 14, 465918, 14, 465919, 14, 466114, 18, 466118, 18, 466130, 18, 466134, 18, 467909, 16, 467911, 16, 467941, 16, 467943, 16, 468160, 13, 468161, 13, 468162, 13, 468163, 13, 468164, 13, 468165, 13, 468166, 13, 468167, 13, 580568, 17, 580570, 17, 580585, 15, 580589, 15, 580598, 14, 580599, 14, 580600, 17, 580601, 15, 580602, 17, 580605, 15, 580606, 14, 580607, 14, 580824, 17, 580826, 17, 580856, 17, 580858, 17, 582633, 15, 582637, 15, 582649, 15, 582653, 15, 582856, 12, 582857, 12, 582872, 12, 582873, 12, 582888, 12, 582889, 12, 582904, 12, 582905, 12, 596982, 14, 596983, 14, 596990, 14, 596991, 14, 597202, 11, 597206, 11, 597210, 11, 597214, 11, 597234, 11, 597238, 11, 597242, 11, 597246, 11, 599013, 10, 599015, 10, 599021, 10, 599023, 10, 599029, 10, 599031, 10, 599037, 10, 599039, 10, 599232, 13, 599233, 13, 599234, 13, 599235, 13, 599236, 13, 599237, 13, 599238, 13, 599239, 13, 599240, 12, 599241, 12, 599250, 11, 599254, 11, 599256, 12, 599257, 12, 599258, 11, 599262, 11, 599269, 10, 599271, 10, 599272, 12, 599273, 12, 599277, 10, 599279, 10, 599282, 11, 599285, 10, 599286, 11, 599287, 10, 599288, 12, 599289, 12, 599290, 11, 599293, 10, 599294, 11, 599295, 10 };
75  private static int[][] base2D = new int[][]
76  {
77  new int[] { 1, 1, 0, 1, 0, 1, 0, 0, 0 },
78  new int[] { 1, 1, 0, 1, 0, 1, 2, 1, 1 }
79  };
80  private static int[][] base3D = new int[][]
81  {
82  new int[] { 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1 },
83  new int[] { 2, 1, 1, 0, 2, 1, 0, 1, 2, 0, 1, 1, 3, 1, 1, 1 },
84  new int[] { 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 2, 1, 1, 0, 2, 1, 0, 1, 2, 0, 1, 1 }
85  };
86  private static int[][] base4D = new int[][]
87  {
88  new int[] { 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 1 },
89  new int[] { 3, 1, 1, 1, 0, 3, 1, 1, 0, 1, 3, 1, 0, 1, 1, 3, 0, 1, 1, 1, 4, 1, 1, 1, 1 },
90  new int[] { 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 1, 2, 1, 1, 0, 0, 2, 1, 0, 1, 0, 2, 1, 0, 0, 1, 2, 0, 1, 1, 0, 2, 0, 1, 0, 1, 2, 0, 0, 1, 1 },
91  new int[] { 3, 1, 1, 1, 0, 3, 1, 1, 0, 1, 3, 1, 0, 1, 1, 3, 0, 1, 1, 1, 2, 1, 1, 0, 0, 2, 1, 0, 1, 0, 2, 1, 0, 0, 1, 2, 0, 1, 1, 0, 2, 0, 1, 0, 1, 2, 0, 0, 1, 1 }
92  };
93 
94  private static Contribution2[] lookup2D;
95  private static Contribution3[] lookup3D;
96  private static Contribution4[] lookup4D;
97 
98  static FastSimplexNoise()
99  {
100  var contributions2D = new Contribution2[p2D.Length / 4];
101  for (int i = 0; i < p2D.Length; i += 4)
102  {
103  var baseSet = base2D[p2D[i]];
104  Contribution2 previous = null, current = null;
105  for (int k = 0; k < baseSet.Length; k += 3)
106  {
107  current = new Contribution2(baseSet[k], baseSet[k + 1], baseSet[k + 2]);
108  if (previous == null)
109  {
110  contributions2D[i / 4] = current;
111  }
112  else
113  {
114  previous.Next = current;
115  }
116  previous = current;
117  }
118  current.Next = new Contribution2(p2D[i + 1], p2D[i + 2], p2D[i + 3]);
119  }
120 
121  lookup2D = new Contribution2[64];
122  for (var i = 0; i < lookupPairs2D.Length; i += 2)
123  {
124  lookup2D[lookupPairs2D[i]] = contributions2D[lookupPairs2D[i + 1]];
125  }
126 
127  var contributions3D = new Contribution3[p3D.Length / 9];
128  for (int i = 0; i < p3D.Length; i += 9)
129  {
130  var baseSet = base3D[p3D[i]];
131  Contribution3 previous = null, current = null;
132  for (int k = 0; k < baseSet.Length; k += 4)
133  {
134  current = new Contribution3(baseSet[k], baseSet[k + 1], baseSet[k + 2], baseSet[k + 3]);
135  if (previous == null)
136  {
137  contributions3D[i / 9] = current;
138  }
139  else
140  {
141  previous.Next = current;
142  }
143  previous = current;
144  }
145  current.Next = new Contribution3(p3D[i + 1], p3D[i + 2], p3D[i + 3], p3D[i + 4]);
146  current.Next.Next = new Contribution3(p3D[i + 5], p3D[i + 6], p3D[i + 7], p3D[i + 8]);
147  }
148 
149  lookup3D = new Contribution3[2048];
150  for (var i = 0; i < lookupPairs3D.Length; i += 2)
151  {
152  lookup3D[lookupPairs3D[i]] = contributions3D[lookupPairs3D[i + 1]];
153  }
154 
155  var contributions4D = new Contribution4[p4D.Length / 16];
156  for (int i = 0; i < p4D.Length; i += 16)
157  {
158  var baseSet = base4D[p4D[i]];
159  Contribution4 previous = null, current = null;
160  for (int k = 0; k < baseSet.Length; k += 5)
161  {
162  current = new Contribution4(baseSet[k], baseSet[k + 1], baseSet[k + 2], baseSet[k + 3], baseSet[k + 4]);
163  if (previous == null)
164  {
165  contributions4D[i / 16] = current;
166  }
167  else
168  {
169  previous.Next = current;
170  }
171  previous = current;
172  }
173  current.Next = new Contribution4(p4D[i + 1], p4D[i + 2], p4D[i + 3], p4D[i + 4], p4D[i + 5]);
174  current.Next.Next = new Contribution4(p4D[i + 6], p4D[i + 7], p4D[i + 8], p4D[i + 9], p4D[i + 10]);
175  current.Next.Next.Next = new Contribution4(p4D[i + 11], p4D[i + 12], p4D[i + 13], p4D[i + 14], p4D[i + 15]);
176  }
177 
178  lookup4D = new Contribution4[1048576];
179  for (var i = 0; i < lookupPairs4D.Length; i += 2)
180  {
181  lookup4D[lookupPairs4D[i]] = contributions4D[lookupPairs4D[i + 1]];
182  }
183  }
184 
185  private static int FastFloor(double x)
186  {
187  var xi = (int)x;
188  return x < xi ? xi - 1 : xi;
189  }
190 
191  public FastSimplexNoise() : this(DateTime.Now.Ticks) { }
192 
193  public FastSimplexNoise(long seed)
194  {
195  perm = new byte[256];
196  perm2D = new byte[256];
197  perm3D = new byte[256];
198  perm4D = new byte[256];
199  var source = new byte[256];
200  for (int i = 0; i < 256; i++)
201  {
202  source[i] = (byte)i;
203  }
204  seed = seed * SEEDVAL_1 + SEEDVAL_2;
205  seed = seed * SEEDVAL_1 + SEEDVAL_2;
206  seed = seed * SEEDVAL_1 + SEEDVAL_2;
207  for (int i = 255; i >= 0; i--)
208  {
209  seed = seed * SEEDVAL_1 + SEEDVAL_2;
210  int r = (int)((seed + 31) % (i + 1));
211  if (r < 0)
212  {
213  r += (i + 1);
214  }
215  perm[i] = source[r];
216  perm2D[i] = (byte)(perm[i] & 0x0E);
217  perm3D[i] = (byte)((perm[i] % 24) * 3);
218  perm4D[i] = (byte)(perm[i] & 0xFC);
219  source[r] = source[i];
220  }
221  }
222 
223  public double Evaluate(double x, double y)
224  {
225  var stretchOffset = (x + y) * STRETCH_2D;
226  var xs = x + stretchOffset;
227  var ys = y + stretchOffset;
228 
229  var xsb = FastFloor(xs);
230  var ysb = FastFloor(ys);
231 
232  var squishOffset = (xsb + ysb) * SQUISH_2D;
233  var dx0 = x - (xsb + squishOffset);
234  var dy0 = y - (ysb + squishOffset);
235 
236  var xins = xs - xsb;
237  var yins = ys - ysb;
238 
239  var inSum = xins + yins;
240 
241  var hash =
242  (int)(xins - yins + 1) |
243  (int)(inSum) << 1 |
244  (int)(inSum + yins) << 2 |
245  (int)(inSum + xins) << 4;
246 
247  var c = lookup2D[hash];
248 
249  var value = 0.0;
250  while (c != null)
251  {
252  var dx = dx0 + c.dx;
253  var dy = dy0 + c.dy;
254  var attn = 2 - dx * dx - dy * dy;
255  if (attn > 0)
256  {
257  var px = xsb + c.xsb;
258  var py = ysb + c.ysb;
259 
260  var i = perm2D[(perm[px & 0xFF] + py) & 0xFF];
261  var valuePart = gradients2D[i] * dx + gradients2D[i + 1] * dy;
262 
263  attn *= attn;
264  value += attn * attn * valuePart;
265  }
266  c = c.Next;
267  }
268  return value * NORM_2D;
269  }
270 
271  public double Evaluate(double x, double y, double z)
272  {
273  var stretchOffset = (x + y + z) * STRETCH_3D;
274  var xs = x + stretchOffset;
275  var ys = y + stretchOffset;
276  var zs = z + stretchOffset;
277 
278  var xsb = FastFloor(xs);
279  var ysb = FastFloor(ys);
280  var zsb = FastFloor(zs);
281 
282  var squishOffset = (xsb + ysb + zsb) * SQUISH_3D;
283  var dx0 = x - (xsb + squishOffset);
284  var dy0 = y - (ysb + squishOffset);
285  var dz0 = z - (zsb + squishOffset);
286 
287  var xins = xs - xsb;
288  var yins = ys - ysb;
289  var zins = zs - zsb;
290 
291  var inSum = xins + yins + zins;
292 
293  var hash =
294  (int)(yins - zins + 1) |
295  (int)(xins - yins + 1) << 1 |
296  (int)(xins - zins + 1) << 2 |
297  (int)inSum << 3 |
298  (int)(inSum + zins) << 5 |
299  (int)(inSum + yins) << 7 |
300  (int)(inSum + xins) << 9;
301 
302  var c = lookup3D[hash];
303 
304  var value = 0.0;
305  while (c != null)
306  {
307  var dx = dx0 + c.dx;
308  var dy = dy0 + c.dy;
309  var dz = dz0 + c.dz;
310  var attn = 2 - dx * dx - dy * dy - dz * dz;
311  if (attn > 0)
312  {
313  var px = xsb + c.xsb;
314  var py = ysb + c.ysb;
315  var pz = zsb + c.zsb;
316 
317  var i = perm3D[(perm[(perm[px & 0xFF] + py) & 0xFF] + pz) & 0xFF];
318  var valuePart = gradients3D[i] * dx + gradients3D[i + 1] * dy + gradients3D[i + 2] * dz;
319 
320  attn *= attn;
321  value += attn * attn * valuePart;
322  }
323 
324  c = c.Next;
325  }
326  return value * NORM_3D;
327  }
328 
329  public double Evaluate(double x, double y, double z, double w)
330  {
331  var stretchOffset = (x + y + z + w) * STRETCH_4D;
332  var xs = x + stretchOffset;
333  var ys = y + stretchOffset;
334  var zs = z + stretchOffset;
335  var ws = w + stretchOffset;
336 
337  var xsb = FastFloor(xs);
338  var ysb = FastFloor(ys);
339  var zsb = FastFloor(zs);
340  var wsb = FastFloor(ws);
341 
342  var squishOffset = (xsb + ysb + zsb + wsb) * SQUISH_4D;
343  var dx0 = x - (xsb + squishOffset);
344  var dy0 = y - (ysb + squishOffset);
345  var dz0 = z - (zsb + squishOffset);
346  var dw0 = w - (wsb + squishOffset);
347 
348  var xins = xs - xsb;
349  var yins = ys - ysb;
350  var zins = zs - zsb;
351  var wins = ws - wsb;
352 
353  var inSum = xins + yins + zins + wins;
354 
355  var hash =
356  (int)(zins - wins + 1) |
357  (int)(yins - zins + 1) << 1 |
358  (int)(yins - wins + 1) << 2 |
359  (int)(xins - yins + 1) << 3 |
360  (int)(xins - zins + 1) << 4 |
361  (int)(xins - wins + 1) << 5 |
362  (int)inSum << 6 |
363  (int)(inSum + wins) << 8 |
364  (int)(inSum + zins) << 11 |
365  (int)(inSum + yins) << 14 |
366  (int)(inSum + xins) << 17;
367 
368  var c = lookup4D[hash];
369 
370  var value = 0.0;
371  while (c != null)
372  {
373  var dx = dx0 + c.dx;
374  var dy = dy0 + c.dy;
375  var dz = dz0 + c.dz;
376  var dw = dw0 + c.dw;
377  var attn = 2 - dx * dx - dy * dy - dz * dz - dw * dw;
378  if (attn > 0)
379  {
380  var px = xsb + c.xsb;
381  var py = ysb + c.ysb;
382  var pz = zsb + c.zsb;
383  var pw = wsb + c.wsb;
384 
385  var i = perm4D[(perm[(perm[(perm[px & 0xFF] + py) & 0xFF] + pz) & 0xFF] + pw) & 0xFF];
386  var valuePart = gradients4D[i] * dx + gradients4D[i + 1] * dy + gradients4D[i + 2] * dz + gradients4D[i + 3] * dw;
387 
388  attn *= attn;
389  value += attn * attn * valuePart;
390  }
391 
392  c = c.Next;
393  }
394  return value * NORM_4D;
395  }
396 
397  private class Contribution2
398  {
399  public double dx, dy;
400  public int xsb, ysb;
401  public Contribution2 Next;
402 
403  public Contribution2(double multiplier, int xsb, int ysb)
404  {
405  dx = -xsb - multiplier * SQUISH_2D;
406  dy = -ysb - multiplier * SQUISH_2D;
407  this.xsb = xsb;
408  this.ysb = ysb;
409  }
410  }
411 
412  private class Contribution3
413  {
414  public double dx, dy, dz;
415  public int xsb, ysb, zsb;
416  public Contribution3 Next;
417 
418  public Contribution3(double multiplier, int xsb, int ysb, int zsb)
419  {
420  dx = -xsb - multiplier * SQUISH_3D;
421  dy = -ysb - multiplier * SQUISH_3D;
422  dz = -zsb - multiplier * SQUISH_3D;
423  this.xsb = xsb;
424  this.ysb = ysb;
425  this.zsb = zsb;
426  }
427  }
428 
429  private class Contribution4
430  {
431  public double dx, dy, dz, dw;
432  public int xsb, ysb, zsb, wsb;
433  public Contribution4 Next;
434 
435  public Contribution4(double multiplier, int xsb, int ysb, int zsb, int wsb)
436  {
437  dx = -xsb - multiplier * SQUISH_4D;
438  dy = -ysb - multiplier * SQUISH_4D;
439  dz = -zsb - multiplier * SQUISH_4D;
440  dw = -wsb - multiplier * SQUISH_4D;
441  this.xsb = xsb;
442  this.ysb = ysb;
443  this.zsb = zsb;
444  this.wsb = wsb;
445  }
446  }
447 }
FastSimplexNoise.Evaluate
double Evaluate(double x, double y, double z, double w)
Definition: FastSimplexNoise.cs:329
FastSimplexNoise
A conglomeration of open-source simplex libraries in C# with an emphasis on performance ...
Definition: FastSimplexNoise.cs:9
FastSimplexNoise.Evaluate
double Evaluate(double x, double y)
Definition: FastSimplexNoise.cs:223
FastSimplexNoise.Evaluate
double Evaluate(double x, double y, double z)
Definition: FastSimplexNoise.cs:271
FastSimplexNoise.FastSimplexNoise
FastSimplexNoise(long seed)
Definition: FastSimplexNoise.cs:193