it is now possible to have multiple map models loaded at once - removed
[divverent/darkplaces.git] / model_brush.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
21 #include "quakedef.h"
22 #include "image.h"
23 #include "r_shadow.h"
24 #include "polygon.h"
25 #include "curves.h"
26 #include "wad.h"
27
28
29 //cvar_t r_subdivide_size = {CVAR_SAVE, "r_subdivide_size", "128", "how large water polygons should be (smaller values produce more polygons which give better warping effects)"};
30 cvar_t r_novis = {0, "r_novis", "0", "draws whole level, see also sv_cullentities_pvs 0"};
31 cvar_t r_picmipworld = {CVAR_SAVE, "r_picmipworld", "1", "whether gl_picmip shall apply to world textures too"};
32 cvar_t r_nosurftextures = {0, "r_nosurftextures", "0", "pretends there was no texture lump found in the q1bsp/hlbsp loading (useful for debugging this rare case)"};
33 cvar_t r_subdivisions_tolerance = {0, "r_subdivisions_tolerance", "4", "maximum error tolerance on curve subdivision for rendering purposes (in other words, the curves will be given as many polygons as necessary to represent curves at this quality)"};
34 cvar_t r_subdivisions_mintess = {0, "r_subdivisions_mintess", "0", "minimum number of subdivisions (values above 0 will smooth curves that don't need it)"};
35 cvar_t r_subdivisions_maxtess = {0, "r_subdivisions_maxtess", "1024", "maximum number of subdivisions (prevents curves beyond a certain detail level, limits smoothing)"};
36 cvar_t r_subdivisions_maxvertices = {0, "r_subdivisions_maxvertices", "65536", "maximum vertices allowed per subdivided curve"};
37 cvar_t r_subdivisions_collision_tolerance = {0, "r_subdivisions_collision_tolerance", "15", "maximum error tolerance on curve subdivision for collision purposes (usually a larger error tolerance than for rendering)"};
38 cvar_t r_subdivisions_collision_mintess = {0, "r_subdivisions_collision_mintess", "0", "minimum number of subdivisions (values above 0 will smooth curves that don't need it)"};
39 cvar_t r_subdivisions_collision_maxtess = {0, "r_subdivisions_collision_maxtess", "1024", "maximum number of subdivisions (prevents curves beyond a certain detail level, limits smoothing)"};
40 cvar_t r_subdivisions_collision_maxvertices = {0, "r_subdivisions_collision_maxvertices", "4225", "maximum vertices allowed per subdivided curve"};
41 cvar_t mod_q3bsp_curves_collisions = {0, "mod_q3bsp_curves_collisions", "1", "enables collisions with curves (SLOW)"};
42 cvar_t mod_q3bsp_optimizedtraceline = {0, "mod_q3bsp_optimizedtraceline", "1", "whether to use optimized traceline code for line traces (as opposed to tracebox code)"};
43 cvar_t mod_q3bsp_debugtracebrush = {0, "mod_q3bsp_debugtracebrush", "0", "selects different tracebrush bsp recursion algorithms (for debugging purposes only)"};
44 cvar_t mod_q3bsp_lightmapmergepower = {CVAR_SAVE, "mod_q3bsp_lightmapmergepower", "4", "merges the quake3 128x128 lightmap textures into larger lightmap group textures to speed up rendering, 1 = 256x256, 2 = 512x512, 3 = 1024x1024, 4 = 2048x2048, 5 = 4096x4096, ..."};
45 cvar_t mod_q3bsp_nolightmaps = {CVAR_SAVE, "mod_q3bsp_nolightmaps", "0", "do not load lightmaps in Q3BSP maps (to save video RAM, but be warned: it looks ugly)"};
46
47 static texture_t mod_q1bsp_texture_solid;
48 static texture_t mod_q1bsp_texture_sky;
49 static texture_t mod_q1bsp_texture_lava;
50 static texture_t mod_q1bsp_texture_slime;
51 static texture_t mod_q1bsp_texture_water;
52
53 void Mod_BrushInit(void)
54 {
55 //      Cvar_RegisterVariable(&r_subdivide_size);
56         Cvar_RegisterVariable(&r_novis);
57         Cvar_RegisterVariable(&r_picmipworld);
58         Cvar_RegisterVariable(&r_nosurftextures);
59         Cvar_RegisterVariable(&r_subdivisions_tolerance);
60         Cvar_RegisterVariable(&r_subdivisions_mintess);
61         Cvar_RegisterVariable(&r_subdivisions_maxtess);
62         Cvar_RegisterVariable(&r_subdivisions_maxvertices);
63         Cvar_RegisterVariable(&r_subdivisions_collision_tolerance);
64         Cvar_RegisterVariable(&r_subdivisions_collision_mintess);
65         Cvar_RegisterVariable(&r_subdivisions_collision_maxtess);
66         Cvar_RegisterVariable(&r_subdivisions_collision_maxvertices);
67         Cvar_RegisterVariable(&mod_q3bsp_curves_collisions);
68         Cvar_RegisterVariable(&mod_q3bsp_optimizedtraceline);
69         Cvar_RegisterVariable(&mod_q3bsp_debugtracebrush);
70         Cvar_RegisterVariable(&mod_q3bsp_lightmapmergepower);
71         Cvar_RegisterVariable(&mod_q3bsp_nolightmaps);
72
73         memset(&mod_q1bsp_texture_solid, 0, sizeof(mod_q1bsp_texture_solid));
74         strlcpy(mod_q1bsp_texture_solid.name, "solid" , sizeof(mod_q1bsp_texture_solid.name));
75         mod_q1bsp_texture_solid.surfaceflags = 0;
76         mod_q1bsp_texture_solid.supercontents = SUPERCONTENTS_SOLID;
77
78         mod_q1bsp_texture_sky = mod_q1bsp_texture_solid;
79         strlcpy(mod_q1bsp_texture_sky.name, "sky", sizeof(mod_q1bsp_texture_sky.name));
80         mod_q1bsp_texture_sky.surfaceflags = Q3SURFACEFLAG_SKY | Q3SURFACEFLAG_NOIMPACT | Q3SURFACEFLAG_NOMARKS | Q3SURFACEFLAG_NODLIGHT | Q3SURFACEFLAG_NOLIGHTMAP;
81         mod_q1bsp_texture_sky.supercontents = SUPERCONTENTS_SKY | SUPERCONTENTS_NODROP;
82
83         mod_q1bsp_texture_lava = mod_q1bsp_texture_solid;
84         strlcpy(mod_q1bsp_texture_lava.name, "*lava", sizeof(mod_q1bsp_texture_lava.name));
85         mod_q1bsp_texture_lava.surfaceflags = Q3SURFACEFLAG_NOMARKS;
86         mod_q1bsp_texture_lava.supercontents = SUPERCONTENTS_LAVA | SUPERCONTENTS_NODROP;
87
88         mod_q1bsp_texture_slime = mod_q1bsp_texture_solid;
89         strlcpy(mod_q1bsp_texture_slime.name, "*slime", sizeof(mod_q1bsp_texture_slime.name));
90         mod_q1bsp_texture_slime.surfaceflags = Q3SURFACEFLAG_NOMARKS;
91         mod_q1bsp_texture_slime.supercontents = SUPERCONTENTS_SLIME;
92
93         mod_q1bsp_texture_water = mod_q1bsp_texture_solid;
94         strlcpy(mod_q1bsp_texture_water.name, "*water", sizeof(mod_q1bsp_texture_water.name));
95         mod_q1bsp_texture_water.surfaceflags = Q3SURFACEFLAG_NOMARKS;
96         mod_q1bsp_texture_water.supercontents = SUPERCONTENTS_WATER;
97 }
98
99 static mleaf_t *Mod_Q1BSP_PointInLeaf(dp_model_t *model, const vec3_t p)
100 {
101         mnode_t *node;
102
103         if (model == NULL)
104                 return NULL;
105
106         // LordHavoc: modified to start at first clip node,
107         // in other words: first node of the (sub)model
108         node = model->brush.data_nodes + model->brushq1.hulls[0].firstclipnode;
109         while (node->plane)
110                 node = node->children[(node->plane->type < 3 ? p[node->plane->type] : DotProduct(p,node->plane->normal)) < node->plane->dist];
111
112         return (mleaf_t *)node;
113 }
114
115 static void Mod_Q1BSP_AmbientSoundLevelsForPoint(dp_model_t *model, const vec3_t p, unsigned char *out, int outsize)
116 {
117         int i;
118         mleaf_t *leaf;
119         leaf = Mod_Q1BSP_PointInLeaf(model, p);
120         if (leaf)
121         {
122                 i = min(outsize, (int)sizeof(leaf->ambient_sound_level));
123                 if (i)
124                 {
125                         memcpy(out, leaf->ambient_sound_level, i);
126                         out += i;
127                         outsize -= i;
128                 }
129         }
130         if (outsize)
131                 memset(out, 0, outsize);
132 }
133
134 static int Mod_Q1BSP_FindBoxClusters(dp_model_t *model, const vec3_t mins, const vec3_t maxs, int maxclusters, int *clusterlist)
135 {
136         int numclusters = 0;
137         int nodestackindex = 0;
138         mnode_t *node, *nodestack[1024];
139         if (!model->brush.num_pvsclusters)
140                 return -1;
141         node = model->brush.data_nodes;
142         for (;;)
143         {
144 #if 1
145                 if (node->plane)
146                 {
147                         // node - recurse down the BSP tree
148                         int sides = BoxOnPlaneSide(mins, maxs, node->plane);
149                         if (sides < 3)
150                         {
151                                 if (sides == 0)
152                                         return -1; // ERROR: NAN bounding box!
153                                 // box is on one side of plane, take that path
154                                 node = node->children[sides-1];
155                         }
156                         else
157                         {
158                                 // box crosses plane, take one path and remember the other
159                                 if (nodestackindex < 1024)
160                                         nodestack[nodestackindex++] = node->children[0];
161                                 node = node->children[1];
162                         }
163                         continue;
164                 }
165                 else
166                 {
167                         // leaf - add clusterindex to list
168                         if (numclusters < maxclusters)
169                                 clusterlist[numclusters] = ((mleaf_t *)node)->clusterindex;
170                         numclusters++;
171                 }
172 #else
173                 if (BoxesOverlap(mins, maxs, node->mins, node->maxs))
174                 {
175                         if (node->plane)
176                         {
177                                 if (nodestackindex < 1024)
178                                         nodestack[nodestackindex++] = node->children[0];
179                                 node = node->children[1];
180                                 continue;
181                         }
182                         else
183                         {
184                                 // leaf - add clusterindex to list
185                                 if (numclusters < maxclusters)
186                                         clusterlist[numclusters] = ((mleaf_t *)node)->clusterindex;
187                                 numclusters++;
188                         }
189                 }
190 #endif
191                 // try another path we didn't take earlier
192                 if (nodestackindex == 0)
193                         break;
194                 node = nodestack[--nodestackindex];
195         }
196         // return number of clusters found (even if more than the maxclusters)
197         return numclusters;
198 }
199
200 static int Mod_Q1BSP_BoxTouchingPVS(dp_model_t *model, const unsigned char *pvs, const vec3_t mins, const vec3_t maxs)
201 {
202         int nodestackindex = 0;
203         mnode_t *node, *nodestack[1024];
204         if (!model->brush.num_pvsclusters)
205                 return true;
206         node = model->brush.data_nodes;
207         for (;;)
208         {
209 #if 1
210                 if (node->plane)
211                 {
212                         // node - recurse down the BSP tree
213                         int sides = BoxOnPlaneSide(mins, maxs, node->plane);
214                         if (sides < 3)
215                         {
216                                 if (sides == 0)
217                                         return -1; // ERROR: NAN bounding box!
218                                 // box is on one side of plane, take that path
219                                 node = node->children[sides-1];
220                         }
221                         else
222                         {
223                                 // box crosses plane, take one path and remember the other
224                                 if (nodestackindex < 1024)
225                                         nodestack[nodestackindex++] = node->children[0];
226                                 node = node->children[1];
227                         }
228                         continue;
229                 }
230                 else
231                 {
232                         // leaf - check cluster bit
233                         int clusterindex = ((mleaf_t *)node)->clusterindex;
234                         if (CHECKPVSBIT(pvs, clusterindex))
235                         {
236                                 // it is visible, return immediately with the news
237                                 return true;
238                         }
239                 }
240 #else
241                 if (BoxesOverlap(mins, maxs, node->mins, node->maxs))
242                 {
243                         if (node->plane)
244                         {
245                                 if (nodestackindex < 1024)
246                                         nodestack[nodestackindex++] = node->children[0];
247                                 node = node->children[1];
248                                 continue;
249                         }
250                         else
251                         {
252                                 // leaf - check cluster bit
253                                 int clusterindex = ((mleaf_t *)node)->clusterindex;
254                                 if (CHECKPVSBIT(pvs, clusterindex))
255                                 {
256                                         // it is visible, return immediately with the news
257                                         return true;
258                                 }
259                         }
260                 }
261 #endif
262                 // nothing to see here, try another path we didn't take earlier
263                 if (nodestackindex == 0)
264                         break;
265                 node = nodestack[--nodestackindex];
266         }
267         // it is not visible
268         return false;
269 }
270
271 static int Mod_Q1BSP_BoxTouchingLeafPVS(dp_model_t *model, const unsigned char *pvs, const vec3_t mins, const vec3_t maxs)
272 {
273         int nodestackindex = 0;
274         mnode_t *node, *nodestack[1024];
275         if (!model->brush.num_leafs)
276                 return true;
277         node = model->brush.data_nodes;
278         for (;;)
279         {
280 #if 1
281                 if (node->plane)
282                 {
283                         // node - recurse down the BSP tree
284                         int sides = BoxOnPlaneSide(mins, maxs, node->plane);
285                         if (sides < 3)
286                         {
287                                 if (sides == 0)
288                                         return -1; // ERROR: NAN bounding box!
289                                 // box is on one side of plane, take that path
290                                 node = node->children[sides-1];
291                         }
292                         else
293                         {
294                                 // box crosses plane, take one path and remember the other
295                                 if (nodestackindex < 1024)
296                                         nodestack[nodestackindex++] = node->children[0];
297                                 node = node->children[1];
298                         }
299                         continue;
300                 }
301                 else
302                 {
303                         // leaf - check cluster bit
304                         int clusterindex = ((mleaf_t *)node) - model->brush.data_leafs;
305                         if (CHECKPVSBIT(pvs, clusterindex))
306                         {
307                                 // it is visible, return immediately with the news
308                                 return true;
309                         }
310                 }
311 #else
312                 if (BoxesOverlap(mins, maxs, node->mins, node->maxs))
313                 {
314                         if (node->plane)
315                         {
316                                 if (nodestackindex < 1024)
317                                         nodestack[nodestackindex++] = node->children[0];
318                                 node = node->children[1];
319                                 continue;
320                         }
321                         else
322                         {
323                                 // leaf - check cluster bit
324                                 int clusterindex = ((mleaf_t *)node) - model->brush.data_leafs;
325                                 if (CHECKPVSBIT(pvs, clusterindex))
326                                 {
327                                         // it is visible, return immediately with the news
328                                         return true;
329                                 }
330                         }
331                 }
332 #endif
333                 // nothing to see here, try another path we didn't take earlier
334                 if (nodestackindex == 0)
335                         break;
336                 node = nodestack[--nodestackindex];
337         }
338         // it is not visible
339         return false;
340 }
341
342 static int Mod_Q1BSP_BoxTouchingVisibleLeafs(dp_model_t *model, const unsigned char *visibleleafs, const vec3_t mins, const vec3_t maxs)
343 {
344         int nodestackindex = 0;
345         mnode_t *node, *nodestack[1024];
346         if (!model->brush.num_leafs)
347                 return true;
348         node = model->brush.data_nodes;
349         for (;;)
350         {
351 #if 1
352                 if (node->plane)
353                 {
354                         // node - recurse down the BSP tree
355                         int sides = BoxOnPlaneSide(mins, maxs, node->plane);
356                         if (sides < 3)
357                         {
358                                 if (sides == 0)
359                                         return -1; // ERROR: NAN bounding box!
360                                 // box is on one side of plane, take that path
361                                 node = node->children[sides-1];
362                         }
363                         else
364                         {
365                                 // box crosses plane, take one path and remember the other
366                                 if (nodestackindex < 1024)
367                                         nodestack[nodestackindex++] = node->children[0];
368                                 node = node->children[1];
369                         }
370                         continue;
371                 }
372                 else
373                 {
374                         // leaf - check if it is visible
375                         if (visibleleafs[(mleaf_t *)node - model->brush.data_leafs])
376                         {
377                                 // it is visible, return immediately with the news
378                                 return true;
379                         }
380                 }
381 #else
382                 if (BoxesOverlap(mins, maxs, node->mins, node->maxs))
383                 {
384                         if (node->plane)
385                         {
386                                 if (nodestackindex < 1024)
387                                         nodestack[nodestackindex++] = node->children[0];
388                                 node = node->children[1];
389                                 continue;
390                         }
391                         else
392                         {
393                                 // leaf - check if it is visible
394                                 if (visibleleafs[(mleaf_t *)node - model->brush.data_leafs])
395                                 {
396                                         // it is visible, return immediately with the news
397                                         return true;
398                                 }
399                         }
400                 }
401 #endif
402                 // nothing to see here, try another path we didn't take earlier
403                 if (nodestackindex == 0)
404                         break;
405                 node = nodestack[--nodestackindex];
406         }
407         // it is not visible
408         return false;
409 }
410
411 typedef struct findnonsolidlocationinfo_s
412 {
413         vec3_t center;
414         vec_t radius;
415         vec3_t nudge;
416         vec_t bestdist;
417         dp_model_t *model;
418 }
419 findnonsolidlocationinfo_t;
420
421 static void Mod_Q1BSP_FindNonSolidLocation_r_Leaf(findnonsolidlocationinfo_t *info, mleaf_t *leaf)
422 {
423         int i, surfacenum, k, *tri, *mark;
424         float dist, f, vert[3][3], edge[3][3], facenormal[3], edgenormal[3][3], point[3];
425         msurface_t *surface;
426         for (surfacenum = 0, mark = leaf->firstleafsurface;surfacenum < leaf->numleafsurfaces;surfacenum++, mark++)
427         {
428                 surface = info->model->data_surfaces + *mark;
429                 if (surface->texture->supercontents & SUPERCONTENTS_SOLID)
430                 {
431                         for (k = 0;k < surface->num_triangles;k++)
432                         {
433                                 tri = (info->model->surfmesh.data_element3i + 3 * surface->num_firsttriangle) + k * 3;
434                                 VectorCopy((info->model->surfmesh.data_vertex3f + tri[0] * 3), vert[0]);
435                                 VectorCopy((info->model->surfmesh.data_vertex3f + tri[1] * 3), vert[1]);
436                                 VectorCopy((info->model->surfmesh.data_vertex3f + tri[2] * 3), vert[2]);
437                                 VectorSubtract(vert[1], vert[0], edge[0]);
438                                 VectorSubtract(vert[2], vert[1], edge[1]);
439                                 CrossProduct(edge[1], edge[0], facenormal);
440                                 if (facenormal[0] || facenormal[1] || facenormal[2])
441                                 {
442                                         VectorNormalize(facenormal);
443                                         f = DotProduct(info->center, facenormal) - DotProduct(vert[0], facenormal);
444                                         if (f <= info->bestdist && f >= -info->bestdist)
445                                         {
446                                                 VectorSubtract(vert[0], vert[2], edge[2]);
447                                                 VectorNormalize(edge[0]);
448                                                 VectorNormalize(edge[1]);
449                                                 VectorNormalize(edge[2]);
450                                                 CrossProduct(facenormal, edge[0], edgenormal[0]);
451                                                 CrossProduct(facenormal, edge[1], edgenormal[1]);
452                                                 CrossProduct(facenormal, edge[2], edgenormal[2]);
453                                                 // face distance
454                                                 if (DotProduct(info->center, edgenormal[0]) < DotProduct(vert[0], edgenormal[0])
455                                                  && DotProduct(info->center, edgenormal[1]) < DotProduct(vert[1], edgenormal[1])
456                                                  && DotProduct(info->center, edgenormal[2]) < DotProduct(vert[2], edgenormal[2]))
457                                                 {
458                                                         // we got lucky, the center is within the face
459                                                         dist = DotProduct(info->center, facenormal) - DotProduct(vert[0], facenormal);
460                                                         if (dist < 0)
461                                                         {
462                                                                 dist = -dist;
463                                                                 if (info->bestdist > dist)
464                                                                 {
465                                                                         info->bestdist = dist;
466                                                                         VectorScale(facenormal, (info->radius - -dist), info->nudge);
467                                                                 }
468                                                         }
469                                                         else
470                                                         {
471                                                                 if (info->bestdist > dist)
472                                                                 {
473                                                                         info->bestdist = dist;
474                                                                         VectorScale(facenormal, (info->radius - dist), info->nudge);
475                                                                 }
476                                                         }
477                                                 }
478                                                 else
479                                                 {
480                                                         // check which edge or vertex the center is nearest
481                                                         for (i = 0;i < 3;i++)
482                                                         {
483                                                                 f = DotProduct(info->center, edge[i]);
484                                                                 if (f >= DotProduct(vert[0], edge[i])
485                                                                  && f <= DotProduct(vert[1], edge[i]))
486                                                                 {
487                                                                         // on edge
488                                                                         VectorMA(info->center, -f, edge[i], point);
489                                                                         dist = sqrt(DotProduct(point, point));
490                                                                         if (info->bestdist > dist)
491                                                                         {
492                                                                                 info->bestdist = dist;
493                                                                                 VectorScale(point, (info->radius / dist), info->nudge);
494                                                                         }
495                                                                         // skip both vertex checks
496                                                                         // (both are further away than this edge)
497                                                                         i++;
498                                                                 }
499                                                                 else
500                                                                 {
501                                                                         // not on edge, check first vertex of edge
502                                                                         VectorSubtract(info->center, vert[i], point);
503                                                                         dist = sqrt(DotProduct(point, point));
504                                                                         if (info->bestdist > dist)
505                                                                         {
506                                                                                 info->bestdist = dist;
507                                                                                 VectorScale(point, (info->radius / dist), info->nudge);
508                                                                         }
509                                                                 }
510                                                         }
511                                                 }
512                                         }
513                                 }
514                         }
515                 }
516         }
517 }
518
519 static void Mod_Q1BSP_FindNonSolidLocation_r(findnonsolidlocationinfo_t *info, mnode_t *node)
520 {
521         if (node->plane)
522         {
523                 float f = PlaneDiff(info->center, node->plane);
524                 if (f >= -info->bestdist)
525                         Mod_Q1BSP_FindNonSolidLocation_r(info, node->children[0]);
526                 if (f <= info->bestdist)
527                         Mod_Q1BSP_FindNonSolidLocation_r(info, node->children[1]);
528         }
529         else
530         {
531                 if (((mleaf_t *)node)->numleafsurfaces)
532                         Mod_Q1BSP_FindNonSolidLocation_r_Leaf(info, (mleaf_t *)node);
533         }
534 }
535
536 static void Mod_Q1BSP_FindNonSolidLocation(dp_model_t *model, const vec3_t in, vec3_t out, float radius)
537 {
538         int i;
539         findnonsolidlocationinfo_t info;
540         if (model == NULL)
541         {
542                 VectorCopy(in, out);
543                 return;
544         }
545         VectorCopy(in, info.center);
546         info.radius = radius;
547         info.model = model;
548         i = 0;
549         do
550         {
551                 VectorClear(info.nudge);
552                 info.bestdist = radius;
553                 Mod_Q1BSP_FindNonSolidLocation_r(&info, model->brush.data_nodes + model->brushq1.hulls[0].firstclipnode);
554                 VectorAdd(info.center, info.nudge, info.center);
555         }
556         while (info.bestdist < radius && ++i < 10);
557         VectorCopy(info.center, out);
558 }
559
560 int Mod_Q1BSP_SuperContentsFromNativeContents(dp_model_t *model, int nativecontents)
561 {
562         switch(nativecontents)
563         {
564                 case CONTENTS_EMPTY:
565                         return 0;
566                 case CONTENTS_SOLID:
567                         return SUPERCONTENTS_SOLID;
568                 case CONTENTS_WATER:
569                         return SUPERCONTENTS_WATER;
570                 case CONTENTS_SLIME:
571                         return SUPERCONTENTS_SLIME;
572                 case CONTENTS_LAVA:
573                         return SUPERCONTENTS_LAVA | SUPERCONTENTS_NODROP;
574                 case CONTENTS_SKY:
575                         return SUPERCONTENTS_SKY | SUPERCONTENTS_NODROP;
576         }
577         return 0;
578 }
579
580 int Mod_Q1BSP_NativeContentsFromSuperContents(dp_model_t *model, int supercontents)
581 {
582         if (supercontents & (SUPERCONTENTS_SOLID | SUPERCONTENTS_BODY))
583                 return CONTENTS_SOLID;
584         if (supercontents & SUPERCONTENTS_SKY)
585                 return CONTENTS_SKY;
586         if (supercontents & SUPERCONTENTS_LAVA)
587                 return CONTENTS_LAVA;
588         if (supercontents & SUPERCONTENTS_SLIME)
589                 return CONTENTS_SLIME;
590         if (supercontents & SUPERCONTENTS_WATER)
591                 return CONTENTS_WATER;
592         return CONTENTS_EMPTY;
593 }
594
595 typedef struct RecursiveHullCheckTraceInfo_s
596 {
597         // the hull we're tracing through
598         const hull_t *hull;
599
600         // the trace structure to fill in
601         trace_t *trace;
602
603         // start, end, and end - start (in model space)
604         double start[3];
605         double end[3];
606         double dist[3];
607 }
608 RecursiveHullCheckTraceInfo_t;
609
610 // 1/32 epsilon to keep floating point happy
611 #define DIST_EPSILON (0.03125)
612
613 #define HULLCHECKSTATE_EMPTY 0
614 #define HULLCHECKSTATE_SOLID 1
615 #define HULLCHECKSTATE_DONE 2
616
617 extern cvar_t collision_prefernudgedfraction;
618 static int Mod_Q1BSP_RecursiveHullCheck(RecursiveHullCheckTraceInfo_t *t, int num, double p1f, double p2f, double p1[3], double p2[3])
619 {
620         // status variables, these don't need to be saved on the stack when
621         // recursing...  but are because this should be thread-safe
622         // (note: tracing against a bbox is not thread-safe, yet)
623         int ret;
624         mplane_t *plane;
625         double t1, t2;
626
627         // variables that need to be stored on the stack when recursing
628         mclipnode_t *node;
629         int side;
630         double midf, mid[3];
631
632         // LordHavoc: a goto!  everyone flee in terror... :)
633 loc0:
634         // check for empty
635         if (num < 0)
636         {
637                 num = Mod_Q1BSP_SuperContentsFromNativeContents(NULL, num);
638                 if (!t->trace->startfound)
639                 {
640                         t->trace->startfound = true;
641                         t->trace->startsupercontents |= num;
642                 }
643                 if (num & SUPERCONTENTS_LIQUIDSMASK)
644                         t->trace->inwater = true;
645                 if (num == 0)
646                         t->trace->inopen = true;
647                 if (num & SUPERCONTENTS_SOLID)
648                         t->trace->hittexture = &mod_q1bsp_texture_solid;
649                 else if (num & SUPERCONTENTS_SKY)
650                         t->trace->hittexture = &mod_q1bsp_texture_sky;
651                 else if (num & SUPERCONTENTS_LAVA)
652                         t->trace->hittexture = &mod_q1bsp_texture_lava;
653                 else if (num & SUPERCONTENTS_SLIME)
654                         t->trace->hittexture = &mod_q1bsp_texture_slime;
655                 else
656                         t->trace->hittexture = &mod_q1bsp_texture_water;
657                 t->trace->hitq3surfaceflags = t->trace->hittexture->surfaceflags;
658                 t->trace->hitsupercontents = num;
659                 if (num & t->trace->hitsupercontentsmask)
660                 {
661                         // if the first leaf is solid, set startsolid
662                         if (t->trace->allsolid)
663                                 t->trace->startsolid = true;
664 #if COLLISIONPARANOID >= 3
665                         Con_Print("S");
666 #endif
667                         return HULLCHECKSTATE_SOLID;
668                 }
669                 else
670                 {
671                         t->trace->allsolid = false;
672 #if COLLISIONPARANOID >= 3
673                         Con_Print("E");
674 #endif
675                         return HULLCHECKSTATE_EMPTY;
676                 }
677         }
678
679         // find the point distances
680         node = t->hull->clipnodes + num;
681
682         plane = t->hull->planes + node->planenum;
683         if (plane->type < 3)
684         {
685                 t1 = p1[plane->type] - plane->dist;
686                 t2 = p2[plane->type] - plane->dist;
687         }
688         else
689         {
690                 t1 = DotProduct (plane->normal, p1) - plane->dist;
691                 t2 = DotProduct (plane->normal, p2) - plane->dist;
692         }
693
694         if (t1 < 0)
695         {
696                 if (t2 < 0)
697                 {
698 #if COLLISIONPARANOID >= 3
699                         Con_Print("<");
700 #endif
701                         num = node->children[1];
702                         goto loc0;
703                 }
704                 side = 1;
705         }
706         else
707         {
708                 if (t2 >= 0)
709                 {
710 #if COLLISIONPARANOID >= 3
711                         Con_Print(">");
712 #endif
713                         num = node->children[0];
714                         goto loc0;
715                 }
716                 side = 0;
717         }
718
719         // the line intersects, find intersection point
720         // LordHavoc: this uses the original trace for maximum accuracy
721 #if COLLISIONPARANOID >= 3
722         Con_Print("M");
723 #endif
724         if (plane->type < 3)
725         {
726                 t1 = t->start[plane->type] - plane->dist;
727                 t2 = t->end[plane->type] - plane->dist;
728         }
729         else
730         {
731                 t1 = DotProduct (plane->normal, t->start) - plane->dist;
732                 t2 = DotProduct (plane->normal, t->end) - plane->dist;
733         }
734
735         midf = t1 / (t1 - t2);
736         midf = bound(p1f, midf, p2f);
737         VectorMA(t->start, midf, t->dist, mid);
738
739         // recurse both sides, front side first
740         ret = Mod_Q1BSP_RecursiveHullCheck(t, node->children[side], p1f, midf, p1, mid);
741         // if this side is not empty, return what it is (solid or done)
742         if (ret != HULLCHECKSTATE_EMPTY)
743                 return ret;
744
745         ret = Mod_Q1BSP_RecursiveHullCheck(t, node->children[side ^ 1], midf, p2f, mid, p2);
746         // if other side is not solid, return what it is (empty or done)
747         if (ret != HULLCHECKSTATE_SOLID)
748                 return ret;
749
750         // front is air and back is solid, this is the impact point...
751         if (side)
752         {
753                 t->trace->plane.dist = -plane->dist;
754                 VectorNegate (plane->normal, t->trace->plane.normal);
755         }
756         else
757         {
758                 t->trace->plane.dist = plane->dist;
759                 VectorCopy (plane->normal, t->trace->plane.normal);
760         }
761
762         // calculate the true fraction
763         t1 = DotProduct(t->trace->plane.normal, t->start) - t->trace->plane.dist;
764         t2 = DotProduct(t->trace->plane.normal, t->end) - t->trace->plane.dist;
765         midf = t1 / (t1 - t2);
766         t->trace->realfraction = bound(0, midf, 1);
767
768         // calculate the return fraction which is nudged off the surface a bit
769         midf = (t1 - DIST_EPSILON) / (t1 - t2);
770         t->trace->fraction = bound(0, midf, 1);
771
772         if (collision_prefernudgedfraction.integer)
773                 t->trace->realfraction = t->trace->fraction;
774
775 #if COLLISIONPARANOID >= 3
776         Con_Print("D");
777 #endif
778         return HULLCHECKSTATE_DONE;
779 }
780
781 //#if COLLISIONPARANOID < 2
782 static int Mod_Q1BSP_RecursiveHullCheckPoint(RecursiveHullCheckTraceInfo_t *t, int num)
783 {
784         mplane_t *plane;
785         mclipnode_t *nodes = t->hull->clipnodes;
786         mplane_t *planes = t->hull->planes;
787         vec3_t point;
788         VectorCopy(t->start, point);
789         while (num >= 0)
790         {
791                 plane = planes + nodes[num].planenum;
792                 num = nodes[num].children[(plane->type < 3 ? point[plane->type] : DotProduct(plane->normal, point)) < plane->dist];
793         }
794         num = Mod_Q1BSP_SuperContentsFromNativeContents(NULL, num);
795         t->trace->startsupercontents |= num;
796         if (num & SUPERCONTENTS_LIQUIDSMASK)
797                 t->trace->inwater = true;
798         if (num == 0)
799                 t->trace->inopen = true;
800         if (num & t->trace->hitsupercontentsmask)
801         {
802                 t->trace->allsolid = t->trace->startsolid = true;
803                 return HULLCHECKSTATE_SOLID;
804         }
805         else
806         {
807                 t->trace->allsolid = t->trace->startsolid = false;
808                 return HULLCHECKSTATE_EMPTY;
809         }
810 }
811 //#endif
812
813 static void Mod_Q1BSP_TraceBox(struct model_s *model, int frame, trace_t *trace, const vec3_t start, const vec3_t boxmins, const vec3_t boxmaxs, const vec3_t end, int hitsupercontentsmask)
814 {
815         // this function currently only supports same size start and end
816         double boxsize[3];
817         RecursiveHullCheckTraceInfo_t rhc;
818
819         memset(&rhc, 0, sizeof(rhc));
820         memset(trace, 0, sizeof(trace_t));
821         rhc.trace = trace;
822         rhc.trace->hitsupercontentsmask = hitsupercontentsmask;
823         rhc.trace->fraction = 1;
824         rhc.trace->realfraction = 1;
825         rhc.trace->allsolid = true;
826         VectorSubtract(boxmaxs, boxmins, boxsize);
827         if (boxsize[0] < 3)
828                 rhc.hull = &model->brushq1.hulls[0]; // 0x0x0
829         else if (model->brush.ishlbsp)
830         {
831                 // LordHavoc: this has to have a minor tolerance (the .1) because of
832                 // minor float precision errors from the box being transformed around
833                 if (boxsize[0] < 32.1)
834                 {
835                         if (boxsize[2] < 54) // pick the nearest of 36 or 72
836                                 rhc.hull = &model->brushq1.hulls[3]; // 32x32x36
837                         else
838                                 rhc.hull = &model->brushq1.hulls[1]; // 32x32x72
839                 }
840                 else
841                         rhc.hull = &model->brushq1.hulls[2]; // 64x64x64
842         }
843         else
844         {
845                 // LordHavoc: this has to have a minor tolerance (the .1) because of
846                 // minor float precision errors from the box being transformed around
847                 if (boxsize[0] < 32.1)
848                         rhc.hull = &model->brushq1.hulls[1]; // 32x32x56
849                 else
850                         rhc.hull = &model->brushq1.hulls[2]; // 64x64x88
851         }
852         VectorMAMAM(1, start, 1, boxmins, -1, rhc.hull->clip_mins, rhc.start);
853         VectorMAMAM(1, end, 1, boxmins, -1, rhc.hull->clip_mins, rhc.end);
854         VectorSubtract(rhc.end, rhc.start, rhc.dist);
855 #if COLLISIONPARANOID >= 2
856         Con_Printf("t(%f %f %f,%f %f %f,%i %f %f %f)", rhc.start[0], rhc.start[1], rhc.start[2], rhc.end[0], rhc.end[1], rhc.end[2], rhc.hull - model->brushq1.hulls, rhc.hull->clip_mins[0], rhc.hull->clip_mins[1], rhc.hull->clip_mins[2]);
857         Mod_Q1BSP_RecursiveHullCheck(&rhc, rhc.hull->firstclipnode, 0, 1, rhc.start, rhc.end);
858         {
859
860                 double test[3];
861                 trace_t testtrace;
862                 VectorLerp(rhc.start, rhc.trace->fraction, rhc.end, test);
863                 memset(&testtrace, 0, sizeof(trace_t));
864                 rhc.trace = &testtrace;
865                 rhc.trace->hitsupercontentsmask = hitsupercontentsmask;
866                 rhc.trace->fraction = 1;
867                 rhc.trace->realfraction = 1;
868                 rhc.trace->allsolid = true;
869                 VectorCopy(test, rhc.start);
870                 VectorCopy(test, rhc.end);
871                 VectorClear(rhc.dist);
872                 Mod_Q1BSP_RecursiveHullCheckPoint(&rhc, rhc.hull->firstclipnode);
873                 //Mod_Q1BSP_RecursiveHullCheck(&rhc, rhc.hull->firstclipnode, 0, 1, test, test);
874                 if (!trace->startsolid && testtrace.startsolid)
875                         Con_Printf(" - ended in solid!\n");
876         }
877         Con_Print("\n");
878 #else
879         if (VectorLength2(rhc.dist))
880                 Mod_Q1BSP_RecursiveHullCheck(&rhc, rhc.hull->firstclipnode, 0, 1, rhc.start, rhc.end);
881         else
882                 Mod_Q1BSP_RecursiveHullCheckPoint(&rhc, rhc.hull->firstclipnode);
883 #endif
884 }
885
886 static int Mod_Q1BSP_PointSuperContents(struct model_s *model, int frame, const vec3_t point)
887 {
888         int num = model->brushq1.hulls[0].firstclipnode;
889         mplane_t *plane;
890         mclipnode_t *nodes = model->brushq1.hulls[0].clipnodes;
891         mplane_t *planes = model->brushq1.hulls[0].planes;
892         while (num >= 0)
893         {
894                 plane = planes + nodes[num].planenum;
895                 num = nodes[num].children[(plane->type < 3 ? point[plane->type] : DotProduct(plane->normal, point)) < plane->dist];
896         }
897         return Mod_Q1BSP_SuperContentsFromNativeContents(NULL, num);
898 }
899
900 void Collision_ClipTrace_Box(trace_t *trace, const vec3_t cmins, const vec3_t cmaxs, const vec3_t start, const vec3_t mins, const vec3_t maxs, const vec3_t end, int hitsupercontentsmask, int boxsupercontents, int boxq3surfaceflags, texture_t *boxtexture)
901 {
902 #if 1
903         colbrushf_t cbox;
904         colplanef_t cbox_planes[6];
905         cbox.supercontents = boxsupercontents;
906         cbox.numplanes = 6;
907         cbox.numpoints = 0;
908         cbox.numtriangles = 0;
909         cbox.planes = cbox_planes;
910         cbox.points = NULL;
911         cbox.elements = NULL;
912         cbox.markframe = 0;
913         cbox.mins[0] = 0;
914         cbox.mins[1] = 0;
915         cbox.mins[2] = 0;
916         cbox.maxs[0] = 0;
917         cbox.maxs[1] = 0;
918         cbox.maxs[2] = 0;
919         cbox_planes[0].normal[0] =  1;cbox_planes[0].normal[1] =  0;cbox_planes[0].normal[2] =  0;cbox_planes[0].dist = cmaxs[0] - mins[0];
920         cbox_planes[1].normal[0] = -1;cbox_planes[1].normal[1] =  0;cbox_planes[1].normal[2] =  0;cbox_planes[1].dist = maxs[0] - cmins[0];
921         cbox_planes[2].normal[0] =  0;cbox_planes[2].normal[1] =  1;cbox_planes[2].normal[2] =  0;cbox_planes[2].dist = cmaxs[1] - mins[1];
922         cbox_planes[3].normal[0] =  0;cbox_planes[3].normal[1] = -1;cbox_planes[3].normal[2] =  0;cbox_planes[3].dist = maxs[1] - cmins[1];
923         cbox_planes[4].normal[0] =  0;cbox_planes[4].normal[1] =  0;cbox_planes[4].normal[2] =  1;cbox_planes[4].dist = cmaxs[2] - mins[2];
924         cbox_planes[5].normal[0] =  0;cbox_planes[5].normal[1] =  0;cbox_planes[5].normal[2] = -1;cbox_planes[5].dist = maxs[2] - cmins[2];
925         cbox_planes[0].q3surfaceflags = boxq3surfaceflags;cbox_planes[0].texture = boxtexture;
926         cbox_planes[1].q3surfaceflags = boxq3surfaceflags;cbox_planes[1].texture = boxtexture;
927         cbox_planes[2].q3surfaceflags = boxq3surfaceflags;cbox_planes[2].texture = boxtexture;
928         cbox_planes[3].q3surfaceflags = boxq3surfaceflags;cbox_planes[3].texture = boxtexture;
929         cbox_planes[4].q3surfaceflags = boxq3surfaceflags;cbox_planes[4].texture = boxtexture;
930         cbox_planes[5].q3surfaceflags = boxq3surfaceflags;cbox_planes[5].texture = boxtexture;
931         memset(trace, 0, sizeof(trace_t));
932         trace->hitsupercontentsmask = hitsupercontentsmask;
933         trace->fraction = 1;
934         trace->realfraction = 1;
935         Collision_TraceLineBrushFloat(trace, start, end, &cbox, &cbox);
936 #else
937         RecursiveHullCheckTraceInfo_t rhc;
938         static hull_t box_hull;
939         static mclipnode_t box_clipnodes[6];
940         static mplane_t box_planes[6];
941         // fill in a default trace
942         memset(&rhc, 0, sizeof(rhc));
943         memset(trace, 0, sizeof(trace_t));
944         //To keep everything totally uniform, bounding boxes are turned into small
945         //BSP trees instead of being compared directly.
946         // create a temp hull from bounding box sizes
947         box_planes[0].dist = cmaxs[0] - mins[0];
948         box_planes[1].dist = cmins[0] - maxs[0];
949         box_planes[2].dist = cmaxs[1] - mins[1];
950         box_planes[3].dist = cmins[1] - maxs[1];
951         box_planes[4].dist = cmaxs[2] - mins[2];
952         box_planes[5].dist = cmins[2] - maxs[2];
953 #if COLLISIONPARANOID >= 3
954         Con_Printf("box_planes %f:%f %f:%f %f:%f\ncbox %f %f %f:%f %f %f\nbox %f %f %f:%f %f %f\n", box_planes[0].dist, box_planes[1].dist, box_planes[2].dist, box_planes[3].dist, box_planes[4].dist, box_planes[5].dist, cmins[0], cmins[1], cmins[2], cmaxs[0], cmaxs[1], cmaxs[2], mins[0], mins[1], mins[2], maxs[0], maxs[1], maxs[2]);
955 #endif
956
957         if (box_hull.clipnodes == NULL)
958         {
959                 int i, side;
960
961                 //Set up the planes and clipnodes so that the six floats of a bounding box
962                 //can just be stored out and get a proper hull_t structure.
963
964                 box_hull.clipnodes = box_clipnodes;
965                 box_hull.planes = box_planes;
966                 box_hull.firstclipnode = 0;
967                 box_hull.lastclipnode = 5;
968
969                 for (i = 0;i < 6;i++)
970                 {
971                         box_clipnodes[i].planenum = i;
972
973                         side = i&1;
974
975                         box_clipnodes[i].children[side] = CONTENTS_EMPTY;
976                         if (i != 5)
977                                 box_clipnodes[i].children[side^1] = i + 1;
978                         else
979                                 box_clipnodes[i].children[side^1] = CONTENTS_SOLID;
980
981                         box_planes[i].type = i>>1;
982                         box_planes[i].normal[i>>1] = 1;
983                 }
984         }
985
986         // trace a line through the generated clipping hull
987         //rhc.boxsupercontents = boxsupercontents;
988         rhc.hull = &box_hull;
989         rhc.trace = trace;
990         rhc.trace->hitsupercontentsmask = hitsupercontentsmask;
991         rhc.trace->fraction = 1;
992         rhc.trace->realfraction = 1;
993         rhc.trace->allsolid = true;
994         VectorCopy(start, rhc.start);
995         VectorCopy(end, rhc.end);
996         VectorSubtract(rhc.end, rhc.start, rhc.dist);
997         Mod_Q1BSP_RecursiveHullCheck(&rhc, rhc.hull->firstclipnode, 0, 1, rhc.start, rhc.end);
998         //VectorMA(rhc.start, rhc.trace->fraction, rhc.dist, rhc.trace->endpos);
999         if (rhc.trace->startsupercontents)
1000                 rhc.trace->startsupercontents = boxsupercontents;
1001 #endif
1002 }
1003
1004 static int Mod_Q1BSP_TraceLineOfSight_RecursiveNodeCheck(mnode_t *node, double p1[3], double p2[3])
1005 {
1006         double t1, t2;
1007         double midf, mid[3];
1008         int ret, side;
1009
1010         // check for empty
1011         while (node->plane)
1012         {
1013                 // find the point distances
1014                 mplane_t *plane = node->plane;
1015                 if (plane->type < 3)
1016                 {
1017                         t1 = p1[plane->type] - plane->dist;
1018                         t2 = p2[plane->type] - plane->dist;
1019                 }
1020                 else
1021                 {
1022                         t1 = DotProduct (plane->normal, p1) - plane->dist;
1023                         t2 = DotProduct (plane->normal, p2) - plane->dist;
1024                 }
1025
1026                 if (t1 < 0)
1027                 {
1028                         if (t2 < 0)
1029                         {
1030                                 node = node->children[1];
1031                                 continue;
1032                         }
1033                         side = 1;
1034                 }
1035                 else
1036                 {
1037                         if (t2 >= 0)
1038                         {
1039                                 node = node->children[0];
1040                                 continue;
1041                         }
1042                         side = 0;
1043                 }
1044
1045                 midf = t1 / (t1 - t2);
1046                 VectorLerp(p1, midf, p2, mid);
1047
1048                 // recurse both sides, front side first
1049                 // return 2 if empty is followed by solid (hit something)
1050                 // do not return 2 if both are solid or both empty,
1051                 // or if start is solid and end is empty
1052                 // as these degenerate cases usually indicate the eye is in solid and
1053                 // should see the target point anyway
1054                 ret = Mod_Q1BSP_TraceLineOfSight_RecursiveNodeCheck(node->children[side    ], p1, mid);
1055                 if (ret != 0)
1056                         return ret;
1057                 ret = Mod_Q1BSP_TraceLineOfSight_RecursiveNodeCheck(node->children[side ^ 1], mid, p2);
1058                 if (ret != 1)
1059                         return ret;
1060                 return 2;
1061         }
1062         return ((mleaf_t *)node)->clusterindex < 0;
1063 }
1064
1065 static qboolean Mod_Q1BSP_TraceLineOfSight(struct model_s *model, const vec3_t start, const vec3_t end)
1066 {
1067         // this function currently only supports same size start and end
1068         double tracestart[3], traceend[3];
1069         VectorCopy(start, tracestart);
1070         VectorCopy(end, traceend);
1071         return Mod_Q1BSP_TraceLineOfSight_RecursiveNodeCheck(model->brush.data_nodes, tracestart, traceend) != 2;
1072 }
1073
1074 static int Mod_Q1BSP_LightPoint_RecursiveBSPNode(dp_model_t *model, vec3_t ambientcolor, vec3_t diffusecolor, vec3_t diffusenormal, const mnode_t *node, float x, float y, float startz, float endz)
1075 {
1076         int side;
1077         float front, back;
1078         float mid, distz = endz - startz;
1079
1080 loc0:
1081         if (!node->plane)
1082                 return false;           // didn't hit anything
1083
1084         switch (node->plane->type)
1085         {
1086         case PLANE_X:
1087                 node = node->children[x < node->plane->dist];
1088                 goto loc0;
1089         case PLANE_Y:
1090                 node = node->children[y < node->plane->dist];
1091                 goto loc0;
1092         case PLANE_Z:
1093                 side = startz < node->plane->dist;
1094                 if ((endz < node->plane->dist) == side)
1095                 {
1096                         node = node->children[side];
1097                         goto loc0;
1098                 }
1099                 // found an intersection
1100                 mid = node->plane->dist;
1101                 break;
1102         default:
1103                 back = front = x * node->plane->normal[0] + y * node->plane->normal[1];
1104                 front += startz * node->plane->normal[2];
1105                 back += endz * node->plane->normal[2];
1106                 side = front < node->plane->dist;
1107                 if ((back < node->plane->dist) == side)
1108                 {
1109                         node = node->children[side];
1110                         goto loc0;
1111                 }
1112                 // found an intersection
1113                 mid = startz + distz * (front - node->plane->dist) / (front - back);
1114                 break;
1115         }
1116
1117         // go down front side
1118         if (node->children[side]->plane && Mod_Q1BSP_LightPoint_RecursiveBSPNode(model, ambientcolor, diffusecolor, diffusenormal, node->children[side], x, y, startz, mid))
1119                 return true;    // hit something
1120         else
1121         {
1122                 // check for impact on this node
1123                 if (node->numsurfaces)
1124                 {
1125                         int i, dsi, dti, lmwidth, lmheight;
1126                         float ds, dt;
1127                         msurface_t *surface;
1128                         unsigned char *lightmap;
1129                         int maps, line3, size3;
1130                         float dsfrac;
1131                         float dtfrac;
1132                         float scale, w, w00, w01, w10, w11;
1133
1134                         surface = model->data_surfaces + node->firstsurface;
1135                         for (i = 0;i < node->numsurfaces;i++, surface++)
1136                         {
1137                                 if (!(surface->texture->basematerialflags & MATERIALFLAG_WALL) || !surface->lightmapinfo->samples)
1138                                         continue;       // no lightmaps
1139
1140                                 // location we want to sample in the lightmap
1141                                 ds = ((x * surface->lightmapinfo->texinfo->vecs[0][0] + y * surface->lightmapinfo->texinfo->vecs[0][1] + mid * surface->lightmapinfo->texinfo->vecs[0][2] + surface->lightmapinfo->texinfo->vecs[0][3]) - surface->lightmapinfo->texturemins[0]) * 0.0625f;
1142                                 dt = ((x * surface->lightmapinfo->texinfo->vecs[1][0] + y * surface->lightmapinfo->texinfo->vecs[1][1] + mid * surface->lightmapinfo->texinfo->vecs[1][2] + surface->lightmapinfo->texinfo->vecs[1][3]) - surface->lightmapinfo->texturemins[1]) * 0.0625f;
1143
1144                                 // check the bounds
1145                                 dsi = (int)ds;
1146                                 dti = (int)dt;
1147                                 lmwidth = ((surface->lightmapinfo->extents[0]>>4)+1);
1148                                 lmheight = ((surface->lightmapinfo->extents[1]>>4)+1);
1149
1150                                 // is it in bounds?
1151                                 if (dsi >= 0 && dsi < lmwidth-1 && dti >= 0 && dti < lmheight-1)
1152                                 {
1153                                         // calculate bilinear interpolation factors
1154                                         // and also multiply by fixedpoint conversion factors
1155                                         dsfrac = ds - dsi;
1156                                         dtfrac = dt - dti;
1157                                         w00 = (1 - dsfrac) * (1 - dtfrac) * (1.0f / 32768.0f);
1158                                         w01 = (    dsfrac) * (1 - dtfrac) * (1.0f / 32768.0f);
1159                                         w10 = (1 - dsfrac) * (    dtfrac) * (1.0f / 32768.0f);
1160                                         w11 = (    dsfrac) * (    dtfrac) * (1.0f / 32768.0f);
1161
1162                                         // values for pointer math
1163                                         line3 = lmwidth * 3; // LordHavoc: *3 for colored lighting
1164                                         size3 = lmwidth * lmheight * 3; // LordHavoc: *3 for colored lighting
1165
1166                                         // look up the pixel
1167                                         lightmap = surface->lightmapinfo->samples + dti * line3 + dsi*3; // LordHavoc: *3 for colored lighting
1168
1169                                         // bilinear filter each lightmap style, and sum them
1170                                         for (maps = 0;maps < MAXLIGHTMAPS && surface->lightmapinfo->styles[maps] != 255;maps++)
1171                                         {
1172                                                 scale = r_refdef.scene.lightstylevalue[surface->lightmapinfo->styles[maps]];
1173                                                 w = w00 * scale;VectorMA(ambientcolor, w, lightmap            , ambientcolor);
1174                                                 w = w01 * scale;VectorMA(ambientcolor, w, lightmap + 3        , ambientcolor);
1175                                                 w = w10 * scale;VectorMA(ambientcolor, w, lightmap + line3    , ambientcolor);
1176                                                 w = w11 * scale;VectorMA(ambientcolor, w, lightmap + line3 + 3, ambientcolor);
1177                                                 lightmap += size3;
1178                                         }
1179
1180                                         return true; // success
1181                                 }
1182                         }
1183                 }
1184
1185                 // go down back side
1186                 node = node->children[side ^ 1];
1187                 startz = mid;
1188                 distz = endz - startz;
1189                 goto loc0;
1190         }
1191 }
1192
1193 void Mod_Q1BSP_LightPoint(dp_model_t *model, const vec3_t p, vec3_t ambientcolor, vec3_t diffusecolor, vec3_t diffusenormal)
1194 {
1195         // pretend lighting is coming down from above (due to lack of a lightgrid to know primary lighting direction)
1196         VectorSet(diffusenormal, 0, 0, 1);
1197
1198         if (!model->brushq1.lightdata)
1199         {
1200                 VectorSet(ambientcolor, 1, 1, 1);
1201                 VectorSet(diffusecolor, 0, 0, 0);
1202                 return;
1203         }
1204
1205         Mod_Q1BSP_LightPoint_RecursiveBSPNode(model, ambientcolor, diffusecolor, diffusenormal, model->brush.data_nodes + model->brushq1.hulls[0].firstclipnode, p[0], p[1], p[2] + 0.125, p[2] - 65536);
1206 }
1207
1208 static void Mod_Q1BSP_DecompressVis(const unsigned char *in, const unsigned char *inend, unsigned char *out, unsigned char *outend)
1209 {
1210         int c;
1211         unsigned char *outstart = out;
1212         while (out < outend)
1213         {
1214                 if (in == inend)
1215                 {
1216                         Con_Printf("Mod_Q1BSP_DecompressVis: input underrun on model \"%s\" (decompressed %i of %i output bytes)\n", loadmodel->name, (int)(out - outstart), (int)(outend - outstart));
1217                         return;
1218                 }
1219                 c = *in++;
1220                 if (c)
1221                         *out++ = c;
1222                 else
1223                 {
1224                         if (in == inend)
1225                         {
1226                                 Con_Printf("Mod_Q1BSP_DecompressVis: input underrun (during zero-run) on model \"%s\" (decompressed %i of %i output bytes)\n", loadmodel->name, (int)(out - outstart), (int)(outend - outstart));
1227                                 return;
1228                         }
1229                         for (c = *in++;c > 0;c--)
1230                         {
1231                                 if (out == outend)
1232                                 {
1233                                         Con_Printf("Mod_Q1BSP_DecompressVis: output overrun on model \"%s\" (decompressed %i of %i output bytes)\n", loadmodel->name, (int)(out - outstart), (int)(outend - outstart));
1234                                         return;
1235                                 }
1236                                 *out++ = 0;
1237                         }
1238                 }
1239         }
1240 }
1241
1242 /*
1243 =============
1244 R_Q1BSP_LoadSplitSky
1245
1246 A sky texture is 256*128, with the right side being a masked overlay
1247 ==============
1248 */
1249 void R_Q1BSP_LoadSplitSky (unsigned char *src, int width, int height, int bytesperpixel)
1250 {
1251         int i, j;
1252         unsigned solidpixels[128*128], alphapixels[128*128];
1253
1254         // allocate a texture pool if we need it
1255         if (loadmodel->texturepool == NULL && cls.state != ca_dedicated)
1256                 loadmodel->texturepool = R_AllocTexturePool();
1257
1258         if (bytesperpixel == 4)
1259         {
1260                 for (i = 0;i < 128;i++)
1261                 {
1262                         for (j = 0;j < 128;j++)
1263                         {
1264                                 solidpixels[(i*128) + j] = ((unsigned *)src)[i*256+j+128];
1265                                 alphapixels[(i*128) + j] = ((unsigned *)src)[i*256+j];
1266                         }
1267                 }
1268         }
1269         else
1270         {
1271                 // make an average value for the back to avoid
1272                 // a fringe on the top level
1273                 int p, r, g, b;
1274                 union
1275                 {
1276                         unsigned int i;
1277                         unsigned char b[4];
1278                 }
1279                 bgra;
1280                 r = g = b = 0;
1281                 for (i = 0;i < 128;i++)
1282                 {
1283                         for (j = 0;j < 128;j++)
1284                         {
1285                                 p = src[i*256 + j + 128];
1286                                 r += palette_rgb[p][0];
1287                                 g += palette_rgb[p][1];
1288                                 b += palette_rgb[p][2];
1289                         }
1290                 }
1291                 bgra.b[2] = r/(128*128);
1292                 bgra.b[1] = g/(128*128);
1293                 bgra.b[0] = b/(128*128);
1294                 bgra.b[3] = 0;
1295                 for (i = 0;i < 128;i++)
1296                 {
1297                         for (j = 0;j < 128;j++)
1298                         {
1299                                 solidpixels[(i*128) + j] = palette_bgra_complete[src[i*256 + j + 128]];
1300                                 p = src[i*256 + j];
1301                                 alphapixels[(i*128) + j] = p ? palette_bgra_complete[p] : bgra.i;
1302                         }
1303                 }
1304         }
1305
1306         loadmodel->brush.solidskytexture = R_LoadTexture2D(loadmodel->texturepool, "sky_solidtexture", 128, 128, (unsigned char *) solidpixels, TEXTYPE_BGRA, TEXF_PRECACHE, NULL);
1307         loadmodel->brush.alphaskytexture = R_LoadTexture2D(loadmodel->texturepool, "sky_alphatexture", 128, 128, (unsigned char *) alphapixels, TEXTYPE_BGRA, TEXF_ALPHA | TEXF_PRECACHE, NULL);
1308 }
1309
1310 static void Mod_Q1BSP_LoadTextures(lump_t *l)
1311 {
1312         int i, j, k, num, max, altmax, mtwidth, mtheight, *dofs, incomplete;
1313         skinframe_t *skinframe;
1314         miptex_t *dmiptex;
1315         texture_t *tx, *tx2, *anims[10], *altanims[10];
1316         dmiptexlump_t *m;
1317         unsigned char *data, *mtdata;
1318         const char *s;
1319         char mapname[MAX_QPATH], name[MAX_QPATH];
1320         unsigned char zero[4];
1321
1322         memset(zero, 0, sizeof(zero));
1323
1324         loadmodel->data_textures = NULL;
1325
1326         // add two slots for notexture walls and notexture liquids
1327         if (l->filelen)
1328         {
1329                 m = (dmiptexlump_t *)(mod_base + l->fileofs);
1330                 m->nummiptex = LittleLong (m->nummiptex);
1331                 loadmodel->num_textures = m->nummiptex + 2;
1332                 loadmodel->num_texturesperskin = loadmodel->num_textures;
1333         }
1334         else
1335         {
1336                 m = NULL;
1337                 loadmodel->num_textures = 2;
1338                 loadmodel->num_texturesperskin = loadmodel->num_textures;
1339         }
1340
1341         loadmodel->data_textures = (texture_t *)Mem_Alloc(loadmodel->mempool, loadmodel->num_textures * sizeof(texture_t));
1342
1343         // fill out all slots with notexture
1344         if (cls.state != ca_dedicated)
1345                 skinframe = R_SkinFrame_LoadMissing();
1346         else
1347                 skinframe = NULL;
1348         for (i = 0, tx = loadmodel->data_textures;i < loadmodel->num_textures;i++, tx++)
1349         {
1350                 strlcpy(tx->name, "NO TEXTURE FOUND", sizeof(tx->name));
1351                 tx->width = 16;
1352                 tx->height = 16;
1353                 if (cls.state != ca_dedicated)
1354                 {
1355                         tx->numskinframes = 1;
1356                         tx->skinframerate = 1;
1357                         tx->skinframes[0] = skinframe;
1358                         tx->currentskinframe = tx->skinframes[0];
1359                 }
1360                 tx->basematerialflags = MATERIALFLAG_WALL;
1361                 if (i == loadmodel->num_textures - 1)
1362                 {
1363                         tx->basematerialflags |= MATERIALFLAG_WATERSCROLL | MATERIALFLAG_LIGHTBOTHSIDES | MATERIALFLAG_NOSHADOW;
1364                         tx->supercontents = mod_q1bsp_texture_water.supercontents;
1365                         tx->surfaceflags = mod_q1bsp_texture_water.surfaceflags;
1366                 }
1367                 else
1368                 {
1369                         tx->supercontents = mod_q1bsp_texture_solid.supercontents;
1370                         tx->surfaceflags = mod_q1bsp_texture_solid.surfaceflags;
1371                 }
1372                 tx->currentframe = tx;
1373
1374                 // clear water settings
1375                 tx->reflectmin = 0;
1376                 tx->reflectmax = 1;
1377                 tx->refractfactor = 1;
1378                 Vector4Set(tx->refractcolor4f, 1, 1, 1, 1);
1379                 tx->reflectfactor = 1;
1380                 Vector4Set(tx->reflectcolor4f, 1, 1, 1, 1);
1381                 tx->r_water_wateralpha = 1;
1382         }
1383
1384         if (!m)
1385         {
1386                 Con_Printf("%s: no miptex lump to load textures from\n", loadmodel->name);
1387                 return;
1388         }
1389
1390         s = loadmodel->name;
1391         if (!strncasecmp(s, "maps/", 5))
1392                 s += 5;
1393         FS_StripExtension(s, mapname, sizeof(mapname));
1394
1395         // just to work around bounds checking when debugging with it (array index out of bounds error thing)
1396         dofs = m->dataofs;
1397         // LordHavoc: mostly rewritten map texture loader
1398         for (i = 0;i < m->nummiptex;i++)
1399         {
1400                 dofs[i] = LittleLong(dofs[i]);
1401                 if (r_nosurftextures.integer)
1402                         continue;
1403                 if (dofs[i] == -1)
1404                 {
1405                         Con_DPrintf("%s: miptex #%i missing\n", loadmodel->name, i);
1406                         continue;
1407                 }
1408                 dmiptex = (miptex_t *)((unsigned char *)m + dofs[i]);
1409
1410                 // copy name, but only up to 16 characters
1411                 // (the output buffer can hold more than this, but the input buffer is
1412                 //  only 16)
1413                 for (j = 0;j < 16 && dmiptex->name[j];j++)
1414                         name[j] = dmiptex->name[j];
1415                 name[j] = 0;
1416
1417                 if (!name[0])
1418                 {
1419                         dpsnprintf(name, sizeof(name), "unnamed%i", i);
1420                         Con_DPrintf("%s: warning: renaming unnamed texture to %s\n", loadmodel->name, name);
1421                 }
1422
1423                 mtwidth = LittleLong(dmiptex->width);
1424                 mtheight = LittleLong(dmiptex->height);
1425                 mtdata = NULL;
1426                 j = LittleLong(dmiptex->offsets[0]);
1427                 if (j)
1428                 {
1429                         // texture included
1430                         if (j < 40 || j + mtwidth * mtheight > l->filelen)
1431                         {
1432                                 Con_Printf("%s: Texture \"%s\" is corrupt or incomplete\n", loadmodel->name, dmiptex->name);
1433                                 continue;
1434                         }
1435                         mtdata = (unsigned char *)dmiptex + j;
1436                 }
1437
1438                 if ((mtwidth & 15) || (mtheight & 15))
1439                         Con_DPrintf("%s: warning: texture \"%s\" is not 16 aligned\n", loadmodel->name, dmiptex->name);
1440
1441                 // LordHavoc: force all names to lowercase
1442                 for (j = 0;name[j];j++)
1443                         if (name[j] >= 'A' && name[j] <= 'Z')
1444                                 name[j] += 'a' - 'A';
1445
1446                 if (dmiptex->name[0] && Mod_LoadTextureFromQ3Shader(loadmodel->data_textures + i, name, false, false, 0))
1447                         continue;
1448
1449                 tx = loadmodel->data_textures + i;
1450                 strlcpy(tx->name, name, sizeof(tx->name));
1451                 tx->width = mtwidth;
1452                 tx->height = mtheight;
1453
1454                 if (tx->name[0] == '*')
1455                 {
1456                         if (!strncmp(tx->name, "*lava", 5))
1457                         {
1458                                 tx->supercontents = mod_q1bsp_texture_lava.supercontents;
1459                                 tx->surfaceflags = mod_q1bsp_texture_lava.surfaceflags;
1460                         }
1461                         else if (!strncmp(tx->name, "*slime", 6))
1462                         {
1463                                 tx->supercontents = mod_q1bsp_texture_slime.supercontents;
1464                                 tx->surfaceflags = mod_q1bsp_texture_slime.surfaceflags;
1465                         }
1466                         else
1467                         {
1468                                 tx->supercontents = mod_q1bsp_texture_water.supercontents;
1469                                 tx->surfaceflags = mod_q1bsp_texture_water.surfaceflags;
1470                         }
1471                 }
1472                 else if (!strncmp(tx->name, "sky", 3))
1473                 {
1474                         tx->supercontents = mod_q1bsp_texture_sky.supercontents;
1475                         tx->surfaceflags = mod_q1bsp_texture_sky.surfaceflags;
1476                 }
1477                 else
1478                 {
1479                         tx->supercontents = mod_q1bsp_texture_solid.supercontents;
1480                         tx->surfaceflags = mod_q1bsp_texture_solid.surfaceflags;
1481                 }
1482
1483                 if (cls.state != ca_dedicated)
1484                 {
1485                         // LordHavoc: HL sky textures are entirely different than quake
1486                         if (!loadmodel->brush.ishlbsp && !strncmp(tx->name, "sky", 3) && mtwidth == 256 && mtheight == 128)
1487                         {
1488                                 data = loadimagepixelsbgra(tx->name, false, false);
1489                                 if (data && image_width == 256 && image_height == 128)
1490                                 {
1491                                         R_Q1BSP_LoadSplitSky(data, image_width, image_height, 4);
1492                                         Mem_Free(data);
1493                                 }
1494                                 else if (mtdata != NULL)
1495                                         R_Q1BSP_LoadSplitSky(mtdata, mtwidth, mtheight, 1);
1496                         }
1497                         else
1498                         {
1499                                 skinframe = R_SkinFrame_LoadExternal(gamemode == GAME_TENEBRAE ? tx->name : va("textures/%s/%s", mapname, tx->name), TEXF_ALPHA | TEXF_MIPMAP | TEXF_PRECACHE | (r_picmipworld.integer ? TEXF_PICMIP : 0) | TEXF_COMPRESS, false);
1500                                 if (!skinframe)
1501                                         skinframe = R_SkinFrame_LoadExternal(gamemode == GAME_TENEBRAE ? tx->name : va("textures/%s", tx->name), TEXF_ALPHA | TEXF_MIPMAP | TEXF_PRECACHE | (r_picmipworld.integer ? TEXF_PICMIP : 0) | TEXF_COMPRESS, false);
1502                                 if (!skinframe)
1503                                 {
1504                                         // did not find external texture, load it from the bsp or wad3
1505                                         if (loadmodel->brush.ishlbsp)
1506                                         {
1507                                                 // internal texture overrides wad
1508                                                 unsigned char *pixels, *freepixels;
1509                                                 pixels = freepixels = NULL;
1510                                                 if (mtdata)
1511                                                         pixels = W_ConvertWAD3TextureBGRA(dmiptex);
1512                                                 if (pixels == NULL)
1513                                                         pixels = freepixels = W_GetTextureBGRA(tx->name);
1514                                                 if (pixels != NULL)
1515                                                 {
1516                                                         tx->width = image_width;
1517                                                         tx->height = image_height;
1518                                                         skinframe = R_SkinFrame_LoadInternalBGRA(tx->name, TEXF_ALPHA | TEXF_MIPMAP | TEXF_PRECACHE | (r_picmipworld.integer ? TEXF_PICMIP : 0), pixels, image_width, image_height);
1519                                                 }
1520                                                 if (freepixels)
1521                                                         Mem_Free(freepixels);
1522                                         }
1523                                         else if (mtdata) // texture included
1524                                                 skinframe = R_SkinFrame_LoadInternalQuake(tx->name, TEXF_ALPHA | TEXF_MIPMAP | TEXF_PRECACHE | (r_picmipworld.integer ? TEXF_PICMIP : 0), false, r_fullbrights.integer, mtdata, tx->width, tx->height);
1525                                 }
1526                                 // if skinframe is still NULL the "missing" texture will be used
1527                                 if (skinframe)
1528                                         tx->skinframes[0] = skinframe;
1529                         }
1530
1531                         tx->basematerialflags = MATERIALFLAG_WALL;
1532                         if (tx->name[0] == '*')
1533                         {
1534                                 // LordHavoc: some turbulent textures should not be affected by wateralpha
1535                                 if (!strncmp(tx->name, "*glassmirror", 12)) // Tenebrae
1536                                 {
1537                                         // replace the texture with transparent black
1538                                         Vector4Set(zero, 128, 128, 128, 128);
1539                                         tx->skinframes[0] = R_SkinFrame_LoadInternalBGRA(tx->name, TEXF_MIPMAP | TEXF_PRECACHE | TEXF_ALPHA, zero, 1, 1);
1540                                         tx->basematerialflags |= MATERIALFLAG_NOSHADOW | MATERIALFLAG_ADD | MATERIALFLAG_BLENDED | MATERIALFLAG_REFLECTION;
1541                                 }
1542                                 else if (!strncmp(tx->name,"*lava",5)
1543                                  || !strncmp(tx->name,"*teleport",9)
1544                                  || !strncmp(tx->name,"*rift",5)) // Scourge of Armagon texture
1545                                         tx->basematerialflags |= MATERIALFLAG_WATERSCROLL | MATERIALFLAG_LIGHTBOTHSIDES | MATERIALFLAG_NOSHADOW;
1546                                 else
1547                                         tx->basematerialflags |= MATERIALFLAG_WATERSCROLL | MATERIALFLAG_LIGHTBOTHSIDES | MATERIALFLAG_NOSHADOW | MATERIALFLAG_WATERALPHA | MATERIALFLAG_WATERSHADER;
1548                                 if (tx->skinframes[0] && tx->skinframes[0]->fog)
1549                                         tx->basematerialflags |= MATERIALFLAG_ALPHA | MATERIALFLAG_BLENDED | MATERIALFLAG_NOSHADOW;
1550                         }
1551                         else if (!strncmp(tx->name, "mirror", 6)) // Tenebrae
1552                         {
1553                                 // replace the texture with black
1554                                 tx->skinframes[0] = R_SkinFrame_LoadInternalBGRA(tx->name, TEXF_PRECACHE, zero, 1, 1);
1555                                 tx->basematerialflags |= MATERIALFLAG_REFLECTION;
1556                         }
1557                         else if (!strncmp(tx->name, "sky", 3))
1558                                 tx->basematerialflags = MATERIALFLAG_SKY | MATERIALFLAG_NOSHADOW;
1559                         else if (!strcmp(tx->name, "caulk"))
1560                                 tx->basematerialflags = MATERIALFLAG_NODRAW | MATERIALFLAG_NOSHADOW;
1561                         else if (tx->skinframes[0] && tx->skinframes[0]->fog)
1562                                 tx->basematerialflags |= MATERIALFLAG_ALPHA | MATERIALFLAG_BLENDED | MATERIALFLAG_NOSHADOW;
1563
1564                         // start out with no animation
1565                         tx->currentframe = tx;
1566                         tx->currentskinframe = tx->skinframes[0];
1567                 }
1568         }
1569
1570         // sequence the animations
1571         for (i = 0;i < m->nummiptex;i++)
1572         {
1573                 tx = loadmodel->data_textures + i;
1574                 if (!tx || tx->name[0] != '+' || tx->name[1] == 0 || tx->name[2] == 0)
1575                         continue;
1576                 if (tx->anim_total[0] || tx->anim_total[1])
1577                         continue;       // already sequenced
1578
1579                 // find the number of frames in the animation
1580                 memset(anims, 0, sizeof(anims));
1581                 memset(altanims, 0, sizeof(altanims));
1582
1583                 for (j = i;j < m->nummiptex;j++)
1584                 {
1585                         tx2 = loadmodel->data_textures + j;
1586                         if (!tx2 || tx2->name[0] != '+' || strcmp(tx2->name+2, tx->name+2))
1587                                 continue;
1588
1589                         num = tx2->name[1];
1590                         if (num >= '0' && num <= '9')
1591                                 anims[num - '0'] = tx2;
1592                         else if (num >= 'a' && num <= 'j')
1593                                 altanims[num - 'a'] = tx2;
1594                         else
1595                                 Con_Printf("Bad animating texture %s\n", tx->name);
1596                 }
1597
1598                 max = altmax = 0;
1599                 for (j = 0;j < 10;j++)
1600                 {
1601                         if (anims[j])
1602                                 max = j + 1;
1603                         if (altanims[j])
1604                                 altmax = j + 1;
1605                 }
1606                 //Con_Printf("linking animation %s (%i:%i frames)\n\n", tx->name, max, altmax);
1607
1608                 incomplete = false;
1609                 for (j = 0;j < max;j++)
1610                 {
1611                         if (!anims[j])
1612                         {
1613                                 Con_Printf("Missing frame %i of %s\n", j, tx->name);
1614                                 incomplete = true;
1615                         }
1616                 }
1617                 for (j = 0;j < altmax;j++)
1618                 {
1619                         if (!altanims[j])
1620                         {
1621                                 Con_Printf("Missing altframe %i of %s\n", j, tx->name);
1622                                 incomplete = true;
1623                         }
1624                 }
1625                 if (incomplete)
1626                         continue;
1627
1628                 if (altmax < 1)
1629                 {
1630                         // if there is no alternate animation, duplicate the primary
1631                         // animation into the alternate
1632                         altmax = max;
1633                         for (k = 0;k < 10;k++)
1634                                 altanims[k] = anims[k];
1635                 }
1636
1637                 // link together the primary animation
1638                 for (j = 0;j < max;j++)
1639                 {
1640                         tx2 = anims[j];
1641                         tx2->animated = true;
1642                         tx2->anim_total[0] = max;
1643                         tx2->anim_total[1] = altmax;
1644                         for (k = 0;k < 10;k++)
1645                         {
1646                                 tx2->anim_frames[0][k] = anims[k];
1647                                 tx2->anim_frames[1][k] = altanims[k];
1648                         }
1649                 }
1650
1651                 // if there really is an alternate anim...
1652                 if (anims[0] != altanims[0])
1653                 {
1654                         // link together the alternate animation
1655                         for (j = 0;j < altmax;j++)
1656                         {
1657                                 tx2 = altanims[j];
1658                                 tx2->animated = true;
1659                                 // the primary/alternate are reversed here
1660                                 tx2->anim_total[0] = altmax;
1661                                 tx2->anim_total[1] = max;
1662                                 for (k = 0;k < 10;k++)
1663                                 {
1664                                         tx2->anim_frames[0][k] = altanims[k];
1665                                         tx2->anim_frames[1][k] = anims[k];
1666                                 }
1667                         }
1668                 }
1669         }
1670 }
1671
1672 static void Mod_Q1BSP_LoadLighting(lump_t *l)
1673 {
1674         int i;
1675         unsigned char *in, *out, *data, d;
1676         char litfilename[MAX_QPATH];
1677         char dlitfilename[MAX_QPATH];
1678         fs_offset_t filesize;
1679         if (loadmodel->brush.ishlbsp) // LordHavoc: load the colored lighting data straight
1680         {
1681                 loadmodel->brushq1.lightdata = (unsigned char *)Mem_Alloc(loadmodel->mempool, l->filelen);
1682                 for (i=0; i<l->filelen; i++)
1683                         loadmodel->brushq1.lightdata[i] = mod_base[l->fileofs+i] >>= 1;
1684         }
1685         else // LordHavoc: bsp version 29 (normal white lighting)
1686         {
1687                 // LordHavoc: hope is not lost yet, check for a .lit file to load
1688                 strlcpy (litfilename, loadmodel->name, sizeof (litfilename));
1689                 FS_StripExtension (litfilename, litfilename, sizeof (litfilename));
1690                 strlcpy (dlitfilename, litfilename, sizeof (dlitfilename));
1691                 strlcat (litfilename, ".lit", sizeof (litfilename));
1692                 strlcat (dlitfilename, ".dlit", sizeof (dlitfilename));
1693                 data = (unsigned char*) FS_LoadFile(litfilename, tempmempool, false, &filesize);
1694                 if (data)
1695                 {
1696                         if (filesize == (fs_offset_t)(8 + l->filelen * 3) && data[0] == 'Q' && data[1] == 'L' && data[2] == 'I' && data[3] == 'T')
1697                         {
1698                                 i = LittleLong(((int *)data)[1]);
1699                                 if (i == 1)
1700                                 {
1701                                         if (developer_loading.integer)
1702                                                 Con_Printf("loaded %s\n", litfilename);
1703                                         loadmodel->brushq1.lightdata = (unsigned char *)Mem_Alloc(loadmodel->mempool, filesize - 8);
1704                                         memcpy(loadmodel->brushq1.lightdata, data + 8, filesize - 8);
1705                                         Mem_Free(data);
1706                                         data = (unsigned char*) FS_LoadFile(dlitfilename, tempmempool, false, &filesize);
1707                                         if (data)
1708                                         {
1709                                                 if (filesize == (fs_offset_t)(8 + l->filelen * 3) && data[0] == 'Q' && data[1] == 'L' && data[2] == 'I' && data[3] == 'T')
1710                                                 {
1711                                                         i = LittleLong(((int *)data)[1]);
1712                                                         if (i == 1)
1713                                                         {
1714                                                                 if (developer_loading.integer)
1715                                                                         Con_Printf("loaded %s\n", dlitfilename);
1716                                                                 loadmodel->brushq1.nmaplightdata = (unsigned char *)Mem_Alloc(loadmodel->mempool, filesize - 8);
1717                                                                 memcpy(loadmodel->brushq1.nmaplightdata, data + 8, filesize - 8);
1718                                                                 loadmodel->brushq3.deluxemapping_modelspace = false;
1719                                                                 loadmodel->brushq3.deluxemapping = true;
1720                                                         }
1721                                                 }
1722                                                 Mem_Free(data);
1723                                                 data = NULL;
1724                                         }
1725                                         return;
1726                                 }
1727                                 else
1728                                         Con_Printf("Unknown .lit file version (%d)\n", i);
1729                         }
1730                         else if (filesize == 8)
1731                                 Con_Print("Empty .lit file, ignoring\n");
1732                         else
1733                                 Con_Printf("Corrupt .lit file (file size %i bytes, should be %i bytes), ignoring\n", (int) filesize, (int) (8 + l->filelen * 3));
1734                         if (data)
1735                         {
1736                                 Mem_Free(data);
1737                                 data = NULL;
1738                         }
1739                 }
1740                 // LordHavoc: oh well, expand the white lighting data
1741                 if (!l->filelen)
1742                         return;
1743                 loadmodel->brushq1.lightdata = (unsigned char *)Mem_Alloc(loadmodel->mempool, l->filelen*3);
1744                 in = mod_base + l->fileofs;
1745                 out = loadmodel->brushq1.lightdata;
1746                 for (i = 0;i < l->filelen;i++)
1747                 {
1748                         d = *in++;
1749                         *out++ = d;
1750                         *out++ = d;
1751                         *out++ = d;
1752                 }
1753         }
1754 }
1755
1756 static void Mod_Q1BSP_LoadVisibility(lump_t *l)
1757 {
1758         loadmodel->brushq1.num_compressedpvs = 0;
1759         loadmodel->brushq1.data_compressedpvs = NULL;
1760         if (!l->filelen)
1761                 return;
1762         loadmodel->brushq1.num_compressedpvs = l->filelen;
1763         loadmodel->brushq1.data_compressedpvs = (unsigned char *)Mem_Alloc(loadmodel->mempool, l->filelen);
1764         memcpy(loadmodel->brushq1.data_compressedpvs, mod_base + l->fileofs, l->filelen);
1765 }
1766
1767 // used only for HalfLife maps
1768 static void Mod_Q1BSP_ParseWadsFromEntityLump(const char *data)
1769 {
1770         char key[128], value[4096];
1771         int i, j, k;
1772         if (!data)
1773                 return;
1774         if (!COM_ParseToken_Simple(&data, false, false))
1775                 return; // error
1776         if (com_token[0] != '{')
1777                 return; // error
1778         while (1)
1779         {
1780                 if (!COM_ParseToken_Simple(&data, false, false))
1781                         return; // error
1782                 if (com_token[0] == '}')
1783                         break; // end of worldspawn
1784                 if (com_token[0] == '_')
1785                         strlcpy(key, com_token + 1, sizeof(key));
1786                 else
1787                         strlcpy(key, com_token, sizeof(key));
1788                 while (key[strlen(key)-1] == ' ') // remove trailing spaces
1789                         key[strlen(key)-1] = 0;
1790                 if (!COM_ParseToken_Simple(&data, false, false))
1791                         return; // error
1792                 dpsnprintf(value, sizeof(value), "%s", com_token);
1793                 if (!strcmp("wad", key)) // for HalfLife maps
1794                 {
1795                         if (loadmodel->brush.ishlbsp)
1796                         {
1797                                 j = 0;
1798                                 for (i = 0;i < (int)sizeof(value);i++)
1799                                         if (value[i] != ';' && value[i] != '\\' && value[i] != '/' && value[i] != ':')
1800                                                 break;
1801                                 if (value[i])
1802                                 {
1803                                         for (;i < (int)sizeof(value);i++)
1804                                         {
1805                                                 // ignore path - the \\ check is for HalfLife... stupid windoze 'programmers'...
1806                                                 if (value[i] == '\\' || value[i] == '/' || value[i] == ':')
1807                                                         j = i+1;
1808                                                 else if (value[i] == ';' || value[i] == 0)
1809                                                 {
1810                                                         k = value[i];
1811                                                         value[i] = 0;
1812                                                         W_LoadTextureWadFile(&value[j], false);
1813                                                         j = i+1;
1814                                                         if (!k)
1815                                                                 break;
1816                                                 }
1817                                         }
1818                                 }
1819                         }
1820                 }
1821         }
1822 }
1823
1824 static void Mod_Q1BSP_LoadEntities(lump_t *l)
1825 {
1826         loadmodel->brush.entities = NULL;
1827         if (!l->filelen)
1828                 return;
1829         loadmodel->brush.entities = (char *)Mem_Alloc(loadmodel->mempool, l->filelen + 1);
1830         memcpy(loadmodel->brush.entities, mod_base + l->fileofs, l->filelen);
1831         loadmodel->brush.entities[l->filelen] = 0;
1832         if (loadmodel->brush.ishlbsp)
1833                 Mod_Q1BSP_ParseWadsFromEntityLump(loadmodel->brush.entities);
1834 }
1835
1836
1837 static void Mod_Q1BSP_LoadVertexes(lump_t *l)
1838 {
1839         dvertex_t       *in;
1840         mvertex_t       *out;
1841         int                     i, count;
1842
1843         in = (dvertex_t *)(mod_base + l->fileofs);
1844         if (l->filelen % sizeof(*in))
1845                 Host_Error("Mod_Q1BSP_LoadVertexes: funny lump size in %s",loadmodel->name);
1846         count = l->filelen / sizeof(*in);
1847         out = (mvertex_t *)Mem_Alloc(loadmodel->mempool, count*sizeof(*out));
1848
1849         loadmodel->brushq1.vertexes = out;
1850         loadmodel->brushq1.numvertexes = count;
1851
1852         for ( i=0 ; i<count ; i++, in++, out++)
1853         {
1854                 out->position[0] = LittleFloat(in->point[0]);
1855                 out->position[1] = LittleFloat(in->point[1]);
1856                 out->position[2] = LittleFloat(in->point[2]);
1857         }
1858 }
1859
1860 // The following two functions should be removed and MSG_* or SZ_* function sets adjusted so they
1861 // can be used for this
1862 // REMOVEME
1863 int SB_ReadInt (unsigned char **buffer)
1864 {
1865         int     i;
1866         i = ((*buffer)[0]) + 256*((*buffer)[1]) + 65536*((*buffer)[2]) + 16777216*((*buffer)[3]);
1867         (*buffer) += 4;
1868         return i;
1869 }
1870
1871 // REMOVEME
1872 float SB_ReadFloat (unsigned char **buffer)
1873 {
1874         union
1875         {
1876                 int             i;
1877                 float   f;
1878         } u;
1879
1880         u.i = SB_ReadInt (buffer);
1881         return u.f;
1882 }
1883
1884 static void Mod_Q1BSP_LoadSubmodels(lump_t *l, hullinfo_t *hullinfo)
1885 {
1886         unsigned char           *index;
1887         dmodel_t        *out;
1888         int                     i, j, count;
1889
1890         index = (unsigned char *)(mod_base + l->fileofs);
1891         if (l->filelen % (48+4*hullinfo->filehulls))
1892                 Host_Error ("Mod_Q1BSP_LoadSubmodels: funny lump size in %s", loadmodel->name);
1893
1894         count = l->filelen / (48+4*hullinfo->filehulls);
1895         out = (dmodel_t *)Mem_Alloc (loadmodel->mempool, count*sizeof(*out));
1896
1897         loadmodel->brushq1.submodels = out;
1898         loadmodel->brush.numsubmodels = count;
1899
1900         for (i = 0; i < count; i++, out++)
1901         {
1902         // spread out the mins / maxs by a pixel
1903                 out->mins[0] = SB_ReadFloat (&index) - 1;
1904                 out->mins[1] = SB_ReadFloat (&index) - 1;
1905                 out->mins[2] = SB_ReadFloat (&index) - 1;
1906                 out->maxs[0] = SB_ReadFloat (&index) + 1;
1907                 out->maxs[1] = SB_ReadFloat (&index) + 1;
1908                 out->maxs[2] = SB_ReadFloat (&index) + 1;
1909                 out->origin[0] = SB_ReadFloat (&index);
1910                 out->origin[1] = SB_ReadFloat (&index);
1911                 out->origin[2] = SB_ReadFloat (&index);
1912                 for (j = 0; j < hullinfo->filehulls; j++)
1913                         out->headnode[j] = SB_ReadInt (&index);
1914                 out->visleafs = SB_ReadInt (&index);
1915                 out->firstface = SB_ReadInt (&index);
1916                 out->numfaces = SB_ReadInt (&index);
1917         }
1918 }
1919
1920 static void Mod_Q1BSP_LoadEdges(lump_t *l)
1921 {
1922         dedge_t *in;
1923         medge_t *out;
1924         int     i, count;
1925
1926         in = (dedge_t *)(mod_base + l->fileofs);
1927         if (l->filelen % sizeof(*in))
1928                 Host_Error("Mod_Q1BSP_LoadEdges: funny lump size in %s",loadmodel->name);
1929         count = l->filelen / sizeof(*in);
1930         out = (medge_t *)Mem_Alloc(loadmodel->mempool, count * sizeof(*out));
1931
1932         loadmodel->brushq1.edges = out;
1933         loadmodel->brushq1.numedges = count;
1934
1935         for ( i=0 ; i<count ; i++, in++, out++)
1936         {
1937                 out->v[0] = (unsigned short)LittleShort(in->v[0]);
1938                 out->v[1] = (unsigned short)LittleShort(in->v[1]);
1939                 if (out->v[0] >= loadmodel->brushq1.numvertexes || out->v[1] >= loadmodel->brushq1.numvertexes)
1940                 {
1941                         Con_Printf("Mod_Q1BSP_LoadEdges: %s has invalid vertex indices in edge %i (vertices %i %i >= numvertices %i)\n", loadmodel->name, i, out->v[0], out->v[1], loadmodel->brushq1.numvertexes);
1942                         out->v[0] = 0;
1943                         out->v[1] = 0;
1944                 }
1945         }
1946 }
1947
1948 static void Mod_Q1BSP_LoadTexinfo(lump_t *l)
1949 {
1950         texinfo_t *in;
1951         mtexinfo_t *out;
1952         int i, j, k, count, miptex;
1953
1954         in = (texinfo_t *)(mod_base + l->fileofs);
1955         if (l->filelen % sizeof(*in))
1956                 Host_Error("Mod_Q1BSP_LoadTexinfo: funny lump size in %s",loadmodel->name);
1957         count = l->filelen / sizeof(*in);
1958         out = (mtexinfo_t *)Mem_Alloc(loadmodel->mempool, count * sizeof(*out));
1959
1960         loadmodel->brushq1.texinfo = out;
1961         loadmodel->brushq1.numtexinfo = count;
1962
1963         for (i = 0;i < count;i++, in++, out++)
1964         {
1965                 for (k = 0;k < 2;k++)
1966                         for (j = 0;j < 4;j++)
1967                                 out->vecs[k][j] = LittleFloat(in->vecs[k][j]);
1968
1969                 miptex = LittleLong(in->miptex);
1970                 out->flags = LittleLong(in->flags);
1971
1972                 out->texture = NULL;
1973                 if (loadmodel->data_textures)
1974                 {
1975                         if ((unsigned int) miptex >= (unsigned int) loadmodel->num_textures)
1976                                 Con_Printf("error in model \"%s\": invalid miptex index %i(of %i)\n", loadmodel->name, miptex, loadmodel->num_textures);
1977                         else
1978                                 out->texture = loadmodel->data_textures + miptex;
1979                 }
1980                 if (out->flags & TEX_SPECIAL)
1981                 {
1982                         // if texture chosen is NULL or the shader needs a lightmap,
1983                         // force to notexture water shader
1984                         if (out->texture == NULL)
1985                                 out->texture = loadmodel->data_textures + (loadmodel->num_textures - 1);
1986                 }
1987                 else
1988                 {
1989                         // if texture chosen is NULL, force to notexture
1990                         if (out->texture == NULL)
1991                                 out->texture = loadmodel->data_textures + (loadmodel->num_textures - 2);
1992                 }
1993         }
1994 }
1995
1996 #if 0
1997 void BoundPoly(int numverts, float *verts, vec3_t mins, vec3_t maxs)
1998 {
1999         int             i, j;
2000         float   *v;
2001
2002         mins[0] = mins[1] = mins[2] = 9999;
2003         maxs[0] = maxs[1] = maxs[2] = -9999;
2004         v = verts;
2005         for (i = 0;i < numverts;i++)
2006         {
2007                 for (j = 0;j < 3;j++, v++)
2008                 {
2009                         if (*v < mins[j])
2010                                 mins[j] = *v;
2011                         if (*v > maxs[j])
2012                                 maxs[j] = *v;
2013                 }
2014         }
2015 }
2016
2017 #define MAX_SUBDIVPOLYTRIANGLES 4096
2018 #define MAX_SUBDIVPOLYVERTS(MAX_SUBDIVPOLYTRIANGLES * 3)
2019
2020 static int subdivpolyverts, subdivpolytriangles;
2021 static int subdivpolyindex[MAX_SUBDIVPOLYTRIANGLES][3];
2022 static float subdivpolyvert[MAX_SUBDIVPOLYVERTS][3];
2023
2024 static int subdivpolylookupvert(vec3_t v)
2025 {
2026         int i;
2027         for (i = 0;i < subdivpolyverts;i++)
2028                 if (subdivpolyvert[i][0] == v[0]
2029                  && subdivpolyvert[i][1] == v[1]
2030                  && subdivpolyvert[i][2] == v[2])
2031                         return i;
2032         if (subdivpolyverts >= MAX_SUBDIVPOLYVERTS)
2033                 Host_Error("SubDividePolygon: ran out of vertices in buffer, please increase your r_subdivide_size");
2034         VectorCopy(v, subdivpolyvert[subdivpolyverts]);
2035         return subdivpolyverts++;
2036 }
2037
2038 static void SubdividePolygon(int numverts, float *verts)
2039 {
2040         int             i, i1, i2, i3, f, b, c, p;
2041         vec3_t  mins, maxs, front[256], back[256];
2042         float   m, *pv, *cv, dist[256], frac;
2043
2044         if (numverts > 250)
2045                 Host_Error("SubdividePolygon: ran out of verts in buffer");
2046
2047         BoundPoly(numverts, verts, mins, maxs);
2048
2049         for (i = 0;i < 3;i++)
2050         {
2051                 m = (mins[i] + maxs[i]) * 0.5;
2052                 m = r_subdivide_size.value * floor(m/r_subdivide_size.value + 0.5);
2053                 if (maxs[i] - m < 8)
2054                         continue;
2055                 if (m - mins[i] < 8)
2056                         continue;
2057
2058                 // cut it
2059                 for (cv = verts, c = 0;c < numverts;c++, cv += 3)
2060                         dist[c] = cv[i] - m;
2061
2062                 f = b = 0;
2063                 for (p = numverts - 1, c = 0, pv = verts + p * 3, cv = verts;c < numverts;p = c, c++, pv = cv, cv += 3)
2064                 {
2065                         if (dist[p] >= 0)
2066                         {
2067                                 VectorCopy(pv, front[f]);
2068                                 f++;
2069                         }
2070                         if (dist[p] <= 0)
2071                         {
2072                                 VectorCopy(pv, back[b]);
2073                                 b++;
2074                         }
2075                         if (dist[p] == 0 || dist[c] == 0)
2076                                 continue;
2077                         if ((dist[p] > 0) != (dist[c] > 0) )
2078                         {
2079                                 // clip point
2080                                 frac = dist[p] / (dist[p] - dist[c]);
2081                                 front[f][0] = back[b][0] = pv[0] + frac * (cv[0] - pv[0]);
2082                                 front[f][1] = back[b][1] = pv[1] + frac * (cv[1] - pv[1]);
2083                                 front[f][2] = back[b][2] = pv[2] + frac * (cv[2] - pv[2]);
2084                                 f++;
2085                                 b++;
2086                         }
2087                 }
2088
2089                 SubdividePolygon(f, front[0]);
2090                 SubdividePolygon(b, back[0]);
2091                 return;
2092         }
2093
2094         i1 = subdivpolylookupvert(verts);
2095         i2 = subdivpolylookupvert(verts + 3);
2096         for (i = 2;i < numverts;i++)
2097         {
2098                 if (subdivpolytriangles >= MAX_SUBDIVPOLYTRIANGLES)
2099                 {
2100                         Con_Print("SubdividePolygon: ran out of triangles in buffer, please increase your r_subdivide_size\n");
2101                         return;
2102                 }
2103
2104                 i3 = subdivpolylookupvert(verts + i * 3);
2105                 subdivpolyindex[subdivpolytriangles][0] = i1;
2106                 subdivpolyindex[subdivpolytriangles][1] = i2;
2107                 subdivpolyindex[subdivpolytriangles][2] = i3;
2108                 i2 = i3;
2109                 subdivpolytriangles++;
2110         }
2111 }
2112
2113 //Breaks a polygon up along axial 64 unit
2114 //boundaries so that turbulent and sky warps
2115 //can be done reasonably.
2116 static void Mod_Q1BSP_GenerateWarpMesh(msurface_t *surface)
2117 {
2118         int i, j;
2119         surfvertex_t *v;
2120         surfmesh_t *mesh;
2121
2122         subdivpolytriangles = 0;
2123         subdivpolyverts = 0;
2124         SubdividePolygon(surface->num_vertices, (surface->mesh->data_vertex3f + 3 * surface->num_firstvertex));
2125         if (subdivpolytriangles < 1)
2126                 Host_Error("Mod_Q1BSP_GenerateWarpMesh: no triangles?");
2127
2128         surface->mesh = mesh = Mem_Alloc(loadmodel->mempool, sizeof(surfmesh_t) + subdivpolytriangles * sizeof(int[3]) + subdivpolyverts * sizeof(surfvertex_t));
2129         mesh->num_vertices = subdivpolyverts;
2130         mesh->num_triangles = subdivpolytriangles;
2131         mesh->vertex = (surfvertex_t *)(mesh + 1);
2132         mesh->index = (int *)(mesh->vertex + mesh->num_vertices);
2133         memset(mesh->vertex, 0, mesh->num_vertices * sizeof(surfvertex_t));
2134
2135         for (i = 0;i < mesh->num_triangles;i++)
2136                 for (j = 0;j < 3;j++)
2137                         mesh->index[i*3+j] = subdivpolyindex[i][j];
2138
2139         for (i = 0, v = mesh->vertex;i < subdivpolyverts;i++, v++)
2140         {
2141                 VectorCopy(subdivpolyvert[i], v->v);
2142                 v->st[0] = DotProduct(v->v, surface->lightmapinfo->texinfo->vecs[0]);
2143                 v->st[1] = DotProduct(v->v, surface->lightmapinfo->texinfo->vecs[1]);
2144         }
2145 }
2146 #endif
2147
2148 /* Maximum size of a single LM */
2149 #define MAX_SINGLE_LM_SIZE    256
2150
2151 struct alloc_lm_row
2152 {
2153         int rowY;
2154         int currentX;
2155 };
2156
2157 struct alloc_lm_state
2158 {
2159         int currentY;
2160         struct alloc_lm_row rows[MAX_SINGLE_LM_SIZE];
2161 };
2162
2163 static void init_alloc_lm_state (struct alloc_lm_state* state)
2164 {
2165         int r;
2166
2167         state->currentY = 0;
2168         for (r = 0; r < MAX_SINGLE_LM_SIZE; r++)
2169         {
2170           state->rows[r].currentX = 0;
2171           state->rows[r].rowY = -1;
2172         }
2173 }
2174
2175 static qboolean Mod_Q1BSP_AllocLightmapBlock(struct alloc_lm_state* state, int totalwidth, int totalheight, int blockwidth, int blockheight, int *outx, int *outy)
2176 {
2177         struct alloc_lm_row* row;
2178         int r;
2179
2180         row = &(state->rows[blockheight]);
2181         if ((row->rowY < 0) || (row->currentX + blockwidth > totalwidth))
2182         {
2183                 if (state->currentY + blockheight <= totalheight)
2184                 {
2185                         row->rowY = state->currentY;
2186                         row->currentX = 0;
2187                         state->currentY += blockheight;
2188                 }
2189                 else
2190                 {
2191                         /* See if we can stuff the block into a higher row */
2192                         row = NULL;
2193                         for (r = blockheight; r < MAX_SINGLE_LM_SIZE; r++)
2194                         {
2195                                 if ((state->rows[r].rowY >= 0)
2196                                   && (state->rows[r].currentX + blockwidth <= totalwidth))
2197                                 {
2198                                         row = &(state->rows[r]);
2199                                         break;
2200                                 }
2201                         }
2202                         if (row == NULL) return false;
2203                 }
2204         }
2205         *outy = row->rowY;
2206         *outx = row->currentX;
2207         row->currentX += blockwidth;
2208
2209         return true;
2210 }
2211
2212 extern cvar_t gl_max_size;
2213 static void Mod_Q1BSP_LoadFaces(lump_t *l)
2214 {
2215         dface_t *in;
2216         msurface_t *surface;
2217         int i, j, count, surfacenum, planenum, smax, tmax, ssize, tsize, firstedge, numedges, totalverts, totaltris, lightmapnumber, lightmapsize, totallightmapsamples;
2218         float texmins[2], texmaxs[2], val;
2219         rtexture_t *lightmaptexture, *deluxemaptexture;
2220
2221         in = (dface_t *)(mod_base + l->fileofs);
2222         if (l->filelen % sizeof(*in))
2223                 Host_Error("Mod_Q1BSP_LoadFaces: funny lump size in %s",loadmodel->name);
2224         count = l->filelen / sizeof(*in);
2225         loadmodel->data_surfaces = (msurface_t *)Mem_Alloc(loadmodel->mempool, count*sizeof(msurface_t));
2226         loadmodel->data_surfaces_lightmapinfo = (msurface_lightmapinfo_t *)Mem_Alloc(loadmodel->mempool, count*sizeof(msurface_lightmapinfo_t));
2227
2228         loadmodel->num_surfaces = count;
2229
2230         loadmodel->brushq1.lightmapupdateflags = (unsigned char *)Mem_Alloc(loadmodel->mempool, count*sizeof(unsigned char));
2231
2232         totalverts = 0;
2233         totaltris = 0;
2234         for (surfacenum = 0, in = (dface_t *)(mod_base + l->fileofs);surfacenum < count;surfacenum++, in++)
2235         {
2236                 numedges = (unsigned short)LittleShort(in->numedges);
2237                 totalverts += numedges;
2238                 totaltris += numedges - 2;
2239         }
2240
2241         Mod_AllocSurfMesh(loadmodel->mempool, totalverts, totaltris, true, false, false);
2242
2243         lightmaptexture = NULL;
2244         deluxemaptexture = r_texture_blanknormalmap;
2245         lightmapnumber = 1;
2246         lightmapsize = max(256, gl_max_size.integer);
2247         totallightmapsamples = 0;
2248
2249         totalverts = 0;
2250         totaltris = 0;
2251         for (surfacenum = 0, in = (dface_t *)(mod_base + l->fileofs), surface = loadmodel->data_surfaces;surfacenum < count;surfacenum++, in++, surface++)
2252         {
2253                 surface->lightmapinfo = loadmodel->data_surfaces_lightmapinfo + surfacenum;
2254
2255                 // FIXME: validate edges, texinfo, etc?
2256                 firstedge = LittleLong(in->firstedge);
2257                 numedges = (unsigned short)LittleShort(in->numedges);
2258                 if ((unsigned int) firstedge > (unsigned int) loadmodel->brushq1.numsurfedges || (unsigned int) numedges > (unsigned int) loadmodel->brushq1.numsurfedges || (unsigned int) firstedge + (unsigned int) numedges > (unsigned int) loadmodel->brushq1.numsurfedges)
2259                         Host_Error("Mod_Q1BSP_LoadFaces: invalid edge range (firstedge %i, numedges %i, model edges %i)", firstedge, numedges, loadmodel->brushq1.numsurfedges);
2260                 i = (unsigned short)LittleShort(in->texinfo);
2261                 if ((unsigned int) i >= (unsigned int) loadmodel->brushq1.numtexinfo)
2262                         Host_Error("Mod_Q1BSP_LoadFaces: invalid texinfo index %i(model has %i texinfos)", i, loadmodel->brushq1.numtexinfo);
2263                 surface->lightmapinfo->texinfo = loadmodel->brushq1.texinfo + i;
2264                 surface->texture = surface->lightmapinfo->texinfo->texture;
2265
2266                 planenum = (unsigned short)LittleShort(in->planenum);
2267                 if ((unsigned int) planenum >= (unsigned int) loadmodel->brush.num_planes)
2268                         Host_Error("Mod_Q1BSP_LoadFaces: invalid plane index %i (model has %i planes)", planenum, loadmodel->brush.num_planes);
2269
2270                 //surface->flags = surface->texture->flags;
2271                 //if (LittleShort(in->side))
2272                 //      surface->flags |= SURF_PLANEBACK;
2273                 //surface->plane = loadmodel->brush.data_planes + planenum;
2274
2275                 surface->num_firstvertex = totalverts;
2276                 surface->num_vertices = numedges;
2277                 surface->num_firsttriangle = totaltris;
2278                 surface->num_triangles = numedges - 2;
2279                 totalverts += numedges;
2280                 totaltris += numedges - 2;
2281
2282                 // convert edges back to a normal polygon
2283                 for (i = 0;i < surface->num_vertices;i++)
2284                 {
2285                         int lindex = loadmodel->brushq1.surfedges[firstedge + i];
2286                         float s, t;
2287                         if (lindex > 0)
2288                                 VectorCopy(loadmodel->brushq1.vertexes[loadmodel->brushq1.edges[lindex].v[0]].position, (loadmodel->surfmesh.data_vertex3f + 3 * surface->num_firstvertex) + i * 3);
2289                         else
2290                                 VectorCopy(loadmodel->brushq1.vertexes[loadmodel->brushq1.edges[-lindex].v[1]].position, (loadmodel->surfmesh.data_vertex3f + 3 * surface->num_firstvertex) + i * 3);
2291                         s = DotProduct(((loadmodel->surfmesh.data_vertex3f + 3 * surface->num_firstvertex) + i * 3), surface->lightmapinfo->texinfo->vecs[0]) + surface->lightmapinfo->texinfo->vecs[0][3];
2292                         t = DotProduct(((loadmodel->surfmesh.data_vertex3f + 3 * surface->num_firstvertex) + i * 3), surface->lightmapinfo->texinfo->vecs[1]) + surface->lightmapinfo->texinfo->vecs[1][3];
2293                         (loadmodel->surfmesh.data_texcoordtexture2f + 2 * surface->num_firstvertex)[i * 2 + 0] = s / surface->texture->width;
2294                         (loadmodel->surfmesh.data_texcoordtexture2f + 2 * surface->num_firstvertex)[i * 2 + 1] = t / surface->texture->height;
2295                         (loadmodel->surfmesh.data_texcoordlightmap2f + 2 * surface->num_firstvertex)[i * 2 + 0] = 0;
2296                         (loadmodel->surfmesh.data_texcoordlightmap2f + 2 * surface->num_firstvertex)[i * 2 + 1] = 0;
2297                         (loadmodel->surfmesh.data_lightmapoffsets + surface->num_firstvertex)[i] = 0;
2298                 }
2299
2300                 for (i = 0;i < surface->num_triangles;i++)
2301                 {
2302                         (loadmodel->surfmesh.data_element3i + 3 * surface->num_firsttriangle)[i * 3 + 0] = 0 + surface->num_firstvertex;
2303                         (loadmodel->surfmesh.data_element3i + 3 * surface->num_firsttriangle)[i * 3 + 1] = i + 1 + surface->num_firstvertex;
2304                         (loadmodel->surfmesh.data_element3i + 3 * surface->num_firsttriangle)[i * 3 + 2] = i + 2 + surface->num_firstvertex;
2305                 }
2306
2307                 // compile additional data about the surface geometry
2308                 Mod_BuildNormals(surface->num_firstvertex, surface->num_vertices, surface->num_triangles, loadmodel->surfmesh.data_vertex3f, (loadmodel->surfmesh.data_element3i + 3 * surface->num_firsttriangle), loadmodel->surfmesh.data_normal3f, true);
2309                 Mod_BuildTextureVectorsFromNormals(surface->num_firstvertex, surface->num_vertices, surface->num_triangles, loadmodel->surfmesh.data_vertex3f, loadmodel->surfmesh.data_texcoordtexture2f, loadmodel->surfmesh.data_normal3f, (loadmodel->surfmesh.data_element3i + 3 * surface->num_firsttriangle), loadmodel->surfmesh.data_svector3f, loadmodel->surfmesh.data_tvector3f, true);
2310                 BoxFromPoints(surface->mins, surface->maxs, surface->num_vertices, (loadmodel->surfmesh.data_vertex3f + 3 * surface->num_firstvertex));
2311
2312                 // generate surface extents information
2313                 texmins[0] = texmaxs[0] = DotProduct((loadmodel->surfmesh.data_vertex3f + 3 * surface->num_firstvertex), surface->lightmapinfo->texinfo->vecs[0]) + surface->lightmapinfo->texinfo->vecs[0][3];
2314                 texmins[1] = texmaxs[1] = DotProduct((loadmodel->surfmesh.data_vertex3f + 3 * surface->num_firstvertex), surface->lightmapinfo->texinfo->vecs[1]) + surface->lightmapinfo->texinfo->vecs[1][3];
2315                 for (i = 1;i < surface->num_vertices;i++)
2316                 {
2317                         for (j = 0;j < 2;j++)
2318                         {
2319                                 val = DotProduct((loadmodel->surfmesh.data_vertex3f + 3 * surface->num_firstvertex) + i * 3, surface->lightmapinfo->texinfo->vecs[j]) + surface->lightmapinfo->texinfo->vecs[j][3];
2320                                 texmins[j] = min(texmins[j], val);
2321                                 texmaxs[j] = max(texmaxs[j], val);
2322                         }
2323                 }
2324                 for (i = 0;i < 2;i++)
2325                 {
2326                         surface->lightmapinfo->texturemins[i] = (int) floor(texmins[i] / 16.0) * 16;
2327                         surface->lightmapinfo->extents[i] = (int) ceil(texmaxs[i] / 16.0) * 16 - surface->lightmapinfo->texturemins[i];
2328                 }
2329
2330                 smax = surface->lightmapinfo->extents[0] >> 4;
2331                 tmax = surface->lightmapinfo->extents[1] >> 4;
2332                 ssize = (surface->lightmapinfo->extents[0] >> 4) + 1;
2333                 tsize = (surface->lightmapinfo->extents[1] >> 4) + 1;
2334
2335                 // lighting info
2336                 for (i = 0;i < MAXLIGHTMAPS;i++)
2337                         surface->lightmapinfo->styles[i] = in->styles[i];
2338                 surface->lightmaptexture = NULL;
2339                 surface->deluxemaptexture = r_texture_blanknormalmap;
2340                 i = LittleLong(in->lightofs);
2341                 if (i == -1)
2342                 {
2343                         surface->lightmapinfo->samples = NULL;
2344                         // give non-lightmapped water a 1x white lightmap
2345                         if (surface->texture->name[0] == '*' && (surface->lightmapinfo->texinfo->flags & TEX_SPECIAL) && ssize <= 256 && tsize <= 256)
2346                         {
2347                                 surface->lightmapinfo->samples = (unsigned char *)Mem_Alloc(loadmodel->mempool, ssize * tsize * 3);
2348                                 surface->lightmapinfo->styles[0] = 0;
2349                                 memset(surface->lightmapinfo->samples, 128, ssize * tsize * 3);
2350                         }
2351                 }
2352                 else if (loadmodel->brush.ishlbsp) // LordHavoc: HalfLife map (bsp version 30)
2353                         surface->lightmapinfo->samples = loadmodel->brushq1.lightdata + i;
2354                 else // LordHavoc: white lighting (bsp version 29)
2355                 {
2356                         surface->lightmapinfo->samples = loadmodel->brushq1.lightdata + (i * 3);
2357                         if (loadmodel->brushq1.nmaplightdata)
2358                                 surface->lightmapinfo->nmapsamples = loadmodel->brushq1.nmaplightdata + (i * 3);
2359                 }
2360
2361                 // check if we should apply a lightmap to this
2362                 if (!(surface->lightmapinfo->texinfo->flags & TEX_SPECIAL) || surface->lightmapinfo->samples)
2363                 {
2364                         if (ssize > 256 || tsize > 256)
2365                                 Host_Error("Bad surface extents");
2366
2367                         if (lightmapsize < ssize)
2368                                 lightmapsize = ssize;
2369                         if (lightmapsize < tsize)
2370                                 lightmapsize = tsize;
2371
2372                         totallightmapsamples += ssize*tsize;
2373
2374                         // force lightmap upload on first time seeing the surface
2375                         //
2376                         // additionally this is used by the later code to see if a
2377                         // lightmap is needed on this surface (rather than duplicating the
2378                         // logic above)
2379                         loadmodel->brushq1.lightmapupdateflags[surfacenum] = true;
2380                 }
2381         }
2382
2383         // small maps (such as ammo boxes especially) don't need big lightmap
2384         // textures, so this code tries to guess a good size based on
2385         // totallightmapsamples (size of the lightmaps lump basically), as well as
2386         // trying to max out the gl_max_size if there is a lot of lightmap data to
2387         // store
2388         // additionally, never choose a lightmapsize that is smaller than the
2389         // largest surface encountered (as it would fail)
2390         i = lightmapsize;
2391         for (lightmapsize = 64; (lightmapsize < i) && (lightmapsize < gl_max_size.integer) && (totallightmapsamples > lightmapsize*lightmapsize); lightmapsize*=2)
2392                 ;
2393
2394         // now that we've decided the lightmap texture size, we can do the rest
2395         if (cls.state != ca_dedicated)
2396         {
2397                 struct alloc_lm_state allocState;
2398
2399                 for (surfacenum = 0, surface = loadmodel->data_surfaces;surfacenum < count;surfacenum++, surface++)
2400                 {
2401                         int i, iu, iv, lightmapx = 0, lightmapy = 0;
2402                         float u, v, ubase, vbase, uscale, vscale;
2403
2404                         smax = surface->lightmapinfo->extents[0] >> 4;
2405                         tmax = surface->lightmapinfo->extents[1] >> 4;
2406                         ssize = (surface->lightmapinfo->extents[0] >> 4) + 1;
2407                         tsize = (surface->lightmapinfo->extents[1] >> 4) + 1;
2408
2409                         if (!lightmaptexture || !Mod_Q1BSP_AllocLightmapBlock(&allocState, lightmapsize, lightmapsize, ssize, tsize, &lightmapx, &lightmapy))
2410                         {
2411                                 // allocate a texture pool if we need it
2412                                 if (loadmodel->texturepool == NULL)
2413                                         loadmodel->texturepool = R_AllocTexturePool();
2414                                 // could not find room, make a new lightmap
2415                                 lightmaptexture = R_LoadTexture2D(loadmodel->texturepool, va("lightmap%i", lightmapnumber), lightmapsize, lightmapsize, NULL, TEXTYPE_BGRA, TEXF_FORCELINEAR | TEXF_PRECACHE, NULL);
2416                                 if (loadmodel->brushq1.nmaplightdata)
2417                                         deluxemaptexture = R_LoadTexture2D(loadmodel->texturepool, va("deluxemap%i", lightmapnumber), lightmapsize, lightmapsize, NULL, TEXTYPE_BGRA, TEXF_FORCELINEAR | TEXF_PRECACHE, NULL);
2418                                 lightmapnumber++;
2419                                 init_alloc_lm_state (&allocState);
2420                                 Mod_Q1BSP_AllocLightmapBlock(&allocState, lightmapsize, lightmapsize, ssize, tsize, &lightmapx, &lightmapy);
2421                         }
2422                         surface->lightmaptexture = lightmaptexture;
2423                         surface->deluxemaptexture = deluxemaptexture;
2424                         surface->lightmapinfo->lightmaporigin[0] = lightmapx;
2425                         surface->lightmapinfo->lightmaporigin[1] = lightmapy;
2426
2427                         uscale = 1.0f / (float)lightmapsize;
2428                         vscale = 1.0f / (float)lightmapsize;
2429                         ubase = lightmapx * uscale;
2430                         vbase = lightmapy * vscale;
2431
2432                         for (i = 0;i < surface->num_vertices;i++)
2433                         {
2434                                 u = ((DotProduct(((loadmodel->surfmesh.data_vertex3f + 3 * surface->num_firstvertex) + i * 3), surface->lightmapinfo->texinfo->vecs[0]) + surface->lightmapinfo->texinfo->vecs[0][3]) + 8 - surface->lightmapinfo->texturemins[0]) * (1.0 / 16.0);
2435                                 v = ((DotProduct(((loadmodel->surfmesh.data_vertex3f + 3 * surface->num_firstvertex) + i * 3), surface->lightmapinfo->texinfo->vecs[1]) + surface->lightmapinfo->texinfo->vecs[1][3]) + 8 - surface->lightmapinfo->texturemins[1]) * (1.0 / 16.0);
2436                                 (loadmodel->surfmesh.data_texcoordlightmap2f + 2 * surface->num_firstvertex)[i * 2 + 0] = u * uscale + ubase;
2437                                 (loadmodel->surfmesh.data_texcoordlightmap2f + 2 * surface->num_firstvertex)[i * 2 + 1] = v * vscale + vbase;
2438                                 // LordHavoc: calc lightmap data offset for vertex lighting to use
2439                                 iu = (int) u;
2440                                 iv = (int) v;
2441                                 (loadmodel->surfmesh.data_lightmapoffsets + surface->num_firstvertex)[i] = (bound(0, iv, tmax) * ssize + bound(0, iu, smax)) * 3;
2442                         }
2443
2444                 }
2445
2446                 if (cl_stainmaps.integer)
2447                 {
2448                         // allocate stainmaps for permanent marks on walls
2449                         int stainmapsize = 0;
2450                         unsigned char *stainsamples = NULL;
2451                         for (surfacenum = 0, surface = loadmodel->data_surfaces;surfacenum < count;surfacenum++, surface++)
2452                         {
2453                                 ssize = (surface->lightmapinfo->extents[0] >> 4) + 1;
2454                                 tsize = (surface->lightmapinfo->extents[1] >> 4) + 1;
2455                                 stainmapsize += ssize * tsize * 3;
2456                         }
2457                         // allocate and clear to white
2458                         stainsamples = (unsigned char *)Mem_Alloc(loadmodel->mempool, stainmapsize);
2459                         memset(stainsamples, 255, stainmapsize);
2460                         // assign pointers
2461                         for (surfacenum = 0, surface = loadmodel->data_surfaces;surfacenum < count;surfacenum++, surface++)
2462                         {
2463                                 ssize = (surface->lightmapinfo->extents[0] >> 4) + 1;
2464                                 tsize = (surface->lightmapinfo->extents[1] >> 4) + 1;
2465                                 surface->lightmapinfo->stainsamples = stainsamples;
2466                                 stainsamples += ssize * tsize * 3;
2467                         }
2468                 }
2469         }
2470
2471         // generate ushort elements array if possible
2472         if (loadmodel->surfmesh.data_element3s)
2473                 for (i = 0;i < loadmodel->surfmesh.num_triangles*3;i++)
2474                         loadmodel->surfmesh.data_element3s[i] = loadmodel->surfmesh.data_element3i[i];
2475 }
2476
2477 static void Mod_Q1BSP_LoadNodes_RecursiveSetParent(mnode_t *node, mnode_t *parent)
2478 {
2479         //if (node->parent)
2480         //      Host_Error("Mod_Q1BSP_LoadNodes_RecursiveSetParent: runaway recursion");
2481         node->parent = parent;
2482         if (node->plane)
2483         {
2484                 // this is a node, recurse to children
2485                 Mod_Q1BSP_LoadNodes_RecursiveSetParent(node->children[0], node);
2486                 Mod_Q1BSP_LoadNodes_RecursiveSetParent(node->children[1], node);
2487                 // combine supercontents of children
2488                 node->combinedsupercontents = node->children[0]->combinedsupercontents | node->children[1]->combinedsupercontents;
2489         }
2490         else
2491         {
2492                 int j;
2493                 mleaf_t *leaf = (mleaf_t *)node;
2494                 // if this is a leaf, calculate supercontents mask from all collidable
2495                 // primitives in the leaf (brushes and collision surfaces)
2496                 // also flag if the leaf contains any collision surfaces
2497                 leaf->combinedsupercontents = 0;
2498                 // combine the supercontents values of all brushes in this leaf
2499                 for (j = 0;j < leaf->numleafbrushes;j++)
2500                         leaf->combinedsupercontents |= loadmodel->brush.data_brushes[leaf->firstleafbrush[j]].texture->supercontents;
2501                 // check if this leaf contains any collision surfaces (q3 patches)
2502                 for (j = 0;j < leaf->numleafsurfaces;j++)
2503                 {
2504                         msurface_t *surface = loadmodel->data_surfaces + leaf->firstleafsurface[j];
2505                         if (surface->num_collisiontriangles)
2506                         {
2507                                 leaf->containscollisionsurfaces = true;
2508                                 leaf->combinedsupercontents |= surface->texture->supercontents;
2509                         }
2510                 }
2511         }
2512 }
2513
2514 static void Mod_Q1BSP_LoadNodes(lump_t *l)
2515 {
2516         int                     i, j, count, p;
2517         dnode_t         *in;
2518         mnode_t         *out;
2519
2520         in = (dnode_t *)(mod_base + l->fileofs);
2521         if (l->filelen % sizeof(*in))
2522                 Host_Error("Mod_Q1BSP_LoadNodes: funny lump size in %s",loadmodel->name);
2523         count = l->filelen / sizeof(*in);
2524         out = (mnode_t *)Mem_Alloc(loadmodel->mempool, count*sizeof(*out));
2525
2526         loadmodel->brush.data_nodes = out;
2527         loadmodel->brush.num_nodes = count;
2528
2529         for ( i=0 ; i<count ; i++, in++, out++)
2530         {
2531                 for (j=0 ; j<3 ; j++)
2532                 {
2533                         out->mins[j] = LittleShort(in->mins[j]);
2534                         out->maxs[j] = LittleShort(in->maxs[j]);
2535                 }
2536
2537                 p = LittleLong(in->planenum);
2538                 out->plane = loadmodel->brush.data_planes + p;
2539
2540                 out->firstsurface = (unsigned short)LittleShort(in->firstface);
2541                 out->numsurfaces = (unsigned short)LittleShort(in->numfaces);
2542
2543                 for (j=0 ; j<2 ; j++)
2544                 {
2545                         // LordHavoc: this code supports broken bsp files produced by
2546                         // arguire qbsp which can produce more than 32768 nodes, any value
2547                         // below count is assumed to be a node number, any other value is
2548                         // assumed to be a leaf number
2549                         p = (unsigned short)LittleShort(in->children[j]);
2550                         if (p < count)
2551                         {
2552                                 if (p < loadmodel->brush.num_nodes)
2553                                         out->children[j] = loadmodel->brush.data_nodes + p;
2554                                 else
2555                                 {
2556                                         Con_Printf("Mod_Q1BSP_LoadNodes: invalid node index %i (file has only %i nodes)\n", p, loadmodel->brush.num_nodes);
2557                                         // map it to the solid leaf
2558                                         out->children[j] = (mnode_t *)loadmodel->brush.data_leafs;
2559                                 }
2560                         }
2561                         else
2562                         {
2563                                 // note this uses 65535 intentionally, -1 is leaf 0
2564                                 p = 65535 - p;
2565                                 if (p < loadmodel->brush.num_leafs)
2566                                         out->children[j] = (mnode_t *)(loadmodel->brush.data_leafs + p);
2567                                 else
2568                                 {
2569                                         Con_Printf("Mod_Q1BSP_LoadNodes: invalid leaf index %i (file has only %i leafs)\n", p, loadmodel->brush.num_leafs);
2570                                         // map it to the solid leaf
2571                                         out->children[j] = (mnode_t *)loadmodel->brush.data_leafs;
2572                                 }
2573                         }
2574                 }
2575         }
2576
2577         Mod_Q1BSP_LoadNodes_RecursiveSetParent(loadmodel->brush.data_nodes, NULL);      // sets nodes and leafs
2578 }
2579
2580 static void Mod_Q1BSP_LoadLeafs(lump_t *l)
2581 {
2582         dleaf_t *in;
2583         mleaf_t *out;
2584         int i, j, count, p;
2585
2586         in = (dleaf_t *)(mod_base + l->fileofs);
2587         if (l->filelen % sizeof(*in))
2588                 Host_Error("Mod_Q1BSP_LoadLeafs: funny lump size in %s",loadmodel->name);
2589         count = l->filelen / sizeof(*in);
2590         out = (mleaf_t *)Mem_Alloc(loadmodel->mempool, count*sizeof(*out));
2591
2592         loadmodel->brush.data_leafs = out;
2593         loadmodel->brush.num_leafs = count;
2594         // get visleafs from the submodel data
2595         loadmodel->brush.num_pvsclusters = loadmodel->brushq1.submodels[0].visleafs;
2596         loadmodel->brush.num_pvsclusterbytes = (loadmodel->brush.num_pvsclusters+7)>>3;
2597         loadmodel->brush.data_pvsclusters = (unsigned char *)Mem_Alloc(loadmodel->mempool, loadmodel->brush.num_pvsclusters * loadmodel->brush.num_pvsclusterbytes);
2598         memset(loadmodel->brush.data_pvsclusters, 0xFF, loadmodel->brush.num_pvsclusters * loadmodel->brush.num_pvsclusterbytes);
2599
2600         for ( i=0 ; i<count ; i++, in++, out++)
2601         {
2602                 for (j=0 ; j<3 ; j++)
2603                 {
2604                         out->mins[j] = LittleShort(in->mins[j]);
2605                         out->maxs[j] = LittleShort(in->maxs[j]);
2606                 }
2607
2608                 // FIXME: this function could really benefit from some error checking
2609
2610                 out->contents = LittleLong(in->contents);
2611
2612                 out->firstleafsurface = loadmodel->brush.data_leafsurfaces + (unsigned short)LittleShort(in->firstmarksurface);
2613                 out->numleafsurfaces = (unsigned short)LittleShort(in->nummarksurfaces);
2614                 if (out->firstleafsurface < 0 || (unsigned short)LittleShort(in->firstmarksurface) + out->numleafsurfaces > loadmodel->brush.num_leafsurfaces)
2615                 {
2616                         Con_Printf("Mod_Q1BSP_LoadLeafs: invalid leafsurface range %i:%i outside range %i:%i\n", (int)(out->firstleafsurface - loadmodel->brush.data_leafsurfaces), (int)(out->firstleafsurface + out->numleafsurfaces - loadmodel->brush.data_leafsurfaces), 0, loadmodel->brush.num_leafsurfaces);
2617                         out->firstleafsurface = NULL;
2618                         out->numleafsurfaces = 0;
2619                 }
2620
2621                 out->clusterindex = i - 1;
2622                 if (out->clusterindex >= loadmodel->brush.num_pvsclusters)
2623                         out->clusterindex = -1;
2624
2625                 p = LittleLong(in->visofs);
2626                 // ignore visofs errors on leaf 0 (solid)
2627                 if (p >= 0 && out->clusterindex >= 0)
2628                 {
2629                         if (p >= loadmodel->brushq1.num_compressedpvs)
2630                                 Con_Print("Mod_Q1BSP_LoadLeafs: invalid visofs\n");
2631                         else
2632                                 Mod_Q1BSP_DecompressVis(loadmodel->brushq1.data_compressedpvs + p, loadmodel->brushq1.data_compressedpvs + loadmodel->brushq1.num_compressedpvs, loadmodel->brush.data_pvsclusters + out->clusterindex * loadmodel->brush.num_pvsclusterbytes, loadmodel->brush.data_pvsclusters + (out->clusterindex + 1) * loadmodel->brush.num_pvsclusterbytes);
2633                 }
2634
2635                 for (j = 0;j < 4;j++)
2636                         out->ambient_sound_level[j] = in->ambient_level[j];
2637
2638                 // FIXME: Insert caustics here
2639         }
2640 }
2641
2642 qboolean Mod_Q1BSP_CheckWaterAlphaSupport(void)
2643 {
2644         int i, j;
2645         mleaf_t *leaf;
2646         const unsigned char *pvs;
2647         // if there's no vis data, assume supported (because everything is visible all the time)
2648         if (!loadmodel->brush.data_pvsclusters)
2649                 return true;
2650         // check all liquid leafs to see if they can see into empty leafs, if any
2651         // can we can assume this map supports r_wateralpha
2652         for (i = 0, leaf = loadmodel->brush.data_leafs;i < loadmodel->brush.num_leafs;i++, leaf++)
2653         {
2654                 if ((leaf->contents == CONTENTS_WATER || leaf->contents == CONTENTS_SLIME) && leaf->clusterindex >= 0)
2655                 {
2656                         pvs = loadmodel->brush.data_pvsclusters + leaf->clusterindex * loadmodel->brush.num_pvsclusterbytes;
2657                         for (j = 0;j < loadmodel->brush.num_leafs;j++)
2658                                 if (CHECKPVSBIT(pvs, loadmodel->brush.data_leafs[j].clusterindex) && loadmodel->brush.data_leafs[j].contents == CONTENTS_EMPTY)
2659                                         return true;
2660                 }
2661         }
2662         return false;
2663 }
2664
2665 static void Mod_Q1BSP_LoadClipnodes(lump_t *l, hullinfo_t *hullinfo)
2666 {
2667         dclipnode_t *in;
2668         mclipnode_t *out;
2669         int                     i, count;
2670         hull_t          *hull;
2671
2672         in = (dclipnode_t *)(mod_base + l->fileofs);
2673         if (l->filelen % sizeof(*in))
2674                 Host_Error("Mod_Q1BSP_LoadClipnodes: funny lump size in %s",loadmodel->name);
2675         count = l->filelen / sizeof(*in);
2676         out = (mclipnode_t *)Mem_Alloc(loadmodel->mempool, count*sizeof(*out));
2677
2678         loadmodel->brushq1.clipnodes = out;
2679         loadmodel->brushq1.numclipnodes = count;
2680
2681         for (i = 1; i < MAX_MAP_HULLS; i++)
2682         {
2683                 hull = &loadmodel->brushq1.hulls[i];
2684                 hull->clipnodes = out;
2685                 hull->firstclipnode = 0;
2686                 hull->lastclipnode = count-1;
2687                 hull->planes = loadmodel->brush.data_planes;
2688                 hull->clip_mins[0] = hullinfo->hullsizes[i][0][0];
2689                 hull->clip_mins[1] = hullinfo->hullsizes[i][0][1];
2690                 hull->clip_mins[2] = hullinfo->hullsizes[i][0][2];
2691                 hull->clip_maxs[0] = hullinfo->hullsizes[i][1][0];
2692                 hull->clip_maxs[1] = hullinfo->hullsizes[i][1][1];
2693                 hull->clip_maxs[2] = hullinfo->hullsizes[i][1][2];
2694                 VectorSubtract(hull->clip_maxs, hull->clip_mins, hull->clip_size);
2695         }
2696
2697         for (i=0 ; i<count ; i++, out++, in++)
2698         {
2699                 out->planenum = LittleLong(in->planenum);
2700                 // LordHavoc: this code supports arguire qbsp's broken clipnodes indices (more than 32768 clipnodes), values above count are assumed to be contents values
2701                 out->children[0] = (unsigned short)LittleShort(in->children[0]);
2702                 out->children[1] = (unsigned short)LittleShort(in->children[1]);
2703                 if (out->children[0] >= count)
2704                         out->children[0] -= 65536;
2705                 if (out->children[1] >= count)
2706                         out->children[1] -= 65536;
2707                 if (out->planenum < 0 || out->planenum >= loadmodel->brush.num_planes)
2708                         Host_Error("Corrupt clipping hull(out of range planenum)");
2709         }
2710 }
2711
2712 //Duplicate the drawing hull structure as a clipping hull
2713 static void Mod_Q1BSP_MakeHull0(void)
2714 {
2715         mnode_t         *in;
2716         mclipnode_t *out;
2717         int                     i;
2718         hull_t          *hull;
2719
2720         hull = &loadmodel->brushq1.hulls[0];
2721
2722         in = loadmodel->brush.data_nodes;
2723         out = (mclipnode_t *)Mem_Alloc(loadmodel->mempool, loadmodel->brush.num_nodes * sizeof(*out));
2724
2725         hull->clipnodes = out;
2726         hull->firstclipnode = 0;
2727         hull->lastclipnode = loadmodel->brush.num_nodes - 1;
2728         hull->planes = loadmodel->brush.data_planes;
2729
2730         for (i = 0;i < loadmodel->brush.num_nodes;i++, out++, in++)
2731         {
2732                 out->planenum = in->plane - loadmodel->brush.data_planes;
2733                 out->children[0] = in->children[0]->plane ? in->children[0] - loadmodel->brush.data_nodes : ((mleaf_t *)in->children[0])->contents;
2734                 out->children[1] = in->children[1]->plane ? in->children[1] - loadmodel->brush.data_nodes : ((mleaf_t *)in->children[1])->contents;
2735         }
2736 }
2737
2738 static void Mod_Q1BSP_LoadLeaffaces(lump_t *l)
2739 {
2740         int i, j;
2741         short *in;
2742
2743         in = (short *)(mod_base + l->fileofs);
2744         if (l->filelen % sizeof(*in))
2745                 Host_Error("Mod_Q1BSP_LoadLeaffaces: funny lump size in %s",loadmodel->name);
2746         loadmodel->brush.num_leafsurfaces = l->filelen / sizeof(*in);
2747         loadmodel->brush.data_leafsurfaces = (int *)Mem_Alloc(loadmodel->mempool, loadmodel->brush.num_leafsurfaces * sizeof(int));
2748
2749         for (i = 0;i < loadmodel->brush.num_leafsurfaces;i++)
2750         {
2751                 j = (unsigned short) LittleShort(in[i]);
2752                 if (j >= loadmodel->num_surfaces)
2753                         Host_Error("Mod_Q1BSP_LoadLeaffaces: bad surface number");
2754                 loadmodel->brush.data_leafsurfaces[i] = j;
2755         }
2756 }
2757
2758 static void Mod_Q1BSP_LoadSurfedges(lump_t *l)
2759 {
2760         int             i;
2761         int             *in;
2762
2763         in = (int *)(mod_base + l->fileofs);
2764         if (l->filelen % sizeof(*in))
2765                 Host_Error("Mod_Q1BSP_LoadSurfedges: funny lump size in %s",loadmodel->name);
2766         loadmodel->brushq1.numsurfedges = l->filelen / sizeof(*in);
2767         loadmodel->brushq1.surfedges = (int *)Mem_Alloc(loadmodel->mempool, loadmodel->brushq1.numsurfedges * sizeof(int));
2768
2769         for (i = 0;i < loadmodel->brushq1.numsurfedges;i++)
2770                 loadmodel->brushq1.surfedges[i] = LittleLong(in[i]);
2771 }
2772
2773
2774 static void Mod_Q1BSP_LoadPlanes(lump_t *l)
2775 {
2776         int                     i;
2777         mplane_t        *out;
2778         dplane_t        *in;
2779
2780         in = (dplane_t *)(mod_base + l->fileofs);
2781         if (l->filelen % sizeof(*in))
2782                 Host_Error("Mod_Q1BSP_LoadPlanes: funny lump size in %s", loadmodel->name);
2783
2784         loadmodel->brush.num_planes = l->filelen / sizeof(*in);
2785         loadmodel->brush.data_planes = out = (mplane_t *)Mem_Alloc(loadmodel->mempool, loadmodel->brush.num_planes * sizeof(*out));
2786
2787         for (i = 0;i < loadmodel->brush.num_planes;i++, in++, out++)
2788         {
2789                 out->normal[0] = LittleFloat(in->normal[0]);
2790                 out->normal[1] = LittleFloat(in->normal[1]);
2791                 out->normal[2] = LittleFloat(in->normal[2]);
2792                 out->dist = LittleFloat(in->dist);
2793
2794                 PlaneClassify(out);
2795         }
2796 }
2797
2798 static void Mod_Q1BSP_LoadMapBrushes(void)
2799 {
2800 #if 0
2801 // unfinished
2802         int submodel, numbrushes;
2803         qboolean firstbrush;
2804         char *text, *maptext;
2805         char mapfilename[MAX_QPATH];
2806         FS_StripExtension (loadmodel->name, mapfilename, sizeof (mapfilename));
2807         strlcat (mapfilename, ".map", sizeof (mapfilename));
2808         maptext = (unsigned char*) FS_LoadFile(mapfilename, tempmempool, false, NULL);
2809         if (!maptext)
2810                 return;
2811         text = maptext;
2812         if (!COM_ParseToken_Simple(&data, false, false))
2813                 return; // error
2814         submodel = 0;
2815         for (;;)
2816         {
2817                 if (!COM_ParseToken_Simple(&data, false, false))
2818                         break;
2819                 if (com_token[0] != '{')
2820                         return; // error
2821                 // entity
2822                 firstbrush = true;
2823                 numbrushes = 0;
2824                 maxbrushes = 256;
2825                 brushes = Mem_Alloc(loadmodel->mempool, maxbrushes * sizeof(mbrush_t));
2826                 for (;;)
2827                 {
2828                         if (!COM_ParseToken_Simple(&data, false, false))
2829                                 return; // error
2830                         if (com_token[0] == '}')
2831                                 break; // end of entity
2832                         if (com_token[0] == '{')
2833                         {
2834                                 // brush
2835                                 if (firstbrush)
2836                                 {
2837                                         if (submodel)
2838                                         {
2839                                                 if (submodel > loadmodel->brush.numsubmodels)
2840                                                 {
2841                                                         Con_Printf("Mod_Q1BSP_LoadMapBrushes: .map has more submodels than .bsp!\n");
2842                                                         model = NULL;
2843                                                 }
2844                                                 else
2845                                                         model = loadmodel->brush.submodels[submodel];
2846                                         }
2847                                         else
2848                                                 model = loadmodel;
2849                                 }
2850                                 for (;;)
2851                                 {
2852                                         if (!COM_ParseToken_Simple(&data, false, false))
2853                                                 return; // error
2854                                         if (com_token[0] == '}')
2855                                                 break; // end of brush
2856                                         // each brush face should be this format:
2857                                         // ( x y z ) ( x y z ) ( x y z ) texture scroll_s scroll_t rotateangle scale_s scale_t
2858                                         // FIXME: support hl .map format
2859                                         for (pointnum = 0;pointnum < 3;pointnum++)
2860                                         {
2861                                                 COM_ParseToken_Simple(&data, false, false);
2862                                                 for (componentnum = 0;componentnum < 3;componentnum++)
2863                                                 {
2864                                                         COM_ParseToken_Simple(&data, false, false);
2865                                                         point[pointnum][componentnum] = atof(com_token);
2866                                                 }
2867                                                 COM_ParseToken_Simple(&data, false, false);
2868                                         }
2869                                         COM_ParseToken_Simple(&data, false, false);
2870                                         strlcpy(facetexture, com_token, sizeof(facetexture));
2871                                         COM_ParseToken_Simple(&data, false, false);
2872                                         //scroll_s = atof(com_token);
2873                                         COM_ParseToken_Simple(&data, false, false);
2874                                         //scroll_t = atof(com_token);
2875                                         COM_ParseToken_Simple(&data, false, false);
2876                                         //rotate = atof(com_token);
2877                                         COM_ParseToken_Simple(&data, false, false);
2878                                         //scale_s = atof(com_token);
2879                                         COM_ParseToken_Simple(&data, false, false);
2880                                         //scale_t = atof(com_token);
2881                                         TriangleNormal(point[0], point[1], point[2], planenormal);
2882                                         VectorNormalizeDouble(planenormal);
2883                                         planedist = DotProduct(point[0], planenormal);
2884                                         //ChooseTexturePlane(planenormal, texturevector[0], texturevector[1]);
2885                                 }
2886                                 continue;
2887                         }
2888                 }
2889         }
2890 #endif
2891 }
2892
2893
2894 #define MAX_PORTALPOINTS 64
2895
2896 typedef struct portal_s
2897 {
2898         mplane_t plane;
2899         mnode_t *nodes[2];              // [0] = front side of plane
2900         struct portal_s *next[2];
2901         int numpoints;
2902         double points[3*MAX_PORTALPOINTS];
2903         struct portal_s *chain; // all portals are linked into a list
2904 }
2905 portal_t;
2906
2907 static portal_t *portalchain;
2908
2909 /*
2910 ===========
2911 AllocPortal
2912 ===========
2913 */
2914 static portal_t *AllocPortal(void)
2915 {
2916         portal_t *p;
2917         p = (portal_t *)Mem_Alloc(loadmodel->mempool, sizeof(portal_t));
2918         p->chain = portalchain;
2919         portalchain = p;
2920         return p;
2921 }
2922
2923 static void FreePortal(portal_t *p)
2924 {
2925         Mem_Free(p);
2926 }
2927
2928 static void Mod_Q1BSP_RecursiveRecalcNodeBBox(mnode_t *node)
2929 {
2930         // process only nodes (leafs already had their box calculated)
2931         if (!node->plane)
2932                 return;
2933
2934         // calculate children first
2935         Mod_Q1BSP_RecursiveRecalcNodeBBox(node->children[0]);
2936         Mod_Q1BSP_RecursiveRecalcNodeBBox(node->children[1]);
2937
2938         // make combined bounding box from children
2939         node->mins[0] = min(node->children[0]->mins[0], node->children[1]->mins[0]);
2940         node->mins[1] = min(node->children[0]->mins[1], node->children[1]->mins[1]);
2941         node->mins[2] = min(node->children[0]->mins[2], node->children[1]->mins[2]);
2942         node->maxs[0] = max(node->children[0]->maxs[0], node->children[1]->maxs[0]);
2943         node->maxs[1] = max(node->children[0]->maxs[1], node->children[1]->maxs[1]);
2944         node->maxs[2] = max(node->children[0]->maxs[2], node->children[1]->maxs[2]);
2945 }
2946
2947 static void Mod_Q1BSP_FinalizePortals(void)
2948 {
2949         int i, j, numportals, numpoints;
2950         portal_t *p, *pnext;
2951         mportal_t *portal;
2952         mvertex_t *point;
2953         mleaf_t *leaf, *endleaf;
2954
2955         // tally up portal and point counts and recalculate bounding boxes for all
2956         // leafs (because qbsp is very sloppy)
2957         leaf = loadmodel->brush.data_leafs;
2958         endleaf = leaf + loadmodel->brush.num_leafs;
2959         for (;leaf < endleaf;leaf++)
2960         {
2961                 VectorSet(leaf->mins,  2000000000,  2000000000,  2000000000);
2962                 VectorSet(leaf->maxs, -2000000000, -2000000000, -2000000000);
2963         }
2964         p = portalchain;
2965         numportals = 0;
2966         numpoints = 0;
2967         while (p)
2968         {
2969                 // note: this check must match the one below or it will usually corrupt memory
2970                 // the nodes[0] != nodes[1] check is because leaf 0 is the shared solid leaf, it can have many portals inside with leaf 0 on both sides
2971                 if (p->numpoints >= 3 && p->nodes[0] != p->nodes[1] && ((mleaf_t *)p->nodes[0])->clusterindex >= 0 && ((mleaf_t *)p->nodes[1])->clusterindex >= 0)
2972                 {
2973                         numportals += 2;
2974                         numpoints += p->numpoints * 2;
2975                 }
2976                 p = p->chain;
2977         }
2978         loadmodel->brush.data_portals = (mportal_t *)Mem_Alloc(loadmodel->mempool, numportals * sizeof(mportal_t) + numpoints * sizeof(mvertex_t));
2979         loadmodel->brush.num_portals = numportals;
2980         loadmodel->brush.data_portalpoints = (mvertex_t *)((unsigned char *) loadmodel->brush.data_portals + numportals * sizeof(mportal_t));
2981         loadmodel->brush.num_portalpoints = numpoints;
2982         // clear all leaf portal chains
2983         for (i = 0;i < loadmodel->brush.num_leafs;i++)
2984                 loadmodel->brush.data_leafs[i].portals = NULL;
2985         // process all portals in the global portal chain, while freeing them
2986         portal = loadmodel->brush.data_portals;
2987         point = loadmodel->brush.data_portalpoints;
2988         p = portalchain;
2989         portalchain = NULL;
2990         while (p)
2991         {
2992                 pnext = p->chain;
2993
2994                 if (p->numpoints >= 3 && p->nodes[0] != p->nodes[1])
2995                 {
2996                         // note: this check must match the one above or it will usually corrupt memory
2997                         // the nodes[0] != nodes[1] check is because leaf 0 is the shared solid leaf, it can have many portals inside with leaf 0 on both sides
2998                         if (((mleaf_t *)p->nodes[0])->clusterindex >= 0 && ((mleaf_t *)p->nodes[1])->clusterindex >= 0)
2999                         {
3000                                 // first make the back to front portal(forward portal)
3001                                 portal->points = point;
3002                                 portal->numpoints = p->numpoints;
3003                                 portal->plane.dist = p->plane.dist;
3004                                 VectorCopy(p->plane.normal, portal->plane.normal);
3005                                 portal->here = (mleaf_t *)p->nodes[1];
3006                                 portal->past = (mleaf_t *)p->nodes[0];
3007                                 // copy points
3008                                 for (j = 0;j < portal->numpoints;j++)
3009                                 {
3010                                         VectorCopy(p->points + j*3, point->position);
3011                                         point++;
3012                                 }
3013                                 BoxFromPoints(portal->mins, portal->maxs, portal->numpoints, portal->points->position);
3014                                 PlaneClassify(&portal->plane);
3015
3016                                 // link into leaf's portal chain
3017                                 portal->next = portal->here->portals;
3018                                 portal->here->portals = portal;
3019
3020                                 // advance to next portal
3021                                 portal++;
3022
3023                                 // then make the front to back portal(backward portal)
3024                                 portal->points = point;
3025                                 portal->numpoints = p->numpoints;
3026                                 portal->plane.dist = -p->plane.dist;
3027                                 VectorNegate(p->plane.normal, portal->plane.normal);
3028                                 portal->here = (mleaf_t *)p->nodes[0];
3029                                 portal->past = (mleaf_t *)p->nodes[1];
3030                                 // copy points
3031                                 for (j = portal->numpoints - 1;j >= 0;j--)
3032                                 {
3033                                         VectorCopy(p->points + j*3, point->position);
3034                                         point++;
3035                                 }
3036                                 BoxFromPoints(portal->mins, portal->maxs, portal->numpoints, portal->points->position);
3037                                 PlaneClassify(&portal->plane);
3038
3039                                 // link into leaf's portal chain
3040                                 portal->next = portal->here->portals;
3041                                 portal->here->portals = portal;
3042
3043                                 // advance to next portal
3044                                 portal++;
3045                         }
3046                         // add the portal's polygon points to the leaf bounding boxes
3047                         for (i = 0;i < 2;i++)
3048                         {
3049                                 leaf = (mleaf_t *)p->nodes[i];
3050                                 for (j = 0;j < p->numpoints;j++)
3051                                 {
3052                                         if (leaf->mins[0] > p->points[j*3+0]) leaf->mins[0] = p->points[j*3+0];
3053                                         if (leaf->mins[1] > p->points[j*3+1]) leaf->mins[1] = p->points[j*3+1];
3054                                         if (leaf->mins[2] > p->points[j*3+2]) leaf->mins[2] = p->points[j*3+2];
3055                                         if (leaf->maxs[0] < p->points[j*3+0]) leaf->maxs[0] = p->points[j*3+0];
3056                                         if (leaf->maxs[1] < p->points[j*3+1]) leaf->maxs[1] = p->points[j*3+1];
3057                                         if (leaf->maxs[2] < p->points[j*3+2]) leaf->maxs[2] = p->points[j*3+2];
3058                                 }
3059                         }
3060                 }
3061                 FreePortal(p);
3062                 p = pnext;
3063         }
3064         // now recalculate the node bounding boxes from the leafs
3065         Mod_Q1BSP_RecursiveRecalcNodeBBox(loadmodel->brush.data_nodes);
3066 }
3067
3068 /*
3069 =============
3070 AddPortalToNodes
3071 =============
3072 */
3073 static void AddPortalToNodes(portal_t *p, mnode_t *front, mnode_t *back)
3074 {
3075         if (!front)
3076                 Host_Error("AddPortalToNodes: NULL front node");
3077         if (!back)
3078                 Host_Error("AddPortalToNodes: NULL back node");
3079         if (p->nodes[0] || p->nodes[1])
3080                 Host_Error("AddPortalToNodes: already included");
3081         // note: front == back is handled gracefully, because leaf 0 is the shared solid leaf, it can often have portals with the same leaf on both sides
3082
3083         p->nodes[0] = front;
3084         p->next[0] = (portal_t *)front->portals;
3085         front->portals = (mportal_t *)p;
3086
3087         p->nodes[1] = back;
3088         p->next[1] = (portal_t *)back->portals;
3089         back->portals = (mportal_t *)p;
3090 }
3091
3092 /*
3093 =============
3094 RemovePortalFromNode
3095 =============
3096 */
3097 static void RemovePortalFromNodes(portal_t *portal)
3098 {
3099         int i;
3100         mnode_t *node;
3101         void **portalpointer;
3102         portal_t *t;
3103         for (i = 0;i < 2;i++)
3104         {
3105                 node = portal->nodes[i];
3106
3107                 portalpointer = (void **) &node->portals;
3108                 while (1)
3109                 {
3110                         t = (portal_t *)*portalpointer;
3111                         if (!t)
3112                                 Host_Error("RemovePortalFromNodes: portal not in leaf");
3113
3114                         if (t == portal)
3115                         {
3116                                 if (portal->nodes[0] == node)
3117                                 {
3118                                         *portalpointer = portal->next[0];
3119                                         portal->nodes[0] = NULL;
3120                                 }
3121                                 else if (portal->nodes[1] == node)
3122                                 {
3123                                         *portalpointer = portal->next[1];
3124                                         portal->nodes[1] = NULL;
3125                                 }
3126                                 else
3127                                         Host_Error("RemovePortalFromNodes: portal not bounding leaf");
3128                                 break;
3129                         }
3130
3131                         if (t->nodes[0] == node)
3132                                 portalpointer = (void **) &t->next[0];
3133                         else if (t->nodes[1] == node)
3134                                 portalpointer = (void **) &t->next[1];
3135                         else
3136                                 Host_Error("RemovePortalFromNodes: portal not bounding leaf");
3137                 }
3138         }
3139 }
3140
3141 #define PORTAL_DIST_EPSILON (1.0 / 32.0)
3142 static void Mod_Q1BSP_RecursiveNodePortals(mnode_t *node)
3143 {
3144         int i, side;
3145         mnode_t *front, *back, *other_node;
3146         mplane_t clipplane, *plane;
3147         portal_t *portal, *nextportal, *nodeportal, *splitportal, *temp;
3148         int numfrontpoints, numbackpoints;
3149         double frontpoints[3*MAX_PORTALPOINTS], backpoints[3*MAX_PORTALPOINTS];
3150
3151         // if a leaf, we're done
3152         if (!node->plane)
3153                 return;
3154
3155         plane = node->plane;
3156
3157         front = node->children[0];
3158         back = node->children[1];
3159         if (front == back)
3160                 Host_Error("Mod_Q1BSP_RecursiveNodePortals: corrupt node hierarchy");
3161
3162         // create the new portal by generating a polygon for the node plane,
3163         // and clipping it by all of the other portals(which came from nodes above this one)
3164         nodeportal = AllocPortal();
3165         nodeportal->plane = *plane;
3166
3167         // TODO: calculate node bounding boxes during recursion and calculate a maximum plane size accordingly to improve precision (as most maps do not need 1 billion unit plane polygons)
3168         PolygonD_QuadForPlane(nodeportal->points, nodeportal->plane.normal[0], nodeportal->plane.normal[1], nodeportal->plane.normal[2], nodeportal->plane.dist, 1024.0*1024.0*1024.0);
3169         nodeportal->numpoints = 4;
3170         side = 0;       // shut up compiler warning
3171         for (portal = (portal_t *)node->portals;portal;portal = portal->next[side])
3172         {
3173                 clipplane = portal->plane;
3174                 if (portal->nodes[0] == portal->nodes[1])
3175                         Host_Error("Mod_Q1BSP_RecursiveNodePortals: portal has same node on both sides(1)");
3176                 if (portal->nodes[0] == node)
3177                         side = 0;
3178                 else if (portal->nodes[1] == node)
3179                 {
3180                         clipplane.dist = -clipplane.dist;
3181                         VectorNegate(clipplane.normal, clipplane.normal);
3182                         side = 1;
3183                 }
3184                 else
3185                         Host_Error("Mod_Q1BSP_RecursiveNodePortals: mislinked portal");
3186
3187                 for (i = 0;i < nodeportal->numpoints*3;i++)
3188                         frontpoints[i] = nodeportal->points[i];
3189                 PolygonD_Divide(nodeportal->numpoints, frontpoints, clipplane.normal[0], clipplane.normal[1], clipplane.normal[2], clipplane.dist, PORTAL_DIST_EPSILON, MAX_PORTALPOINTS, nodeportal->points, &nodeportal->numpoints, 0, NULL, NULL, NULL);
3190                 if (nodeportal->numpoints <= 0 || nodeportal->numpoints >= MAX_PORTALPOINTS)
3191                         break;
3192         }
3193
3194         if (nodeportal->numpoints < 3)
3195         {
3196                 Con_Print("Mod_Q1BSP_RecursiveNodePortals: WARNING: new portal was clipped away\n");
3197                 nodeportal->numpoints = 0;
3198         }
3199         else if (nodeportal->numpoints >= MAX_PORTALPOINTS)
3200         {
3201                 Con_Print("Mod_Q1BSP_RecursiveNodePortals: WARNING: new portal has too many points\n");
3202                 nodeportal->numpoints = 0;
3203         }
3204
3205         AddPortalToNodes(nodeportal, front, back);
3206
3207         // split the portals of this node along this node's plane and assign them to the children of this node
3208         // (migrating the portals downward through the tree)
3209         for (portal = (portal_t *)node->portals;portal;portal = nextportal)
3210         {
3211                 if (portal->nodes[0] == portal->nodes[1])
3212                         Host_Error("Mod_Q1BSP_RecursiveNodePortals: portal has same node on both sides(2)");
3213                 if (portal->nodes[0] == node)
3214                         side = 0;
3215                 else if (portal->nodes[1] == node)
3216                         side = 1;
3217                 else
3218                         Host_Error("Mod_Q1BSP_RecursiveNodePortals: mislinked portal");
3219                 nextportal = portal->next[side];
3220                 if (!portal->numpoints)
3221                         continue;
3222
3223                 other_node = portal->nodes[!side];
3224                 RemovePortalFromNodes(portal);
3225
3226                 // cut the portal into two portals, one on each side of the node plane
3227                 PolygonD_Divide(portal->numpoints, portal->points, plane->normal[0], plane->normal[1], plane->normal[2], plane->dist, PORTAL_DIST_EPSILON, MAX_PORTALPOINTS, frontpoints, &numfrontpoints, MAX_PORTALPOINTS, backpoints, &numbackpoints, NULL);
3228
3229                 if (!numfrontpoints)
3230                 {
3231                         if (side == 0)
3232                                 AddPortalToNodes(portal, back, other_node);
3233                         else
3234                                 AddPortalToNodes(portal, other_node, back);
3235                         continue;
3236                 }
3237                 if (!numbackpoints)
3238                 {
3239                         if (side == 0)
3240                                 AddPortalToNodes(portal, front, other_node);
3241                         else
3242                                 AddPortalToNodes(portal, other_node, front);
3243                         continue;
3244                 }
3245
3246                 // the portal is split
3247                 splitportal = AllocPortal();
3248                 temp = splitportal->chain;
3249                 *splitportal = *portal;
3250                 splitportal->chain = temp;
3251                 for (i = 0;i < numbackpoints*3;i++)
3252                         splitportal->points[i] = backpoints[i];
3253                 splitportal->numpoints = numbackpoints;
3254                 for (i = 0;i < numfrontpoints*3;i++)
3255                         portal->points[i] = frontpoints[i];
3256                 portal->numpoints = numfrontpoints;
3257
3258                 if (side == 0)
3259                 {
3260                         AddPortalToNodes(portal, front, other_node);
3261                         AddPortalToNodes(splitportal, back, other_node);
3262                 }
3263                 else
3264                 {
3265                         AddPortalToNodes(portal, other_node, front);
3266                         AddPortalToNodes(splitportal, other_node, back);
3267                 }
3268         }
3269
3270         Mod_Q1BSP_RecursiveNodePortals(front);
3271         Mod_Q1BSP_RecursiveNodePortals(back);
3272 }
3273
3274 static void Mod_Q1BSP_MakePortals(void)
3275 {
3276         portalchain = NULL;
3277         Mod_Q1BSP_RecursiveNodePortals(loadmodel->brush.data_nodes);
3278         Mod_Q1BSP_FinalizePortals();
3279 }
3280
3281 //Returns PVS data for a given point
3282 //(note: can return NULL)
3283 static unsigned char *Mod_Q1BSP_GetPVS(dp_model_t *model, const vec3_t p)
3284 {
3285         mnode_t *node;
3286         node = model->brush.data_nodes;
3287         while (node->plane)
3288                 node = node->children[(node->plane->type < 3 ? p[node->plane->type] : DotProduct(p,node->plane->normal)) < node->plane->dist];
3289         if (((mleaf_t *)node)->clusterindex >= 0)
3290                 return model->brush.data_pvsclusters + ((mleaf_t *)node)->clusterindex * model->brush.num_pvsclusterbytes;
3291         else
3292                 return NULL;
3293 }
3294
3295 static void Mod_Q1BSP_FatPVS_RecursiveBSPNode(dp_model_t *model, const vec3_t org, vec_t radius, unsigned char *pvsbuffer, int pvsbytes, mnode_t *node)
3296 {
3297         while (node->plane)
3298         {
3299                 float d = PlaneDiff(org, node->plane);
3300                 if (d > radius)
3301                         node = node->children[0];
3302                 else if (d < -radius)
3303                         node = node->children[1];
3304                 else
3305                 {
3306                         // go down both sides
3307                         Mod_Q1BSP_FatPVS_RecursiveBSPNode(model, org, radius, pvsbuffer, pvsbytes, node->children[0]);
3308                         node = node->children[1];
3309                 }
3310         }
3311         // if this leaf is in a cluster, accumulate the pvs bits
3312         if (((mleaf_t *)node)->clusterindex >= 0)
3313         {
3314                 int i;
3315                 unsigned char *pvs = model->brush.data_pvsclusters + ((mleaf_t *)node)->clusterindex * model->brush.num_pvsclusterbytes;
3316                 for (i = 0;i < pvsbytes;i++)