ReactOS 0.4.16-dev-125-g798ea90
mmath.c
Go to the documentation of this file.
1/* $Id: mmath.c,v 1.3 1997/07/24 01:23:16 brianp Exp $ */
2
3/*
4 * Mesa 3-D graphics library
5 * Version: 2.4
6 * Copyright (C) 1995-1997 Brian Paul
7 *
8 * This library is free software; you can redistribute it and/or
9 * modify it under the terms of the GNU Library General Public
10 * License as published by the Free Software Foundation; either
11 * version 2 of the License, or (at your option) any later version.
12 *
13 * This library is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 * Library General Public License for more details.
17 *
18 * You should have received a copy of the GNU Library General Public
19 * License along with this library; if not, write to the Free
20 * Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
21 */
22
23
24/*
25 * $Log: mmath.c,v $
26 * Revision 1.3 1997/07/24 01:23:16 brianp
27 * changed precompiled header symbol from PCH to PC_HEADER
28 *
29 * Revision 1.2 1997/05/28 03:25:43 brianp
30 * added precompiled header (PCH) support
31 *
32 * Revision 1.1 1997/05/01 01:41:54 brianp
33 * Initial revision
34 *
35 */
36
37
38#ifdef PC_HEADER
39#include "all.h"
40#else
41#include "GL/gl.h"
42#include "mmath.h"
43#endif
44
45
46
47/*
48 * A High Speed, Low Precision Square Root
49 * by Paul Lalonde and Robert Dawson
50 * from "Graphics Gems", Academic Press, 1990
51 */
52
53/*
54 * SPARC implementation of a fast square root by table
55 * lookup.
56 * SPARC floating point format is as follows:
57 *
58 * BIT 31 30 23 22 0
59 * sign exponent mantissa
60 */
61static short sqrttab[0x100]; /* declare table of square roots */
62
63static void init_sqrt(void)
64{
65#ifdef FAST_MATH
66 unsigned short i;
67 float f;
68 unsigned int *fi = (unsigned int *)&f;
69 /* to access the bits of a float in */
70 /* C quickly we must misuse pointers */
71
72 for(i=0; i<= 0x7f; i++) {
73 *fi = 0;
74
75 /*
76 * Build a float with the bit pattern i as mantissa
77 * and an exponent of 0, stored as 127
78 */
79
80 *fi = (i << 16) | (127 << 23);
81 f = sqrt(f);
82
83 /*
84 * Take the square root then strip the first 7 bits of
85 * the mantissa into the table
86 */
87
88 sqrttab[i] = (*fi & 0x7fffff) >> 16;
89
90 /*
91 * Repeat the process, this time with an exponent of
92 * 1, stored as 128
93 */
94
95 *fi = 0;
96 *fi = (i << 16) | (128 << 23);
97 f = sqrt(f);
98 sqrttab[i+0x80] = (*fi & 0x7fffff) >> 16;
99 }
100#endif /*FAST_MATH*/
101}
102
103
104float gl_sqrt( float x )
105{
106#ifdef FAST_MATH
107 unsigned int *num = (unsigned int *)&x;
108 /* to access the bits of a float in C
109 * we must misuse pointers */
110
111 short e; /* the exponent */
112 if (x == 0.0F) return 0.0F; /* check for square root of 0 */
113 e = (*num >> 23) - 127; /* get the exponent - on a SPARC the */
114 /* exponent is stored with 127 added */
115 *num &= 0x7fffff; /* leave only the mantissa */
116 if (e & 0x01) *num |= 0x800000;
117 /* the exponent is odd so we have to */
118 /* look it up in the second half of */
119 /* the lookup table, so we set the */
120 /* high bit */
121 e >>= 1; /* divide the exponent by two */
122 /* note that in C the shift */
123 /* operators are sign preserving */
124 /* for signed operands */
125 /* Do the table lookup, based on the quaternary mantissa,
126 * then reconstruct the result back into a float
127 */
128 *num = ((sqrttab[*num >> 16]) << 16) | ((e + 127) << 23);
129 return x;
130#else
131 return sqrt(x);
132#endif
133}
134
135
136
137/*
138 * Initialize tables, etc for fast math functions.
139 */
140void gl_init_math(void)
141{
143
144 if (!initialized) {
145 init_sqrt();
146
147
149 }
150}
_STLP_DECLSPEC complex< float > _STLP_CALL sqrt(const complex< float > &)
Definition: complex.cpp:188
#define GL_TRUE
Definition: gl.h:174
GLint GLint GLint GLint GLint x
Definition: gl.h:1548
#define GL_FALSE
Definition: gl.h:173
unsigned char GLboolean
Definition: gl.h:151
GLfloat f
Definition: glext.h:7540
GLuint GLuint num
Definition: glext.h:9618
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
Definition: glfuncs.h:248
#define e
Definition: ke_i.h:82
#define f
Definition: ke_i.h:83
float gl_sqrt(float x)
Definition: mmath.c:104
void gl_init_math(void)
Definition: mmath.c:140
static short sqrttab[0x100]
Definition: mmath.c:61
static void init_sqrt(void)
Definition: mmath.c:63
static BOOL initialized
Definition: syslog.c:39