]> de.git.xonotic.org Git - xonotic/netradiant.git/blobdiff - tools/quake3/q3map2/vis.c
Merge branch 'q3map2help' into 'master'
[xonotic/netradiant.git] / tools / quake3 / q3map2 / vis.c
index f3b2f4ae51d4ee97da956ba5e9fc32c6697e80ff..b6340473daea2945925fbb5cfb25a417239e0a6a 100644 (file)
-/*\r
-Copyright (C) 1999-2007 id Software, Inc. and contributors.\r
-For a list of contributors, see the accompanying CONTRIBUTORS file.\r
-\r
-This file is part of GtkRadiant.\r
-\r
-GtkRadiant is free software; you can redistribute it and/or modify\r
-it under the terms of the GNU General Public License as published by\r
-the Free Software Foundation; either version 2 of the License, or\r
-(at your option) any later version.\r
-\r
-GtkRadiant is distributed in the hope that it will be useful,\r
-but WITHOUT ANY WARRANTY; without even the implied warranty of\r
-MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the\r
-GNU General Public License for more details.\r
-\r
-You should have received a copy of the GNU General Public License\r
-along with GtkRadiant; if not, write to the Free Software\r
-Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA\r
-\r
-----------------------------------------------------------------------------------\r
-\r
-This code has been altered significantly from its original form, to support\r
-several games based on the Quake III Arena engine, in the form of "Q3Map2."\r
-\r
-------------------------------------------------------------------------------- */\r
-\r
-\r
-\r
-/* marker */\r
-#define VIS_C\r
-\r
-\r
-\r
-/* dependencies */\r
-#include "q3map2.h"\r
-\r
-\r
-\r
-\r
-void PlaneFromWinding (fixedWinding_t *w, visPlane_t *plane)\r
-{\r
-       vec3_t          v1, v2;\r
-\r
-// calc plane\r
-       VectorSubtract (w->points[2], w->points[1], v1);\r
-       VectorSubtract (w->points[0], w->points[1], v2);\r
-       CrossProduct (v2, v1, plane->normal);\r
-       VectorNormalize (plane->normal, plane->normal);\r
-       plane->dist = DotProduct (w->points[0], plane->normal);\r
-}\r
-\r
-\r
-/*\r
-NewFixedWinding()\r
-returns a new fixed winding\r
-ydnar: altered this a bit to reconcile multiply-defined winding_t\r
-*/\r
-\r
-fixedWinding_t *NewFixedWinding( int points )\r
-{\r
-       fixedWinding_t  *w;\r
-       int                     size;\r
-       \r
-       if (points > MAX_POINTS_ON_WINDING)\r
-               Error ("NewWinding: %i points", points);\r
-       \r
-       size = (int)((fixedWinding_t *)0)->points[points];\r
-       w = safe_malloc (size);\r
-       memset (w, 0, size);\r
-       \r
-       return w;\r
-}\r
-\r
-\r
-\r
-void prl(leaf_t *l)\r
-{\r
-       int                     i;\r
-       vportal_t       *p;\r
-       visPlane_t      pl;\r
-       \r
-       for (i=0 ; i<l->numportals ; i++)\r
-       {\r
-               p = l->portals[i];\r
-               pl = p->plane;\r
-               Sys_Printf ("portal %4i to leaf %4i : %7.1f : (%4.1f, %4.1f, %4.1f)\n",(int)(p-portals),p->leaf,pl.dist, pl.normal[0], pl.normal[1], pl.normal[2]);\r
-       }\r
-}\r
-\r
-\r
-//=============================================================================\r
-\r
-/*\r
-=============\r
-SortPortals\r
-\r
-Sorts the portals from the least complex, so the later ones can reuse\r
-the earlier information.\r
-=============\r
-*/\r
-int PComp (const void *a, const void *b)\r
-{\r
-       if ( (*(vportal_t **)a)->nummightsee == (*(vportal_t **)b)->nummightsee)\r
-               return 0;\r
-       if ( (*(vportal_t **)a)->nummightsee < (*(vportal_t **)b)->nummightsee)\r
-               return -1;\r
-       return 1;\r
-}\r
-void SortPortals (void)\r
-{\r
-       int             i;\r
-       \r
-       for (i=0 ; i<numportals*2 ; i++)\r
-               sorted_portals[i] = &portals[i];\r
-\r
-       if (nosort)\r
-               return;\r
-       qsort (sorted_portals, numportals*2, sizeof(sorted_portals[0]), PComp);\r
-}\r
-\r
-\r
-/*\r
-==============\r
-LeafVectorFromPortalVector\r
-==============\r
-*/\r
-int LeafVectorFromPortalVector (byte *portalbits, byte *leafbits)\r
-{\r
-       int                     i, j, leafnum;\r
-       vportal_t       *p;\r
-       int                     c_leafs;\r
-\r
-\r
-       for (i=0 ; i<numportals*2 ; i++)\r
-       {\r
-               if (portalbits[i>>3] & (1<<(i&7)) )\r
-               {\r
-                       p = portals+i;\r
-                       leafbits[p->leaf>>3] |= (1<<(p->leaf&7));\r
-               }\r
-       }\r
-\r
-       for (j = 0; j < portalclusters; j++)\r
-       {\r
-               leafnum = j;\r
-               while (leafs[leafnum].merged >= 0)\r
-                       leafnum = leafs[leafnum].merged;\r
-               //if the merged leaf is visible then the original leaf is visible\r
-               if (leafbits[leafnum>>3] & (1<<(leafnum&7)))\r
-               {\r
-                       leafbits[j>>3] |= (1<<(j&7));\r
-               }\r
-       }\r
-\r
-       c_leafs = CountBits (leafbits, portalclusters);\r
-\r
-       return c_leafs;\r
-}\r
-\r
-\r
-/*\r
-===============\r
-ClusterMerge\r
-\r
-Merges the portal visibility for a leaf\r
-===============\r
-*/\r
-void ClusterMerge (int leafnum)\r
-{\r
-       leaf_t          *leaf;\r
-       byte            portalvector[MAX_PORTALS/8];\r
-       byte            uncompressed[MAX_MAP_LEAFS/8];\r
-       int                     i, j;\r
-       int                     numvis, mergedleafnum;\r
-       vportal_t       *p;\r
-       int                     pnum;\r
-\r
-       // OR together all the portalvis bits\r
-\r
-       mergedleafnum = leafnum;\r
-       while(leafs[mergedleafnum].merged >= 0)\r
-               mergedleafnum = leafs[mergedleafnum].merged;\r
-\r
-       memset (portalvector, 0, portalbytes);\r
-       leaf = &leafs[mergedleafnum];\r
-       for (i = 0; i < leaf->numportals; i++)\r
-       {\r
-               p = leaf->portals[i];\r
-               if (p->removed)\r
-                       continue;\r
-\r
-               if (p->status != stat_done)\r
-                       Error ("portal not done");\r
-               for (j=0 ; j<portallongs ; j++)\r
-                       ((long *)portalvector)[j] |= ((long *)p->portalvis)[j];\r
-               pnum = p - portals;\r
-               portalvector[pnum>>3] |= 1<<(pnum&7);\r
-       }\r
-\r
-       memset (uncompressed, 0, leafbytes);\r
-\r
-       uncompressed[mergedleafnum>>3] |= (1<<(mergedleafnum&7));\r
-       // convert portal bits to leaf bits\r
-       numvis = LeafVectorFromPortalVector (portalvector, uncompressed);\r
-\r
-//     if (uncompressed[leafnum>>3] & (1<<(leafnum&7)))\r
-//             Sys_Printf ("WARNING: Leaf portals saw into leaf\n");\r
-               \r
-//     uncompressed[leafnum>>3] |= (1<<(leafnum&7));\r
-\r
-       numvis++;               // count the leaf itself\r
-\r
-       totalvis += numvis;\r
-\r
-       Sys_FPrintf (SYS_VRB,"cluster %4i : %4i visible\n", leafnum, numvis);\r
-\r
-       memcpy (bspVisBytes + VIS_HEADER_SIZE + leafnum*leafbytes, uncompressed, leafbytes);\r
-}\r
-\r
-/*\r
-==================\r
-CalcPortalVis\r
-==================\r
-*/\r
-void CalcPortalVis (void)\r
-{\r
-#ifdef MREDEBUG\r
-       Sys_Printf("%6d portals out of %d", 0, numportals*2);\r
-       //get rid of the counter\r
-       RunThreadsOnIndividual (numportals*2, qfalse, PortalFlow);\r
-#else\r
-       RunThreadsOnIndividual (numportals*2, qtrue, PortalFlow);\r
-#endif\r
-\r
-}\r
-\r
-/*\r
-==================\r
-CalcPassageVis\r
-==================\r
-*/\r
-void CalcPassageVis(void)\r
-{\r
-       PassageMemory();\r
-\r
-#ifdef MREDEBUG\r
-       _printf("%6d portals out of %d", 0, numportals*2);\r
-       RunThreadsOnIndividual (numportals*2, qfalse, CreatePassages);\r
-       _printf("\n");\r
-       _printf("%6d portals out of %d", 0, numportals*2);\r
-       RunThreadsOnIndividual (numportals*2, qfalse, PassageFlow);\r
-       _printf("\n");\r
-#else\r
-       Sys_Printf( "\n--- CreatePassages (%d) ---\n", numportals * 2 );\r
-       RunThreadsOnIndividual( numportals*2, qtrue, CreatePassages );\r
-       \r
-       Sys_Printf( "\n--- PassageFlow (%d) ---\n", numportals * 2 );\r
-       RunThreadsOnIndividual( numportals * 2, qtrue, PassageFlow );\r
-#endif\r
-}\r
-\r
-/*\r
-==================\r
-CalcPassagePortalVis\r
-==================\r
-*/\r
-void CalcPassagePortalVis(void)\r
-{\r
-       PassageMemory();\r
-\r
-#ifdef MREDEBUG\r
-       Sys_Printf("%6d portals out of %d", 0, numportals*2);\r
-       RunThreadsOnIndividual (numportals*2, qfalse, CreatePassages);\r
-       Sys_Printf("\n");\r
-       Sys_Printf("%6d portals out of %d", 0, numportals*2);\r
-       RunThreadsOnIndividual (numportals*2, qfalse, PassagePortalFlow);\r
-       Sys_Printf("\n");\r
-#else\r
-       Sys_Printf( "\n--- CreatePassages (%d) ---\n", numportals * 2 );\r
-       RunThreadsOnIndividual( numportals * 2, qtrue, CreatePassages);\r
-       \r
-       Sys_Printf( "\n--- PassagePortalFlow (%d) ---\n", numportals * 2 );\r
-       RunThreadsOnIndividual( numportals * 2, qtrue, PassagePortalFlow );\r
-#endif\r
-}\r
-\r
-/*\r
-==================\r
-CalcFastVis\r
-==================\r
-*/\r
-void CalcFastVis(void)\r
-{\r
-       int             i;\r
-\r
-       // fastvis just uses mightsee for a very loose bound\r
-       for (i=0 ; i<numportals*2 ; i++)\r
-       {\r
-               portals[i].portalvis = portals[i].portalflood;\r
-               portals[i].status = stat_done;\r
-       }\r
-}\r
-\r
-/*\r
-==================\r
-CalcVis\r
-==================\r
-*/\r
-void CalcVis (void)\r
-{\r
-       int                     i;\r
-       const char      *value;\r
-       \r
-       \r
-       /* ydnar: rr2do2's farplane code */\r
-       farPlaneDist = 0.0f;\r
-       value = ValueForKey( &entities[ 0 ], "_farplanedist" );         /* proper '_' prefixed key */\r
-       if( value[ 0 ] == '\0' )\r
-               value = ValueForKey( &entities[ 0 ], "fogclip" );               /* wolf compatibility */\r
-       if( value[ 0 ] == '\0' )\r
-               value = ValueForKey( &entities[ 0 ], "distancecull" );  /* sof2 compatibility */\r
-       if( value[ 0 ] != '\0' )\r
-       {\r
-               farPlaneDist = atof( value );\r
-               if( farPlaneDist > 0.0f )\r
-                       Sys_Printf( "farplane distance = %.1f\n", farPlaneDist );\r
-               else\r
-                       farPlaneDist = 0.0f;\r
-       }\r
-       \r
-       \r
-       \r
-       Sys_Printf( "\n--- BasePortalVis (%d) ---\n", numportals * 2 );\r
-       RunThreadsOnIndividual( numportals * 2, qtrue, BasePortalVis );\r
-       \r
-//     RunThreadsOnIndividual (numportals*2, qtrue, BetterPortalVis);\r
-\r
-       SortPortals ();\r
-\r
-  if (fastvis) {\r
-    CalcFastVis();\r
-  }\r
-  else if ( noPassageVis ) {\r
-    CalcPortalVis();\r
-  }\r
-  else if ( passageVisOnly ) {\r
-    CalcPassageVis();\r
-  }\r
-  else {\r
-    CalcPassagePortalVis();\r
-  }\r
-       //\r
-       // assemble the leaf vis lists by oring and compressing the portal lists\r
-       //\r
-       Sys_Printf("creating leaf vis...\n");\r
-       for (i=0 ; i<portalclusters ; i++)\r
-               ClusterMerge (i);\r
-\r
-  Sys_Printf( "Total visible clusters: %i\n", totalvis );\r
-  Sys_Printf( "Average clusters visible: %i\n", totalvis / portalclusters );\r
-}\r
-\r
-/*\r
-==================\r
-SetPortalSphere\r
-==================\r
-*/\r
-void SetPortalSphere (vportal_t *p)\r
-{\r
-       int             i;\r
-       vec3_t  total, dist;\r
-       fixedWinding_t  *w;\r
-       float   r, bestr;\r
-\r
-       w = p->winding;\r
-       VectorCopy (vec3_origin, total);\r
-       for (i=0 ; i<w->numpoints ; i++)\r
-       {\r
-               VectorAdd (total, w->points[i], total);\r
-       }\r
-       \r
-       for (i=0 ; i<3 ; i++)\r
-               total[i] /= w->numpoints;\r
-\r
-       bestr = 0;              \r
-       for (i=0 ; i<w->numpoints ; i++)\r
-       {\r
-               VectorSubtract (w->points[i], total, dist);\r
-               r = VectorLength (dist);\r
-               if (r > bestr)\r
-                       bestr = r;\r
-       }\r
-       VectorCopy (total, p->origin);\r
-       p->radius = bestr;\r
-}\r
-\r
-/*\r
-=============\r
-Winding_PlanesConcave\r
-=============\r
-*/\r
-#define WCONVEX_EPSILON                0.2\r
-\r
-int Winding_PlanesConcave(fixedWinding_t *w1, fixedWinding_t *w2,\r
-                                                        vec3_t normal1, vec3_t normal2,\r
-                                                        float dist1, float dist2)\r
-{\r
-       int i;\r
-\r
-       if (!w1 || !w2) return qfalse;\r
-\r
-       // check if one of the points of winding 1 is at the front of the plane of winding 2\r
-       for (i = 0; i < w1->numpoints; i++)\r
-       {\r
-               if (DotProduct(normal2, w1->points[i]) - dist2 > WCONVEX_EPSILON) return qtrue;\r
-       }\r
-       // check if one of the points of winding 2 is at the front of the plane of winding 1\r
-       for (i = 0; i < w2->numpoints; i++)\r
-       {\r
-               if (DotProduct(normal1, w2->points[i]) - dist1 > WCONVEX_EPSILON) return qtrue;\r
-       }\r
-\r
-       return qfalse;\r
-}\r
-\r
-/*\r
-============\r
-TryMergeLeaves\r
-============\r
-*/\r
-int TryMergeLeaves(int l1num, int l2num)\r
-{\r
-       int i, j, k, n, numportals;\r
-       visPlane_t plane1, plane2;\r
-       leaf_t *l1, *l2;\r
-       vportal_t *p1, *p2;\r
-       vportal_t *portals[MAX_PORTALS_ON_LEAF];\r
-\r
-       for (k = 0; k < 2; k++)\r
-       {\r
-               if (k) l1 = &leafs[l1num];\r
-               else l1 = &faceleafs[l1num];\r
-               for (i = 0; i < l1->numportals; i++)\r
-               {\r
-                       p1 = l1->portals[i];\r
-                       if (p1->leaf == l2num) continue;\r
-                       for (n = 0; n < 2; n++)\r
-                       {\r
-                               if (n) l2 = &leafs[l2num];\r
-                               else l2 = &faceleafs[l2num];\r
-                               for (j = 0; j < l2->numportals; j++)\r
-                               {\r
-                                       p2 = l2->portals[j];\r
-                                       if (p2->leaf == l1num) continue;\r
-                                       //\r
-                                       plane1 = p1->plane;\r
-                                       plane2 = p2->plane;\r
-                                       if (Winding_PlanesConcave(p1->winding, p2->winding, plane1.normal, plane2.normal, plane1.dist, plane2.dist))\r
-                                               return qfalse;\r
-                               }\r
-                       }\r
-               }\r
-       }\r
-       for (k = 0; k < 2; k++)\r
-       {\r
-               if (k)\r
-               {\r
-                       l1 = &leafs[l1num];\r
-                       l2 = &leafs[l2num];\r
-               }\r
-               else\r
-               {\r
-                       l1 = &faceleafs[l1num];\r
-                       l2 = &faceleafs[l2num];\r
-               }\r
-               numportals = 0;\r
-               //the leaves can be merged now\r
-               for (i = 0; i < l1->numportals; i++)\r
-               {\r
-                       p1 = l1->portals[i];\r
-                       if (p1->leaf == l2num)\r
-                       {\r
-                               p1->removed = qtrue;\r
-                               continue;\r
-                       }\r
-                       portals[numportals++] = p1;\r
-               }\r
-               for (j = 0; j < l2->numportals; j++)\r
-               {\r
-                       p2 = l2->portals[j];\r
-                       if (p2->leaf == l1num)\r
-                       {\r
-                               p2->removed = qtrue;\r
-                               continue;\r
-                       }\r
-                       portals[numportals++] = p2;\r
-               }\r
-               for (i = 0; i < numportals; i++)\r
-               {\r
-                       l2->portals[i] = portals[i];\r
-               }\r
-               l2->numportals = numportals;\r
-               l1->merged = l2num;\r
-       }\r
-       return qtrue;\r
-}\r
-\r
-/*\r
-============\r
-UpdatePortals\r
-============\r
-*/\r
-void UpdatePortals(void)\r
-{\r
-       int i;\r
-       vportal_t *p;\r
-\r
-       for (i = 0; i < numportals * 2; i++)\r
-       {\r
-               p = &portals[i];\r
-               if (p->removed)\r
-                       continue;\r
-               while(leafs[p->leaf].merged >= 0)\r
-                       p->leaf = leafs[p->leaf].merged;\r
-       }\r
-}\r
-\r
-/*\r
-============\r
-MergeLeaves\r
-\r
-try to merge leaves but don't merge through hint splitters\r
-============\r
-*/\r
-void MergeLeaves(void)\r
-{\r
-       int i, j, nummerges, totalnummerges;\r
-       leaf_t *leaf;\r
-       vportal_t *p;\r
-\r
-       totalnummerges = 0;\r
-       do\r
-       {\r
-               nummerges = 0;\r
-               for (i = 0; i < portalclusters; i++)\r
-               {\r
-                       leaf = &leafs[i];\r
-                       //if this leaf is merged already\r
-\r
-                       /* ydnar: vmods: merge all non-hint portals */\r
-                       if( leaf->merged >= 0 && hint == qfalse )\r
-                               continue;\r
-\r
-\r
-                       for (j = 0; j < leaf->numportals; j++)\r
-                       {\r
-                               p = leaf->portals[j];\r
-                               //\r
-                               if (p->removed)\r
-                                       continue;\r
-                               //never merge through hint portals\r
-                               if (p->hint)\r
-                                       continue;\r
-                               if (TryMergeLeaves(i, p->leaf))\r
-                               {\r
-                                       UpdatePortals();\r
-                                       nummerges++;\r
-                                       break;\r
-                               }\r
-                       }\r
-               }\r
-               totalnummerges += nummerges;\r
-       } while (nummerges);\r
-       Sys_Printf("%6d leaves merged\n", totalnummerges);\r
-}\r
-\r
-/*\r
-============\r
-TryMergeWinding\r
-============\r
-*/\r
-#define        CONTINUOUS_EPSILON      0.005\r
-\r
-fixedWinding_t *TryMergeWinding (fixedWinding_t *f1, fixedWinding_t *f2, vec3_t planenormal)\r
-{\r
-       vec_t           *p1, *p2, *p3, *p4, *back;\r
-       fixedWinding_t  *newf;\r
-       int                     i, j, k, l;\r
-       vec3_t          normal, delta;\r
-       vec_t           dot;\r
-       qboolean        keep1, keep2;\r
-       \r
-\r
-       //\r
-       // find a common edge\r
-       //      \r
-       p1 = p2 = NULL; // stop compiler warning\r
-       j = 0;                  // \r
-       \r
-       for (i = 0; i < f1->numpoints; i++)\r
-       {\r
-               p1 = f1->points[i];\r
-               p2 = f1->points[(i+1) % f1->numpoints];\r
-               for (j = 0; j < f2->numpoints; j++)\r
-               {\r
-                       p3 = f2->points[j];\r
-                       p4 = f2->points[(j+1) % f2->numpoints];\r
-                       for (k = 0; k < 3; k++)\r
-                       {\r
-                               if (fabs(p1[k] - p4[k]) > 0.1)//EQUAL_EPSILON) //ME\r
-                                       break;\r
-                               if (fabs(p2[k] - p3[k]) > 0.1)//EQUAL_EPSILON) //ME\r
-                                       break;\r
-                       } //end for\r
-                       if (k==3)\r
-                               break;\r
-               } //end for\r
-               if (j < f2->numpoints)\r
-                       break;\r
-       } //end for\r
-       \r
-       if (i == f1->numpoints)\r
-               return NULL;                    // no matching edges\r
-\r
-       //\r
-       // check slope of connected lines\r
-       // if the slopes are colinear, the point can be removed\r
-       //\r
-       back = f1->points[(i+f1->numpoints-1)%f1->numpoints];\r
-       VectorSubtract (p1, back, delta);\r
-       CrossProduct (planenormal, delta, normal);\r
-       VectorNormalize (normal, normal);\r
-       \r
-       back = f2->points[(j+2)%f2->numpoints];\r
-       VectorSubtract (back, p1, delta);\r
-       dot = DotProduct (delta, normal);\r
-       if (dot > CONTINUOUS_EPSILON)\r
-               return NULL;                    // not a convex polygon\r
-       keep1 = (qboolean)(dot < -CONTINUOUS_EPSILON);\r
-       \r
-       back = f1->points[(i+2)%f1->numpoints];\r
-       VectorSubtract (back, p2, delta);\r
-       CrossProduct (planenormal, delta, normal);\r
-       VectorNormalize (normal, normal);\r
-\r
-       back = f2->points[(j+f2->numpoints-1)%f2->numpoints];\r
-       VectorSubtract (back, p2, delta);\r
-       dot = DotProduct (delta, normal);\r
-       if (dot > CONTINUOUS_EPSILON)\r
-               return NULL;                    // not a convex polygon\r
-       keep2 = (qboolean)(dot < -CONTINUOUS_EPSILON);\r
-\r
-       //\r
-       // build the new polygon\r
-       //\r
-       newf = NewFixedWinding (f1->numpoints + f2->numpoints);\r
-       \r
-       // copy first polygon\r
-       for (k=(i+1)%f1->numpoints ; k != i ; k=(k+1)%f1->numpoints)\r
-       {\r
-               if (k==(i+1)%f1->numpoints && !keep2)\r
-                       continue;\r
-               \r
-               VectorCopy (f1->points[k], newf->points[newf->numpoints]);\r
-               newf->numpoints++;\r
-       }\r
-       \r
-       // copy second polygon\r
-       for (l= (j+1)%f2->numpoints ; l != j ; l=(l+1)%f2->numpoints)\r
-       {\r
-               if (l==(j+1)%f2->numpoints && !keep1)\r
-                       continue;\r
-               VectorCopy (f2->points[l], newf->points[newf->numpoints]);\r
-               newf->numpoints++;\r
-       }\r
-\r
-       return newf;\r
-}\r
-\r
-/*\r
-============\r
-MergeLeafPortals\r
-============\r
-*/\r
-void MergeLeafPortals(void)\r
-{\r
-       int i, j, k, nummerges, hintsmerged;\r
-       leaf_t *leaf;\r
-       vportal_t *p1, *p2;\r
-       fixedWinding_t *w;\r
-\r
-       nummerges = 0;\r
-       hintsmerged = 0;\r
-       for (i = 0; i < portalclusters; i++)\r
-       {\r
-               leaf = &leafs[i];\r
-               if (leaf->merged >= 0) continue;\r
-               for (j = 0; j < leaf->numportals; j++)\r
-               {\r
-                       p1 = leaf->portals[j];\r
-                       if (p1->removed)\r
-                               continue;\r
-                       for (k = j+1; k < leaf->numportals; k++)\r
-                       {\r
-                               p2 = leaf->portals[k];\r
-                               if (p2->removed)\r
-                                       continue;\r
-                               if (p1->leaf == p2->leaf)\r
-                               {\r
-                                       w = TryMergeWinding(p1->winding, p2->winding, p1->plane.normal);\r
-                                       if (w)\r
-                                       {\r
-                                               free( p1->winding );    //% FreeWinding(p1->winding);\r
-                                               p1->winding = w;\r
-                                               if (p1->hint && p2->hint)\r
-                                                       hintsmerged++;\r
-                                               p1->hint |= p2->hint;\r
-                                               SetPortalSphere(p1);\r
-                                               p2->removed = qtrue;\r
-                                               nummerges++;\r
-                                               i--;\r
-                                               break;\r
-                                       }\r
-                               }\r
-                       }\r
-                       if (k < leaf->numportals)\r
-                               break;\r
-               }\r
-       }\r
-       Sys_Printf("%6d portals merged\n", nummerges);\r
-       Sys_Printf("%6d hint portals merged\n", hintsmerged);\r
-}\r
-\r
-\r
-/*\r
-============\r
-WritePortals\r
-============\r
-*/\r
-int CountActivePortals(void)\r
-{\r
-       int num, hints, j;\r
-       vportal_t *p;\r
-\r
-       num = 0;\r
-       hints = 0;\r
-       for (j = 0; j < numportals * 2; j++)\r
-       {\r
-               p = portals + j;\r
-               if (p->removed)\r
-                       continue;\r
-               if (p->hint)\r
-                       hints++;\r
-               num++;\r
-       }\r
-       Sys_Printf("%6d active portals\n", num);\r
-       Sys_Printf("%6d hint portals\n", hints);\r
-       return num;\r
-}\r
-\r
-/*\r
-============\r
-WritePortals\r
-============\r
-*/\r
-void WriteFloat (FILE *f, vec_t v);\r
-\r
-void WritePortals(char *filename)\r
-{\r
-       int i, j, num;\r
-       FILE *pf;\r
-       vportal_t *p;\r
-       fixedWinding_t *w;\r
-\r
-       // write the file\r
-       pf = fopen (filename, "w");\r
-       if (!pf)\r
-               Error ("Error opening %s", filename);\r
-\r
-       num = 0;\r
-       for (j = 0; j < numportals * 2; j++)\r
-       {\r
-               p = portals + j;\r
-               if (p->removed)\r
-                       continue;\r
-//             if (!p->hint)\r
-//                     continue;\r
-               num++;\r
-       }\r
-\r
-       fprintf (pf, "%s\n", PORTALFILE);\r
-       fprintf (pf, "%i\n", 0);\r
-       fprintf (pf, "%i\n", num);// + numfaces);\r
-       fprintf (pf, "%i\n", 0);\r
-\r
-       for (j = 0; j < numportals * 2; j++)\r
-       {\r
-               p = portals + j;\r
-               if (p->removed)\r
-                       continue;\r
-//             if (!p->hint)\r
-//                     continue;\r
-               w = p->winding;\r
-               fprintf (pf,"%i %i %i ",w->numpoints, 0, 0);\r
-               fprintf (pf, "%d ", p->hint);\r
-               for (i=0 ; i<w->numpoints ; i++)\r
-               {\r
-                       fprintf (pf,"(");\r
-                       WriteFloat (pf, w->points[i][0]);\r
-                       WriteFloat (pf, w->points[i][1]);\r
-                       WriteFloat (pf, w->points[i][2]);\r
-                       fprintf (pf,") ");\r
-               }\r
-               fprintf (pf,"\n");\r
-       }\r
-\r
-       /*\r
-       for (j = 0; j < numfaces; j++)\r
-       {\r
-               p = faces + j;\r
-               w = p->winding;\r
-               fprintf (pf,"%i %i %i ",w->numpoints, 0, 0);\r
-               fprintf (pf, "0 ");\r
-               for (i=0 ; i<w->numpoints ; i++)\r
-               {\r
-                       fprintf (pf,"(");\r
-                       WriteFloat (pf, w->points[i][0]);\r
-                       WriteFloat (pf, w->points[i][1]);\r
-                       WriteFloat (pf, w->points[i][2]);\r
-                       fprintf (pf,") ");\r
-               }\r
-               fprintf (pf,"\n");\r
-       }*/\r
-\r
-       fclose (pf);\r
-}\r
-\r
-/*\r
-============\r
-LoadPortals\r
-============\r
-*/\r
-void LoadPortals (char *name)\r
-{\r
-       int                     i, j, hint;\r
-       vportal_t       *p;\r
-       leaf_t          *l;\r
-       char            magic[80];\r
-       FILE            *f;\r
-       int                     numpoints;\r
-       fixedWinding_t  *w;\r
-       int                     leafnums[2];\r
-       visPlane_t      plane;\r
-       \r
-       if (!strcmp(name,"-"))\r
-               f = stdin;\r
-       else\r
-       {\r
-               f = fopen(name, "r");\r
-               if (!f)\r
-                       Error ("LoadPortals: couldn't read %s\n",name);\r
-       }\r
-\r
-       if (fscanf (f,"%79s\n%i\n%i\n%i\n",magic, &portalclusters, &numportals, &numfaces) != 4)\r
-               Error ("LoadPortals: failed to read header");\r
-       if (strcmp(magic,PORTALFILE))\r
-               Error ("LoadPortals: not a portal file");\r
-\r
-       Sys_Printf ("%6i portalclusters\n", portalclusters);\r
-       Sys_Printf ("%6i numportals\n", numportals);\r
-       Sys_Printf ("%6i numfaces\n", numfaces);\r
-\r
-       // these counts should take advantage of 64 bit systems automatically\r
-       leafbytes = ((portalclusters+63)&~63)>>3;\r
-       leaflongs = leafbytes/sizeof(long);\r
-       \r
-       portalbytes = ((numportals*2+63)&~63)>>3;\r
-       portallongs = portalbytes/sizeof(long);\r
-\r
-       // each file portal is split into two memory portals\r
-       portals = safe_malloc(2*numportals*sizeof(vportal_t));\r
-       memset (portals, 0, 2*numportals*sizeof(vportal_t));\r
-       \r
-       leafs = safe_malloc(portalclusters*sizeof(leaf_t));\r
-       memset (leafs, 0, portalclusters*sizeof(leaf_t));\r
-\r
-       for (i = 0; i < portalclusters; i++)\r
-               leafs[i].merged = -1;\r
-\r
-       numBSPVisBytes = VIS_HEADER_SIZE + portalclusters*leafbytes;\r
-\r
-       if (numBSPVisBytes > MAX_MAP_VISIBILITY)\r
-         Error("MAX_MAP_VISIBILITY exceeded");\r
-\r
-       ((int *)bspVisBytes)[0] = portalclusters;\r
-       ((int *)bspVisBytes)[1] = leafbytes;\r
-               \r
-       for (i=0, p=portals ; i<numportals ; i++)\r
-       {\r
-               if (fscanf (f, "%i %i %i ", &numpoints, &leafnums[0], &leafnums[1]) != 3)\r
-                       Error ("LoadPortals: reading portal %i", i);\r
-               if (numpoints > MAX_POINTS_ON_WINDING)\r
-                       Error ("LoadPortals: portal %i has too many points", i);\r
-               if ( (unsigned)leafnums[0] > portalclusters\r
-               || (unsigned)leafnums[1] > portalclusters)\r
-                       Error ("LoadPortals: reading portal %i", i);\r
-               if (fscanf (f, "%i ", &hint) != 1)\r
-                       Error ("LoadPortals: reading hint state");\r
-               \r
-               w = p->winding = NewFixedWinding (numpoints);\r
-               w->numpoints = numpoints;\r
-               \r
-               for (j=0 ; j<numpoints ; j++)\r
-               {\r
-                       double  v[3];\r
-                       int             k;\r
-\r
-                       // scanf into double, then assign to vec_t\r
-                       // so we don't care what size vec_t is\r
-                       if (fscanf (f, "(%lf %lf %lf ) "\r
-                       , &v[0], &v[1], &v[2]) != 3)\r
-                               Error ("LoadPortals: reading portal %i", i);\r
-                       for (k=0 ; k<3 ; k++)\r
-                               w->points[j][k] = v[k];\r
-               }\r
-               fscanf (f, "\n");\r
-               \r
-               // calc plane\r
-               PlaneFromWinding (w, &plane);\r
-\r
-               // create forward portal\r
-               l = &leafs[leafnums[0]];\r
-               if (l->numportals == MAX_PORTALS_ON_LEAF)\r
-                       Error ("Leaf with too many portals");\r
-               l->portals[l->numportals] = p;\r
-               l->numportals++;\r
-               \r
-               p->num = i+1;\r
-               p->hint = hint;\r
-               p->winding = w;\r
-               VectorSubtract (vec3_origin, plane.normal, p->plane.normal);\r
-               p->plane.dist = -plane.dist;\r
-               p->leaf = leafnums[1];\r
-               SetPortalSphere (p);\r
-               p++;\r
-               \r
-               // create backwards portal\r
-               l = &leafs[leafnums[1]];\r
-               if (l->numportals == MAX_PORTALS_ON_LEAF)\r
-                       Error ("Leaf with too many portals");\r
-               l->portals[l->numportals] = p;\r
-               l->numportals++;\r
-               \r
-               p->num = i+1;\r
-               p->hint = hint;\r
-               p->winding = NewFixedWinding(w->numpoints);\r
-               p->winding->numpoints = w->numpoints;\r
-               for (j=0 ; j<w->numpoints ; j++)\r
-               {\r
-                       VectorCopy (w->points[w->numpoints-1-j], p->winding->points[j]);\r
-               }\r
-\r
-               p->plane = plane;\r
-               p->leaf = leafnums[0];\r
-               SetPortalSphere (p);\r
-               p++;\r
-\r
-       }\r
-\r
-       faces = safe_malloc(2*numfaces*sizeof(vportal_t));\r
-       memset (faces, 0, 2*numfaces*sizeof(vportal_t));\r
-\r
-       faceleafs = safe_malloc(portalclusters*sizeof(leaf_t));\r
-       memset(faceleafs, 0, portalclusters*sizeof(leaf_t));\r
-\r
-       for (i = 0, p = faces; i < numfaces; i++)\r
-       {\r
-               if (fscanf (f, "%i %i ", &numpoints, &leafnums[0]) != 2)\r
-                       Error ("LoadPortals: reading portal %i", i);\r
-\r
-               w = p->winding = NewFixedWinding (numpoints);\r
-               w->numpoints = numpoints;\r
-               \r
-               for (j=0 ; j<numpoints ; j++)\r
-               {\r
-                       double  v[3];\r
-                       int             k;\r
-\r
-                       // scanf into double, then assign to vec_t\r
-                       // so we don't care what size vec_t is\r
-                       if (fscanf (f, "(%lf %lf %lf ) "\r
-                       , &v[0], &v[1], &v[2]) != 3)\r
-                               Error ("LoadPortals: reading portal %i", i);\r
-                       for (k=0 ; k<3 ; k++)\r
-                               w->points[j][k] = v[k];\r
-               }\r
-               fscanf (f, "\n");\r
-               \r
-               // calc plane\r
-               PlaneFromWinding (w, &plane);\r
-\r
-               l = &faceleafs[leafnums[0]];\r
-               l->merged = -1;\r
-               if (l->numportals == MAX_PORTALS_ON_LEAF)\r
-                       Error ("Leaf with too many faces");\r
-               l->portals[l->numportals] = p;\r
-               l->numportals++;\r
-               \r
-               p->num = i+1;\r
-               p->winding = w;\r
-               // normal pointing out of the leaf\r
-               VectorSubtract (vec3_origin, plane.normal, p->plane.normal);\r
-               p->plane.dist = -plane.dist;\r
-               p->leaf = -1;\r
-               SetPortalSphere (p);\r
-               p++;\r
-       }\r
-       \r
-       fclose (f);\r
-}\r
-\r
-\r
-\r
-/*\r
-===========\r
-VisMain\r
-===========\r
-*/\r
-int VisMain (int argc, char **argv)\r
-{\r
-       char            portalfile[1024];\r
-       int                     i;\r
-       \r
-       \r
-       /* note it */\r
-       Sys_Printf( "--- Vis ---\n" );\r
-       \r
-       /* process arguments */\r
-       for (i=1 ; i < (argc - 1) ; i++)\r
-       {\r
-               if (!strcmp(argv[i], "-fast")) {\r
-                       Sys_Printf ("fastvis = true\n");\r
-                       fastvis = qtrue;\r
-               } else if (!strcmp(argv[i], "-merge")) {\r
-                       Sys_Printf ("merge = true\n");\r
-                       mergevis = qtrue;\r
-               } else if (!strcmp(argv[i], "-nopassage")) {\r
-                       Sys_Printf ("nopassage = true\n");\r
-                       noPassageVis = qtrue;\r
-               } else if (!strcmp(argv[i], "-passageOnly")) {\r
-                       Sys_Printf ("passageOnly = true\n");\r
-                       passageVisOnly = qtrue;\r
-               } else if (!strcmp (argv[i],"-nosort")) {\r
-                       Sys_Printf ("nosort = true\n");\r
-                       nosort = qtrue;\r
-               } else if (!strcmp (argv[i],"-saveprt")) {\r
-                       Sys_Printf ("saveprt = true\n");\r
-                       saveprt = qtrue;\r
-               } else if (!strcmp (argv[i],"-tmpin")) {\r
-                       strcpy (inbase, "/tmp");\r
-               } else if (!strcmp (argv[i],"-tmpout")) {\r
-                       strcpy (outbase, "/tmp");\r
-               }\r
-               \r
-       \r
-               /* ydnar: -hint to merge all but hint portals */\r
-               else if( !strcmp( argv[ i ], "-hint" ) )\r
-               {\r
-                       Sys_Printf( "hint = true\n" );\r
-                       hint = qtrue;\r
-                       mergevis = qtrue;\r
-               }\r
-               \r
-               else\r
-                       Sys_Printf( "WARNING: Unknown option \"%s\"\n", argv[ i ] );\r
-       }\r
-\r
-       if( i != argc - 1 )\r
-               Error( "usage: vis [-threads #] [-level 0-4] [-fast] [-v] bspfile" );\r
-       \r
-\r
-       /* load the bsp */\r
-       sprintf( source, "%s%s", inbase, ExpandArg( argv[ i ] ) );\r
-       StripExtension( source );\r
-       strcat( source, ".bsp" );\r
-       Sys_Printf( "Loading %s\n", source );\r
-       LoadBSPFile( source );\r
-       \r
-       /* load the portal file */\r
-       sprintf( portalfile, "%s%s", inbase, ExpandArg( argv[ i ] ) );\r
-       StripExtension( portalfile );\r
-       strcat( portalfile, ".prt" );\r
-       Sys_Printf( "Loading %s\n", portalfile );\r
-       LoadPortals( portalfile );\r
-       \r
-       /* ydnar: for getting far plane */\r
-       ParseEntities();\r
-       \r
-       if( mergevis )\r
-       {\r
-               MergeLeaves();\r
-               MergeLeafPortals();\r
-       }\r
-       \r
-       CountActivePortals();\r
-       /* WritePortals( "maps/hints.prs" );*/\r
-       \r
-       Sys_Printf( "visdatasize:%i\n", numBSPVisBytes );\r
-       \r
-       CalcVis();\r
-       \r
-       /* delete the prt file */\r
-       if( !saveprt )\r
-               remove( portalfile );\r
-\r
-       /* write the bsp file */\r
-       Sys_Printf( "Writing %s\n", source );\r
-       WriteBSPFile( source );\r
-\r
-       return 0;\r
-}\r
+/* -------------------------------------------------------------------------------
+
+   Copyright (C) 1999-2007 id Software, Inc. and contributors.
+   For a list of contributors, see the accompanying CONTRIBUTORS file.
+
+   This file is part of GtkRadiant.
+
+   GtkRadiant is free software; you can redistribute it and/or modify
+   it under the terms of the GNU General Public License as published by
+   the Free Software Foundation; either version 2 of the License, or
+   (at your option) any later version.
+
+   GtkRadiant is distributed in the hope that it will be useful,
+   but WITHOUT ANY WARRANTY; without even the implied warranty of
+   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+   GNU General Public License for more details.
+
+   You should have received a copy of the GNU General Public License
+   along with GtkRadiant; if not, write to the Free Software
+   Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
+
+   ----------------------------------------------------------------------------------
+
+   This code has been altered significantly from its original form, to support
+   several games based on the Quake III Arena engine, in the form of "Q3Map2."
+
+   ------------------------------------------------------------------------------- */
+
+
+
+/* marker */
+#define VIS_C
+
+
+
+/* dependencies */
+#include "q3map2.h"
+
+
+
+
+void PlaneFromWinding( fixedWinding_t *w, visPlane_t *plane ){
+       vec3_t v1, v2;
+
+// calc plane
+       VectorSubtract( w->points[2], w->points[1], v1 );
+       VectorSubtract( w->points[0], w->points[1], v2 );
+       CrossProduct( v2, v1, plane->normal );
+       VectorNormalize( plane->normal, plane->normal );
+       plane->dist = DotProduct( w->points[0], plane->normal );
+}
+
+
+/*
+   NewFixedWinding()
+   returns a new fixed winding
+   ydnar: altered this a bit to reconcile multiply-defined winding_t
+ */
+
+fixedWinding_t *NewFixedWinding( int points ){
+       fixedWinding_t  *w;
+       int size;
+
+       if ( points > MAX_POINTS_ON_WINDING ) {
+               Error( "NewWinding: %i points", points );
+       }
+
+       size = (int)( (size_t)( (fixedWinding_t *)0 )->points[points] );
+       w = safe_malloc( size );
+       memset( w, 0, size );
+
+       return w;
+}
+
+
+
+void prl( leaf_t *l ){
+       int i;
+       vportal_t   *p;
+       visPlane_t pl;
+
+       for ( i = 0 ; i < l->numportals ; i++ )
+       {
+               p = l->portals[i];
+               pl = p->plane;
+               Sys_Printf( "portal %4i to leaf %4i : %7.1f : (%4.1f, %4.1f, %4.1f)\n",(int)( p - portals ),p->leaf,pl.dist, pl.normal[0], pl.normal[1], pl.normal[2] );
+       }
+}
+
+
+//=============================================================================
+
+/*
+   =============
+   SortPortals
+
+   Sorts the portals from the least complex, so the later ones can reuse
+   the earlier information.
+   =============
+ */
+int PComp( const void *a, const void *b ){
+       if ( ( *(const vportal_t *const *)a )->nummightsee == ( *(const vportal_t *const *)b )->nummightsee ) {
+               return 0;
+       }
+       if ( ( *(const vportal_t *const *)a )->nummightsee < ( *(const vportal_t *const *)b )->nummightsee ) {
+               return -1;
+       }
+       return 1;
+}
+void SortPortals( void ){
+       int i;
+
+       for ( i = 0 ; i < numportals * 2 ; i++ )
+               sorted_portals[i] = &portals[i];
+
+       if ( nosort ) {
+               return;
+       }
+       qsort( sorted_portals, numportals * 2, sizeof( sorted_portals[0] ), PComp );
+}
+
+
+/*
+   ==============
+   LeafVectorFromPortalVector
+   ==============
+ */
+int LeafVectorFromPortalVector( byte *portalbits, byte *leafbits ){
+       int i, j, leafnum;
+       vportal_t   *p;
+       int c_leafs;
+
+
+       for ( i = 0 ; i < numportals * 2 ; i++ )
+       {
+               if ( portalbits[i >> 3] & ( 1 << ( i & 7 ) ) ) {
+                       p = portals + i;
+                       leafbits[p->leaf >> 3] |= ( 1 << ( p->leaf & 7 ) );
+               }
+       }
+
+       for ( j = 0; j < portalclusters; j++ )
+       {
+               leafnum = j;
+               while ( leafs[leafnum].merged >= 0 )
+                       leafnum = leafs[leafnum].merged;
+               //if the merged leaf is visible then the original leaf is visible
+               if ( leafbits[leafnum >> 3] & ( 1 << ( leafnum & 7 ) ) ) {
+                       leafbits[j >> 3] |= ( 1 << ( j & 7 ) );
+               }
+       }
+
+       c_leafs = CountBits( leafbits, portalclusters );
+
+       return c_leafs;
+}
+
+
+/*
+   ===============
+   ClusterMerge
+
+   Merges the portal visibility for a leaf
+   ===============
+ */
+static int clustersizehistogram[MAX_MAP_LEAFS] = {0};
+void ClusterMerge( int leafnum ){
+       leaf_t      *leaf;
+       byte portalvector[MAX_PORTALS / 8];
+       byte uncompressed[MAX_MAP_LEAFS / 8];
+       int i, j;
+       int numvis, mergedleafnum;
+       vportal_t   *p;
+       int pnum;
+
+       // OR together all the portalvis bits
+
+       mergedleafnum = leafnum;
+       while ( leafs[mergedleafnum].merged >= 0 )
+               mergedleafnum = leafs[mergedleafnum].merged;
+
+       memset( portalvector, 0, portalbytes );
+       leaf = &leafs[mergedleafnum];
+       for ( i = 0; i < leaf->numportals; i++ )
+       {
+               p = leaf->portals[i];
+               if ( p->removed ) {
+                       continue;
+               }
+
+               if ( p->status != stat_done ) {
+                       Error( "portal not done" );
+               }
+               for ( j = 0 ; j < portallongs ; j++ )
+                       ( (long *)portalvector )[j] |= ( (long *)p->portalvis )[j];
+               pnum = p - portals;
+               portalvector[pnum >> 3] |= 1 << ( pnum & 7 );
+       }
+
+       memset( uncompressed, 0, leafbytes );
+
+       uncompressed[mergedleafnum >> 3] |= ( 1 << ( mergedleafnum & 7 ) );
+       // convert portal bits to leaf bits
+       numvis = LeafVectorFromPortalVector( portalvector, uncompressed );
+
+//     if (uncompressed[leafnum>>3] & (1<<(leafnum&7)))
+//             Sys_Printf ("WARNING: Leaf portals saw into leaf\n");
+
+//     uncompressed[leafnum>>3] |= (1<<(leafnum&7));
+
+       numvis++;       // count the leaf itself
+
+       //Sys_FPrintf( SYS_VRB,"cluster %4i : %4i visible\n", leafnum, numvis );
+       ++clustersizehistogram[numvis];
+
+       memcpy( bspVisBytes + VIS_HEADER_SIZE + leafnum * leafbytes, uncompressed, leafbytes );
+}
+
+/*
+   ==================
+   CalcPortalVis
+   ==================
+ */
+void CalcPortalVis( void ){
+#ifdef MREDEBUG
+       Sys_Printf( "%6d portals out of %d", 0, numportals * 2 );
+       //get rid of the counter
+       RunThreadsOnIndividual( numportals * 2, qfalse, PortalFlow );
+#else
+       RunThreadsOnIndividual( numportals * 2, qtrue, PortalFlow );
+#endif
+
+}
+
+/*
+   ==================
+   CalcPassageVis
+   ==================
+ */
+void CalcPassageVis( void ){
+       PassageMemory();
+
+#ifdef MREDEBUG
+       _printf( "%6d portals out of %d", 0, numportals * 2 );
+       RunThreadsOnIndividual( numportals * 2, qfalse, CreatePassages );
+       _printf( "\n" );
+       _printf( "%6d portals out of %d", 0, numportals * 2 );
+       RunThreadsOnIndividual( numportals * 2, qfalse, PassageFlow );
+       _printf( "\n" );
+#else
+       Sys_Printf( "\n--- CreatePassages (%d) ---\n", numportals * 2 );
+       RunThreadsOnIndividual( numportals * 2, qtrue, CreatePassages );
+
+       Sys_Printf( "\n--- PassageFlow (%d) ---\n", numportals * 2 );
+       RunThreadsOnIndividual( numportals * 2, qtrue, PassageFlow );
+#endif
+}
+
+/*
+   ==================
+   CalcPassagePortalVis
+   ==================
+ */
+void CalcPassagePortalVis( void ){
+       PassageMemory();
+
+#ifdef MREDEBUG
+       Sys_Printf( "%6d portals out of %d", 0, numportals * 2 );
+       RunThreadsOnIndividual( numportals * 2, qfalse, CreatePassages );
+       Sys_Printf( "\n" );
+       Sys_Printf( "%6d portals out of %d", 0, numportals * 2 );
+       RunThreadsOnIndividual( numportals * 2, qfalse, PassagePortalFlow );
+       Sys_Printf( "\n" );
+#else
+       Sys_Printf( "\n--- CreatePassages (%d) ---\n", numportals * 2 );
+       RunThreadsOnIndividual( numportals * 2, qtrue, CreatePassages );
+
+       Sys_Printf( "\n--- PassagePortalFlow (%d) ---\n", numportals * 2 );
+       RunThreadsOnIndividual( numportals * 2, qtrue, PassagePortalFlow );
+#endif
+}
+
+/*
+   ==================
+   CalcFastVis
+   ==================
+ */
+void CalcFastVis( void ){
+       int i;
+
+       // fastvis just uses mightsee for a very loose bound
+       for ( i = 0 ; i < numportals * 2 ; i++ )
+       {
+               portals[i].portalvis = portals[i].portalflood;
+               portals[i].status = stat_done;
+       }
+}
+
+/*
+   ==================
+   CalcVis
+   ==================
+ */
+void CalcVis( void ){
+       int i, minvis, maxvis;
+       const char  *value;
+       double mu, sigma, totalvis, totalvis2;
+
+
+       /* ydnar: rr2do2's farplane code */
+       farPlaneDist = 0.0f;
+       value = ValueForKey( &entities[ 0 ], "_farplanedist" );     /* proper '_' prefixed key */
+       if ( value[ 0 ] == '\0' ) {
+               value = ValueForKey( &entities[ 0 ], "fogclip" );       /* wolf compatibility */
+       }
+       if ( value[ 0 ] == '\0' ) {
+               value = ValueForKey( &entities[ 0 ], "distancecull" );  /* sof2 compatibility */
+       }
+       if ( value[ 0 ] != '\0' ) {
+               farPlaneDist = atof( value );
+               if ( farPlaneDist > 0.0f ) {
+                       Sys_Printf( "farplane distance = %.1f\n", farPlaneDist );
+               }
+               else{
+                       farPlaneDist = 0.0f;
+               }
+       }
+
+
+
+       Sys_Printf( "\n--- BasePortalVis (%d) ---\n", numportals * 2 );
+       RunThreadsOnIndividual( numportals * 2, qtrue, BasePortalVis );
+
+//     RunThreadsOnIndividual (numportals*2, qtrue, BetterPortalVis);
+
+       SortPortals();
+
+       if ( fastvis ) {
+               CalcFastVis();
+       }
+       else if ( noPassageVis ) {
+               CalcPortalVis();
+       }
+       else if ( passageVisOnly ) {
+               CalcPassageVis();
+       }
+       else {
+               CalcPassagePortalVis();
+       }
+       //
+       // assemble the leaf vis lists by oring and compressing the portal lists
+       //
+       Sys_Printf( "creating leaf vis...\n" );
+       for ( i = 0 ; i < portalclusters ; i++ )
+               ClusterMerge( i );
+
+       totalvis = 0;
+       totalvis2 = 0;
+       minvis = -1;
+       maxvis = -1;
+       for ( i = 0; i < MAX_MAP_LEAFS; ++i )
+               if ( clustersizehistogram[i] ) {
+                       if ( debugCluster ) {
+                               Sys_FPrintf( SYS_VRB, "%4i clusters have exactly %4i visible clusters\n", clustersizehistogram[i], i );
+                       }
+                       /* cast is to prevent integer overflow */
+                       totalvis  += ( (double) i )                * ( (double) clustersizehistogram[i] );
+                       totalvis2 += ( (double) i ) * ( (double) i ) * ( (double) clustersizehistogram[i] );
+
+                       if ( minvis < 0 ) {
+                               minvis = i;
+                       }
+                       maxvis = i;
+               }
+
+       mu = totalvis / portalclusters;
+       sigma = sqrt( totalvis2 / portalclusters - mu * mu );
+
+       Sys_Printf( "Total clusters: %i\n", portalclusters );
+       Sys_Printf( "Total visible clusters: %.0f\n", totalvis );
+       Sys_Printf( "Average clusters visible: %.2f (%.3f%%/total)\n", mu, mu / portalclusters * 100.0 );
+       Sys_Printf( "  Standard deviation: %.2f (%.3f%%/total, %.3f%%/avg)\n", sigma, sigma / portalclusters * 100.0, sigma / mu * 100.0 );
+       Sys_Printf( "  Minimum: %i (%.3f%%/total, %.3f%%/avg)\n", minvis, minvis / (double) portalclusters * 100.0, minvis / mu * 100.0 );
+       Sys_Printf( "  Maximum: %i (%.3f%%/total, %.3f%%/avg)\n", maxvis, maxvis / (double) portalclusters * 100.0, maxvis / mu * 100.0 );
+}
+
+/*
+   ==================
+   SetPortalSphere
+   ==================
+ */
+void SetPortalSphere( vportal_t *p ){
+       int i;
+       vec3_t total, dist;
+       fixedWinding_t  *w;
+       float r, bestr;
+
+       w = p->winding;
+       VectorCopy( vec3_origin, total );
+       for ( i = 0 ; i < w->numpoints ; i++ )
+       {
+               VectorAdd( total, w->points[i], total );
+       }
+
+       for ( i = 0 ; i < 3 ; i++ )
+               total[i] /= w->numpoints;
+
+       bestr = 0;
+       for ( i = 0 ; i < w->numpoints ; i++ )
+       {
+               VectorSubtract( w->points[i], total, dist );
+               r = VectorLength( dist );
+               if ( r > bestr ) {
+                       bestr = r;
+               }
+       }
+       VectorCopy( total, p->origin );
+       p->radius = bestr;
+}
+
+/*
+   =============
+   Winding_PlanesConcave
+   =============
+ */
+#define WCONVEX_EPSILON     0.2
+
+int Winding_PlanesConcave( fixedWinding_t *w1, fixedWinding_t *w2,
+                                                  vec3_t normal1, vec3_t normal2,
+                                                  float dist1, float dist2 ){
+       int i;
+
+       if ( !w1 || !w2 ) {
+               return qfalse;
+       }
+
+       // check if one of the points of winding 1 is at the front of the plane of winding 2
+       for ( i = 0; i < w1->numpoints; i++ )
+       {
+               if ( DotProduct( normal2, w1->points[i] ) - dist2 > WCONVEX_EPSILON ) {
+                       return qtrue;
+               }
+       }
+       // check if one of the points of winding 2 is at the front of the plane of winding 1
+       for ( i = 0; i < w2->numpoints; i++ )
+       {
+               if ( DotProduct( normal1, w2->points[i] ) - dist1 > WCONVEX_EPSILON ) {
+                       return qtrue;
+               }
+       }
+
+       return qfalse;
+}
+
+/*
+   ============
+   TryMergeLeaves
+   ============
+ */
+int TryMergeLeaves( int l1num, int l2num ){
+       int i, j, k, n, numportals;
+       visPlane_t plane1, plane2;
+       leaf_t *l1, *l2;
+       vportal_t *p1, *p2;
+       vportal_t *portals[MAX_PORTALS_ON_LEAF];
+
+       for ( k = 0; k < 2; k++ )
+       {
+               if ( k ) {
+                       l1 = &leafs[l1num];
+               }
+               else{l1 = &faceleafs[l1num]; }
+               for ( i = 0; i < l1->numportals; i++ )
+               {
+                       p1 = l1->portals[i];
+                       if ( p1->leaf == l2num ) {
+                               continue;
+                       }
+                       for ( n = 0; n < 2; n++ )
+                       {
+                               if ( n ) {
+                                       l2 = &leafs[l2num];
+                               }
+                               else{l2 = &faceleafs[l2num]; }
+                               for ( j = 0; j < l2->numportals; j++ )
+                               {
+                                       p2 = l2->portals[j];
+                                       if ( p2->leaf == l1num ) {
+                                               continue;
+                                       }
+                                       //
+                                       plane1 = p1->plane;
+                                       plane2 = p2->plane;
+                                       if ( Winding_PlanesConcave( p1->winding, p2->winding, plane1.normal, plane2.normal, plane1.dist, plane2.dist ) ) {
+                                               return qfalse;
+                                       }
+                               }
+                       }
+               }
+       }
+       for ( k = 0; k < 2; k++ )
+       {
+               if ( k ) {
+                       l1 = &leafs[l1num];
+                       l2 = &leafs[l2num];
+               }
+               else
+               {
+                       l1 = &faceleafs[l1num];
+                       l2 = &faceleafs[l2num];
+               }
+               numportals = 0;
+               //the leaves can be merged now
+               for ( i = 0; i < l1->numportals; i++ )
+               {
+                       p1 = l1->portals[i];
+                       if ( p1->leaf == l2num ) {
+                               p1->removed = qtrue;
+                               continue;
+                       }
+                       portals[numportals++] = p1;
+               }
+               for ( j = 0; j < l2->numportals; j++ )
+               {
+                       p2 = l2->portals[j];
+                       if ( p2->leaf == l1num ) {
+                               p2->removed = qtrue;
+                               continue;
+                       }
+                       portals[numportals++] = p2;
+               }
+               for ( i = 0; i < numportals; i++ )
+               {
+                       l2->portals[i] = portals[i];
+               }
+               l2->numportals = numportals;
+               l1->merged = l2num;
+       }
+       return qtrue;
+}
+
+/*
+   ============
+   UpdatePortals
+   ============
+ */
+void UpdatePortals( void ){
+       int i;
+       vportal_t *p;
+
+       for ( i = 0; i < numportals * 2; i++ )
+       {
+               p = &portals[i];
+               if ( p->removed ) {
+                       continue;
+               }
+               while ( leafs[p->leaf].merged >= 0 )
+                       p->leaf = leafs[p->leaf].merged;
+       }
+}
+
+/*
+   ============
+   MergeLeaves
+
+   try to merge leaves but don't merge through hint splitters
+   ============
+ */
+void MergeLeaves( void ){
+       int i, j, nummerges, totalnummerges;
+       leaf_t *leaf;
+       vportal_t *p;
+
+       totalnummerges = 0;
+       do
+       {
+               nummerges = 0;
+               for ( i = 0; i < portalclusters; i++ )
+               {
+                       leaf = &leafs[i];
+                       //if this leaf is merged already
+
+                       /* ydnar: vmods: merge all non-hint portals */
+                       if ( leaf->merged >= 0 && hint == qfalse ) {
+                               continue;
+                       }
+
+
+                       for ( j = 0; j < leaf->numportals; j++ )
+                       {
+                               p = leaf->portals[j];
+                               //
+                               if ( p->removed ) {
+                                       continue;
+                               }
+                               //never merge through hint portals
+                               if ( p->hint ) {
+                                       continue;
+                               }
+                               if ( TryMergeLeaves( i, p->leaf ) ) {
+                                       UpdatePortals();
+                                       nummerges++;
+                                       break;
+                               }
+                       }
+               }
+               totalnummerges += nummerges;
+       } while ( nummerges );
+       Sys_Printf( "%6d leaves merged\n", totalnummerges );
+}
+
+/*
+   ============
+   TryMergeWinding
+   ============
+ */
+#define CONTINUOUS_EPSILON  0.005
+
+fixedWinding_t *TryMergeWinding( fixedWinding_t *f1, fixedWinding_t *f2, vec3_t planenormal ){
+       vec_t       *p1, *p2, *p3, *p4, *back;
+       fixedWinding_t  *newf;
+       int i, j, k, l;
+       vec3_t normal, delta;
+       vec_t dot;
+       qboolean keep1, keep2;
+
+
+       //
+       // find a common edge
+       //
+       p1 = p2 = NULL; // stop compiler warning
+       j = 0;          //
+
+       for ( i = 0; i < f1->numpoints; i++ )
+       {
+               p1 = f1->points[i];
+               p2 = f1->points[( i + 1 ) % f1->numpoints];
+               for ( j = 0; j < f2->numpoints; j++ )
+               {
+                       p3 = f2->points[j];
+                       p4 = f2->points[( j + 1 ) % f2->numpoints];
+                       for ( k = 0; k < 3; k++ )
+                       {
+                               if ( fabs( p1[k] - p4[k] ) > 0.1 ) { //EQUAL_EPSILON) //ME
+                                       break;
+                               }
+                               if ( fabs( p2[k] - p3[k] ) > 0.1 ) { //EQUAL_EPSILON) //ME
+                                       break;
+                               }
+                       } //end for
+                       if ( k == 3 ) {
+                               break;
+                       }
+               } //end for
+               if ( j < f2->numpoints ) {
+                       break;
+               }
+       } //end for
+
+       if ( i == f1->numpoints ) {
+               return NULL;            // no matching edges
+
+       }
+       //
+       // check slope of connected lines
+       // if the slopes are colinear, the point can be removed
+       //
+       back = f1->points[( i + f1->numpoints - 1 ) % f1->numpoints];
+       VectorSubtract( p1, back, delta );
+       CrossProduct( planenormal, delta, normal );
+       VectorNormalize( normal, normal );
+
+       back = f2->points[( j + 2 ) % f2->numpoints];
+       VectorSubtract( back, p1, delta );
+       dot = DotProduct( delta, normal );
+       if ( dot > CONTINUOUS_EPSILON ) {
+               return NULL;            // not a convex polygon
+       }
+       keep1 = (qboolean)( dot < -CONTINUOUS_EPSILON );
+
+       back = f1->points[( i + 2 ) % f1->numpoints];
+       VectorSubtract( back, p2, delta );
+       CrossProduct( planenormal, delta, normal );
+       VectorNormalize( normal, normal );
+
+       back = f2->points[( j + f2->numpoints - 1 ) % f2->numpoints];
+       VectorSubtract( back, p2, delta );
+       dot = DotProduct( delta, normal );
+       if ( dot > CONTINUOUS_EPSILON ) {
+               return NULL;            // not a convex polygon
+       }
+       keep2 = (qboolean)( dot < -CONTINUOUS_EPSILON );
+
+       //
+       // build the new polygon
+       //
+       newf = NewFixedWinding( f1->numpoints + f2->numpoints );
+
+       // copy first polygon
+       for ( k = ( i + 1 ) % f1->numpoints ; k != i ; k = ( k + 1 ) % f1->numpoints )
+       {
+               if ( k == ( i + 1 ) % f1->numpoints && !keep2 ) {
+                       continue;
+               }
+
+               VectorCopy( f1->points[k], newf->points[newf->numpoints] );
+               newf->numpoints++;
+       }
+
+       // copy second polygon
+       for ( l = ( j + 1 ) % f2->numpoints ; l != j ; l = ( l + 1 ) % f2->numpoints )
+       {
+               if ( l == ( j + 1 ) % f2->numpoints && !keep1 ) {
+                       continue;
+               }
+               VectorCopy( f2->points[l], newf->points[newf->numpoints] );
+               newf->numpoints++;
+       }
+
+       return newf;
+}
+
+/*
+   ============
+   MergeLeafPortals
+   ============
+ */
+void MergeLeafPortals( void ){
+       int i, j, k, nummerges, hintsmerged;
+       leaf_t *leaf;
+       vportal_t *p1, *p2;
+       fixedWinding_t *w;
+
+       nummerges = 0;
+       hintsmerged = 0;
+       for ( i = 0; i < portalclusters; i++ )
+       {
+               leaf = &leafs[i];
+               if ( leaf->merged >= 0 ) {
+                       continue;
+               }
+               for ( j = 0; j < leaf->numportals; j++ )
+               {
+                       p1 = leaf->portals[j];
+                       if ( p1->removed ) {
+                               continue;
+                       }
+                       for ( k = j + 1; k < leaf->numportals; k++ )
+                       {
+                               p2 = leaf->portals[k];
+                               if ( p2->removed ) {
+                                       continue;
+                               }
+                               if ( p1->leaf == p2->leaf ) {
+                                       w = TryMergeWinding( p1->winding, p2->winding, p1->plane.normal );
+                                       if ( w ) {
+                                               free( p1->winding );    //% FreeWinding(p1->winding);
+                                               p1->winding = w;
+                                               if ( p1->hint && p2->hint ) {
+                                                       hintsmerged++;
+                                               }
+                                               p1->hint |= p2->hint;
+                                               SetPortalSphere( p1 );
+                                               p2->removed = qtrue;
+                                               nummerges++;
+                                               i--;
+                                               break;
+                                       }
+                               }
+                       }
+                       if ( k < leaf->numportals ) {
+                               break;
+                       }
+               }
+       }
+       Sys_Printf( "%6d portals merged\n", nummerges );
+       Sys_Printf( "%6d hint portals merged\n", hintsmerged );
+}
+
+
+/*
+   ============
+   WritePortals
+   ============
+ */
+int CountActivePortals( void ){
+       int num, hints, j;
+       vportal_t *p;
+
+       num = 0;
+       hints = 0;
+       for ( j = 0; j < numportals * 2; j++ )
+       {
+               p = portals + j;
+               if ( p->removed ) {
+                       continue;
+               }
+               if ( p->hint ) {
+                       hints++;
+               }
+               num++;
+       }
+       Sys_Printf( "%6d active portals\n", num );
+       Sys_Printf( "%6d hint portals\n", hints );
+       return num;
+}
+
+/*
+   ============
+   WritePortals
+   ============
+ */
+void WriteFloat( FILE *f, vec_t v );
+
+void WritePortals( char *portalpathfile ){
+       int i, j, num;
+       FILE *pf;
+       vportal_t *p;
+       fixedWinding_t *w;
+
+       // write the file
+       pf = fopen( portalpathfile, "w" );
+       if ( !pf ) {
+               Error( "Error opening %s", portalpathfile );
+       }
+
+       num = 0;
+       for ( j = 0; j < numportals * 2; j++ )
+       {
+               p = portals + j;
+               if ( p->removed ) {
+                       continue;
+               }
+//             if (!p->hint)
+//                     continue;
+               num++;
+       }
+
+       fprintf( pf, "%s\n", PORTALFILE );
+       fprintf( pf, "%i\n", 0 );
+       fprintf( pf, "%i\n", num ); // + numfaces);
+       fprintf( pf, "%i\n", 0 );
+
+       for ( j = 0; j < numportals * 2; j++ )
+       {
+               p = portals + j;
+               if ( p->removed ) {
+                       continue;
+               }
+//             if (!p->hint)
+//                     continue;
+               w = p->winding;
+               fprintf( pf,"%i %i %i ",w->numpoints, 0, 0 );
+               fprintf( pf, "%d ", p->hint );
+               for ( i = 0 ; i < w->numpoints ; i++ )
+               {
+                       fprintf( pf,"(" );
+                       WriteFloat( pf, w->points[i][0] );
+                       WriteFloat( pf, w->points[i][1] );
+                       WriteFloat( pf, w->points[i][2] );
+                       fprintf( pf,") " );
+               }
+               fprintf( pf,"\n" );
+       }
+
+       /*
+          for (j = 0; j < numfaces; j++)
+          {
+           p = faces + j;
+           w = p->winding;
+           fprintf (pf,"%i %i %i ",w->numpoints, 0, 0);
+           fprintf (pf, "0 ");
+           for (i=0 ; i<w->numpoints ; i++)
+           {
+               fprintf (pf,"(");
+               WriteFloat (pf, w->points[i][0]);
+               WriteFloat (pf, w->points[i][1]);
+               WriteFloat (pf, w->points[i][2]);
+               fprintf (pf,") ");
+           }
+           fprintf (pf,"\n");
+          }*/
+
+       fclose( pf );
+}
+
+/*
+   ============
+   LoadPortals
+   ============
+ */
+void LoadPortals( char *name ){
+       int i, j, flags;
+       vportal_t   *p;
+       leaf_t      *l;
+       char magic[80];
+       FILE        *f;
+       int numpoints;
+       fixedWinding_t  *w;
+       int leafnums[2];
+       visPlane_t plane;
+
+       if ( !strcmp( name,"-" ) ) {
+               f = stdin;
+       }
+       else
+       {
+               f = fopen( name, "r" );
+               if ( !f ) {
+                       Error( "LoadPortals: couldn't read %s\n",name );
+               }
+       }
+
+       if ( fscanf( f,"%79s\n%i\n%i\n%i\n",magic, &portalclusters, &numportals, &numfaces ) != 4 ) {
+               Error( "LoadPortals: failed to read header" );
+       }
+       if ( strcmp( magic,PORTALFILE ) ) {
+               Error( "LoadPortals: not a portal file" );
+       }
+
+       Sys_Printf( "%6i portalclusters\n", portalclusters );
+       Sys_Printf( "%6i numportals\n", numportals );
+       Sys_Printf( "%6i numfaces\n", numfaces );
+
+       if ( numportals > MAX_PORTALS ) {
+               Error( "MAX_PORTALS" );
+       }
+
+       // these counts should take advantage of 64 bit systems automatically
+       leafbytes = ( ( portalclusters + 63 ) & ~63 ) >> 3;
+       leaflongs = leafbytes / sizeof( long );
+
+       portalbytes = ( ( numportals * 2 + 63 ) & ~63 ) >> 3;
+       portallongs = portalbytes / sizeof( long );
+
+       // each file portal is split into two memory portals
+       portals = safe_malloc( 2 * numportals * sizeof( vportal_t ) );
+       memset( portals, 0, 2 * numportals * sizeof( vportal_t ) );
+
+       leafs = safe_malloc( portalclusters * sizeof( leaf_t ) );
+       memset( leafs, 0, portalclusters * sizeof( leaf_t ) );
+
+       for ( i = 0; i < portalclusters; i++ )
+               leafs[i].merged = -1;
+
+       numBSPVisBytes = VIS_HEADER_SIZE + portalclusters * leafbytes;
+
+       if ( numBSPVisBytes > MAX_MAP_VISIBILITY ) {
+               Error( "MAX_MAP_VISIBILITY exceeded" );
+       }
+
+       ( (int *)bspVisBytes )[0] = portalclusters;
+       ( (int *)bspVisBytes )[1] = leafbytes;
+
+       for ( i = 0, p = portals ; i < numportals ; i++ )
+       {
+               if ( fscanf( f, "%i %i %i ", &numpoints, &leafnums[0], &leafnums[1] ) != 3 ) {
+                       Error( "LoadPortals: reading portal %i", i );
+               }
+               if ( numpoints > MAX_POINTS_ON_WINDING ) {
+                       Error( "LoadPortals: portal %i has too many points", i );
+               }
+               if ( leafnums[0] > portalclusters
+                        || leafnums[1] > portalclusters ) {
+                       Error( "LoadPortals: reading portal %i", i );
+               }
+               if ( fscanf( f, "%i ", &flags ) != 1 ) {
+                       Error( "LoadPortals: reading flags" );
+               }
+
+               w = p->winding = NewFixedWinding( numpoints );
+               w->numpoints = numpoints;
+
+               for ( j = 0 ; j < numpoints ; j++ )
+               {
+                       double v[3];
+                       int k;
+
+                       // scanf into double, then assign to vec_t
+                       // so we don't care what size vec_t is
+                       if ( fscanf( f, "(%lf %lf %lf ) "
+                                                , &v[0], &v[1], &v[2] ) != 3 ) {
+                               Error( "LoadPortals: reading portal %i", i );
+                       }
+                       for ( k = 0 ; k < 3 ; k++ )
+                               w->points[j][k] = v[k];
+               }
+               if ( fscanf( f, "\n" ) != 0 ) {
+                       // silence gcc warning
+               }
+
+               // calc plane
+               PlaneFromWinding( w, &plane );
+
+               // create forward portal
+               l = &leafs[leafnums[0]];
+               if ( l->numportals == MAX_PORTALS_ON_LEAF ) {
+                       Error( "Leaf with too many portals" );
+               }
+               l->portals[l->numportals] = p;
+               l->numportals++;
+
+               p->num = i + 1;
+               p->hint = ((flags & 1) != 0);
+               p->sky = ((flags & 2) != 0);
+               p->winding = w;
+               VectorSubtract( vec3_origin, plane.normal, p->plane.normal );
+               p->plane.dist = -plane.dist;
+               p->leaf = leafnums[1];
+               SetPortalSphere( p );
+               p++;
+
+               // create backwards portal
+               l = &leafs[leafnums[1]];
+               if ( l->numportals == MAX_PORTALS_ON_LEAF ) {
+                       Error( "Leaf with too many portals" );
+               }
+               l->portals[l->numportals] = p;
+               l->numportals++;
+
+               p->num = i + 1;
+               p->hint = hint;
+               p->winding = NewFixedWinding( w->numpoints );
+               p->winding->numpoints = w->numpoints;
+               for ( j = 0 ; j < w->numpoints ; j++ )
+               {
+                       VectorCopy( w->points[w->numpoints - 1 - j], p->winding->points[j] );
+               }
+
+               p->plane = plane;
+               p->leaf = leafnums[0];
+               SetPortalSphere( p );
+               p++;
+
+       }
+
+       faces = safe_malloc( 2 * numfaces * sizeof( vportal_t ) );
+       memset( faces, 0, 2 * numfaces * sizeof( vportal_t ) );
+
+       faceleafs = safe_malloc( portalclusters * sizeof( leaf_t ) );
+       memset( faceleafs, 0, portalclusters * sizeof( leaf_t ) );
+
+       for ( i = 0, p = faces; i < numfaces; i++ )
+       {
+               if ( fscanf( f, "%i %i ", &numpoints, &leafnums[0] ) != 2 ) {
+                       Error( "LoadPortals: reading portal %i", i );
+               }
+
+               w = p->winding = NewFixedWinding( numpoints );
+               w->numpoints = numpoints;
+
+               for ( j = 0 ; j < numpoints ; j++ )
+               {
+                       double v[3];
+                       int k;
+
+                       // scanf into double, then assign to vec_t
+                       // so we don't care what size vec_t is
+                       if ( fscanf( f, "(%lf %lf %lf ) "
+                                                , &v[0], &v[1], &v[2] ) != 3 ) {
+                               Error( "LoadPortals: reading portal %i", i );
+                       }
+                       for ( k = 0 ; k < 3 ; k++ )
+                               w->points[j][k] = v[k];
+               }
+               if ( fscanf( f, "\n" ) != 0 ) {
+                       // silence gcc warning
+               }
+
+               // calc plane
+               PlaneFromWinding( w, &plane );
+
+               l = &faceleafs[leafnums[0]];
+               l->merged = -1;
+               if ( l->numportals == MAX_PORTALS_ON_LEAF ) {
+                       Error( "Leaf with too many faces" );
+               }
+               l->portals[l->numportals] = p;
+               l->numportals++;
+
+               p->num = i + 1;
+               p->winding = w;
+               // normal pointing out of the leaf
+               VectorSubtract( vec3_origin, plane.normal, p->plane.normal );
+               p->plane.dist = -plane.dist;
+               p->leaf = -1;
+               SetPortalSphere( p );
+               p++;
+       }
+
+       fclose( f );
+}
+
+
+
+/*
+   ===========
+   VisMain
+   ===========
+ */
+int VisMain( int argc, char **argv ){
+       int i;
+       char portalFilePath[ 1024 ];
+       portalFilePath[0] = 0;
+
+
+       /* note it */
+       Sys_Printf( "--- Vis ---\n" );
+
+       /* process arguments */
+       for ( i = 1 ; i < ( argc - 1 ) ; i++ )
+       {
+               if ( !strcmp( argv[i], "-fast" ) ) {
+                       Sys_Printf( "fastvis = true\n" );
+                       fastvis = qtrue;
+               }
+               else if ( !strcmp( argv[i], "-merge" ) ) {
+                       Sys_Printf( "merge = true\n" );
+                       mergevis = qtrue;
+               }
+               else if ( !strcmp( argv[i], "-mergeportals" ) ) {
+                       Sys_Printf( "mergeportals = true\n" );
+                       mergevisportals = qtrue;
+               }
+               else if ( !strcmp( argv[i], "-nopassage" ) ) {
+                       Sys_Printf( "nopassage = true\n" );
+                       noPassageVis = qtrue;
+               }
+               else if ( !strcmp( argv[i], "-passageOnly" ) ) {
+                       Sys_Printf( "passageOnly = true\n" );
+                       passageVisOnly = qtrue;
+               }
+               else if ( !strcmp( argv[i],"-nosort" ) ) {
+                       Sys_Printf( "nosort = true\n" );
+                       nosort = qtrue;
+               }
+               else if ( !strcmp( argv[i],"-saveprt" ) ) {
+                       Sys_Printf( "saveprt = true\n" );
+                       saveprt = qtrue;
+               }
+               else if ( !strcmp( argv[ i ], "-v" ) ) {
+                       debugCluster = qtrue;
+                       Sys_Printf( "Extra verbous mode enabled\n" );
+               }
+               else if ( !strcmp( argv[i],"-tmpin" ) ) {
+                       strcpy( inbase, "/tmp" );
+               }
+               else if ( !strcmp( argv[i],"-tmpout" ) ) {
+                       strcpy( outbase, "/tmp" );
+               }
+
+               /* ydnar: -hint to merge all but hint portals */
+               else if ( !strcmp( argv[ i ], "-hint" ) ) {
+                       Sys_Printf( "hint = true\n" );
+                       hint = qtrue;
+                       mergevis = qtrue;
+               }
+               else if ( !strcmp( argv[ i ], "-prtfile" ) )
+               {
+                       strcpy( portalFilePath, argv[i + 1] );
+                       i++;
+                       Sys_Printf( "Use %s as portal file\n", portalFilePath );
+               }
+
+               else{
+                       Sys_FPrintf( SYS_WRN, "WARNING: Unknown option \"%s\"\n", argv[ i ] );
+               }
+       }
+
+       if ( i != argc - 1 ) {
+               Error( "usage: vis [-threads #] [-fast] [-v] BSPFilePath" );
+       }
+
+
+       /* load the bsp */
+       sprintf( source, "%s%s", inbase, ExpandArg( argv[ i ] ) );
+       StripExtension( source );
+       strcat( source, ".bsp" );
+       Sys_Printf( "Loading %s\n", source );
+       LoadBSPFile( source );
+
+       /* load the portal file */
+       if (!portalFilePath[0]) {
+               sprintf( portalFilePath, "%s%s", inbase, ExpandArg( argv[ i ] ) );
+               StripExtension( portalFilePath );
+               strcat( portalFilePath, ".prt" );
+       }
+       Sys_Printf( "Loading %s\n", portalFilePath );
+       LoadPortals( portalFilePath );
+
+       /* ydnar: exit if no portals, hence no vis */
+       if ( numportals == 0 ) {
+               Sys_Printf( "No portals means no vis, exiting.\n" );
+               return 0;
+       }
+
+       /* ydnar: for getting far plane */
+       ParseEntities();
+
+       /* inject command line parameters */
+       InjectCommandLine( argv, 0, argc - 1 );
+       UnparseEntities();
+
+       if ( mergevis ) {
+               MergeLeaves();
+       }
+
+       if ( mergevis || mergevisportals ) {
+               MergeLeafPortals();
+       }
+
+       CountActivePortals();
+       /* WritePortals( "maps/hints.prs" );*/
+
+       Sys_Printf( "visdatasize:%i\n", numBSPVisBytes );
+
+       CalcVis();
+
+       /* delete the prt file */
+       if ( !saveprt ) {
+               remove( portalFilePath );
+       }
+
+       /* write the bsp file */
+       Sys_Printf( "Writing %s\n", source );
+       WriteBSPFile( source );
+
+       return 0;
+}