Files
Raytracer/perlin.h
2024-07-09 17:31:25 -04:00

98 lines
2.6 KiB
C++

#ifndef PERLIN_H
#define PERLIN_H
#include "rayTracer.h"
class perlin {
public:
perlin() {
randFloat = new double[pointCount];
for (int i = 0; i < pointCount; i++) {
randFloat[i] = randomDouble();
}
permX = perlinGeneratePerm();
permY = perlinGeneratePerm();
permZ = perlinGeneratePerm();
}
~perlin() {
delete[] randFloat;
delete[] permX;
delete[] permY;
delete[] permZ;
}
double noise(const point3& p) const {
auto u = p.x() - floor(p.x());
auto v = p.y() - floor(p.y());
auto w = p.z() - floor(p.z());
u = u*u*(3-2*u);
v = v*v*(3-2*v);
w = w*w*(3-2*w);
auto i = int(floor(p.x()));
auto j = int(floor(p.y()));
auto k = int(floor(p.z()));
double c[2][2][2];
for (int di = 0; di < 2; di++) {
for (int dj = 0; dj < 2; dj++) {
for (int dk = 0; dk < 2; dk++) {
c[di][dj][dk] = randFloat[
permX[(i+di) & 255] ^
permY[(j+dj) & 255] ^
permZ[(k+dk) & 255]
];
}
}
}
return trileanerInterpolation(c, u, v, w);
}
private:
static const int pointCount = 256;
double* randFloat;
int* permX;
int* permY;
int* permZ;
static int* perlinGeneratePerm(void) {
auto p = new int[pointCount];
for (int i = 0; i < pointCount; i++) {
p[i] = i;
}
permute(p, pointCount);
return p;
}
static void permute(int* p, int n) {
for (int i = n - 1; i > 0; i--) {
int target = randomInt(0, i);
int tmp = p[i];
p[i] = p[target];
p[target] = tmp;
}
}
static double trileanerInterpolation(double c[2][2][2], double u, double v, double w) {
auto accum = 0.0;
for (int i =0; i < 2; i++) {
for (int j = 0; j < 2; j++) {
for (int k = 0; k < 2; k++) {
accum += (i*u + (1-i)*(1-u))
* (j*v + (1-j)*(1-v))
* (k*w + (1-k)*(1-w))
* c[i][j][k];
}
}
}
return accum;
}
};
#endif