00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015 #ifndef NOISES_H
00016 #define NOISES_H 1
00017
00018
00019 #ifndef FILTERWIDTH_H
00020 #include "filterwidth.h"
00021 #endif
00022
00023 #ifndef PATTERNS_H
00024 #include "patterns.h"
00025 #endif
00026
00027
00028
00029 #ifndef snoise
00030
00031
00032
00033
00034
00035
00036
00037
00038 #define snoise(p) (2 * (float noise(p)) - 1)
00039 #define snoiset(p,t) (2 * (float noise(p,t)) - 1)
00040 #define snoisexy(x,y) (2 * (float noise(x,y)) - 1)
00041 #define vsnoise(p) (2 * (vector noise(p)) - 1)
00042 #define vsnoiset(p,t) (2 * (vector noise(p,t)) - 1)
00043 #endif
00044
00045
00046
00047
00048
00049 #define filteredsnoise(p,width) (snoise(p) * (1-smoothstep (0.2,0.75,width)))
00050 #define filteredvsnoise(p,width) (vsnoise(p) * (1-smoothstep (0.2,0.75,width)))
00051 #define filteredsnoiset(p,width,t) \
00052 (snoiset(p,t) * (1-smoothstep (0.2,0.75,width)))
00053 #define filteredvsnoiset(p,width,t) \
00054 (vsnoiset(p,t) * (1-smoothstep (0.2,0.75,width)))
00055
00056
00057
00058
00059
00060
00061
00062
00063 float fBm (point p; float filtwidth;
00064 uniform float octaves, lacunarity, gain, fourthdimension)
00065 {
00066 uniform float amp = 1;
00067 varying point pp = p;
00068 varying float sum = 0, fw = filtwidth;
00069 uniform float i;
00070
00071 for (i = 0; i < octaves; i += 1) {
00072 #pragma nolint
00073 sum += amp * filteredsnoiset (pp, fw, fourthdimension);
00074 amp *= gain; pp *= lacunarity; fw *= lacunarity;
00075 }
00076 return sum;
00077 }
00078
00079
00080
00081 #define fBm_default(p) fBm (p, filterwidthp(p), 4, 2, 0.5)
00082
00083
00084
00085
00086
00087
00088 vector
00089 vfBm (point p; float filtwidth;
00090 uniform float octaves, lacunarity, gain)
00091 {
00092 uniform float amp = 1;
00093 varying point pp = p;
00094 varying vector sum = 0;
00095 varying float fw = filtwidth;
00096 uniform float i;
00097
00098 for (i = 0; i < octaves; i += 1) {
00099 #pragma nolint
00100 sum += amp * filteredvsnoise (pp, fw);
00101 amp *= gain; pp *= lacunarity; fw *= lacunarity;
00102 }
00103 return sum;
00104 }
00105
00106
00107
00108 #define vfBm_default(p) vfBm (p, filterwidthp(p), 4, 2, 0.5)
00109
00110
00111
00112
00113 #define VLNoise(Pt,scale) (snoise(vsnoise(Pt)*scale+Pt))
00114 #define filteredVLNoise(Pt,fwidth,scale) \
00115 (filteredsnoise(filteredvsnoise(Pt,fwidth)*scale+Pt,fwidth))
00116
00117
00118 float
00119 VLfBm (point p; float filtwidth;
00120 uniform float octaves, lacunarity, gain, scale)
00121 {
00122 uniform float amp = 1;
00123 varying point pp = p;
00124 varying float sum = 0;
00125 varying float fw = filtwidth;
00126 uniform float i;
00127
00128 for (i = 0; i < octaves; i += 1) {
00129 #pragma nolint
00130 sum += amp * filteredVLNoise (pp, fw, scale);
00131 amp *= gain; pp *= lacunarity; fw *= lacunarity;
00132 }
00133 return sum;
00134 }
00135
00136
00137
00138 #define VLfBm_default(p) VLfBm (p, filterwidthp(p), 4, 2, 0.5, 1.0)
00139
00140
00141
00142
00143
00144
00145
00146 float turbulence (point p; float filtwidth;
00147 uniform float octaves, lacunarity, gain)
00148 {
00149 extern float du, dv;
00150 uniform float amp = 1;
00151 varying point pp = p;
00152 varying float sum = 0, fw = filtwidth;
00153 uniform float i;
00154
00155 for (i = 0; i < octaves; i += 1) {
00156 #pragma nolint
00157 float n = filteredsnoise (pp, fw);
00158 sum += amp * filteredabs (n, fw);
00159 amp *= gain; pp *= lacunarity; fw *= lacunarity;
00160 }
00161 return sum;
00162 }
00163
00164
00165
00166 #define turbulence_default(p) turbulence (p, filterwidthp(p), 4, 2, 0.5)
00167
00168
00169
00170
00171
00172
00173
00174
00175
00176
00177
00178
00179
00180
00181
00182
00183
00184
00185
00186
00187
00188
00189
00190
00191
00192
00193
00194 void
00195 voronoi_f1_3d (point P;
00196 float jitter;
00197 output float f1;
00198 output point pos1;
00199 )
00200 {
00201 point thiscell = point (floor(xcomp(P))+0.5, floor(ycomp(P))+0.5,
00202 floor(zcomp(P))+0.5);
00203 f1 = 1000;
00204 uniform float i, j, k;
00205 for (i = -1; i <= 1; i += 1) {
00206 for (j = -1; j <= 1; j += 1) {
00207 for (k = -1; k <= 1; k += 1) {
00208 point testcell = thiscell + vector(i,j,k);
00209 point pos = testcell +
00210 jitter * (vector cellnoise (testcell) - 0.5);
00211 vector offset = pos - P;
00212 float dist = offset . offset;
00213 if (dist < f1) {
00214 f1 = dist; pos1 = pos;
00215 }
00216 }
00217 }
00218 }
00219 f1 = sqrt(f1);
00220 }
00221
00222
00223
00224 void
00225 voronoi_f1f2_3d (point P;
00226 float jitter;
00227 output float f1; output point pos1;
00228 output float f2; output point pos2;
00229 )
00230 {
00231 point thiscell = point (floor(xcomp(P))+0.5, floor(ycomp(P))+0.5,
00232 floor(zcomp(P))+0.5);
00233 f1 = f2 = 1000;
00234 uniform float i, j, k;
00235 for (i = -1; i <= 1; i += 1) {
00236 for (j = -1; j <= 1; j += 1) {
00237 for (k = -1; k <= 1; k += 1) {
00238 point testcell = thiscell + vector(i,j,k);
00239 point pos = testcell +
00240 jitter * (vector cellnoise (testcell) - 0.5);
00241 vector offset = pos - P;
00242 float dist = offset . offset;
00243 if (dist < f1) {
00244 f2 = f1; pos2 = pos1;
00245 f1 = dist; pos1 = pos;
00246 } else if (dist < f2) {
00247 f2 = dist; pos2 = pos;
00248 }
00249 }
00250 }
00251 }
00252 f1 = sqrt(f1); f2 = sqrt(f2);
00253 }
00254
00255
00256
00257 void
00258 voronoi_f1_2d (float ss, tt;
00259 float jitter;
00260 output float f1;
00261 output float spos1, tpos1;
00262 )
00263 {
00264 float sthiscell = floor(ss)+0.5, tthiscell = floor(tt)+0.5;
00265 f1 = 1000;
00266 uniform float i, j;
00267 for (i = -1; i <= 1; i += 1) {
00268 float stestcell = sthiscell + i;
00269 for (j = -1; j <= 1; j += 1) {
00270 float ttestcell = tthiscell + j;
00271 float spos = stestcell +
00272 jitter * (float cellnoise(stestcell, ttestcell) - 0.5);
00273 float tpos = ttestcell +
00274 jitter * (float cellnoise(stestcell+23, ttestcell-87) - 0.5);
00275 float soffset = spos - ss;
00276 float toffset = tpos - tt;
00277 float dist = soffset*soffset + toffset*toffset;
00278 if (dist < f1) {
00279 f1 = dist;
00280 spos1 = spos; tpos1 = tpos;
00281 }
00282 }
00283 }
00284 f1 = sqrt(f1);
00285 }
00286
00287
00288
00289 void
00290 voronoi_f1f2_2d (float ss, tt;
00291 float jitter;
00292 output float f1;
00293 output float spos1, tpos1;
00294 output float f2;
00295 output float spos2, tpos2;
00296 )
00297 {
00298 float sthiscell = floor(ss)+0.5, tthiscell = floor(tt)+0.5;
00299 f1 = f2 = 1000;
00300 uniform float i, j;
00301 for (i = -1; i <= 1; i += 1) {
00302 float stestcell = sthiscell + i;
00303 for (j = -1; j <= 1; j += 1) {
00304 float ttestcell = tthiscell + j;
00305 float spos = stestcell +
00306 jitter * (cellnoise(stestcell, ttestcell) - 0.5);
00307 float tpos = ttestcell +
00308 jitter * (cellnoise(stestcell+23, ttestcell-87) - 0.5);
00309 float soffset = spos - ss;
00310 float toffset = tpos - tt;
00311 float dist = soffset*soffset + toffset*toffset;
00312 if (dist < f1) {
00313 f2 = f1; spos2 = spos1; tpos2 = tpos1;
00314 f1 = dist; spos1 = spos; tpos1 = tpos;
00315 } else if (dist < f2) {
00316 f2 = dist;
00317 spos2 = spos; tpos2 = tpos;
00318 }
00319 }
00320 }
00321 f1 = sqrt(f1); f2 = sqrt(f2);
00322 }
00323
00324
00325 #endif