This commit was generated by cvs2svn to compensate for changes in r2,
[btb/d2x.git] / 3d / points.c
1 /*
2 THE COMPUTER CODE CONTAINED HEREIN IS THE SOLE PROPERTY OF PARALLAX
3 SOFTWARE CORPORATION ("PARALLAX").  PARALLAX, IN DISTRIBUTING THE CODE TO
4 END-USERS, AND SUBJECT TO ALL OF THE TERMS AND CONDITIONS HEREIN, GRANTS A
5 ROYALTY-FREE, PERPETUAL LICENSE TO SUCH END-USERS FOR USE BY SUCH END-USERS
6 IN USING, DISPLAYING,  AND CREATING DERIVATIVE WORKS THEREOF, SO LONG AS
7 SUCH USE, DISPLAY OR CREATION IS FOR NON-COMMERCIAL, ROYALTY OR REVENUE
8 FREE PURPOSES.  IN NO EVENT SHALL THE END-USER USE THE COMPUTER CODE
9 CONTAINED HEREIN FOR REVENUE-BEARING PURPOSES.  THE END-USER UNDERSTANDS
10 AND AGREES TO THE TERMS HEREIN AND ACCEPTS THE SAME BY USE OF THIS FILE.  
11 COPYRIGHT 1993-1998 PARALLAX SOFTWARE CORPORATION.  ALL RIGHTS RESERVED.
12 */
13 /*
14  * $Source: /cvs/cvsroot/d2x/3d/points.c,v $
15  * $Revision: 1.1.1.1 $
16  * $Author: bradleyb $
17  * $Date: 2001-01-19 03:29:58 $
18  * 
19  * Routines for point definition, rotation, etc.
20  * 
21  * $Log: not supported by cvs2svn $
22  * Revision 1.1.1.1  1999/06/14 21:57:49  donut
23  * Import of d1x 1.37 source.
24  *
25  * Revision 1.3  1995/09/21  17:29:40  allender
26  * changed project_point to overflow if z <= 0
27  *
28  * Revision 1.2  1995/09/13  11:31:28  allender
29  * removed checkmuldiv from g3_project_point
30  *
31  * Revision 1.1  1995/05/05  08:52:35  allender
32  * Initial revision
33  *
34  * Revision 1.1  1995/04/17  04:32:25  matt
35  * Initial revision
36  * 
37  * 
38  */
39
40 #ifdef RCS
41 static char rcsid[] = "$Id: points.c,v 1.1.1.1 2001-01-19 03:29:58 bradleyb Exp $";
42 #endif
43
44
45 #include <conf.h>
46 #include "fix.h"
47 #include "vecmat.h"
48 #include "gr.h"
49 #include "3d.h"
50 #include "globvars.h"
51
52
53 //code a point.  fills in the p3_codes field of the point, and returns the codes
54 ubyte g3_code_point(g3s_point *p)
55 {
56         ubyte cc=0;
57
58         if (p->p3_x > p->p3_z)
59                 cc |= CC_OFF_RIGHT;
60
61         if (p->p3_y > p->p3_z)
62                 cc |= CC_OFF_TOP;
63
64         if (p->p3_x < -p->p3_z)
65                 cc |= CC_OFF_LEFT;
66
67         if (p->p3_y < -p->p3_z)
68                 cc |= CC_OFF_BOT;
69
70         if (p->p3_z < 0)
71                 cc |= CC_BEHIND;
72
73         return p->p3_codes = cc;
74
75 }
76
77 //rotates a point. returns codes.  does not check if already rotated
78 ubyte g3_rotate_point(g3s_point *dest,vms_vector *src)
79 {
80         vms_vector tempv;
81
82 #ifdef D1XD3D
83         dest->p3_orig = *src;
84 #endif
85
86         vm_vec_sub(&tempv,src,&View_position);
87
88         vm_vec_rotate(&dest->p3_vec,&tempv,&View_matrix);
89
90         dest->p3_flags = 0;     //no projected
91
92         return g3_code_point(dest);
93
94 }
95
96 //checks for overflow & divides if ok, fillig in r
97 //returns true if div is ok, else false
98 int checkmuldiv(fix *r,fix a,fix b,fix c)
99 {
100         quad q,qt;
101
102         q.low=q.high=0;
103         fixmulaccum(&q,a,b);
104
105         qt = q;
106         if (qt.high < 0)
107                 fixquadnegate(&qt);
108
109         qt.high *= 2;
110         if (qt.low > 0x7fff)
111                 qt.high++;
112
113         if (qt.high >= c)
114                 return 0;
115         else {
116                 *r = fixdivquadlong(q.low,q.high,c);
117                 return 1;
118         }
119 }
120
121 //projects a point
122 void g3_project_point(g3s_point *p)
123 {
124 #ifndef __powerc
125         fix tx,ty;
126
127         if (p->p3_flags & PF_PROJECTED || p->p3_codes & CC_BEHIND)
128                 return;
129
130         if (checkmuldiv(&tx,p->p3_x,Canv_w2,p->p3_z) && checkmuldiv(&ty,p->p3_y,Canv_h2,p->p3_z)) {
131                 p->p3_sx = Canv_w2 + tx;
132                 p->p3_sy = Canv_h2 - ty;
133                 p->p3_flags |= PF_PROJECTED;
134         }
135         else
136                 p->p3_flags |= PF_OVERFLOW;
137 #else
138         double fz;
139         
140         if ((p->p3_flags & PF_PROJECTED) || (p->p3_codes & CC_BEHIND))
141                 return;
142         
143         if ( p->p3_z <= 0 )     {
144                 p->p3_flags |= PF_OVERFLOW;
145                 return;
146         }
147
148         fz = f2fl(p->p3_z);
149         p->p3_sx = fl2f(fCanv_w2 + (f2fl(p->p3_x)*fCanv_w2 / fz));
150         p->p3_sy = fl2f(fCanv_h2 - (f2fl(p->p3_y)*fCanv_h2 / fz));
151
152         p->p3_flags |= PF_PROJECTED;
153 #endif
154 }
155
156 //from a 2d point, compute the vector through that point
157 void g3_point_2_vec(vms_vector *v,short sx,short sy)
158 {
159         vms_vector tempv;
160         vms_matrix tempm;
161
162         tempv.x =  fixmuldiv(fixdiv((sx<<16) - Canv_w2,Canv_w2),Matrix_scale.z,Matrix_scale.x);
163         tempv.y = -fixmuldiv(fixdiv((sy<<16) - Canv_h2,Canv_h2),Matrix_scale.z,Matrix_scale.y);
164         tempv.z = f1_0;
165
166         vm_vec_normalize(&tempv);
167
168         vm_copy_transpose_matrix(&tempm,&Unscaled_matrix);
169
170         vm_vec_rotate(v,&tempv,&tempm);
171
172 }
173
174 //delta rotation functions
175 vms_vector *g3_rotate_delta_x(vms_vector *dest,fix dx)
176 {
177         dest->x = fixmul(View_matrix.rvec.x,dx);
178         dest->y = fixmul(View_matrix.uvec.x,dx);
179         dest->z = fixmul(View_matrix.fvec.x,dx);
180
181         return dest;
182 }
183
184 vms_vector *g3_rotate_delta_y(vms_vector *dest,fix dy)
185 {
186         dest->x = fixmul(View_matrix.rvec.y,dy);
187         dest->y = fixmul(View_matrix.uvec.y,dy);
188         dest->z = fixmul(View_matrix.fvec.y,dy);
189
190         return dest;
191 }
192
193 vms_vector *g3_rotate_delta_z(vms_vector *dest,fix dz)
194 {
195         dest->x = fixmul(View_matrix.rvec.z,dz);
196         dest->y = fixmul(View_matrix.uvec.z,dz);
197         dest->z = fixmul(View_matrix.fvec.z,dz);
198
199         return dest;
200 }
201
202
203 vms_vector *g3_rotate_delta_vec(vms_vector *dest,vms_vector *src)
204 {
205         return vm_vec_rotate(dest,src,&View_matrix);
206 }
207
208 ubyte g3_add_delta_vec(g3s_point *dest,g3s_point *src,vms_vector *deltav)
209 {
210         vm_vec_add(&dest->p3_vec,&src->p3_vec,deltav);
211
212         dest->p3_flags = 0;             //not projected
213
214         return g3_code_point(dest);
215 }
216
217 //calculate the depth of a point - returns the z coord of the rotated point
218 fix g3_calc_point_depth(vms_vector *pnt)
219 {
220         quad q;
221
222         q.low=q.high=0;
223         fixmulaccum(&q,(pnt->x - View_position.x),View_matrix.fvec.x);
224         fixmulaccum(&q,(pnt->y - View_position.y),View_matrix.fvec.y);
225         fixmulaccum(&q,(pnt->z - View_position.z),View_matrix.fvec.z);
226
227         return fixquadadjust(&q);
228 }
229
230
231