nuclear@10: /*
nuclear@10: libvmath - a vector math library
nuclear@10: Copyright (C) 2004-2015 John Tsiombikas <nuclear@member.fsf.org>
nuclear@10: 
nuclear@10: This program is free software: you can redistribute it and/or modify
nuclear@10: it under the terms of the GNU Lesser General Public License as published
nuclear@10: by the Free Software Foundation, either version 3 of the License, or
nuclear@10: (at your option) any later version.
nuclear@10: 
nuclear@10: This program is distributed in the hope that it will be useful,
nuclear@10: but WITHOUT ANY WARRANTY; without even the implied warranty of
nuclear@10: MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
nuclear@10: GNU Lesser General Public License for more details.
nuclear@10: 
nuclear@10: You should have received a copy of the GNU Lesser General Public License
nuclear@10: along with this program.  If not, see <http://www.gnu.org/licenses/>.
nuclear@10: */
nuclear@10: 
nuclear@10: #ifndef VMATH_H_
nuclear@10: #define VMATH_H_
nuclear@10: 
nuclear@10: #include <math.h>
nuclear@10: #include "vmath_types.h"
nuclear@10: 
nuclear@10: #ifndef M_PI
nuclear@10: #define M_PI	PI
nuclear@10: #endif
nuclear@10: 
nuclear@10: #ifndef M_E
nuclear@10: #define M_E				2.718281828459045
nuclear@10: #endif
nuclear@10: 
nuclear@10: #define PI				3.141592653589793
nuclear@10: #define HALF_PI			1.570796326794897
nuclear@10: #define QUARTER_PI		0.785398163397448
nuclear@10: #define TWO_PI			6.283185307179586
nuclear@10: 
nuclear@10: 
nuclear@10: #define RAD_TO_DEG(a) ((((scalar_t)a) * 360.0) / TWO_PI)
nuclear@10: #define DEG_TO_RAD(a) (((scalar_t)a) * (PI / 180.0))
nuclear@10: 
nuclear@10: #define SQ(x) ((x) * (x))
nuclear@10: 
nuclear@10: #define MIN(a, b)	((a) < (b) ? (a) : (b))
nuclear@10: #define MAX(a, b)	((a) > (b) ? (a) : (b))
nuclear@10: 
nuclear@10: #if (__STDC_VERSION__ < 199901L) && (__cplusplus < 201103L) && !defined(__GNUC__) && (_MSC_VER < 1800)
nuclear@10: #define round(x)	((x) >= 0 ? (x) + 0.5 : (x) - 0.5)
nuclear@10: #endif
nuclear@10: 
nuclear@10: #ifdef __cplusplus
nuclear@10: extern "C" {
nuclear@10: #endif	/* __cplusplus */
nuclear@10: 
nuclear@10: void enable_fpexcept(void);
nuclear@10: void disable_fpexcept(void);
nuclear@10: 
nuclear@10: static inline scalar_t smoothstep(float a, float b, float x);
nuclear@10: 
nuclear@10: static inline scalar_t frand(scalar_t range);
nuclear@10: static inline vec3_t sphrand(scalar_t rad);
nuclear@10: 
nuclear@10: scalar_t integral(scalar_t (*f)(scalar_t), scalar_t low, scalar_t high, int samples);
nuclear@10: scalar_t gaussian(scalar_t x, scalar_t mean, scalar_t sdev);
nuclear@10: 
nuclear@10: static inline scalar_t lerp(scalar_t a, scalar_t b, scalar_t t);
nuclear@10: 
nuclear@10: scalar_t bspline(scalar_t a, scalar_t b, scalar_t c, scalar_t d, scalar_t t);
nuclear@10: scalar_t spline(scalar_t a, scalar_t b, scalar_t c, scalar_t d, scalar_t t);
nuclear@10: scalar_t bezier(scalar_t a, scalar_t b, scalar_t c, scalar_t d, scalar_t t);
nuclear@10: 
nuclear@10: scalar_t noise1(scalar_t x);
nuclear@10: scalar_t noise2(scalar_t x, scalar_t y);
nuclear@10: scalar_t noise3(scalar_t x, scalar_t y, scalar_t z);
nuclear@10: 
nuclear@10: scalar_t fbm1(scalar_t x, int octaves);
nuclear@10: scalar_t fbm2(scalar_t x, scalar_t y, int octaves);
nuclear@10: scalar_t fbm3(scalar_t x, scalar_t y, scalar_t z, int octaves);
nuclear@10: 
nuclear@10: scalar_t turbulence1(scalar_t x, int octaves);
nuclear@10: scalar_t turbulence2(scalar_t x, scalar_t y, int octaves);
nuclear@10: scalar_t turbulence3(scalar_t x, scalar_t y, scalar_t z, int octaves);
nuclear@10: 
nuclear@10: #ifdef __cplusplus
nuclear@10: }
nuclear@10: #endif	/* __cplusplus */
nuclear@10: 
nuclear@10: #include "vmath.inl"
nuclear@10: 
nuclear@10: #include "vector.h"
nuclear@10: #include "matrix.h"
nuclear@10: #include "quat.h"
nuclear@10: #include "sphvec.h"
nuclear@10: #include "ray.h"
nuclear@10: #include "geom.h"
nuclear@10: 
nuclear@10: #endif	/* VMATH_H_ */