48 #define TOLERANCE 0.0001 55 #define MAX_DIMENSION 4 61 static void bezierCurveEvalfast(
float u0,
float u1,
int order,
float *ctlpoints,
int stride,
int dimension,
float u,
float retpoint[]);
77 float uprime = (
u-u0)/(
u1-u0);
78 float *ctlptr = ctlpoints;
79 float oneMinusX = 1.0f-uprime;
83 for(
k=0;
k<dimension;
k++)
84 retpoint[
k] = (*(ctlptr +
k));
89 for(
k=0;
k<dimension;
k++) {
99 void bezierCurveEvalfast(
float u0,
float u1,
int order,
float *ctlpoints,
int stride,
int dimension,
float u,
float retpoint[])
101 float uprime = (
u-u0)/(
u1-u0);
103 float* ctlptr = ctlpoints;
106 for(
j=0;
j<dimension;
j++)
112 for(
j=0;
j<dimension;
j++)
117 for(
j=0;
j<dimension;
j++)
129 float *ctlptr = ctlpoints;
133 for(
k=0;
k<dimension;
k++)
137 for(
k=0;
k<dimension;
k++) {
149 float *ctlptr = ctlpoints;
154 for(
k=0;
k<dimension;
k++){
161 for(
r=1;
r<=der;
r++){
163 for(
k=0;
k<dimension;
k++){
176 void bezierSurfEvalDerGen(
int uder,
int vder,
float u0,
float u1,
int uorder,
float v0,
float v1,
int vorder,
int dimension,
float *ctlpoints,
int ustride,
int vstride,
float u,
float v,
float ret[])
192 void bezierSurfEval(
float u0,
float u1,
int uorder,
float v0,
float v1,
int vorder,
int dimension,
float *ctlpoints,
int ustride,
int vstride,
float u,
float v,
float ret[])
194 bezierSurfEvalDerGen(0, 0, u0,
u1,
uorder,
v0,
v1,
vorder, dimension, ctlpoints,
ustride,
vstride,
u,
v,
ret);
202 void bezierSurfEvalNormal(
float u0,
float u1,
int uorder,
float v0,
float v1,
int vorder,
int dimension,
float *ctlpoints,
int ustride,
int vstride,
float u,
float v,
float retNormal[])
206 assert(dimension>=3 && dimension <=4);
207 bezierSurfEvalDerGen(1,0, u0,
u1,
uorder,
v0,
v1,
vorder, dimension, ctlpoints,
ustride,
vstride,
u,
v, partialU);
208 bezierSurfEvalDerGen(0,1, u0,
u1,
uorder,
v0,
v1,
vorder, dimension, ctlpoints,
ustride,
vstride,
u,
v, partialV);
222 bezierSurfEvalDerGen(0,0, u0,
u1,
uorder,
v0,
v1,
vorder, dimension, ctlpoints,
ustride,
vstride,
u,
v,
val);
225 newPartialU[
i] = partialU[
i] *
val[3] -
val[
i] * partialU[3];
226 newPartialV[
i] = partialV[
i] *
val[3] -
val[
i] * partialV[3];
255 ret[0] =
x[1]*
y[2] -
y[1]*
x[2];
256 ret[1] =
x[2]*
y[0] -
y[2]*
x[0];
257 ret[2] =
x[0]*
y[1] -
y[0]*
x[1];
GLsizei GLenum const GLvoid GLsizei GLenum GLbyte GLbyte GLbyte GLdouble GLdouble GLdouble GLfloat GLfloat GLfloat GLint GLint GLint GLshort GLshort GLshort GLubyte GLubyte GLubyte GLuint GLuint GLuint GLushort GLushort GLushort GLbyte GLbyte GLbyte GLbyte GLdouble GLdouble GLdouble GLdouble GLfloat GLfloat GLfloat GLfloat GLint GLint GLint GLint GLshort GLshort GLshort GLshort GLubyte GLubyte GLubyte GLubyte GLuint GLuint GLuint GLuint GLushort GLushort GLushort GLushort GLboolean const GLdouble const GLfloat const GLint const GLshort const GLbyte const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLdouble const GLfloat const GLfloat const GLint const GLint const GLshort const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort GLenum GLenum GLenum GLfloat GLenum GLint GLenum GLenum GLenum GLfloat GLenum GLenum GLint GLenum GLfloat GLenum GLint GLint GLushort GLenum GLenum GLfloat GLenum GLenum GLint GLfloat const GLubyte GLenum GLenum GLenum const GLfloat GLenum GLenum const GLint GLenum GLint GLint GLsizei GLsizei GLint GLenum GLenum const GLvoid GLenum GLenum const GLfloat GLenum GLenum const GLint GLenum GLenum const GLdouble GLenum GLenum const GLfloat GLenum GLenum const GLint GLsizei GLuint GLfloat GLuint GLbitfield GLfloat GLint GLuint GLboolean GLenum GLfloat GLenum GLbitfield GLenum GLfloat GLfloat GLint GLint const GLfloat GLenum GLfloat GLfloat GLint GLint GLfloat GLfloat GLint GLint const GLfloat GLint GLfloat GLfloat GLint GLfloat GLfloat GLint GLfloat GLfloat const GLdouble * u
_STLP_DECLSPEC complex< float > _STLP_CALL sqrt(const complex< float > &)
GLdouble GLdouble GLint GLint GLdouble GLdouble GLint vstride
GLint GLint GLsizei width
GLdouble GLdouble GLint ustride
GLuint GLdouble GLdouble GLint GLint order
void bezierSurfEvalDerGen(int uder, int vder, float u0, float u1, int uorder, float v0, float v1, int vorder, int dimension, float *ctlpoints, int ustride, int vstride, float u, float v, float ret[])
GLdouble GLdouble GLdouble r
static void normalize(float vec[3])
GLint GLint GLint GLint GLint x
void bezierCurveEval(float u0, float u1, int order, float *ctlpoints, int stride, int dimension, float u, float retpoint[])
GLdouble GLdouble GLint GLint GLdouble GLdouble GLint GLint vorder
GLenum GLuint GLenum GLsizei const GLchar * buf
_Check_return_opt_ _CRTIMP int __cdecl fprintf(_Inout_ FILE *_File, _In_z_ _Printf_format_string_ const char *_Format,...)
static float binomialCoefficients[8][8]
static void crossProduct(float x[3], float y[3], float ret[3])
void bezierSurfEvalNormal(float u0, float u1, int uorder, float v0, float v1, int vorder, int dimension, float *ctlpoints, int ustride, int vstride, float u, float v, float retNormal[])
GLsizei GLenum const GLvoid GLsizei GLenum GLbyte GLbyte GLbyte GLdouble GLdouble GLdouble GLfloat GLfloat GLfloat GLint GLint GLint GLshort GLshort GLshort GLubyte GLubyte GLubyte GLuint GLuint GLuint GLushort GLushort GLushort GLbyte GLbyte GLbyte GLbyte GLdouble GLdouble GLdouble GLdouble GLfloat GLfloat GLfloat GLfloat GLint GLint GLint GLint GLshort GLshort GLshort GLshort GLubyte GLubyte GLubyte GLubyte GLuint GLuint GLuint GLuint GLushort GLushort GLushort GLushort GLboolean const GLdouble const GLfloat const GLint const GLshort const GLbyte const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLdouble const GLfloat const GLfloat const GLint const GLint const GLshort const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort GLenum GLenum GLenum GLfloat GLenum GLint GLenum GLenum GLenum GLfloat GLenum GLenum GLint GLenum GLfloat GLenum GLint GLint GLushort GLenum GLenum GLfloat GLenum GLenum GLint GLfloat const GLubyte GLenum GLenum GLenum const GLfloat GLenum GLenum const GLint GLenum GLint GLint GLsizei GLsizei GLint GLenum GLenum const GLvoid GLenum GLenum const GLfloat GLenum GLenum const GLint GLenum GLenum const GLdouble GLenum GLenum const GLfloat GLenum GLenum const GLint GLsizei GLuint GLfloat GLuint GLbitfield GLfloat GLint GLuint GLboolean GLenum GLfloat GLenum GLbitfield GLenum GLfloat GLfloat GLint GLint const GLfloat GLenum GLfloat GLfloat GLint GLint GLfloat GLfloat GLint GLint const GLfloat GLint GLfloat GLfloat GLint GLfloat GLfloat GLint GLfloat GLfloat const GLdouble const GLfloat const GLdouble const GLfloat GLint GLint GLint j
GLdouble GLdouble GLint GLint uorder
void bezierSurfEval(float u0, float u1, int uorder, float v0, float v1, int vorder, int dimension, float *ctlpoints, int ustride, int vstride, float u, float v, float ret[])
GLsizei GLenum const GLvoid GLsizei GLenum GLbyte GLbyte GLbyte GLdouble GLdouble GLdouble GLfloat GLfloat GLfloat GLint GLint GLint GLshort GLshort GLshort GLubyte GLubyte GLubyte GLuint GLuint GLuint GLushort GLushort GLushort GLbyte GLbyte GLbyte GLbyte GLdouble GLdouble GLdouble GLdouble GLfloat GLfloat GLfloat GLfloat GLint GLint GLint GLint GLshort GLshort GLshort GLshort GLubyte GLubyte GLubyte GLubyte GLuint GLuint GLuint GLuint GLushort GLushort GLushort GLushort GLboolean const GLdouble const GLfloat const GLint const GLshort const GLbyte const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLdouble const GLfloat const GLfloat const GLint const GLint const GLshort const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort GLenum GLenum GLenum GLfloat GLenum GLint GLenum GLenum GLenum GLfloat GLenum GLenum GLint GLenum GLfloat GLenum GLint GLint GLushort GLenum GLenum GLfloat GLenum GLenum GLint GLfloat const GLubyte GLenum GLenum GLenum const GLfloat GLenum GLenum const GLint GLenum GLint GLint GLsizei GLsizei GLint GLenum GLenum const GLvoid GLenum GLenum const GLfloat GLenum GLenum const GLint GLenum GLenum const GLdouble GLenum GLenum const GLfloat GLenum GLenum const GLint GLsizei GLuint GLfloat GLuint GLbitfield GLfloat GLint GLuint GLboolean GLenum GLfloat GLenum GLbitfield GLenum GLfloat GLfloat GLint GLint const GLfloat GLenum GLfloat GLfloat GLint GLint GLfloat GLfloat GLint GLint const GLfloat GLint GLfloat GLfloat GLint GLfloat GLfloat GLint GLfloat GLfloat const GLdouble const GLfloat const GLdouble const GLfloat GLint i
void bezierCurveEvalDer(float u0, float u1, int order, float *ctlpoints, int stride, int dimension, float u, float retDer[])
static float(__cdecl *square_half_float)(float x
GLint GLint GLint GLint GLint GLint y
void bezierCurveEvalDerGen(int der, float u0, float u1, int order, float *ctlpoints, int stride, int dimension, float u, float retDer[])