1 /* -------------------------------------------------------------------------------
3 Copyright (C) 1999-2007 id Software, Inc. and contributors.
4 For a list of contributors, see the accompanying CONTRIBUTORS file.
6 This file is part of GtkRadiant.
8 GtkRadiant is free software; you can redistribute it and/or modify
9 it under the terms of the GNU General Public License as published by
10 the Free Software Foundation; either version 2 of the License, or
11 (at your option) any later version.
13 GtkRadiant is distributed in the hope that it will be useful,
14 but WITHOUT ANY WARRANTY; without even the implied warranty of
15 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 GNU General Public License for more details.
18 You should have received a copy of the GNU General Public License
19 along with GtkRadiant; if not, write to the Free Software
20 Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
22 ----------------------------------------------------------------------------------
24 This code has been altered significantly from its original form, to support
25 several games based on the Quake III Arena engine, in the form of "Q3Map2."
27 ------------------------------------------------------------------------------- */
32 #define SURFACE_META_C
41 #define LIGHTMAP_EXCEEDED -1
44 #define ST_EXCEEDED -4
45 #define UNSUITABLE_TRIANGLE -10
46 #define VERTS_EXCEEDED -1000
47 #define INDEXES_EXCEEDED -2000
49 #define GROW_META_VERTS 1024
50 #define GROW_META_TRIANGLES 1024
52 static int numMetaSurfaces, numPatchMetaSurfaces;
54 static int maxMetaVerts = 0;
55 static int numMetaVerts = 0;
56 static int firstSearchMetaVert = 0;
57 static bspDrawVert_t *metaVerts = NULL;
59 static int maxMetaTriangles = 0;
60 static int numMetaTriangles = 0;
61 static metaTriangle_t *metaTriangles = NULL;
67 called before staring a new entity to clear out the triangle list
70 void ClearMetaTriangles( void )
80 finds a matching metavertex in the global list, returning its index
83 static int FindMetaVertex( bspDrawVert_t *src )
86 bspDrawVert_t *v, *temp;
89 /* try to find an existing drawvert */
90 for( i = firstSearchMetaVert, v = &metaVerts[ i ]; i < numMetaVerts; i++, v++ )
92 if( memcmp( src, v, sizeof( bspDrawVert_t ) ) == 0 )
97 if( numMetaVerts >= maxMetaVerts )
99 /* reallocate more room */
100 maxMetaVerts += GROW_META_VERTS;
101 temp = safe_malloc( maxMetaVerts * sizeof( bspDrawVert_t ) );
102 if( metaVerts != NULL )
104 memcpy( temp, metaVerts, numMetaVerts * sizeof( bspDrawVert_t ) );
110 /* add the triangle */
111 memcpy( &metaVerts[ numMetaVerts ], src, sizeof( bspDrawVert_t ) );
114 /* return the count */
115 return (numMetaVerts - 1);
122 adds a new meta triangle, allocating more memory if necessary
125 static int AddMetaTriangle( void )
127 metaTriangle_t *temp;
131 if( numMetaTriangles >= maxMetaTriangles )
133 /* reallocate more room */
134 maxMetaTriangles += GROW_META_TRIANGLES;
135 temp = safe_malloc( maxMetaTriangles * sizeof( metaTriangle_t ) );
136 if( metaTriangles != NULL )
138 memcpy( temp, metaTriangles, numMetaTriangles * sizeof( metaTriangle_t ) );
139 free( metaTriangles );
141 metaTriangles = temp;
144 /* increment and return */
146 return numMetaTriangles - 1;
153 finds a matching metatriangle in the global list,
154 otherwise adds it and returns the index to the metatriangle
157 int FindMetaTriangle( metaTriangle_t *src, bspDrawVert_t *a, bspDrawVert_t *b, bspDrawVert_t *c, int planeNum )
164 /* detect degenerate triangles fixme: do something proper here */
165 VectorSubtract( a->xyz, b->xyz, dir );
166 if( VectorLength( dir ) < 0.125f )
168 VectorSubtract( b->xyz, c->xyz, dir );
169 if( VectorLength( dir ) < 0.125f )
171 VectorSubtract( c->xyz, a->xyz, dir );
172 if( VectorLength( dir ) < 0.125f )
178 /* because of precision issues with small triangles, try to use the specified plane */
179 src->planeNum = planeNum;
180 VectorCopy( mapplanes[ planeNum ].normal, src->plane );
181 src->plane[ 3 ] = mapplanes[ planeNum ].dist;
185 /* calculate a plane from the triangle's points (and bail if a plane can't be constructed) */
187 if( PlaneFromPoints( src->plane, a->xyz, b->xyz, c->xyz ) == qfalse )
191 /* ydnar 2002-10-03: repair any bogus normals (busted ase import kludge) */
192 if( VectorLength( a->normal ) <= 0.0f )
193 VectorCopy( src->plane, a->normal );
194 if( VectorLength( b->normal ) <= 0.0f )
195 VectorCopy( src->plane, b->normal );
196 if( VectorLength( c->normal ) <= 0.0f )
197 VectorCopy( src->plane, c->normal );
199 /* ydnar 2002-10-04: set lightmap axis if not already set */
200 if( !(src->si->compileFlags & C_VERTEXLIT) &&
201 src->lightmapAxis[ 0 ] == 0.0f && src->lightmapAxis[ 1 ] == 0.0f && src->lightmapAxis[ 2 ] == 0.0f )
203 /* the shader can specify an explicit lightmap axis */
204 if( src->si->lightmapAxis[ 0 ] || src->si->lightmapAxis[ 1 ] || src->si->lightmapAxis[ 2 ] )
205 VectorCopy( src->si->lightmapAxis, src->lightmapAxis );
207 /* new axis-finding code */
209 CalcLightmapAxis( src->plane, src->lightmapAxis );
212 /* fill out the src triangle */
213 src->indexes[ 0 ] = FindMetaVertex( a );
214 src->indexes[ 1 ] = FindMetaVertex( b );
215 src->indexes[ 2 ] = FindMetaVertex( c );
217 /* try to find an existing triangle */
218 #ifdef USE_EXHAUSTIVE_SEARCH
224 for( i = 0, tri = metaTriangles; i < numMetaTriangles; i++, tri++ )
226 if( memcmp( src, tri, sizeof( metaTriangle_t ) ) == 0 )
232 /* get a new triangle */
233 triIndex = AddMetaTriangle();
235 /* add the triangle */
236 memcpy( &metaTriangles[ triIndex ], src, sizeof( metaTriangle_t ) );
238 /* return the triangle index */
245 SurfaceToMetaTriangles()
246 converts a classified surface to metatriangles
249 static void SurfaceToMetaTriangles( mapDrawSurface_t *ds )
253 bspDrawVert_t a, b, c;
256 /* only handle certain types of surfaces */
257 if( ds->type != SURFACE_FACE &&
258 ds->type != SURFACE_META &&
259 ds->type != SURFACE_FORCED_META &&
260 ds->type != SURFACE_DECAL )
263 /* speed at the expense of memory */
264 firstSearchMetaVert = numMetaVerts;
266 /* only handle valid surfaces */
267 if( ds->type != SURFACE_BAD && ds->numVerts >= 3 && ds->numIndexes >= 3 )
269 /* walk the indexes and create triangles */
270 for( i = 0; i < ds->numIndexes; i += 3 )
272 /* sanity check the indexes */
273 if( ds->indexes[ i ] == ds->indexes[ i + 1 ] ||
274 ds->indexes[ i ] == ds->indexes[ i + 2 ] ||
275 ds->indexes[ i + 1 ] == ds->indexes[ i + 2 ] )
277 //% Sys_Printf( "%d! ", ds->numVerts );
281 /* build a metatriangle */
282 src.si = ds->shaderInfo;
283 src.side = (ds->sideRef != NULL ? ds->sideRef->side : NULL);
284 src.entityNum = ds->entityNum;
285 src.surfaceNum = ds->surfaceNum;
286 src.planeNum = ds->planeNum;
287 src.castShadows = ds->castShadows;
288 src.recvShadows = ds->recvShadows;
289 src.fogNum = ds->fogNum;
290 src.sampleSize = ds->sampleSize;
291 src.shadeAngleDegrees = ds->shadeAngleDegrees;
292 VectorCopy( ds->lightmapAxis, src.lightmapAxis );
295 memcpy( &a, &ds->verts[ ds->indexes[ i ] ], sizeof( a ) );
296 memcpy( &b, &ds->verts[ ds->indexes[ i + 1 ] ], sizeof( b ) );
297 memcpy( &c, &ds->verts[ ds->indexes[ i + 2 ] ], sizeof( c ) );
298 FindMetaTriangle( &src, &a, &b, &c, ds->planeNum );
305 /* clear the surface (free verts and indexes, sets it to SURFACE_BAD) */
312 TriangulatePatchSurface()
313 creates triangles from a patch
316 void TriangulatePatchSurface( entity_t *e , mapDrawSurface_t *ds )
318 int iterations, x, y, pw[ 5 ], r;
319 mapDrawSurface_t *dsNew;
320 mesh_t src, *subdivided, *mesh;
323 int patchSubdivision;
325 /* vortex: _patchMeta, _patchQuality, _patchSubdivide support */
326 forcePatchMeta = IntForKey(e, "_patchMeta" );
328 forcePatchMeta = IntForKey(e, "patchMeta" );
329 patchQuality = IntForKey(e, "_patchQuality" );
331 patchQuality = IntForKey(e, "patchQuality" );
334 patchSubdivision = IntForKey(e, "_patchSubdivide" );
335 if (!patchSubdivision)
336 patchSubdivision = IntForKey(e, "patchSubdivide" );
338 /* try to early out */
339 if(ds->numVerts == 0 || ds->type != SURFACE_PATCH || ( patchMeta == qfalse && !forcePatchMeta) )
341 /* make a mesh from the drawsurf */
342 src.width = ds->patchWidth;
343 src.height = ds->patchHeight;
344 src.verts = ds->verts;
345 //% subdivided = SubdivideMesh( src, 8, 999 );
346 if (patchSubdivision)
347 iterations = IterationsForCurve( ds->longestCurve, patchSubdivision );
349 iterations = IterationsForCurve( ds->longestCurve, patchSubdivisions / patchQuality );
351 subdivided = SubdivideMesh2( src, iterations ); //% ds->maxIterations
353 /* fit it to the curve and remove colinear verts on rows/columns */
354 PutMeshOnCurve( *subdivided );
355 mesh = RemoveLinearMeshColumnsRows( subdivided );
356 FreeMesh( subdivided );
357 //% MakeMeshNormals( mesh );
359 /* make a copy of the drawsurface */
360 dsNew = AllocDrawSurface( SURFACE_META );
361 memcpy( dsNew, ds, sizeof( *ds ) );
363 /* if the patch is nonsolid, then discard it */
364 if( !(ds->shaderInfo->compileFlags & C_SOLID) )
367 /* set new pointer */
370 /* basic transmogrification */
371 ds->type = SURFACE_META;
373 ds->indexes = safe_malloc( mesh->width * mesh->height * 6 * sizeof( int ) );
375 /* copy the verts in */
376 ds->numVerts = (mesh->width * mesh->height);
377 ds->verts = mesh->verts;
379 /* iterate through the mesh quads */
380 for( y = 0; y < (mesh->height - 1); y++ )
382 for( x = 0; x < (mesh->width - 1); x++ )
385 pw[ 0 ] = x + (y * mesh->width);
386 pw[ 1 ] = x + ((y + 1) * mesh->width);
387 pw[ 2 ] = x + 1 + ((y + 1) * mesh->width);
388 pw[ 3 ] = x + 1 + (y * mesh->width);
389 pw[ 4 ] = x + (y * mesh->width); /* same as pw[ 0 ] */
394 /* make first triangle */
395 ds->indexes[ ds->numIndexes++ ] = pw[ r + 0 ];
396 ds->indexes[ ds->numIndexes++ ] = pw[ r + 1 ];
397 ds->indexes[ ds->numIndexes++ ] = pw[ r + 2 ];
399 /* make second triangle */
400 ds->indexes[ ds->numIndexes++ ] = pw[ r + 0 ];
401 ds->indexes[ ds->numIndexes++ ] = pw[ r + 2 ];
402 ds->indexes[ ds->numIndexes++ ] = pw[ r + 3 ];
406 /* free the mesh, but not the verts */
410 numPatchMetaSurfaces++;
413 ClassifySurfaces( 1, ds );
416 #define TINY_AREA 1.0f
417 int MaxAreaIndexes(bspDrawVert_t *vert, int cnt, int *indexes)
419 int r, s, t, bestR = 0, bestS = 1, bestT = 2;
421 double A, bestA = -1;
422 vec3_t ab, ac, cross;
428 /* find the triangle with highest area */
429 for(r = 0; r+2 < cnt; ++r)
430 for(s = r+1; s+1 < cnt; ++s)
431 for(t = s+1; t < cnt; ++t)
433 VectorSubtract(vert[s].xyz, vert[r].xyz, ab);
434 VectorSubtract(vert[t].xyz, vert[r].xyz, ac);
435 CrossProduct(ab, ac, cross);
436 A = VectorLength(cross);
446 if(bestA < TINY_AREA)
447 /* the biggest triangle is degenerate - then every other is too, and the other algorithms wouldn't generate anything useful either */
451 indexes[i++] = bestR;
452 indexes[i++] = bestS;
453 indexes[i++] = bestT;
456 /* identify the other fragments */
458 /* full polygon without triangle (bestR,bestS,bestT) = three new polygons:
464 j = i + MaxAreaIndexes(vert + bestR, bestS - bestR + 1, indexes + i);
467 /* uses 3*(bestS-bestR+1)-6 */
468 j = i + MaxAreaIndexes(vert + bestS, bestT - bestS + 1, indexes + i);
471 /* uses 3*(bestT-bestS+1)-6 */
473 /* can'bestT recurse this one directly... therefore, buffering */
474 if(cnt + bestR - bestT + 1 >= 3)
476 buf = safe_malloc(sizeof(*vert) * (cnt + bestR - bestT + 1));
477 memcpy(buf, vert + bestT, sizeof(*vert) * (cnt - bestT));
478 memcpy(buf + (cnt - bestT), vert, sizeof(*vert) * (bestR + 1));
479 j = i + MaxAreaIndexes(buf, cnt + bestR - bestT + 1, indexes + i);
481 indexes[i] = (indexes[i] + bestT) % cnt;
482 /* uses 3*(cnt+bestR-bestT+1)-6 */
486 /* together 3 + 3*(cnt+3) - 18 = 3*cnt-6 q.e.d. */
492 MaxAreaFaceSurface() - divVerent
493 creates a triangle list using max area indexes
496 void MaxAreaFaceSurface(mapDrawSurface_t *ds)
498 /* try to early out */
499 if( !ds->numVerts || (ds->type != SURFACE_FACE && ds->type != SURFACE_DECAL) )
502 /* is this a simple triangle? */
503 if( ds->numVerts == 3 )
506 ds->indexes = safe_malloc( ds->numIndexes * sizeof( int ) );
507 VectorSet( ds->indexes, 0, 1, 2 );
508 numMaxAreaSurfaces++;
513 ds->numIndexes = 3 * ds->numVerts - 6;
514 ds->indexes = safe_malloc( ds->numIndexes * sizeof( int ) );
515 ds->numIndexes = MaxAreaIndexes(ds->verts, ds->numVerts, ds->indexes);
518 numMaxAreaSurfaces++;
521 ClassifySurfaces( 1, ds );
526 FanFaceSurface() - ydnar
527 creates a tri-fan from a brush face winding
528 loosely based on SurfaceAsTriFan()
531 void FanFaceSurface( mapDrawSurface_t *ds )
533 int i, j, k, a, b, c, color[ MAX_LIGHTMAPS ][ 4 ];
534 bspDrawVert_t *verts, *centroid, *dv;
538 /* try to early out */
539 if( !ds->numVerts || (ds->type != SURFACE_FACE && ds->type != SURFACE_DECAL) )
542 /* add a new vertex at the beginning of the surface */
543 verts = safe_malloc( (ds->numVerts + 1) * sizeof( bspDrawVert_t ) );
544 memset( verts, 0, sizeof( bspDrawVert_t ) );
545 memcpy( &verts[ 1 ], ds->verts, ds->numVerts * sizeof( bspDrawVert_t ) );
549 /* add up the drawverts to create a centroid */
550 centroid = &verts[ 0 ];
551 memset( color, 0, 4 * MAX_LIGHTMAPS * sizeof( int ) );
552 for( i = 1, dv = &verts[ 1 ]; i < (ds->numVerts + 1); i++, dv++ )
554 VectorAdd( centroid->xyz, dv->xyz, centroid->xyz );
555 VectorAdd( centroid->normal, dv->normal, centroid->normal );
556 for( j = 0; j < 4; j++ )
558 for( k = 0; k < MAX_LIGHTMAPS; k++ )
559 color[ k ][ j ] += dv->color[ k ][ j ];
562 centroid->st[ j ] += dv->st[ j ];
563 for( k = 0; k < MAX_LIGHTMAPS; k++ )
564 centroid->lightmap[ k ][ j ] += dv->lightmap[ k ][ j ];
569 /* average the centroid */
570 iv = 1.0f / ds->numVerts;
571 VectorScale( centroid->xyz, iv, centroid->xyz );
572 if( VectorNormalize( centroid->normal, centroid->normal ) <= 0 )
573 VectorCopy( verts[ 1 ].normal, centroid->normal );
574 for( j = 0; j < 4; j++ )
576 for( k = 0; k < MAX_LIGHTMAPS; k++ )
578 color[ k ][ j ] /= ds->numVerts;
579 centroid->color[ k ][ j ] = (color[ k ][ j ] < 255.0f ? color[ k ][ j ] : 255);
583 centroid->st[ j ] *= iv;
584 for( k = 0; k < MAX_LIGHTMAPS; k++ )
585 centroid->lightmap[ k ][ j ] *= iv;
589 /* add to vert count */
592 /* fill indexes in triangle fan order */
594 ds->indexes = safe_malloc( ds->numVerts * 3 * sizeof( int ) );
595 for( i = 1; i < ds->numVerts; i++ )
599 c = (i + 1) % ds->numVerts;
601 ds->indexes[ ds->numIndexes++ ] = a;
602 ds->indexes[ ds->numIndexes++ ] = b;
603 ds->indexes[ ds->numIndexes++ ] = c;
610 ClassifySurfaces( 1, ds );
616 StripFaceSurface() - ydnar
617 attempts to create a valid tri-strip w/o degenerate triangles from a brush face winding
618 based on SurfaceAsTriStrip()
621 #define MAX_INDEXES 1024
623 void StripFaceSurface( mapDrawSurface_t *ds )
625 int i, r, least, rotate, numIndexes, ni, a, b, c, indexes[ MAX_INDEXES ];
629 /* try to early out */
630 if( !ds->numVerts || (ds->type != SURFACE_FACE && ds->type != SURFACE_DECAL) )
633 /* is this a simple triangle? */
634 if( ds->numVerts == 3 )
637 VectorSet( indexes, 0, 1, 2 );
641 /* ydnar: find smallest coordinate */
643 if( ds->shaderInfo != NULL && ds->shaderInfo->autosprite == qfalse )
645 for( i = 0; i < ds->numVerts; i++ )
648 v1 = ds->verts[ i ].xyz;
649 v2 = ds->verts[ least ].xyz;
652 if( v1[ 0 ] < v2[ 0 ] ||
653 (v1[ 0 ] == v2[ 0 ] && v1[ 1 ] < v2[ 1 ]) ||
654 (v1[ 0 ] == v2[ 0 ] && v1[ 1 ] == v2[ 1 ] && v1[ 2 ] < v2[ 2 ]) )
659 /* determine the triangle strip order */
660 numIndexes = (ds->numVerts - 2) * 3;
661 if( numIndexes > MAX_INDEXES )
662 Error( "MAX_INDEXES exceeded for surface (%d > %d) (%d verts)", numIndexes, MAX_INDEXES, ds->numVerts );
664 /* try all possible orderings of the points looking for a non-degenerate strip order */
665 for( r = 0; r < ds->numVerts; r++ )
668 rotate = (r + least) % ds->numVerts;
670 /* walk the winding in both directions */
671 for( ni = 0, i = 0; i < ds->numVerts - 2 - i; i++ )
674 a = (ds->numVerts - 1 - i + rotate) % ds->numVerts;
675 b = (i + rotate ) % ds->numVerts;
676 c = (ds->numVerts - 2 - i + rotate) % ds->numVerts;
678 /* test this triangle */
679 if( ds->numVerts > 4 && IsTriangleDegenerate( ds->verts, a, b, c ) )
685 /* handle end case */
686 if( i + 1 != ds->numVerts - 1 - i )
689 a = (ds->numVerts - 2 - i + rotate ) % ds->numVerts;
690 b = (i + rotate ) % ds->numVerts;
691 c = (i + 1 + rotate ) % ds->numVerts;
694 if( ds->numVerts > 4 && IsTriangleDegenerate( ds->verts, a, b, c ) )
703 if( ni == numIndexes )
707 /* if any triangle in the strip is degenerate, render from a centered fan point instead */
708 if( ni < numIndexes )
710 FanFaceSurface( ds );
715 /* copy strip triangle indexes */
716 ds->numIndexes = numIndexes;
717 ds->indexes = safe_malloc( ds->numIndexes * sizeof( int ) );
718 memcpy( ds->indexes, indexes, ds->numIndexes * sizeof( int ) );
724 ClassifySurfaces( 1, ds );
730 vortex: prints meta statistics in general output
735 Sys_Printf( "--- EmitMetaStats ---\n" );
736 Sys_Printf( "%9d total meta surfaces\n", numMetaSurfaces );
737 Sys_Printf( "%9d stripped surfaces\n", numStripSurfaces );
738 Sys_Printf( "%9d fanned surfaces\n", numFanSurfaces );
739 Sys_Printf( "%9d maxarea'd surfaces\n", numMaxAreaSurfaces );
740 Sys_Printf( "%9d patch meta surfaces\n", numPatchMetaSurfaces );
741 Sys_Printf( "%9d meta verts\n", numMetaVerts );
742 Sys_Printf( "%9d meta triangles\n", numMetaTriangles );
746 MakeEntityMetaTriangles()
747 builds meta triangles from brush faces (tristrips and fans)
750 void MakeEntityMetaTriangles( entity_t *e )
752 int i, f, fOld, start;
753 mapDrawSurface_t *ds;
757 Sys_FPrintf( SYS_VRB, "--- MakeEntityMetaTriangles ---\n" );
761 start = I_FloatTime();
763 /* walk the list of surfaces in the entity */
764 for( i = e->firstDrawSurf; i < numMapDrawSurfs; i++ )
767 f = 10 * (i - e->firstDrawSurf) / (numMapDrawSurfs - e->firstDrawSurf);
771 Sys_FPrintf( SYS_VRB, "%d...", f );
775 ds = &mapDrawSurfs[ i ];
776 if( ds->numVerts <= 0 )
779 /* ignore autosprite surfaces */
780 if( ds->shaderInfo->autosprite )
783 /* meta this surface? */
784 if( meta == qfalse && ds->shaderInfo->forceMeta == qfalse )
792 if(maxAreaFaceSurface)
793 MaxAreaFaceSurface( ds );
795 StripFaceSurface( ds );
796 SurfaceToMetaTriangles( ds );
800 TriangulatePatchSurface(e, ds );
803 case SURFACE_TRIANGLES:
806 case SURFACE_FORCED_META:
808 SurfaceToMetaTriangles( ds );
817 if( (numMapDrawSurfs - e->firstDrawSurf) )
818 Sys_FPrintf( SYS_VRB, " (%d)\n", (int) (I_FloatTime() - start) );
820 /* emit some stats */
821 Sys_FPrintf( SYS_VRB, "%9d total meta surfaces\n", numMetaSurfaces );
822 Sys_FPrintf( SYS_VRB, "%9d stripped surfaces\n", numStripSurfaces );
823 Sys_FPrintf( SYS_VRB, "%9d fanned surfaces\n", numFanSurfaces );
824 Sys_FPrintf( SYS_VRB, "%9d maxarea'd surfaces\n", numMaxAreaSurfaces );
825 Sys_FPrintf( SYS_VRB, "%9d patch meta surfaces\n", numPatchMetaSurfaces );
826 Sys_FPrintf( SYS_VRB, "%9d meta verts\n", numMetaVerts );
827 Sys_FPrintf( SYS_VRB, "%9d meta triangles\n", numMetaTriangles );
830 TidyEntitySurfaces( e );
836 PointTriangleIntersect()
837 assuming that all points lie in plane, determine if pt
838 is inside the triangle abc
839 code originally (c) 2001 softSurfer (www.softsurfer.com)
842 #define MIN_OUTSIDE_EPSILON -0.01f
843 #define MAX_OUTSIDE_EPSILON 1.01f
845 static qboolean PointTriangleIntersect( vec3_t pt, vec4_t plane, vec3_t a, vec3_t b, vec3_t c, vec3_t bary )
848 float uu, uv, vv, wu, wv, d;
852 VectorSubtract( b, a, u );
853 VectorSubtract( c, a, v );
854 VectorSubtract( pt, a, w );
857 uu = DotProduct( u, u );
858 uv = DotProduct( u, v );
859 vv = DotProduct( v, v );
860 wu = DotProduct( w, u );
861 wv = DotProduct( w, v );
862 d = uv * uv - uu * vv;
864 /* calculate barycentric coordinates */
865 bary[ 1 ] = (uv * wv - vv * wu) / d;
866 if( bary[ 1 ] < MIN_OUTSIDE_EPSILON || bary[ 1 ] > MAX_OUTSIDE_EPSILON )
868 bary[ 2 ] = (uv * wv - uu * wv) / d;
869 if( bary[ 2 ] < MIN_OUTSIDE_EPSILON || bary[ 2 ] > MAX_OUTSIDE_EPSILON )
871 bary[ 0 ] = 1.0f - (bary[ 1 ] + bary[ 2 ]);
873 /* point is in triangle */
881 sets up an edge structure from a plane and 2 points that the edge ab falls lies in
884 typedef struct edge_s
887 vec_t length, kingpinLength;
893 void CreateEdge( vec4_t plane, vec3_t a, vec3_t b, edge_t *edge )
895 /* copy edge origin */
896 VectorCopy( a, edge->origin );
898 /* create vector aligned with winding direction of edge */
899 VectorSubtract( b, a, edge->edge );
901 if( fabs( edge->edge[ 0 ] ) > fabs( edge->edge[ 1 ] ) && fabs( edge->edge[ 0 ] ) > fabs( edge->edge[ 2 ] ) )
903 else if( fabs( edge->edge[ 1 ] ) > fabs( edge->edge[ 0 ] ) && fabs( edge->edge[ 1 ] ) > fabs( edge->edge[ 2 ] ) )
907 edge->kingpinLength = edge->edge[ edge->kingpin ];
909 VectorNormalize( edge->edge, edge->edge );
910 edge->edge[ 3 ] = DotProduct( a, edge->edge );
911 edge->length = DotProduct( b, edge->edge ) - edge->edge[ 3 ];
913 /* create perpendicular plane that edge lies in */
914 CrossProduct( plane, edge->edge, edge->plane );
915 edge->plane[ 3 ] = DotProduct( a, edge->plane );
922 fixes t-junctions on meta triangles
925 #define TJ_PLANE_EPSILON (1.0f / 8.0f)
926 #define TJ_EDGE_EPSILON (1.0f / 8.0f)
927 #define TJ_POINT_EPSILON (1.0f / 8.0f)
929 void FixMetaTJunctions( void )
931 int i, j, k, f, fOld, start, vertIndex, triIndex, numTJuncs;
932 metaTriangle_t *tri, *newTri;
934 bspDrawVert_t *a, *b, *c, junc;
941 /* this code is crap; revisit later */
945 Sys_FPrintf( SYS_VRB, "--- FixMetaTJunctions ---\n" );
949 start = I_FloatTime();
951 /* walk triangle list */
953 for( i = 0; i < numMetaTriangles; i++ )
956 tri = &metaTriangles[ i ];
959 f = 10 * i / numMetaTriangles;
963 Sys_FPrintf( SYS_VRB, "%d...", f );
966 /* attempt to early out */
968 if( (si->compileFlags & C_NODRAW) || si->autosprite || si->notjunc )
971 /* calculate planes */
972 VectorCopy( tri->plane, plane );
973 plane[ 3 ] = tri->plane[ 3 ];
974 CreateEdge( plane, metaVerts[ tri->indexes[ 0 ] ].xyz, metaVerts[ tri->indexes[ 1 ] ].xyz, &edges[ 0 ] );
975 CreateEdge( plane, metaVerts[ tri->indexes[ 1 ] ].xyz, metaVerts[ tri->indexes[ 2 ] ].xyz, &edges[ 1 ] );
976 CreateEdge( plane, metaVerts[ tri->indexes[ 2 ] ].xyz, metaVerts[ tri->indexes[ 0 ] ].xyz, &edges[ 2 ] );
978 /* walk meta vert list */
979 for( j = 0; j < numMetaVerts; j++ )
982 VectorCopy( metaVerts[ j ].xyz, pt );
984 /* debug code: darken verts */
986 VectorSet( metaVerts[ j ].color[ 0 ], 8, 8, 8 );
988 /* determine if point lies in the triangle's plane */
989 dist = DotProduct( pt, plane ) - plane[ 3 ];
990 if( fabs( dist ) > TJ_PLANE_EPSILON )
993 /* skip this point if it already exists in the triangle */
994 for( k = 0; k < 3; k++ )
996 if( fabs( pt[ 0 ] - metaVerts[ tri->indexes[ k ] ].xyz[ 0 ] ) <= TJ_POINT_EPSILON &&
997 fabs( pt[ 1 ] - metaVerts[ tri->indexes[ k ] ].xyz[ 1 ] ) <= TJ_POINT_EPSILON &&
998 fabs( pt[ 2 ] - metaVerts[ tri->indexes[ k ] ].xyz[ 2 ] ) <= TJ_POINT_EPSILON )
1005 for( k = 0; k < 3; k++ )
1007 /* ignore bogus edges */
1008 if( fabs( edges[ k ].kingpinLength ) < TJ_EDGE_EPSILON )
1011 /* determine if point lies on the edge */
1012 dist = DotProduct( pt, edges[ k ].plane ) - edges[ k ].plane[ 3 ];
1013 if( fabs( dist ) > TJ_EDGE_EPSILON )
1016 /* determine how far along the edge the point lies */
1017 amount = (pt[ edges[ k ].kingpin ] - edges[ k ].origin[ edges[ k ].kingpin ]) / edges[ k ].kingpinLength;
1018 if( amount <= 0.0f || amount >= 1.0f )
1022 dist = DotProduct( pt, edges[ k ].edge ) - edges[ k ].edge[ 3 ];
1023 if( dist <= -0.0f || dist >= edges[ k ].length )
1025 amount = dist / edges[ k ].length;
1028 /* debug code: brighten this point */
1029 //% metaVerts[ j ].color[ 0 ][ 0 ] += 5;
1030 //% metaVerts[ j ].color[ 0 ][ 1 ] += 4;
1031 VectorSet( metaVerts[ tri->indexes[ k ] ].color[ 0 ], 255, 204, 0 );
1032 VectorSet( metaVerts[ tri->indexes[ (k + 1) % 3 ] ].color[ 0 ], 255, 204, 0 );
1035 /* the edge opposite the zero-weighted vertex was hit, so use that as an amount */
1036 a = &metaVerts[ tri->indexes[ k % 3 ] ];
1037 b = &metaVerts[ tri->indexes[ (k + 1) % 3 ] ];
1038 c = &metaVerts[ tri->indexes[ (k + 2) % 3 ] ];
1041 LerpDrawVertAmount( a, b, amount, &junc );
1042 VectorCopy( pt, junc.xyz );
1044 /* compare against existing verts */
1045 if( VectorCompare( junc.xyz, a->xyz ) || VectorCompare( junc.xyz, b->xyz ) || VectorCompare( junc.xyz, c->xyz ) )
1048 /* see if we can just re-use the existing vert */
1049 if( !memcmp( &metaVerts[ j ], &junc, sizeof( junc ) ) )
1053 /* find new vertex (note: a and b are invalid pointers after this) */
1054 firstSearchMetaVert = numMetaVerts;
1055 vertIndex = FindMetaVertex( &junc );
1060 /* make new triangle */
1061 triIndex = AddMetaTriangle();
1066 tri = &metaTriangles[ i ];
1067 newTri = &metaTriangles[ triIndex ];
1069 /* copy the triangle */
1070 memcpy( newTri, tri, sizeof( *tri ) );
1073 tri->indexes[ (k + 1) % 3 ] = vertIndex;
1074 newTri->indexes[ k ] = vertIndex;
1076 /* recalculate edges */
1077 CreateEdge( plane, metaVerts[ tri->indexes[ 0 ] ].xyz, metaVerts[ tri->indexes[ 1 ] ].xyz, &edges[ 0 ] );
1078 CreateEdge( plane, metaVerts[ tri->indexes[ 1 ] ].xyz, metaVerts[ tri->indexes[ 2 ] ].xyz, &edges[ 1 ] );
1079 CreateEdge( plane, metaVerts[ tri->indexes[ 2 ] ].xyz, metaVerts[ tri->indexes[ 0 ] ].xyz, &edges[ 2 ] );
1082 metaVerts[ vertIndex ].color[ 0 ][ 0 ] = 255;
1083 metaVerts[ vertIndex ].color[ 0 ][ 1 ] = 204;
1084 metaVerts[ vertIndex ].color[ 0 ][ 2 ] = 0;
1086 /* add to counter and end processing of this vert */
1094 Sys_FPrintf( SYS_VRB, " (%d)\n", (int) (I_FloatTime() - start) );
1096 /* emit some stats */
1097 Sys_FPrintf( SYS_VRB, "%9d T-junctions added\n", numTJuncs );
1103 SmoothMetaTriangles()
1104 averages coincident vertex normals in the meta triangles
1107 #define MAX_SAMPLES 256
1108 #define THETA_EPSILON 0.000001
1109 #define EQUAL_NORMAL_EPSILON 0.01
1111 void SmoothMetaTriangles( void )
1113 int i, j, k, f, fOld, start, cs, numVerts, numVotes, numSmoothed;
1114 float shadeAngle, defaultShadeAngle, maxShadeAngle, dot, testAngle;
1115 metaTriangle_t *tri;
1118 vec3_t average, diff;
1119 int indexes[ MAX_SAMPLES ];
1120 vec3_t votes[ MAX_SAMPLES ];
1123 Sys_FPrintf( SYS_VRB, "--- SmoothMetaTriangles ---\n" );
1125 /* allocate shade angle table */
1126 shadeAngles = safe_malloc( numMetaVerts * sizeof( float ) );
1127 memset( shadeAngles, 0, numMetaVerts * sizeof( float ) );
1129 /* allocate smoothed table */
1130 cs = (numMetaVerts / 8) + 1;
1131 smoothed = safe_malloc( cs );
1132 memset( smoothed, 0, cs );
1134 /* set default shade angle */
1135 defaultShadeAngle = DEG2RAD( npDegrees );
1136 maxShadeAngle = 0.0f;
1138 /* run through every surface and flag verts belonging to non-lightmapped surfaces
1139 and set per-vertex smoothing angle */
1140 for( i = 0, tri = &metaTriangles[ i ]; i < numMetaTriangles; i++, tri++ )
1142 shadeAngle = defaultShadeAngle;
1144 /* get shade angle from shader */
1145 if( tri->si->shadeAngleDegrees > 0.0f )
1146 shadeAngle = DEG2RAD( tri->si->shadeAngleDegrees );
1147 /* get shade angle from entity */
1148 else if( tri->shadeAngleDegrees > 0.0f )
1149 shadeAngle = DEG2RAD( tri->shadeAngleDegrees );
1151 if( shadeAngle <= 0.0f )
1152 shadeAngle = defaultShadeAngle;
1154 if( shadeAngle > maxShadeAngle )
1155 maxShadeAngle = shadeAngle;
1157 /* flag its verts */
1158 for( j = 0; j < 3; j++ )
1160 shadeAngles[ tri->indexes[ j ] ] = shadeAngle;
1161 if( shadeAngle <= 0 )
1162 smoothed[ tri->indexes[ j ] >> 3 ] |= (1 << (tri->indexes[ j ] & 7));
1166 /* bail if no surfaces have a shade angle */
1167 if( maxShadeAngle <= 0 )
1169 Sys_FPrintf( SYS_VRB, "No smoothing angles specified, aborting\n" );
1170 free( shadeAngles );
1177 start = I_FloatTime();
1179 /* go through the list of vertexes */
1181 for( i = 0; i < numMetaVerts; i++ )
1183 /* print pacifier */
1184 f = 10 * i / numMetaVerts;
1188 Sys_FPrintf( SYS_VRB, "%d...", f );
1191 /* already smoothed? */
1192 if( smoothed[ i >> 3 ] & (1 << (i & 7)) )
1196 VectorClear( average );
1200 /* build a table of coincident vertexes */
1201 for( j = i; j < numMetaVerts && numVerts < MAX_SAMPLES; j++ )
1203 /* already smoothed? */
1204 if( smoothed[ j >> 3 ] & (1 << (j & 7)) )
1208 if( VectorCompare( metaVerts[ i ].xyz, metaVerts[ j ].xyz ) == qfalse )
1211 /* use smallest shade angle */
1212 shadeAngle = (shadeAngles[ i ] < shadeAngles[ j ] ? shadeAngles[ i ] : shadeAngles[ j ]);
1214 /* check shade angle */
1215 dot = DotProduct( metaVerts[ i ].normal, metaVerts[ j ].normal );
1218 else if( dot < -1.0 )
1220 testAngle = acos( dot ) + THETA_EPSILON;
1221 if( testAngle >= shadeAngle )
1224 /* add to the list */
1225 indexes[ numVerts++ ] = j;
1228 smoothed[ j >> 3 ] |= (1 << (j & 7));
1230 /* see if this normal has already been voted */
1231 for( k = 0; k < numVotes; k++ )
1233 VectorSubtract( metaVerts[ j ].normal, votes[ k ], diff );
1234 if( fabs( diff[ 0 ] ) < EQUAL_NORMAL_EPSILON &&
1235 fabs( diff[ 1 ] ) < EQUAL_NORMAL_EPSILON &&
1236 fabs( diff[ 2 ] ) < EQUAL_NORMAL_EPSILON )
1240 /* add a new vote? */
1241 if( k == numVotes && numVotes < MAX_SAMPLES )
1243 VectorAdd( average, metaVerts[ j ].normal, average );
1244 VectorCopy( metaVerts[ j ].normal, votes[ numVotes ] );
1249 /* don't average for less than 2 verts */
1253 /* average normal */
1254 if( VectorNormalize( average, average ) > 0 )
1257 for( j = 0; j < numVerts; j++ )
1258 VectorCopy( average, metaVerts[ indexes[ j ] ].normal );
1263 /* free the tables */
1264 free( shadeAngles );
1268 Sys_FPrintf( SYS_VRB, " (%d)\n", (int) (I_FloatTime() - start) );
1270 /* emit some stats */
1271 Sys_FPrintf( SYS_VRB, "%9d smoothed vertexes\n", numSmoothed );
1277 AddMetaVertToSurface()
1278 adds a drawvert to a surface unless an existing vert matching already exists
1279 returns the index of that vert (or < 0 on failure)
1282 int AddMetaVertToSurface( mapDrawSurface_t *ds, bspDrawVert_t *dv1, int *coincident )
1288 /* go through the verts and find a suitable candidate */
1289 for( i = 0; i < ds->numVerts; i++ )
1292 dv2 = &ds->verts[ i ];
1294 /* compare xyz and normal */
1295 if( VectorCompare( dv1->xyz, dv2->xyz ) == qfalse )
1297 if( VectorCompare( dv1->normal, dv2->normal ) == qfalse )
1300 /* good enough at this point */
1303 /* compare texture coordinates and color */
1304 if( dv1->st[ 0 ] != dv2->st[ 0 ] || dv1->st[ 1 ] != dv2->st[ 1 ] )
1306 if( dv1->color[ 0 ][ 3 ] != dv2->color[ 0 ][ 3 ] )
1309 /* found a winner */
1314 /* overflow check */
1315 if( ds->numVerts >= ((ds->shaderInfo->compileFlags & C_VERTEXLIT) ? maxSurfaceVerts : maxLMSurfaceVerts) )
1316 return VERTS_EXCEEDED;
1318 /* made it this far, add the vert and return */
1319 dv2 = &ds->verts[ ds->numVerts++ ];
1321 return (ds->numVerts - 1);
1328 AddMetaTriangleToSurface()
1329 attempts to add a metatriangle to a surface
1330 returns the score of the triangle added
1333 #define AXIS_SCORE 100000
1334 #define AXIS_MIN 100000
1335 #define VERT_SCORE 10000
1336 #define SURFACE_SCORE 1000
1338 #define ST_SCORE2 (2 * (ST_SCORE))
1340 #define ADEQUATE_SCORE ((AXIS_MIN) + 1 * (VERT_SCORE))
1341 #define GOOD_SCORE ((AXIS_MIN) + 2 * (VERT_SCORE) + 4 * (ST_SCORE))
1342 #define PERFECT_SCORE ((AXIS_MIN) + 3 * (VERT_SCORE) + (SURFACE_SCORE) + 4 * (ST_SCORE))
1343 //#define MAX_BBOX_DISTANCE 16
1345 static int AddMetaTriangleToSurface( mapDrawSurface_t *ds, metaTriangle_t *tri, qboolean testAdd )
1347 int i, score, coincident, ai, bi, ci, oldTexRange[ 2 ];
1349 vec3_t mins, maxs, p;
1350 qboolean inTexRange, es, et;
1351 mapDrawSurface_t old;
1354 /* overflow check */
1355 if( ds->numIndexes >= maxSurfaceIndexes )
1358 /* test the triangle */
1359 if( ds->entityNum != tri->entityNum ) /* ydnar: added 2002-07-06 */
1361 if( ds->castShadows != tri->castShadows || ds->recvShadows != tri->recvShadows )
1363 if( ds->shaderInfo != tri->si || ds->fogNum != tri->fogNum || ds->sampleSize != tri->sampleSize )
1366 if( !(ds->shaderInfo->compileFlags & C_VERTEXLIT) &&
1367 //% VectorCompare( ds->lightmapAxis, tri->lightmapAxis ) == qfalse )
1368 DotProduct( ds->lightmapAxis, tri->plane ) < 0.25f )
1372 /* planar surfaces will only merge with triangles in the same plane */
1373 if( npDegrees == 0.0f && ds->shaderInfo->nonplanar == qfalse && ds->planeNum >= 0 )
1375 if( VectorCompare( mapplanes[ ds->planeNum ].normal, tri->plane ) == qfalse || mapplanes[ ds->planeNum ].dist != tri->plane[ 3 ] )
1377 if( tri->planeNum >= 0 && tri->planeNum != ds->planeNum )
1381 #if MAX_BBOX_DISTANCE > 0
1382 VectorCopy( ds->mins, mins );
1383 VectorCopy( ds->maxs, maxs );
1384 mins[0] -= MAX_BBOX_DISTANCE;
1385 mins[1] -= MAX_BBOX_DISTANCE;
1386 mins[2] -= MAX_BBOX_DISTANCE;
1387 maxs[0] += MAX_BBOX_DISTANCE;
1388 maxs[1] += MAX_BBOX_DISTANCE;
1389 maxs[2] += MAX_BBOX_DISTANCE;
1390 #define CHECK_1D(mins, v, maxs) ((mins) <= (v) && (v) <= (maxs))
1391 #define CHECK_3D(mins, v, maxs) (CHECK_1D((mins)[0], (v)[0], (maxs)[0]) && CHECK_1D((mins)[1], (v)[1], (maxs)[1]) && CHECK_1D((mins)[2], (v)[2], (maxs)[2]))
1392 VectorCopy(metaVerts[ tri->indexes[ 0 ] ].xyz, p);
1393 if(!CHECK_3D(mins, p, maxs))
1395 VectorCopy(metaVerts[ tri->indexes[ 1 ] ].xyz, p);
1396 if(!CHECK_3D(mins, p, maxs))
1398 VectorCopy(metaVerts[ tri->indexes[ 2 ] ].xyz, p);
1399 if(!CHECK_3D(mins, p, maxs))
1407 /* set initial score */
1408 score = tri->surfaceNum == ds->surfaceNum ? SURFACE_SCORE : 0;
1410 /* score the the dot product of lightmap axis to plane */
1411 if( (ds->shaderInfo->compileFlags & C_VERTEXLIT) || VectorCompare( ds->lightmapAxis, tri->lightmapAxis ) )
1412 score += AXIS_SCORE;
1414 score += AXIS_SCORE * DotProduct( ds->lightmapAxis, tri->plane );
1416 /* preserve old drawsurface if this fails */
1417 memcpy( &old, ds, sizeof( *ds ) );
1419 /* attempt to add the verts */
1421 ai = AddMetaVertToSurface( ds, &metaVerts[ tri->indexes[ 0 ] ], &coincident );
1422 bi = AddMetaVertToSurface( ds, &metaVerts[ tri->indexes[ 1 ] ], &coincident );
1423 ci = AddMetaVertToSurface( ds, &metaVerts[ tri->indexes[ 2 ] ], &coincident );
1425 /* check vertex underflow */
1426 if( ai < 0 || bi < 0 || ci < 0 )
1428 memcpy( ds, &old, sizeof( *ds ) );
1432 /* score coincident vertex count (2003-02-14: changed so this only matters on planar surfaces) */
1433 score += (coincident * VERT_SCORE);
1435 /* add new vertex bounds to mins/maxs */
1436 VectorCopy( ds->mins, mins );
1437 VectorCopy( ds->maxs, maxs );
1438 AddPointToBounds( metaVerts[ tri->indexes[ 0 ] ].xyz, mins, maxs );
1439 AddPointToBounds( metaVerts[ tri->indexes[ 1 ] ].xyz, mins, maxs );
1440 AddPointToBounds( metaVerts[ tri->indexes[ 2 ] ].xyz, mins, maxs );
1442 /* check lightmap bounds overflow (after at least 1 triangle has been added) */
1443 if( !(ds->shaderInfo->compileFlags & C_VERTEXLIT) &&
1444 ds->numIndexes > 0 && VectorLength( ds->lightmapAxis ) > 0.0f &&
1445 (VectorCompare( ds->mins, mins ) == qfalse || VectorCompare( ds->maxs, maxs ) == qfalse) )
1447 /* set maximum size before lightmap scaling (normally 2032 units) */
1448 /* 2004-02-24: scale lightmap test size by 2 to catch larger brush faces */
1449 /* 2004-04-11: reverting to actual lightmap size */
1450 lmMax = (ds->sampleSize * (ds->shaderInfo->lmCustomWidth - 1));
1451 for( i = 0; i < 3; i++ )
1453 if( (maxs[ i ] - mins[ i ]) > lmMax )
1455 memcpy( ds, &old, sizeof( *ds ) );
1461 /* check texture range overflow */
1462 oldTexRange[ 0 ] = ds->texRange[ 0 ];
1463 oldTexRange[ 1 ] = ds->texRange[ 1 ];
1464 inTexRange = CalcSurfaceTextureRange( ds );
1466 es = (ds->texRange[ 0 ] > oldTexRange[ 0 ]) ? qtrue : qfalse;
1467 et = (ds->texRange[ 1 ] > oldTexRange[ 1 ]) ? qtrue : qfalse;
1469 if( inTexRange == qfalse && ds->numIndexes > 0 )
1471 memcpy( ds, &old, sizeof( *ds ) );
1472 return UNSUITABLE_TRIANGLE;
1475 /* score texture range */
1476 if( ds->texRange[ 0 ] <= oldTexRange[ 0 ] )
1478 else if( ds->texRange[ 0 ] > oldTexRange[ 0 ] && oldTexRange[ 1 ] > oldTexRange[ 0 ] )
1481 if( ds->texRange[ 1 ] <= oldTexRange[ 1 ] )
1483 else if( ds->texRange[ 1 ] > oldTexRange[ 1 ] && oldTexRange[ 0 ] > oldTexRange[ 1 ] )
1487 /* go through the indexes and try to find an existing triangle that matches abc */
1488 for( i = 0; i < ds->numIndexes; i += 3 )
1490 /* 2002-03-11 (birthday!): rotate the triangle 3x to find an existing triangle */
1491 if( (ai == ds->indexes[ i ] && bi == ds->indexes[ i + 1 ] && ci == ds->indexes[ i + 2 ]) ||
1492 (bi == ds->indexes[ i ] && ci == ds->indexes[ i + 1 ] && ai == ds->indexes[ i + 2 ]) ||
1493 (ci == ds->indexes[ i ] && ai == ds->indexes[ i + 1 ] && bi == ds->indexes[ i + 2 ]) )
1495 /* triangle already present */
1496 memcpy( ds, &old, sizeof( *ds ) );
1501 /* rotate the triangle 3x to find an inverse triangle (error case) */
1502 if( (ai == ds->indexes[ i ] && bi == ds->indexes[ i + 2 ] && ci == ds->indexes[ i + 1 ]) ||
1503 (bi == ds->indexes[ i ] && ci == ds->indexes[ i + 2 ] && ai == ds->indexes[ i + 1 ]) ||
1504 (ci == ds->indexes[ i ] && ai == ds->indexes[ i + 2 ] && bi == ds->indexes[ i + 1 ]) )
1507 Sys_Printf( "WARNING: Flipped triangle: (%6.0f %6.0f %6.0f) (%6.0f %6.0f %6.0f) (%6.0f %6.0f %6.0f)\n",
1508 ds->verts[ ai ].xyz[ 0 ], ds->verts[ ai ].xyz[ 1 ], ds->verts[ ai ].xyz[ 2 ],
1509 ds->verts[ bi ].xyz[ 0 ], ds->verts[ bi ].xyz[ 1 ], ds->verts[ bi ].xyz[ 2 ],
1510 ds->verts[ ci ].xyz[ 0 ], ds->verts[ ci ].xyz[ 1 ], ds->verts[ ci ].xyz[ 2 ] );
1512 /* reverse triangle already present */
1513 memcpy( ds, &old, sizeof( *ds ) );
1519 /* add the triangle indexes */
1520 if( ds->numIndexes < maxSurfaceIndexes )
1521 ds->indexes[ ds->numIndexes++ ] = ai;
1522 if( ds->numIndexes < maxSurfaceIndexes )
1523 ds->indexes[ ds->numIndexes++ ] = bi;
1524 if( ds->numIndexes < maxSurfaceIndexes )
1525 ds->indexes[ ds->numIndexes++ ] = ci;
1527 /* check index overflow */
1528 if( ds->numIndexes >= maxSurfaceIndexes )
1530 memcpy( ds, &old, sizeof( *ds ) );
1534 /* sanity check the indexes */
1535 if( ds->numIndexes >= 3 &&
1536 (ds->indexes[ ds->numIndexes - 3 ] == ds->indexes[ ds->numIndexes - 2 ] ||
1537 ds->indexes[ ds->numIndexes - 3 ] == ds->indexes[ ds->numIndexes - 1 ] ||
1538 ds->indexes[ ds->numIndexes - 2 ] == ds->indexes[ ds->numIndexes - 1 ]) )
1539 Sys_Printf( "DEG:%d! ", ds->numVerts );
1543 memcpy( ds, &old, sizeof( *ds ) );
1546 /* copy bounds back to surface */
1547 VectorCopy( mins, ds->mins );
1548 VectorCopy( maxs, ds->maxs );
1550 /* mark triangle as used */
1554 /* add a side reference */
1555 ds->sideRef = AllocSideRef( tri->side, ds->sideRef );
1557 /* return to sender */
1564 MetaTrianglesToSurface()
1565 creates map drawsurface(s) from the list of possibles
1568 static void MetaTrianglesToSurface( int numPossibles, metaTriangle_t *possibles, int *fOld, int *numAdded )
1570 int i, j, f, best, score, bestScore;
1571 metaTriangle_t *seed, *test;
1572 mapDrawSurface_t *ds;
1573 bspDrawVert_t *verts;
1578 /* allocate arrays */
1579 verts = safe_malloc( sizeof( *verts ) * maxSurfaceVerts );
1580 indexes = safe_malloc( sizeof( *indexes ) * maxSurfaceIndexes );
1582 /* walk the list of triangles */
1583 for( i = 0, seed = possibles; i < numPossibles; i++, seed++ )
1585 /* skip this triangle if it has already been merged */
1586 if( seed->si == NULL )
1589 /* -----------------------------------------------------------------
1590 initial drawsurf construction
1591 ----------------------------------------------------------------- */
1593 /* start a new drawsurface */
1594 ds = AllocDrawSurface( SURFACE_META );
1595 ds->entityNum = seed->entityNum;
1596 ds->surfaceNum = seed->surfaceNum;
1597 ds->castShadows = seed->castShadows;
1598 ds->recvShadows = seed->recvShadows;
1600 ds->shaderInfo = seed->si;
1601 ds->planeNum = seed->planeNum;
1602 ds->fogNum = seed->fogNum;
1603 ds->sampleSize = seed->sampleSize;
1604 ds->shadeAngleDegrees = seed->shadeAngleDegrees;
1606 ds->indexes = indexes;
1607 VectorCopy( seed->lightmapAxis, ds->lightmapAxis );
1608 ds->sideRef = AllocSideRef( seed->side, NULL );
1610 ClearBounds( ds->mins, ds->maxs );
1612 /* clear verts/indexes */
1613 memset( verts, 0, sizeof( verts ) );
1614 memset( indexes, 0, sizeof( indexes ) );
1616 /* add the first triangle */
1617 if( AddMetaTriangleToSurface( ds, seed, qfalse ) )
1620 /* -----------------------------------------------------------------
1622 ----------------------------------------------------------------- */
1624 /* progressively walk the list until no more triangles can be added */
1628 /* print pacifier */
1629 f = 10 * *numAdded / numMetaTriangles;
1633 Sys_FPrintf( SYS_VRB, "%d...", f );
1636 /* reset best score */
1641 /* walk the list of possible candidates for merging */
1642 for( j = i + 1, test = &possibles[ j ]; j < numPossibles; j++, test++ )
1644 /* skip this triangle if it has already been merged */
1645 if( test->si == NULL )
1648 /* score this triangle */
1649 score = AddMetaTriangleToSurface( ds, test, qtrue );
1650 if( score > bestScore )
1655 /* if we have a score over a certain threshold, just use it */
1656 if( bestScore >= GOOD_SCORE )
1658 if( AddMetaTriangleToSurface( ds, &possibles[ best ], qfalse ) )
1669 /* add best candidate */
1670 if( best >= 0 && bestScore > ADEQUATE_SCORE )
1672 if( AddMetaTriangleToSurface( ds, &possibles[ best ], qfalse ) )
1680 /* copy the verts and indexes to the new surface */
1681 ds->verts = safe_malloc( ds->numVerts * sizeof( bspDrawVert_t ) );
1682 memcpy( ds->verts, verts, ds->numVerts * sizeof( bspDrawVert_t ) );
1683 ds->indexes = safe_malloc( ds->numIndexes * sizeof( int ) );
1684 memcpy( ds->indexes, indexes, ds->numIndexes * sizeof( int ) );
1686 /* classify the surface */
1687 ClassifySurfaces( 1, ds );
1690 numMergedSurfaces++;
1701 CompareMetaTriangles()
1702 compare function for qsort()
1705 static int CompareMetaTriangles( const void *a, const void *b )
1708 vec3_t aMins, bMins;
1712 if( ((metaTriangle_t*) a)->si < ((metaTriangle_t*) b)->si )
1714 else if( ((metaTriangle_t*) a)->si > ((metaTriangle_t*) b)->si )
1718 else if( ((metaTriangle_t*) a)->fogNum < ((metaTriangle_t*) b)->fogNum )
1720 else if( ((metaTriangle_t*) a)->fogNum > ((metaTriangle_t*) b)->fogNum )
1725 else if( npDegrees == 0.0f && ((metaTriangle_t*) a)->si->nonplanar == qfalse &&
1726 ((metaTriangle_t*) a)->planeNum >= 0 && ((metaTriangle_t*) a)->planeNum >= 0 )
1728 if( ((metaTriangle_t*) a)->plane[ 3 ] < ((metaTriangle_t*) b)->plane[ 3 ] )
1730 else if( ((metaTriangle_t*) a)->plane[ 3 ] > ((metaTriangle_t*) b)->plane[ 3 ] )
1732 else if( ((metaTriangle_t*) a)->plane[ 0 ] < ((metaTriangle_t*) b)->plane[ 0 ] )
1734 else if( ((metaTriangle_t*) a)->plane[ 0 ] > ((metaTriangle_t*) b)->plane[ 0 ] )
1736 else if( ((metaTriangle_t*) a)->plane[ 1 ] < ((metaTriangle_t*) b)->plane[ 1 ] )
1738 else if( ((metaTriangle_t*) a)->plane[ 1 ] > ((metaTriangle_t*) b)->plane[ 1 ] )
1740 else if( ((metaTriangle_t*) a)->plane[ 2 ] < ((metaTriangle_t*) b)->plane[ 2 ] )
1742 else if( ((metaTriangle_t*) a)->plane[ 2 ] > ((metaTriangle_t*) b)->plane[ 2 ] )
1747 /* then position in world */
1750 VectorSet( aMins, 999999, 999999, 999999 );
1751 VectorSet( bMins, 999999, 999999, 999999 );
1752 for( i = 0; i < 3; i++ )
1754 av = ((metaTriangle_t*) a)->indexes[ i ];
1755 bv = ((metaTriangle_t*) b)->indexes[ i ];
1756 for( j = 0; j < 3; j++ )
1758 if( metaVerts[ av ].xyz[ j ] < aMins[ j ] )
1759 aMins[ j ] = metaVerts[ av ].xyz[ j ];
1760 if( metaVerts[ bv ].xyz[ j ] < bMins[ j ] )
1761 bMins[ j ] = metaVerts[ bv ].xyz[ j ];
1766 for( i = 0; i < 3; i++ )
1768 if( aMins[ i ] < bMins[ i ] )
1770 else if( aMins[ i ] > bMins[ i ] )
1774 /* functionally equivalent */
1781 MergeMetaTriangles()
1782 merges meta triangles into drawsurfaces
1785 void MergeMetaTriangles( void )
1787 int i, j, fOld, start, numAdded;
1788 metaTriangle_t *head, *end;
1791 /* only do this if there are meta triangles */
1792 if( numMetaTriangles <= 0 )
1796 Sys_FPrintf( SYS_VRB, "--- MergeMetaTriangles ---\n" );
1798 /* sort the triangles by shader major, fognum minor */
1799 qsort( metaTriangles, numMetaTriangles, sizeof( metaTriangle_t ), CompareMetaTriangles );
1803 start = I_FloatTime();
1807 for( i = 0, j = 0; i < numMetaTriangles; i = j )
1809 /* get head of list */
1810 head = &metaTriangles[ i ];
1812 /* skip this triangle if it has already been merged */
1813 if( head->si == NULL )
1819 for( j = i + 1; j < numMetaTriangles; j++ )
1821 /* get end of list */
1822 end = &metaTriangles[ j ];
1823 if( head->si != end->si || head->fogNum != end->fogNum )
1828 /* try to merge this list of possible merge candidates */
1829 MetaTrianglesToSurface( (j - i), head, &fOld, &numAdded );
1832 /* clear meta triangle list */
1833 ClearMetaTriangles();
1837 Sys_FPrintf( SYS_VRB, " (%d)\n", (int) (I_FloatTime() - start) );
1839 /* emit some stats */
1840 Sys_FPrintf( SYS_VRB, "%9d surfaces merged\n", numMergedSurfaces );
1841 Sys_FPrintf( SYS_VRB, "%9d vertexes merged\n", numMergedVerts );