23564583a4
This patch extends perlin noise to operate in 1D, 2D, 3D, and 4D space. The noise code has also been refactored to be more readable. The Color output and distortion patterns changed, so this patch breaks backward compatibility. This is due to the fact that we now use random offsets as noise seeds, as opposed to swizzling and constants offsets. Reviewers: brecht, JacquesLucke Differential Revision: https://developer.blender.org/D5560
132 lines
3.9 KiB
C
132 lines
3.9 KiB
C
/*
|
|
* Copyright 2011-2013 Blender Foundation
|
|
*
|
|
* Licensed under the Apache License, Version 2.0 (the "License");
|
|
* you may not use this file except in compliance with the License.
|
|
* You may obtain a copy of the License at
|
|
*
|
|
* http://www.apache.org/licenses/LICENSE-2.0
|
|
*
|
|
* Unless required by applicable law or agreed to in writing, software
|
|
* distributed under the License is distributed on an "AS IS" BASIS,
|
|
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
* See the License for the specific language governing permissions and
|
|
* limitations under the License.
|
|
*/
|
|
|
|
CCL_NAMESPACE_BEGIN
|
|
|
|
/* The fractal_noise_[1-4] functions are all exactly the same except for the input type. */
|
|
ccl_device_noinline float fractal_noise_1d(float p, float octaves)
|
|
{
|
|
float fscale = 1.0f;
|
|
float amp = 1.0f;
|
|
float sum = 0.0f;
|
|
octaves = clamp(octaves, 0.0f, 16.0f);
|
|
int n = float_to_int(octaves);
|
|
for (int i = 0; i <= n; i++) {
|
|
float t = noise_1d(fscale * p);
|
|
sum += t * amp;
|
|
amp *= 0.5f;
|
|
fscale *= 2.0f;
|
|
}
|
|
float rmd = octaves - floorf(octaves);
|
|
if (rmd != 0.0f) {
|
|
float t = noise_1d(fscale * p);
|
|
float sum2 = sum + t * amp;
|
|
sum *= ((float)(1 << n) / (float)((1 << (n + 1)) - 1));
|
|
sum2 *= ((float)(1 << (n + 1)) / (float)((1 << (n + 2)) - 1));
|
|
return (1.0f - rmd) * sum + rmd * sum2;
|
|
}
|
|
else {
|
|
sum *= ((float)(1 << n) / (float)((1 << (n + 1)) - 1));
|
|
return sum;
|
|
}
|
|
}
|
|
|
|
/* The fractal_noise_[1-4] functions are all exactly the same except for the input type. */
|
|
ccl_device_noinline float fractal_noise_2d(float2 p, float octaves)
|
|
{
|
|
float fscale = 1.0f;
|
|
float amp = 1.0f;
|
|
float sum = 0.0f;
|
|
octaves = clamp(octaves, 0.0f, 16.0f);
|
|
int n = float_to_int(octaves);
|
|
for (int i = 0; i <= n; i++) {
|
|
float t = noise_2d(fscale * p);
|
|
sum += t * amp;
|
|
amp *= 0.5f;
|
|
fscale *= 2.0f;
|
|
}
|
|
float rmd = octaves - floorf(octaves);
|
|
if (rmd != 0.0f) {
|
|
float t = noise_2d(fscale * p);
|
|
float sum2 = sum + t * amp;
|
|
sum *= ((float)(1 << n) / (float)((1 << (n + 1)) - 1));
|
|
sum2 *= ((float)(1 << (n + 1)) / (float)((1 << (n + 2)) - 1));
|
|
return (1.0f - rmd) * sum + rmd * sum2;
|
|
}
|
|
else {
|
|
sum *= ((float)(1 << n) / (float)((1 << (n + 1)) - 1));
|
|
return sum;
|
|
}
|
|
}
|
|
|
|
/* The fractal_noise_[1-4] functions are all exactly the same except for the input type. */
|
|
ccl_device_noinline float fractal_noise_3d(float3 p, float octaves)
|
|
{
|
|
float fscale = 1.0f;
|
|
float amp = 1.0f;
|
|
float sum = 0.0f;
|
|
octaves = clamp(octaves, 0.0f, 16.0f);
|
|
int n = float_to_int(octaves);
|
|
for (int i = 0; i <= n; i++) {
|
|
float t = noise_3d(fscale * p);
|
|
sum += t * amp;
|
|
amp *= 0.5f;
|
|
fscale *= 2.0f;
|
|
}
|
|
float rmd = octaves - floorf(octaves);
|
|
if (rmd != 0.0f) {
|
|
float t = noise_3d(fscale * p);
|
|
float sum2 = sum + t * amp;
|
|
sum *= ((float)(1 << n) / (float)((1 << (n + 1)) - 1));
|
|
sum2 *= ((float)(1 << (n + 1)) / (float)((1 << (n + 2)) - 1));
|
|
return (1.0f - rmd) * sum + rmd * sum2;
|
|
}
|
|
else {
|
|
sum *= ((float)(1 << n) / (float)((1 << (n + 1)) - 1));
|
|
return sum;
|
|
}
|
|
}
|
|
|
|
/* The fractal_noise_[1-4] functions are all exactly the same except for the input type. */
|
|
ccl_device_noinline float fractal_noise_4d(float4 p, float octaves)
|
|
{
|
|
float fscale = 1.0f;
|
|
float amp = 1.0f;
|
|
float sum = 0.0f;
|
|
octaves = clamp(octaves, 0.0f, 16.0f);
|
|
int n = float_to_int(octaves);
|
|
for (int i = 0; i <= n; i++) {
|
|
float t = noise_4d(fscale * p);
|
|
sum += t * amp;
|
|
amp *= 0.5f;
|
|
fscale *= 2.0f;
|
|
}
|
|
float rmd = octaves - floorf(octaves);
|
|
if (rmd != 0.0f) {
|
|
float t = noise_4d(fscale * p);
|
|
float sum2 = sum + t * amp;
|
|
sum *= ((float)(1 << n) / (float)((1 << (n + 1)) - 1));
|
|
sum2 *= ((float)(1 << (n + 1)) / (float)((1 << (n + 2)) - 1));
|
|
return (1.0f - rmd) * sum + rmd * sum2;
|
|
}
|
|
else {
|
|
sum *= ((float)(1 << n) / (float)((1 << (n + 1)) - 1));
|
|
return sum;
|
|
}
|
|
}
|
|
|
|
CCL_NAMESPACE_END
|