--- /dev/null
+/*
+GenSurf plugin for GtkRadiant
+Copyright (C) 2001 David Hyde, Loki software and qeradiant.com
+
+This library is free software; you can redistribute it and/or
+modify it under the terms of the GNU Lesser General Public
+License as published by the Free Software Foundation; either
+version 2.1 of the License, or (at your option) any later version.
+
+This library 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
+Lesser General Public License for more details.
+
+You should have received a copy of the GNU Lesser General Public
+License along with this library; if not, write to the Free Software
+Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+*/
+
+#include <stdlib.h>
+#include <math.h>
+#include "gensurf.h"
+
+#define MAX_FACES 128 // Maximum number of faces on a brush
+#define MAX_POINTS_ON_WINDING 64
+#define SIDE_FRONT 0
+#define SIDE_ON 2
+#define SIDE_BACK 1
+#define SIDE_CROSS -2
+
+//vec3 vec3_origin = {0,0,0};
+
+void PlaneFromPoints (float *p0, float *p1, float *p2, PLANE *plane)
+{
+ vec3 t1, t2;
+ vec length;
+
+ VectorSubtract (p0, p1, t1);
+ VectorSubtract (p2, p1, t2);
+ plane->normal[0] = t1[1]*t2[2] - t1[2]*t2[1];
+ plane->normal[1] = t1[2]*t2[0] - t1[0]*t2[2];
+ plane->normal[2] = t1[0]*t2[1] - t1[1]*t2[0];
+
+ length = (vec)(sqrt(plane->normal[0]*plane->normal[0] +
+ plane->normal[1]*plane->normal[1] +
+ plane->normal[2]*plane->normal[2] ));
+ if (length == 0)
+ {
+ VectorClear(plane->normal);
+ }
+ else
+ {
+ plane->normal[0] /= length;
+ plane->normal[1] /= length;
+ plane->normal[2] /= length;
+ }
+ plane->dist = DotProduct (p0, plane->normal);
+}
+/*
+void VectorMA (vec3 va, vec scale, vec3 vb, vec3 vc)
+{
+ vc[0] = va[0] + scale*vb[0];
+ vc[1] = va[1] + scale*vb[1];
+ vc[2] = va[2] + scale*vb[2];
+}
+
+void CrossProduct (vec3 v1, vec3 v2, vec3 cross)
+{
+ cross[0] = v1[1]*v2[2] - v1[2]*v2[1];
+ cross[1] = v1[2]*v2[0] - v1[0]*v2[2];
+ cross[2] = v1[0]*v2[1] - v1[1]*v2[0];
+}
+*/
+/*
+=============
+AllocWinding
+=============
+*/
+MY_WINDING *AllocWinding (int points)
+{
+ MY_WINDING *w;
+ int s;
+
+ s = sizeof(vec)*3*points + sizeof(int);
+ w = (MY_WINDING*)malloc (s);
+ memset (w, 0, s);
+ return w;
+}
+/*
+vec VectorNormalize (vec3 in, vec3 out)
+{
+ vec length, ilength;
+
+ length = (vec)(sqrt (in[0]*in[0] + in[1]*in[1] + in[2]*in[2]));
+ if (length == 0)
+ {
+ VectorClear (out);
+ return 0;
+ }
+
+ ilength = (vec)1.0/length;
+ out[0] = in[0]*ilength;
+ out[1] = in[1]*ilength;
+ out[2] = in[2]*ilength;
+
+ return length;
+}
+*/
+
+/*
+=================
+BaseWindingForPlane
+=================
+*/
+MY_WINDING *BaseWindingForPlane (vec3 normal, vec dist)
+{
+ int i, x;
+ vec max, v;
+ vec3 org, vright, vup;
+ MY_WINDING *w;
+
+// find the major axis
+
+ max = -BOGUS_RANGE;
+ x = -1;
+ for (i=0 ; i<3; i++)
+ {
+ v = (vec)(fabs(normal[i]));
+ if (v > max)
+ {
+ x = i;
+ max = v;
+ }
+ }
+ if (x==-1) x = 2;
+
+ VectorCopy(vec3_origin,vup);
+ switch (x)
+ {
+ case 0:
+ case 1:
+ vup[2] = 1;
+ break;
+ case 2:
+ vup[0] = 1;
+ break;
+ }
+
+ v = DotProduct (vup, normal);
+ VectorMA (vup, -v, normal, vup);
+ VectorNormalize (vup, vup);
+
+ VectorScale (normal, dist, org);
+
+ CrossProduct (vup, normal, vright);
+
+ VectorScale (vup, 65536, vup);
+ VectorScale (vright, 65536, vright);
+
+// project a really big axis aligned box onto the plane
+ w = AllocWinding (4);
+
+ VectorSubtract (org, vright, w->p[0]);
+ VectorAdd (w->p[0], vup, w->p[0]);
+
+ VectorAdd (org, vright, w->p[1]);
+ VectorAdd (w->p[1], vup, w->p[1]);
+
+ VectorAdd (org, vright, w->p[2]);
+ VectorSubtract (w->p[2], vup, w->p[2]);
+
+ VectorSubtract (org, vright, w->p[3]);
+ VectorSubtract (w->p[3], vup, w->p[3]);
+
+ w->numpoints = 4;
+
+ return w;
+}
+
+void FreeWinding (MY_WINDING *w)
+{
+ if (*(unsigned *)w == 0xdeaddead)
+// Error ("FreeWinding: freed a freed winding");
+ return;
+ *(unsigned *)w = 0xdeaddead;
+
+ free (w);
+}
+
+/*
+=============
+ChopWindingInPlace
+=============
+*/
+void ChopWindingInPlace (MY_WINDING **inout, vec3 normal, vec dist, vec epsilon)
+{
+ MY_WINDING *in;
+ vec dists[MAX_POINTS_ON_WINDING+4];
+ int sides[MAX_POINTS_ON_WINDING+4];
+ int counts[3];
+ static vec dot; // VC 4.2 optimizer bug if not static
+ int i, j;
+ vec *p1, *p2;
+ vec3 mid;
+ MY_WINDING *f;
+ int maxpts;
+
+ in = *inout;
+ counts[0] = counts[1] = counts[2] = 0;
+
+// determine sides for each point
+ for (i=0 ; i<in->numpoints ; i++)
+ {
+ dot = DotProduct (in->p[i], normal);
+ dot -= dist;
+ dists[i] = dot;
+ if (dot > epsilon)
+ sides[i] = SIDE_FRONT;
+ else if (dot < -epsilon)
+ sides[i] = SIDE_BACK;
+ else
+ {
+ sides[i] = SIDE_ON;
+ }
+ counts[sides[i]]++;
+ }
+ sides[i] = sides[0];
+ dists[i] = dists[0];
+
+ if (!counts[0])
+ {
+ FreeWinding(in);
+ *inout = NULL;
+ return;
+ }
+ if (!counts[1])
+ return; // inout stays the same
+
+ maxpts = in->numpoints+4; // cant use counts[0]+2 because
+ // of fp grouping errors
+
+ f = AllocWinding (maxpts);
+
+ for (i=0 ; i<in->numpoints ; i++)
+ {
+ p1 = in->p[i];
+
+ if (sides[i] == SIDE_ON)
+ {
+ VectorCopy (p1, f->p[f->numpoints]);
+ f->numpoints++;
+ continue;
+ }
+
+ if (sides[i] == SIDE_FRONT)
+ {
+ VectorCopy (p1, f->p[f->numpoints]);
+ f->numpoints++;
+ }
+
+ if (sides[i+1] == SIDE_ON || sides[i+1] == sides[i])
+ continue;
+
+ // generate a split point
+ p2 = in->p[(i+1)%in->numpoints];
+
+ dot = dists[i] / (dists[i]-dists[i+1]);
+ for (j=0 ; j<3 ; j++)
+ { // avoid round off error when possible
+ if (normal[j] == 1)
+ mid[j] = dist;
+ else if (normal[j] == -1)
+ mid[j] = -dist;
+ else
+ mid[j] = p1[j] + dot*(p2[j]-p1[j]);
+ }
+
+ VectorCopy (mid, f->p[f->numpoints]);
+ f->numpoints++;
+ }
+
+// if (f->numpoints > maxpts)
+// Error ("ClipWinding: points exceeded estimate");
+// if (f->numpoints > MAX_POINTS_ON_WINDING)
+// Error ("ClipWinding: MAX_POINTS_ON_WINDING");
+
+ FreeWinding(in);
+ *inout = f;
+}
+
+void UseFaceBounds()
+{
+ LPVOID vp;
+ float Dot, BestDot;
+ float planepts[3][3];
+ int BestFace;
+ int i, j;
+ int NumFaces;
+ vec3 SurfNormal;
+ vec3 vmin,vmax;
+ PLANE plane[MAX_FACES*2];
+ PLANE pface;
+ MY_WINDING *w;
+
+ switch(Plane)
+ {
+ case PLANE_XY1:
+ SurfNormal[0] = 0.0;
+ SurfNormal[1] = 0.0;
+ SurfNormal[2] =-1.0;
+ break;
+ case PLANE_XZ0:
+ SurfNormal[0] = 0.0;
+ SurfNormal[1] = 1.0;
+ SurfNormal[2] = 0.0;
+ break;
+ case PLANE_XZ1:
+ SurfNormal[0] = 0.0;
+ SurfNormal[1] =-1.0;
+ SurfNormal[2] = 0.0;
+ break;
+ case PLANE_YZ0:
+ SurfNormal[0] = 1.0;
+ SurfNormal[1] = 0.0;
+ SurfNormal[2] = 0.0;
+ break;
+ case PLANE_YZ1:
+ SurfNormal[0] =-1.0;
+ SurfNormal[1] = 0.0;
+ SurfNormal[2] = 0.0;
+ break;
+ default:
+ SurfNormal[0] = 0.0;
+ SurfNormal[1] = 0.0;
+ SurfNormal[2] = 1.0;
+ }
+
+#if 0
+ i = g_FuncTable.m_pfnAllocateSelectedBrushHandles();
+ vp = g_FuncTable.m_pfnGetSelectedBrushHandle(0);
+ NumFaces = g_FuncTable.m_pfnGetFaceCount(vp);
+
+ BestFace = -1;
+ BestDot = 0.0;
+
+ for(i=0; i<NumFaces; i++)
+ {
+ _QERFaceData* QERFaceData = g_FuncTable.m_pfnGetFaceData(vp,i);
+ planepts[0][0] = QERFaceData->m_v1[0];
+ planepts[0][1] = QERFaceData->m_v1[1];
+ planepts[0][2] = QERFaceData->m_v1[2];
+ planepts[1][0] = QERFaceData->m_v2[0];
+ planepts[1][1] = QERFaceData->m_v2[1];
+ planepts[1][2] = QERFaceData->m_v2[2];
+ planepts[2][0] = QERFaceData->m_v3[0];
+ planepts[2][1] = QERFaceData->m_v3[1];
+ planepts[2][2] = QERFaceData->m_v3[2];
+
+ PlaneFromPoints (planepts[0], planepts[1], planepts[2], &plane[2*i]);
+ VectorSubtract (vec3_origin, plane[2*i].normal, plane[2*i+1].normal);
+ plane[2*i+1].dist = -plane[2*i].dist;
+
+ Dot = DotProduct(plane[2*i].normal,SurfNormal);
+ if(Dot > BestDot)
+ {
+ BestDot = Dot;
+ BestFace = i;
+ if(strlen(QERFaceData->m_TextureName))
+ strcpy(Texture[Game][0],QERFaceData->m_TextureName);
+ }
+ }
+ for(i=0; i<NumFaces; i++)
+ {
+ if(i==BestFace) continue;
+ _QERFaceData* QERFaceData = g_FuncTable.m_pfnGetFaceData(vp,i);
+ if(strlen(QERFaceData->m_TextureName))
+ {
+ if(strcmp(Texture[Game][0],QERFaceData->m_TextureName))
+ strcpy(Texture[Game][1],QERFaceData->m_TextureName);
+ }
+ }
+
+
+ g_FuncTable.m_pfnReleaseSelectedBrushHandles();
+
+ w = BaseWindingForPlane (plane[BestFace*2].normal, plane[BestFace*2].dist);
+
+ for (i=0 ; i<NumFaces && w; i++)
+ {
+ if (BestFace == i)
+ continue;
+ ChopWindingInPlace (&w, plane[i*2+1].normal, plane[i*2+1].dist, 0);
+ }
+ if(!w) return;
+
+ // Get bounding box for this face
+ vmin[0] = vmax[0] = w->p[0][0];
+ vmin[1] = vmax[1] = w->p[0][1];
+ vmin[2] = vmax[2] = w->p[0][2];
+ for(j=1; j<w->numpoints; j++)
+ {
+ vmin[0] = min(vmin[0],w->p[j][0]);
+ vmin[1] = min(vmin[1],w->p[j][1]);
+ vmin[2] = min(vmin[2],w->p[j][2]);
+ vmax[0] = max(vmax[0],w->p[j][0]);
+ vmax[1] = max(vmax[1],w->p[j][1]);
+ vmax[2] = max(vmax[2],w->p[j][2]);
+ }
+
+ FreeWinding(w);
+
+ VectorCopy(plane[BestFace*2].normal,pface.normal);
+ pface.dist = plane[BestFace*2].dist;
+ switch(Plane)
+ {
+ case PLANE_XZ0:
+ case PLANE_XZ1:
+ if(pface.normal[1] == 0.) return;
+ Hll = vmin[0];
+ Hur = vmax[0];
+ Vll = vmin[2];
+ Vur = vmax[2];
+ Z00 = (pface.dist - pface.normal[0]*Hll - pface.normal[2]*Vll)/pface.normal[1];
+ Z01 = (pface.dist - pface.normal[0]*Hll - pface.normal[2]*Vur)/pface.normal[1];
+ Z10 = (pface.dist - pface.normal[0]*Hur - pface.normal[2]*Vll)/pface.normal[1];
+ Z11 = (pface.dist - pface.normal[0]*Hur - pface.normal[2]*Vur)/pface.normal[1];
+ break;
+ case PLANE_YZ0:
+ case PLANE_YZ1:
+ if(pface.normal[0] == 0.) return;
+ Hll = vmin[1];
+ Hur = vmax[1];
+ Vll = vmin[2];
+ Vur = vmax[2];
+ Z00 = (pface.dist - pface.normal[1]*Hll - pface.normal[2]*Vll)/pface.normal[0];
+ Z01 = (pface.dist - pface.normal[1]*Hll - pface.normal[2]*Vur)/pface.normal[0];
+ Z10 = (pface.dist - pface.normal[1]*Hur - pface.normal[2]*Vll)/pface.normal[0];
+ Z11 = (pface.dist - pface.normal[1]*Hur - pface.normal[2]*Vur)/pface.normal[0];
+ break;
+ default:
+ if(pface.normal[2] == 0.) return;
+ Hll = vmin[0];
+ Hur = vmax[0];
+ Vll = vmin[1];
+ Vur = vmax[1];
+ Z00 = (pface.dist - pface.normal[0]*Hll - pface.normal[1]*Vll)/pface.normal[2];
+ Z01 = (pface.dist - pface.normal[0]*Hll - pface.normal[1]*Vur)/pface.normal[2];
+ Z10 = (pface.dist - pface.normal[0]*Hur - pface.normal[1]*Vll)/pface.normal[2];
+ Z11 = (pface.dist - pface.normal[0]*Hur - pface.normal[1]*Vur)/pface.normal[2];
+ }
+#endif
+}