Alternative to ci::Perlin class.
This code is originally intended for use with the Cinder C++ library but only depends on GLM.
Simplex noise is better looking, faster and has proper derivatives without the grid artifacts of Perlin noise derivatives. This code is mostly ported from Stefan Gustavson public domain implementation for the simplex noise, simplex flow noise and their analytical derivatives. The curl noise functions are adapted from Robert Bridson papers to use derivatives instead of finite differences. There's also a series of noise sums and other multi-fractals noises...
//! Returns a 1D simplex noisefloatnoise( float x ); //! Returns a 2D simplex noisefloatnoise( const glm::vec2 &v ); //! Returns a 3D simplex noisefloatnoise( const glm::vec3 &v ); //! Returns a 4D simplex noisefloatnoise( const glm::vec4 &v ); //! Returns a 1D simplex ridged noisefloatridgedNoise( float x ); //! Returns a 2D simplex ridged noisefloatridgedNoise( const glm::vec2 &v ); //! Returns a 3D simplex ridged noisefloatridgedNoise( const glm::vec3 &v ); //! Returns a 4D simplex ridged noisefloatridgedNoise( const glm::vec4 &v ); //! Returns a 1D simplex noise with analytical derivative. glm::vec2 dnoise( float x ); //! Returns a 2D simplex noise with analytical derivatives. glm::vec3 dnoise( const glm::vec2 &v ); //! Returns a 3D simplex noise with analytical derivatives. glm::vec4 dnoise( const glm::vec3 &v ); // not optimal but easiest way to return 5 floatstypedef std::array<float,5> vec5; //! Returns a 4D simplex noise with analytical derivatives vec5 dnoise( const glm::vec4 &v ); //! Returns a 2D simplex cellular/worley noisefloatworleyNoise( const glm::vec2 &v ); //! Returns a 3D simplex cellular/worley noisefloatworleyNoise( const glm::vec3 &v ); //! Returns a 2D simplex smooth cellular/worley noisefloatworleyNoise( const glm::vec2 &v, float falloff ); //! Returns a 3D simplex smooth cellular/worley noisefloatworleyNoise( const glm::vec3 &v, float falloff ); //! Returns a 2D simplex noise with rotating gradientsfloatflowNoise( const glm::vec2 &v, float angle ); //! Returns a 3D simplex noise with rotating gradientsfloatflowNoise( const glm::vec3 &v, float angle ); //! Returns a 2D simplex noise with rotating gradients and analytical derivatives glm::vec3 dFlowNoise( const glm::vec2 &v, float angle ); //! Returns a 3D simplex noise with rotating gradients and analytical derivatives glm::vec4 dFlowNoise( const glm::vec3 &v, float angle ); //! Returns the curl of a 2D simplex noise glm::vec2 curlNoise( const glm::vec2 &v ); //! Returns the curl of a 2D simplex flow noise glm::vec2 curlNoise( const glm::vec2 &v, float t ); //! Returns the curl of a 2D simplex noise fractal brownian motion sum glm::vec2 curlNoise( const glm::vec2 &v, uint8_t octaves, float lacunarity, float gain ); //! Returns the curl of a 3D simplex noise glm::vec3 curlNoise( const glm::vec3 &v ); //! Returns the curl of a 3D simplex flow noise glm::vec3 curlNoise( const glm::vec3 &v, float t ); //! Returns the curl approximation of a 3D simplex noise fractal brownian motion sum glm::vec3 curlNoise( const glm::vec3 &v, uint8_t octaves, float lacunarity, float gain ); //! Returns the curl of a custom 2D potential using finite difference approximation glm::vec2 curl( const glm::vec2 &v, const std::function<float(const glm::vec2&)> &potential, float delta = 1e-4f ); //! Returns the curl of a custom 3D potential using finite difference approximation glm::vec3 curl( const glm::vec3 &v, const std::function<glm::vec3(const glm::vec3&)> &potential, float delta = 1e-4f ); //! Returns a 1D simplex noise fractal brownian motion sumfloatfBm( float x, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns a 2D simplex noise fractal brownian motion sumfloatfBm( const glm::vec2 &v, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns a 3D simplex noise fractal brownian motion sumfloatfBm( const glm::vec3 &v, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns a 4D simplex noise fractal brownian motion sumfloatfBm( const glm::vec4 &v, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns a 1D simplex noise fractal brownian motion sum with analytical derivatives glm::vec2 dfBm( float x, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns a 2D simplex noise fractal brownian motion sum with analytical derivatives glm::vec3 dfBm( const glm::vec2 &v, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns a 3D simplex noise fractal brownian motion sum with analytical derivatives glm::vec4 dfBm( const glm::vec3 &v, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns a 4D simplex noise fractal brownian motion sum with analytical derivatives vec5 dfBm( const glm::vec4 &v, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns a 1D simplex ridged multi-fractal noise sumfloatridgedMF( float x, float ridgeOffset = 1.0f, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns a 2D simplex ridged multi-fractal noise sumfloatridgedMF( const glm::vec2 &v, float ridgeOffset = 1.0f, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns a 3D simplex ridged multi-fractal noise sumfloatridgedMF( const glm::vec3 &v, float ridgeOffset = 1.0f, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns a 4D simplex ridged multi-fractal noise sumfloatridgedMF( const glm::vec4 &v, float ridgeOffset = 1.0f, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns the 2D simplex noise fractal brownian motion sum variation by Iñigo QuilezfloatiqfBm( const glm::vec2 &v, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns the 2D simplex noise fractal brownian motion sum variation by Iñigo QuilezfloatiqfBm( const glm::vec3 &v, uint8_t octaves = 4, float lacunarity = 2.0f, float gain = 0.5f ); //! Returns the 2D simplex noise fractal brownian motion sum variation by Iñigo Quilez that use a mat2 to transform each octavefloatiqMatfBm( const glm::vec2 &v, uint8_t octaves = 4, const glm::mat2 &mat = glm::mat2( 1.6, -1.2, 1.2, 1.6 ), float gain = 0.5f ); //! Seeds the permutation table with new random valuesvoidseed( uint32_t s );