2 * Mesa 3-D graphics library
5 * Copyright (C) 1999-2003 Brian Paul All Rights Reserved.
7 * Permission is hereby granted, free of charge, to any person obtaining a
8 * copy of this software and associated documentation files (the "Software"),
9 * to deal in the Software without restriction, including without limitation
10 * the rights to use, copy, modify, merge, publish, distribute, sublicense,
11 * and/or sell copies of the Software, and to permit persons to whom the
12 * Software is furnished to do so, subject to the following conditions:
14 * The above copyright notice and this permission notice shall be included
15 * in all copies or substantial portions of the Software.
17 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
18 * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
19 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
20 * BRIAN PAUL BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN
21 * AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN
22 * CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
27 * Matrix/vertex/vector transformation stuff
31 * 1. 4x4 transformation matrices are stored in memory in column major order.
32 * 2. Points/vertices are to be thought of as column vectors.
33 * 3. Transformation of a point p by a matrix M is: p' = M * p
36 #include "main/glheader.h"
37 #include "main/macros.h"
41 #include "m_translate.h"
51 #include "x86/common_x86_asm.h"
55 #include "x86-64/x86-64.h"
59 #include "sparc/sparc.h"
63 #include "ppc/common_ppc_features.h"
66 clip_func _mesa_clip_tab
[5];
67 clip_func _mesa_clip_np_tab
[5];
68 dotprod_func _mesa_dotprod_tab
[5];
69 vec_copy_func _mesa_copy_tab
[0x10];
70 normal_func _mesa_normal_tab
[0xf];
71 transform_func
*_mesa_transform_tab
[5];
74 /* Raw data format used for:
75 * - Object-to-eye transform prior to culling, although this too
76 * could be culled under some circumstances.
77 * - Eye-to-clip transform (via the function above).
79 * - And everything else too, if culling happens to be disabled.
81 * GH: It's used for everything now, as clipping/culling is done
82 * elsewhere (most often by the driver itself).
85 #define TAG2(x,y) x##y
86 #define STRIDE_LOOP for ( i = 0 ; i < count ; i++, STRIDE_F(from, stride) )
87 #define LOOP for ( i = 0 ; i < n ; i++ )
89 #include "m_xform_tmp.h"
90 #include "m_clip_tmp.h"
91 #include "m_norm_tmp.h"
92 #include "m_dotprod_tmp.h"
93 #include "m_copy_tmp.h"
102 GLvector4f
*_mesa_project_points( GLvector4f
*proj_vec
,
103 const GLvector4f
*clip_vec
)
105 const GLuint stride
= clip_vec
->stride
;
106 const GLfloat
*from
= (GLfloat
*)clip_vec
->start
;
107 const GLuint count
= clip_vec
->count
;
108 GLfloat (*vProj
)[4] = (GLfloat (*)[4])proj_vec
->start
;
111 for (i
= 0 ; i
< count
; i
++, STRIDE_F(from
, stride
))
113 GLfloat oow
= 1.0F
/ from
[3];
115 vProj
[i
][0] = from
[0] * oow
;
116 vProj
[i
][1] = from
[1] * oow
;
117 vProj
[i
][2] = from
[2] * oow
;
120 proj_vec
->flags
|= VEC_SIZE_4
;
122 proj_vec
->count
= clip_vec
->count
;
132 * Transform a 4-element row vector (1x4 matrix) by a 4x4 matrix. This
133 * function is used for transforming clipping plane equations and spotlight
135 * Mathematically, u = v * m.
136 * Input: v - input vector
137 * m - transformation matrix
138 * Output: u - transformed vector
140 void _mesa_transform_vector( GLfloat u
[4], const GLfloat v
[4], const GLfloat m
[16] )
142 GLfloat v0
=v
[0], v1
=v
[1], v2
=v
[2], v3
=v
[3];
143 #define M(row,col) m[row + col*4]
144 u
[0] = v0
* M(0,0) + v1
* M(1,0) + v2
* M(2,0) + v3
* M(3,0);
145 u
[1] = v0
* M(0,1) + v1
* M(1,1) + v2
* M(2,1) + v3
* M(3,1);
146 u
[2] = v0
* M(0,2) + v1
* M(1,2) + v2
* M(2,2) + v3
* M(3,2);
147 u
[3] = v0
* M(0,3) + v1
* M(1,3) + v2
* M(2,3) + v3
* M(3,3);
152 /* Useful for one-off point transformations, as in clipping.
153 * Note that because the matrix isn't analysed we do too many
154 * multiplies, and that the result is always 4-clean.
156 void _mesa_transform_point_sz( GLfloat Q
[4], const GLfloat M
[16],
157 const GLfloat P
[4], GLuint sz
)
164 Q
[0] = M
[0] * P
[0] + M
[4] * P
[1] + M
[8] * P
[2] + M
[12] * P
[3];
165 Q
[1] = M
[1] * P
[0] + M
[5] * P
[1] + M
[9] * P
[2] + M
[13] * P
[3];
166 Q
[2] = M
[2] * P
[0] + M
[6] * P
[1] + M
[10] * P
[2] + M
[14] * P
[3];
167 Q
[3] = M
[3] * P
[0] + M
[7] * P
[1] + M
[11] * P
[2] + M
[15] * P
[3];
171 Q
[0] = M
[0] * P
[0] + M
[4] * P
[1] + M
[8] * P
[2] + M
[12];
172 Q
[1] = M
[1] * P
[0] + M
[5] * P
[1] + M
[9] * P
[2] + M
[13];
173 Q
[2] = M
[2] * P
[0] + M
[6] * P
[1] + M
[10] * P
[2] + M
[14];
174 Q
[3] = M
[3] * P
[0] + M
[7] * P
[1] + M
[11] * P
[2] + M
[15];
178 Q
[0] = M
[0] * P
[0] + M
[4] * P
[1] + M
[12];
179 Q
[1] = M
[1] * P
[0] + M
[5] * P
[1] + M
[13];
180 Q
[2] = M
[2] * P
[0] + M
[6] * P
[1] + M
[14];
181 Q
[3] = M
[3] * P
[0] + M
[7] * P
[1] + M
[15];
185 Q
[0] = M
[0] * P
[0] + M
[12];
186 Q
[1] = M
[1] * P
[0] + M
[13];
187 Q
[2] = M
[2] * P
[0] + M
[14];
188 Q
[3] = M
[3] * P
[0] + M
[15];
194 * This is called only once. It initializes several tables with pointers
195 * to optimized transformation functions. This is where we can test for
196 * AMD 3Dnow! capability, Intel SSE, etc. and hook in the right code.
199 _math_init_transformation( void )
201 init_c_transformations();
202 init_c_norm_transform();
208 _math_test_all_transform_functions( "default" );
209 _math_test_all_normal_transform_functions( "default" );
210 _math_test_all_cliptest_functions( "default" );
214 _mesa_init_all_x86_transform_asm();
215 #elif defined( USE_SPARC_ASM )
216 _mesa_init_all_sparc_transform_asm();
217 #elif defined( USE_PPC_ASM )
218 _mesa_init_all_ppc_transform_asm();
219 #elif defined( USE_X86_64_ASM )
220 _mesa_init_all_x86_64_transform_asm();
227 _math_init_transformation();
228 _math_init_translate();