got rid of leafnums array in edict structure for pvs checking, pvs is now checked...
[divverent/darkplaces.git] / world.c
1 /*
2 Copyright (C) 1996-1997 Id Software, Inc.
3
4 This program is free software; you can redistribute it and/or
5 modify it under the terms of the GNU General Public License
6 as published by the Free Software Foundation; either version 2
7 of the License, or (at your option) any later version.
8
9 This program is distributed in the hope that it will be useful,
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
12
13 See the GNU General Public License for more details.
14
15 You should have received a copy of the GNU General Public License
16 along with this program; if not, write to the Free Software
17 Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
18
19 */
20 // world.c -- world query functions
21
22 #include "quakedef.h"
23
24 /*
25
26 entities never clip against themselves, or their owner
27
28 line of sight checks trace->crosscontent, but bullets don't
29
30 */
31
32 /*
33 typedef struct link_s
34 {
35         struct link_s   *prev, *next;
36 } link_t;
37 */
38
39
40 void ClearLink (link_t *l);
41 void RemoveLink (link_t *l);
42 void InsertLinkBefore (link_t *l, link_t *before);
43 void InsertLinkAfter (link_t *l, link_t *after);
44
45 // (type *)STRUCT_FROM_LINK(link_t *link, type, member)
46 // ent = STRUCT_FROM_LINK(link,entity_t,order)
47 // FIXME: remove this mess!
48 //#define       STRUCT_FROM_LINK(l,t,m) ((t *)((byte *)l - (int)&(((t *)0)->m)))
49
50 #define EDICT_FROM_AREA(l) ((edict_t *)((byte *)l - (int)&(((edict_t *)0)->area)))
51
52 //============================================================================
53
54 // ClearLink is used for new headnodes
55 void ClearLink (link_t *l)
56 {
57         l->prev = l->next = l;
58 }
59
60 void RemoveLink (link_t *l)
61 {
62         l->next->prev = l->prev;
63         l->prev->next = l->next;
64 }
65
66 void InsertLinkBefore (link_t *l, link_t *before)
67 {
68         l->next = before;
69         l->prev = before->prev;
70         l->prev->next = l;
71         l->next->prev = l;
72 }
73 void InsertLinkAfter (link_t *l, link_t *after)
74 {
75         l->next = after->next;
76         l->prev = after;
77         l->prev->next = l;
78         l->next->prev = l;
79 }
80
81
82 typedef struct
83 {
84         vec3_t          boxmins, boxmaxs;// enclose the test object along entire move
85         float           *mins, *maxs;   // size of the moving object
86         vec3_t          mins2, maxs2;   // size when clipping against mosnters
87         float           *start, *end;
88         trace_t         trace;
89         int                     type;
90         edict_t         *passedict;
91 } moveclip_t;
92
93
94 /*
95 ===============================================================================
96
97 HULL BOXES
98
99 ===============================================================================
100 */
101
102
103 static  hull_t          box_hull;
104 static  dclipnode_t     box_clipnodes[6];
105 static  mplane_t        box_planes[6];
106
107 /*
108 ===================
109 SV_InitBoxHull
110
111 Set up the planes and clipnodes so that the six floats of a bounding box
112 can just be stored out and get a proper hull_t structure.
113 ===================
114 */
115 void SV_InitBoxHull (void)
116 {
117         int             i;
118         int             side;
119
120         box_hull.clipnodes = box_clipnodes;
121         box_hull.planes = box_planes;
122         box_hull.firstclipnode = 0;
123         box_hull.lastclipnode = 5;
124
125         for (i=0 ; i<6 ; i++)
126         {
127                 box_clipnodes[i].planenum = i;
128                 
129                 side = i&1;
130                 
131                 box_clipnodes[i].children[side] = CONTENTS_EMPTY;
132                 if (i != 5)
133                         box_clipnodes[i].children[side^1] = i + 1;
134                 else
135                         box_clipnodes[i].children[side^1] = CONTENTS_SOLID;
136                 
137                 box_planes[i].type = i>>1;
138                 box_planes[i].normal[i>>1] = 1;
139         }
140         
141 }
142
143
144 /*
145 ===================
146 SV_HullForBox
147
148 To keep everything totally uniform, bounding boxes are turned into small
149 BSP trees instead of being compared directly.
150 ===================
151 */
152 hull_t  *SV_HullForBox (vec3_t mins, vec3_t maxs)
153 {
154         box_planes[0].dist = maxs[0];
155         box_planes[1].dist = mins[0];
156         box_planes[2].dist = maxs[1];
157         box_planes[3].dist = mins[1];
158         box_planes[4].dist = maxs[2];
159         box_planes[5].dist = mins[2];
160
161         return &box_hull;
162 }
163
164
165
166 /*
167 ================
168 SV_HullForEntity
169
170 Returns a hull that can be used for testing or clipping an object of mins/maxs
171 size.
172 Offset is filled in to contain the adjustment that must be added to the
173 testing object's origin to get a point to use with the returned hull.
174 ================
175 */
176 hull_t *SV_HullForEntity (edict_t *ent, vec3_t mins, vec3_t maxs, vec3_t offset)
177 {
178         model_t         *model;
179         vec3_t          size;
180         vec3_t          hullmins, hullmaxs;
181         hull_t          *hull;
182
183 // decide which clipping hull to use, based on the size
184         if (ent->v.solid == SOLID_BSP)
185         {       // explicit hulls in the BSP model
186                 if (ent->v.movetype != MOVETYPE_PUSH)
187                         Host_Error ("SOLID_BSP without MOVETYPE_PUSH");
188
189                 model = sv.models[ (int)ent->v.modelindex ];
190
191                 // LordHavoc: fixed SOLID_BSP error message
192                 if (!model || model->type != mod_brush)
193                 {
194                         Con_Printf ("SOLID_BSP with a non bsp model, entity dump:\n");
195                         ED_Print (ent);
196                         Host_Error ("SOLID_BSP with a non bsp model\n");
197                 }
198
199                 VectorSubtract (maxs, mins, size);
200                 // LordHavoc: FIXME!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
201                 if (hlbsp)
202                 {
203                         if (size[0] < 3)
204                                 hull = &model->hulls[0]; // 0x0x0
205                         else if (size[0] <= 32)
206                         {
207                                 if (size[2] < 54) // pick the nearest of 36 or 72
208                                         hull = &model->hulls[3]; // 32x32x36
209                                 else
210                                         hull = &model->hulls[1]; // 32x32x72
211                         }
212                         else
213                                 hull = &model->hulls[2]; // 64x64x64
214                 }
215                 else
216                 {
217                         if (size[0] < 3)
218                                 hull = &model->hulls[0]; // 0x0x0
219                         else if (size[0] <= 32)
220                                 hull = &model->hulls[1]; // 32x32x56
221                         else
222                                 hull = &model->hulls[2]; // 64x64x88
223                 }
224
225 // calculate an offset value to center the origin
226                 VectorSubtract (hull->clip_mins, mins, offset);
227                 VectorAdd (offset, ent->v.origin, offset);
228         }
229         else
230         {       // create a temp hull from bounding box sizes
231
232                 VectorSubtract (ent->v.mins, maxs, hullmins);
233                 VectorSubtract (ent->v.maxs, mins, hullmaxs);
234                 hull = SV_HullForBox (hullmins, hullmaxs);
235
236                 VectorCopy (ent->v.origin, offset);
237         }
238
239
240         return hull;
241 }
242
243 /*
244 ===============================================================================
245
246 ENTITY AREA CHECKING
247
248 ===============================================================================
249 */
250
251 typedef struct areanode_s
252 {
253         int             axis;           // -1 = leaf node
254         float   dist;
255         struct areanode_s       *children[2];
256         link_t  trigger_edicts;
257         link_t  solid_edicts;
258 } areanode_t;
259
260 #define AREA_DEPTH      4
261 #define AREA_NODES      32
262
263 static  areanode_t      sv_areanodes[AREA_NODES];
264 static  int                     sv_numareanodes;
265
266 /*
267 ===============
268 SV_CreateAreaNode
269
270 ===============
271 */
272 areanode_t *SV_CreateAreaNode (int depth, vec3_t mins, vec3_t maxs)
273 {
274         areanode_t      *anode;
275         vec3_t          size;
276         vec3_t          mins1, maxs1, mins2, maxs2;
277
278         anode = &sv_areanodes[sv_numareanodes];
279         sv_numareanodes++;
280
281         ClearLink (&anode->trigger_edicts);
282         ClearLink (&anode->solid_edicts);
283
284         if (depth == AREA_DEPTH)
285         {
286                 anode->axis = -1;
287                 anode->children[0] = anode->children[1] = NULL;
288                 return anode;
289         }
290
291         VectorSubtract (maxs, mins, size);
292         if (size[0] > size[1])
293                 anode->axis = 0;
294         else
295                 anode->axis = 1;
296
297         anode->dist = 0.5 * (maxs[anode->axis] + mins[anode->axis]);
298         VectorCopy (mins, mins1);
299         VectorCopy (mins, mins2);
300         VectorCopy (maxs, maxs1);
301         VectorCopy (maxs, maxs2);
302
303         maxs1[anode->axis] = mins2[anode->axis] = anode->dist;
304
305         anode->children[0] = SV_CreateAreaNode (depth+1, mins2, maxs2);
306         anode->children[1] = SV_CreateAreaNode (depth+1, mins1, maxs1);
307
308         return anode;
309 }
310
311 /*
312 ===============
313 SV_ClearWorld
314
315 ===============
316 */
317 void SV_ClearWorld (void)
318 {
319         SV_InitBoxHull ();
320
321         memset (sv_areanodes, 0, sizeof(sv_areanodes));
322         sv_numareanodes = 0;
323         SV_CreateAreaNode (0, sv.worldmodel->mins, sv.worldmodel->maxs);
324 }
325
326
327 /*
328 ===============
329 SV_UnlinkEdict
330
331 ===============
332 */
333 void SV_UnlinkEdict (edict_t *ent)
334 {
335         if (!ent->area.prev)
336                 return;         // not linked in anywhere
337         RemoveLink (&ent->area);
338         ent->area.prev = ent->area.next = NULL;
339 }
340
341
342 /*
343 ====================
344 SV_TouchLinks
345 ====================
346 */
347 void SV_TouchLinks ( edict_t *ent, areanode_t *node )
348 {
349         link_t          *l, *next;
350         edict_t         *touch;
351         int                     old_self, old_other;
352
353 loc0:
354 // touch linked edicts
355         for (l = node->trigger_edicts.next ; l != &node->trigger_edicts ; l = next)
356         {
357                 next = l->next;
358                 touch = EDICT_FROM_AREA(l);
359                 if (touch == ent)
360                         continue;
361                 if (!touch->v.touch || touch->v.solid != SOLID_TRIGGER)
362                         continue;
363                 if (ent->v.absmin[0] > touch->v.absmax[0]
364                  || ent->v.absmin[1] > touch->v.absmax[1]
365                  || ent->v.absmin[2] > touch->v.absmax[2]
366                  || ent->v.absmax[0] < touch->v.absmin[0]
367                  || ent->v.absmax[1] < touch->v.absmin[1]
368                  || ent->v.absmax[2] < touch->v.absmin[2])
369                         continue;
370                 old_self = pr_global_struct->self;
371                 old_other = pr_global_struct->other;
372
373                 pr_global_struct->self = EDICT_TO_PROG(touch);
374                 pr_global_struct->other = EDICT_TO_PROG(ent);
375                 pr_global_struct->time = sv.time;
376                 PR_ExecuteProgram (touch->v.touch, "");
377
378                 pr_global_struct->self = old_self;
379                 pr_global_struct->other = old_other;
380         }
381
382 // recurse down both sides
383         if (node->axis == -1)
384                 return;
385
386         // LordHavoc: optimized recursion
387 //      if (ent->v.absmax[node->axis] > node->dist) SV_TouchLinks (ent, node->children[0]);
388 //      if (ent->v.absmin[node->axis] < node->dist) SV_TouchLinks (ent, node->children[1]);
389         if (ent->v.absmax[node->axis] > node->dist)
390         {
391                 if (ent->v.absmin[node->axis] < node->dist)
392                         SV_TouchLinks(ent, node->children[1]); // order reversed to reduce code
393                 node = node->children[0];
394                 goto loc0;
395         }
396         else
397         {
398                 if (ent->v.absmin[node->axis] < node->dist)
399                 {
400                         node = node->children[1];
401                         goto loc0;
402                 }
403         }
404 }
405
406
407 /*
408 ===============
409 SV_LinkEdict
410
411 ===============
412 */
413 void SV_LinkEdict (edict_t *ent, qboolean touch_triggers)
414 {
415         areanode_t      *node;
416
417         if (ent->area.prev)
418                 SV_UnlinkEdict (ent);   // unlink from old position
419                 
420         if (ent == sv.edicts)
421                 return;         // don't add the world
422
423         if (ent->free)
424                 return;
425
426 // set the abs box
427
428 // LordHavoc: enabling rotating bmodels
429         if (ent->v.solid == SOLID_BSP && (ent->v.angles[0] || ent->v.angles[1] || ent->v.angles[2]))
430         {
431                 // expand for rotation
432                 float           max, v;
433                 int                     i;
434
435                 max = DotProduct(ent->v.mins, ent->v.mins);
436                 v = DotProduct(ent->v.maxs, ent->v.maxs);
437                 if (max < v)
438                         max = v;
439                 max = sqrt(max);
440                 /*
441                 max = 0;
442                 for (i=0 ; i<3 ; i++)
443                 {
444                         v = fabs(ent->v.mins[i]);
445                         if (max < v)
446                                 max = v;
447                         v = fabs(ent->v.maxs[i]);
448                         if (max < v)
449                                 max = v;
450                 }
451                 */
452                 for (i=0 ; i<3 ; i++)
453                 {
454                         ent->v.absmin[i] = ent->v.origin[i] - max;
455                         ent->v.absmax[i] = ent->v.origin[i] + max;
456                 }
457         }
458         else
459         {
460                 VectorAdd (ent->v.origin, ent->v.mins, ent->v.absmin);
461                 VectorAdd (ent->v.origin, ent->v.maxs, ent->v.absmax);
462         }
463
464 //
465 // to make items easier to pick up and allow them to be grabbed off
466 // of shelves, the abs sizes are expanded
467 //
468         if ((int)ent->v.flags & FL_ITEM)
469         {
470                 ent->v.absmin[0] -= 15;
471                 ent->v.absmin[1] -= 15;
472                 ent->v.absmax[0] += 15;
473                 ent->v.absmax[1] += 15;
474         }
475         else
476         {
477                 // because movement is clipped an epsilon away from an actual edge,
478                 // we must fully check even when bounding boxes don't quite touch
479                 ent->v.absmin[0] -= 1;
480                 ent->v.absmin[1] -= 1;
481                 ent->v.absmin[2] -= 1;
482                 ent->v.absmax[0] += 1;
483                 ent->v.absmax[1] += 1;
484                 ent->v.absmax[2] += 1;
485         }
486
487         if (ent->v.solid == SOLID_NOT)
488                 return;
489
490 // find the first node that the ent's box crosses
491         node = sv_areanodes;
492         while (1)
493         {
494                 if (node->axis == -1)
495                         break;
496                 if (ent->v.absmin[node->axis] > node->dist)
497                         node = node->children[0];
498                 else if (ent->v.absmax[node->axis] < node->dist)
499                         node = node->children[1];
500                 else
501                         break;          // crosses the node
502         }
503
504 // link it in
505
506         if (ent->v.solid == SOLID_TRIGGER)
507                 InsertLinkBefore (&ent->area, &node->trigger_edicts);
508         else
509                 InsertLinkBefore (&ent->area, &node->solid_edicts);
510
511 // if touch_triggers, touch all entities at this node and descend for more
512         if (touch_triggers)
513                 SV_TouchLinks ( ent, sv_areanodes );
514 }
515
516
517
518 /*
519 ===============================================================================
520
521 POINT TESTING IN HULLS
522
523 ===============================================================================
524 */
525
526 /*
527 ==================
528 SV_HullPointContents
529
530 ==================
531 */
532 int SV_HullPointContents (hull_t *hull, int num, vec3_t p)
533 {
534         while (num >= 0)
535                 num = hull->clipnodes[num].children[(hull->planes[hull->clipnodes[num].planenum].type < 3 ? p[hull->planes[hull->clipnodes[num].planenum].type] : DotProduct (hull->planes[hull->clipnodes[num].planenum].normal, p)) < hull->planes[hull->clipnodes[num].planenum].dist];
536
537         return num;
538 }
539
540 /*
541 ============
542 SV_TestEntityPosition
543
544 This could be a lot more efficient...
545 ============
546 */
547 edict_t *SV_TestEntityPosition (edict_t *ent)
548 {
549         trace_t trace;
550
551         trace = SV_Move (ent->v.origin, ent->v.mins, ent->v.maxs, ent->v.origin, MOVE_NORMAL, ent);
552
553         if (trace.startsolid)
554                 return sv.edicts;
555                 
556         return NULL;
557 }
558
559
560 /*
561 ===============================================================================
562
563 LINE TESTING IN HULLS
564
565 ===============================================================================
566 */
567
568 // 1/32 epsilon to keep floating point happy
569 #define DIST_EPSILON    (0.03125)
570
571 /*
572 ==================
573 SV_RecursiveHullCheck
574
575 ==================
576 */
577 /*
578 qboolean SV_RecursiveHullCheck (hull_t *hull, int num, float p1f, float p2f, vec3_t p1, vec3_t p2, trace_t *trace)
579 {
580         dclipnode_t     *node;
581         mplane_t        *plane;
582         float           t1, t2;
583         float           frac;
584         int                     i;
585         vec3_t          mid;
586         int                     side;
587         float           midf;
588
589 loc0:
590 // check for empty
591         if (num < 0)
592         {
593                 if (num != CONTENTS_SOLID)
594                 {
595                         trace->allsolid = false;
596                         if (num == CONTENTS_EMPTY)
597                                 trace->inopen = true;
598                         else
599                                 trace->inwater = true;
600                 }
601                 else
602                         trace->startsolid = true;
603                 return true;            // empty
604         }
605
606         if (num < hull->firstclipnode || num > hull->lastclipnode)
607                 Sys_Error ("SV_RecursiveHullCheck: bad node number");
608
609 //
610 // find the point distances
611 //
612         node = hull->clipnodes + num;
613         plane = hull->planes + node->planenum;
614
615         t1 = PlaneDiff(p1, plane);
616         t2 = PlaneDiff(p2, plane);
617         
618 #if 1
619         if (t1 >= 0 && t2 >= 0)
620         // LordHavoc: optimized recursion
621 //              return SV_RecursiveHullCheck (hull, node->children[0], p1f, p2f, p1, p2, trace);
622         {
623                 num = node->children[0];
624                 goto loc0;
625         }
626         if (t1 < 0 && t2 < 0)
627 //              return SV_RecursiveHullCheck (hull, node->children[1], p1f, p2f, p1, p2, trace);
628         {
629                 num = node->children[1];
630                 goto loc0;
631         }
632 #else
633         if ( (t1 >= DIST_EPSILON && t2 >= DIST_EPSILON) || (t2 > t1 && t1 >= 0) )
634                 return SV_RecursiveHullCheck (hull, node->children[0], p1f, p2f, p1, p2, trace);
635         if ( (t1 <= -DIST_EPSILON && t2 <= -DIST_EPSILON) || (t2 < t1 && t1 <= 0) )
636                 return SV_RecursiveHullCheck (hull, node->children[1], p1f, p2f, p1, p2, trace);
637 #endif
638
639 // put the crosspoint DIST_EPSILON pixels on the near side
640         if (t1 < 0)
641                 frac = bound(0, (t1 + DIST_EPSILON)/(t1-t2), 1);
642         else
643                 frac = bound(0, (t1 - DIST_EPSILON)/(t1-t2), 1);
644                 
645         midf = p1f + (p2f - p1f)*frac;
646         mid[0] = p1[0] + frac*(p2[0] - p1[0]);
647         mid[1] = p1[1] + frac*(p2[1] - p1[1]);
648         mid[2] = p1[2] + frac*(p2[2] - p1[2]);
649
650         side = (t1 < 0);
651
652 // move up to the node
653         if (!SV_RecursiveHullCheck (hull, node->children[side], p1f, midf, p1, mid, trace) )
654                 return false;
655
656 #ifdef PARANOID
657         if (SV_HullPointContents (hull, node->children[side], mid) == CONTENTS_SOLID)
658         {
659                 Con_Printf ("mid PointInHullSolid\n");
660                 return false;
661         }
662 #endif
663
664         if (SV_HullPointContents (hull, node->children[side^1], mid) != CONTENTS_SOLID)
665 // go past the node
666                 return SV_RecursiveHullCheck (hull, node->children[side^1], midf, p2f, mid, p2, trace);
667         // mid would need to be duplicated during recursion...
668 */
669         /*
670         {
671                 p1f = midf;
672                 p1 = mid;
673                 num = node->children[side^1];
674                 goto loc0;
675         }
676         */
677 /*
678
679         if (trace->allsolid)
680                 return false;           // never got out of the solid area
681         
682 //==================
683 // the other side of the node is solid, this is the impact point
684 //==================
685         if (!side)
686         {
687                 VectorCopy (plane->normal, trace->plane.normal);
688                 trace->plane.dist = plane->dist;
689         }
690         else
691         {
692                 VectorNegate (plane->normal, trace->plane.normal);
693                 trace->plane.dist = -plane->dist;
694         }
695
696         while (SV_HullPointContents (hull, hull->firstclipnode, mid) == CONTENTS_SOLID)
697         { // shouldn't really happen, but does occasionally
698                 frac -= 0.1;
699                 if (frac < 0)
700                 {
701                         trace->fraction = midf;
702                         VectorCopy (mid, trace->endpos);
703                         Con_DPrintf ("backup past 0\n");
704                         return false;
705                 }
706                 midf = p1f + (p2f - p1f)*frac;
707                 for (i=0 ; i<3 ; i++)
708                         mid[i] = p1[i] + frac*(p2[i] - p1[i]);
709         }
710
711         trace->fraction = midf;
712         VectorCopy (mid, trace->endpos);
713
714         return false;
715 }
716 */
717
718 // LordHavoc: backported from my optimizations to PM_RecursiveHullCheck in QuakeForge newtree (QW)
719 qboolean SV_RecursiveHullCheck (hull_t *hull, int num, float p1f, float p2f, vec3_t p1, vec3_t p2, trace_t *trace)
720 {
721         dclipnode_t     *node;
722         mplane_t        *plane;
723         float           t1, t2;
724         float           frac;
725         int                     i;
726         vec3_t          mid;
727         int                     side;
728         float           midf;
729
730         // LordHavoc: a goto!  everyone flee in terror... :)
731 loc0:
732 // check for empty
733         if (num < 0)
734         {
735                 if (num != CONTENTS_SOLID)
736                 {
737                         trace->allsolid = false;
738                         if (num == CONTENTS_EMPTY)
739                                 trace->inopen = true;
740                         else
741                                 trace->inwater = true;
742                 }
743                 else
744                         trace->startsolid = true;
745                 return true;            // empty
746         }
747
748 // find the point distances
749         node = hull->clipnodes + num;
750         plane = hull->planes + node->planenum;
751
752         if (plane->type < 3)
753         {
754                 t1 = p1[plane->type] - plane->dist;
755                 t2 = p2[plane->type] - plane->dist;
756         }
757         else
758         {
759                 t1 = DotProduct (plane->normal, p1) - plane->dist;
760                 t2 = DotProduct (plane->normal, p2) - plane->dist;
761         }
762
763         // LordHavoc: recursion optimization
764         if (t1 >= 0 && t2 >= 0)
765         {
766                 num = node->children[0];
767                 goto loc0;
768         }
769         if (t1 < 0 && t2 < 0)
770         {
771                 num = node->children[1];
772                 goto loc0;
773         }
774
775 // put the crosspoint DIST_EPSILON pixels on the near side
776         side = (t1 < 0);
777         if (side)
778                 frac = bound(0, (t1 + DIST_EPSILON) / (t1 - t2), 1);
779         else
780                 frac = bound(0, (t1 - DIST_EPSILON) / (t1 - t2), 1);
781                 
782         midf = p1f + (p2f - p1f)*frac;
783         for (i=0 ; i<3 ; i++)
784                 mid[i] = p1[i] + frac*(p2[i] - p1[i]);
785
786 // move up to the node
787         if (!SV_RecursiveHullCheck (hull, node->children[side], p1f, midf, p1, mid, trace) )
788                 return false;
789
790 #ifdef PARANOID
791         if (SV_HullPointContents (pm_hullmodel, mid, node->children[side]) == CONTENTS_SOLID)
792         {
793                 Con_Printf ("mid PointInHullSolid\n");
794                 return false;
795         }
796 #endif
797
798         // LordHavoc: warning to the clumsy, this recursion can not be optimized because mid would need to be duplicated on a stack
799         if (SV_HullPointContents (hull, node->children[side^1], mid) != CONTENTS_SOLID)
800 // go past the node
801                 return SV_RecursiveHullCheck (hull, node->children[side^1], midf, p2f, mid, p2, trace);
802         
803         if (trace->allsolid)
804                 return false;           // never got out of the solid area
805                 
806 //==================
807 // the other side of the node is solid, this is the impact point
808 //==================
809         if (!side)
810         {
811                 VectorCopy (plane->normal, trace->plane.normal);
812                 trace->plane.dist = plane->dist;
813         }
814         else
815         {
816                 VectorNegate (plane->normal, trace->plane.normal);
817                 trace->plane.dist = -plane->dist;
818         }
819
820         while (SV_HullPointContents (hull, hull->firstclipnode, mid) == CONTENTS_SOLID)
821         { // shouldn't really happen, but does occasionally
822                 frac -= 0.1;
823                 if (frac < 0)
824                 {
825                         trace->fraction = midf;
826                         VectorCopy (mid, trace->endpos);
827                         Con_DPrintf ("backup past 0\n");
828                         return false;
829                 }
830                 midf = p1f + (p2f - p1f)*frac;
831                 for (i=0 ; i<3 ; i++)
832                         mid[i] = p1[i] + frac*(p2[i] - p1[i]);
833         }
834
835         trace->fraction = midf;
836         VectorCopy (mid, trace->endpos);
837
838         return false;
839 }
840
841 qboolean SV_TestLine (hull_t *hull, int num, vec3_t p1, vec3_t p2)
842 {
843         dclipnode_t     *node;
844         mplane_t        *plane;
845         float           t1, t2, frac;
846         vec3_t          mid;
847         int                     side;
848
849 loc0:
850 // check for empty
851         if (num < 0)
852                 return num != CONTENTS_SOLID;
853
854         if (num < hull->firstclipnode || num > hull->lastclipnode)
855                 Sys_Error ("SV_RecursiveHullCheck: bad node number");
856
857 //
858 // find the point distances
859 //
860         node = hull->clipnodes + num;
861         plane = hull->planes + node->planenum;
862
863         t1 = PlaneDiff(p1, plane);
864         t2 = PlaneDiff(p2, plane);
865         
866         if (t1 >= 0 && t2 >= 0)
867         {
868                 num = node->children[0];
869                 goto loc0;
870         }
871         if (t1 < 0 && t2 < 0)
872         {
873                 num = node->children[1];
874                 goto loc0;
875         }
876
877 // put the crosspoint DIST_EPSILON pixels on the near side
878         side = (t1 < 0);
879
880         if (side)
881                 frac = bound(0, (t1 + DIST_EPSILON)/(t1-t2), 1);
882         else
883                 frac = bound(0, (t1 - DIST_EPSILON)/(t1-t2), 1);
884                 
885         mid[0] = p1[0] + frac*(p2[0] - p1[0]);
886         mid[1] = p1[1] + frac*(p2[1] - p1[1]);
887         mid[2] = p1[2] + frac*(p2[2] - p1[2]);
888
889         if (node->children[side] < 0)
890         {
891                 if (node->children[side] == CONTENTS_SOLID)
892                         return false;
893                 return SV_TestLine(hull, node->children[!side], mid, p2);
894         }
895         else if (SV_TestLine(hull, node->children[side], p1, mid))
896                 return SV_TestLine(hull, node->children[!side], mid, p2);
897         else
898                 return false;
899 }
900
901
902 /*
903 ==================
904 SV_ClipMoveToEntity
905
906 Handles selection or creation of a clipping hull, and offseting (and
907 eventually rotation) of the end points
908 ==================
909 */
910 trace_t SV_ClipMoveToEntity (edict_t *ent, vec3_t start, vec3_t mins, vec3_t maxs, vec3_t end)
911 {
912         trace_t         trace;
913         vec3_t          offset;
914         vec3_t          start_l, end_l;
915         hull_t          *hull;
916
917 // fill in a default trace
918         memset (&trace, 0, sizeof(trace_t));
919         trace.fraction = 1;
920         trace.allsolid = true;
921         VectorCopy (end, trace.endpos);
922
923 // get the clipping hull
924         hull = SV_HullForEntity (ent, mins, maxs, offset);
925
926         VectorSubtract (start, offset, start_l);
927         VectorSubtract (end, offset, end_l);
928
929 // LordHavoc: enabling rotating bmodels
930         // rotate start and end into the models frame of reference
931         if (ent->v.solid == SOLID_BSP && (ent->v.angles[0] || ent->v.angles[1] || ent->v.angles[2]))
932         {
933                 vec3_t  forward, right, up;
934                 vec3_t  temp;
935
936                 AngleVectors (ent->v.angles, forward, right, up);
937
938                 VectorCopy (start_l, temp);
939                 start_l[0] = DotProduct (temp, forward);
940                 start_l[1] = -DotProduct (temp, right);
941                 start_l[2] = DotProduct (temp, up);
942
943                 VectorCopy (end_l, temp);
944                 end_l[0] = DotProduct (temp, forward);
945                 end_l[1] = -DotProduct (temp, right);
946                 end_l[2] = DotProduct (temp, up);
947         }
948
949 // trace a line through the apropriate clipping hull
950         SV_RecursiveHullCheck (hull, hull->firstclipnode, 0, 1, start_l, end_l, &trace);
951
952 // LordHavoc: enabling rotating bmodels
953         // rotate endpos back to world frame of reference
954         if (ent->v.solid == SOLID_BSP && (ent->v.angles[0] || ent->v.angles[1] || ent->v.angles[2]))
955         {
956                 vec3_t  a;
957                 vec3_t  forward, right, up;
958                 vec3_t  temp;
959
960                 if (trace.fraction != 1)
961                 {
962                         VectorNegate (ent->v.angles, a);
963                         AngleVectors (a, forward, right, up);
964
965                         VectorCopy (trace.endpos, temp);
966                         trace.endpos[0] = DotProduct (temp, forward);
967                         trace.endpos[1] = -DotProduct (temp, right);
968                         trace.endpos[2] = DotProduct (temp, up);
969
970                         VectorCopy (trace.plane.normal, temp);
971                         trace.plane.normal[0] = DotProduct (temp, forward);
972                         trace.plane.normal[1] = -DotProduct (temp, right);
973                         trace.plane.normal[2] = DotProduct (temp, up);
974                 }
975         }
976
977 // fix trace up by the offset
978         if (trace.fraction != 1)
979                 VectorAdd (trace.endpos, offset, trace.endpos);
980
981 // did we clip the move?
982         if (trace.fraction < 1 || trace.startsolid  )
983                 trace.ent = ent;
984
985         return trace;
986 }
987
988 //===========================================================================
989
990 /*
991 ====================
992 SV_ClipToLinks
993
994 Mins and maxs enclose the entire area swept by the move
995 ====================
996 */
997 void SV_ClipToLinks ( areanode_t *node, moveclip_t *clip )
998 {
999         link_t          *l, *next;
1000         edict_t         *touch;
1001         trace_t         trace;
1002
1003 loc0:
1004 // touch linked edicts
1005         for (l = node->solid_edicts.next ; l != &node->solid_edicts ; l = next)
1006         {
1007                 next = l->next;
1008                 touch = EDICT_FROM_AREA(l);
1009                 if (touch->v.solid == SOLID_NOT)
1010                         continue;
1011                 if (touch == clip->passedict)
1012                         continue;
1013                 if (touch->v.solid == SOLID_TRIGGER)
1014                         Sys_Error ("Trigger in clipping list");
1015
1016                 if (clip->type == MOVE_NOMONSTERS && touch->v.solid != SOLID_BSP)
1017                         continue;
1018
1019                 if (clip->boxmins[0] > touch->v.absmax[0]
1020                  || clip->boxmins[1] > touch->v.absmax[1]
1021                  || clip->boxmins[2] > touch->v.absmax[2]
1022                  || clip->boxmaxs[0] < touch->v.absmin[0]
1023                  || clip->boxmaxs[1] < touch->v.absmin[1]
1024                  || clip->boxmaxs[2] < touch->v.absmin[2])
1025                         continue;
1026
1027                 if (clip->passedict != NULL && clip->passedict->v.size[0] && !touch->v.size[0])
1028                         continue;       // points never interact
1029
1030         // might intersect, so do an exact clip
1031                 if (clip->trace.allsolid)
1032                         return;
1033                 if (clip->passedict)
1034                 {
1035                         if (PROG_TO_EDICT(touch->v.owner) == clip->passedict)
1036                                 continue;       // don't clip against own missiles
1037                         if (PROG_TO_EDICT(clip->passedict->v.owner) == touch)
1038                                 continue;       // don't clip against owner
1039                         // LordHavoc: corpse code
1040                         if (clip->passedict->v.solid == SOLID_CORPSE && touch->v.solid == SOLID_SLIDEBOX)
1041                                 continue;
1042                         if (clip->passedict->v.solid == SOLID_SLIDEBOX && touch->v.solid == SOLID_CORPSE)
1043                                 continue;
1044                 }
1045
1046                 if ((int)touch->v.flags & FL_MONSTER)
1047                         trace = SV_ClipMoveToEntity (touch, clip->start, clip->mins2, clip->maxs2, clip->end);
1048                 else
1049                         trace = SV_ClipMoveToEntity (touch, clip->start, clip->mins, clip->maxs, clip->end);
1050                 if (trace.allsolid || trace.startsolid || trace.fraction < clip->trace.fraction)
1051                 {
1052                         trace.ent = touch;
1053                         if (clip->trace.startsolid)
1054                         {
1055                                 clip->trace = trace;
1056                                 clip->trace.startsolid = true;
1057                         }
1058                         else
1059                                 clip->trace = trace;
1060                 }
1061                 else if (trace.startsolid)
1062                         clip->trace.startsolid = true;
1063         }
1064
1065 // recurse down both sides
1066         if (node->axis == -1)
1067                 return;
1068
1069         // LordHavoc: optimized recursion
1070 //      if (clip->boxmaxs[node->axis] > node->dist) SV_ClipToLinks(node->children[0], clip);
1071 //      if (clip->boxmins[node->axis] < node->dist) SV_ClipToLinks(node->children[1], clip);
1072         if (clip->boxmaxs[node->axis] > node->dist)
1073         {
1074                 if (clip->boxmins[node->axis] < node->dist)
1075                         SV_ClipToLinks(node->children[1], clip);
1076                 node = node->children[0];
1077                 goto loc0;
1078         }
1079         else if (clip->boxmins[node->axis] < node->dist)
1080         {
1081                 node = node->children[1];
1082                 goto loc0;
1083         }
1084 }
1085
1086
1087 /*
1088 ==================
1089 SV_MoveBounds
1090 ==================
1091 */
1092 void SV_MoveBounds (vec3_t start, vec3_t mins, vec3_t maxs, vec3_t end, vec3_t boxmins, vec3_t boxmaxs)
1093 {
1094 #if 0
1095 // debug to test against everything
1096 boxmins[0] = boxmins[1] = boxmins[2] = -9999;
1097 boxmaxs[0] = boxmaxs[1] = boxmaxs[2] = 9999;
1098 #else
1099         int             i;
1100         
1101         for (i=0 ; i<3 ; i++)
1102         {
1103                 if (end[i] > start[i])
1104                 {
1105                         boxmins[i] = start[i] + mins[i] - 1;
1106                         boxmaxs[i] = end[i] + maxs[i] + 1;
1107                 }
1108                 else
1109                 {
1110                         boxmins[i] = end[i] + mins[i] - 1;
1111                         boxmaxs[i] = start[i] + maxs[i] + 1;
1112                 }
1113         }
1114 #endif
1115 }
1116
1117 /*
1118 ==================
1119 SV_Move
1120 ==================
1121 */
1122 trace_t SV_Move (vec3_t start, vec3_t mins, vec3_t maxs, vec3_t end, int type, edict_t *passedict)
1123 {
1124         moveclip_t      clip;
1125         int                     i;
1126
1127         memset ( &clip, 0, sizeof ( moveclip_t ) );
1128
1129 // clip to world
1130         clip.trace = SV_ClipMoveToEntity ( sv.edicts, start, mins, maxs, end );
1131
1132         clip.start = start;
1133         clip.end = end;
1134         clip.mins = mins;
1135         clip.maxs = maxs;
1136         clip.type = type;
1137         clip.passedict = passedict;
1138
1139         if (type == MOVE_MISSILE)
1140         {
1141                 for (i=0 ; i<3 ; i++)
1142                 {
1143                         clip.mins2[i] = -15;
1144                         clip.maxs2[i] = 15;
1145                 }
1146         }
1147         else
1148         {
1149                 VectorCopy (mins, clip.mins2);
1150                 VectorCopy (maxs, clip.maxs2);
1151         }
1152         
1153 // create the bounding box of the entire move
1154         SV_MoveBounds ( start, clip.mins2, clip.maxs2, end, clip.boxmins, clip.boxmaxs );
1155
1156 // clip to entities
1157         SV_ClipToLinks ( sv_areanodes, &clip );
1158
1159         return clip.trace;
1160 }