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