]> de.git.xonotic.org Git - xonotic/darkplaces.git/blob - gl_models.c
cleaned up shadow volume culling somewhat
[xonotic/darkplaces.git] / gl_models.c
1
2 #include "quakedef.h"
3 #include "cl_collision.h"
4 #include "r_shadow.h"
5
6 typedef struct
7 {
8         float m[3][4];
9 } zymbonematrix;
10
11 // LordHavoc: vertex arrays
12
13 float *aliasvertcolorbuf;
14 float *aliasvertcolor; // this may point at aliasvertcolorbuf or at vertex arrays in the mesh backend
15 float *aliasvert_svectors;
16 float *aliasvert_tvectors;
17 float *aliasvert_normals;
18
19 float *aliasvertcolor2;
20 int *aliasvertusage;
21 zymbonematrix *zymbonepose;
22
23 mempool_t *gl_models_mempool;
24
25 void gl_models_start(void)
26 {
27         // allocate vertex processing arrays
28         gl_models_mempool = Mem_AllocPool("GL_Models");
29         aliasvertcolor = aliasvertcolorbuf = Mem_Alloc(gl_models_mempool, sizeof(float[MD2MAX_VERTS][4]));
30         aliasvert_svectors = Mem_Alloc(gl_models_mempool, sizeof(float[MD2MAX_VERTS][4]));
31         aliasvert_tvectors = Mem_Alloc(gl_models_mempool, sizeof(float[MD2MAX_VERTS][4]));
32         aliasvert_normals = Mem_Alloc(gl_models_mempool, sizeof(float[MD2MAX_VERTS][4]));
33         aliasvertcolor2 = Mem_Alloc(gl_models_mempool, sizeof(float[MD2MAX_VERTS][4])); // used temporarily for tinted coloring
34         zymbonepose = Mem_Alloc(gl_models_mempool, sizeof(zymbonematrix[256]));
35         aliasvertusage = Mem_Alloc(gl_models_mempool, sizeof(int[MD2MAX_VERTS]));
36 }
37
38 void gl_models_shutdown(void)
39 {
40         Mem_FreePool(&gl_models_mempool);
41 }
42
43 void gl_models_newmap(void)
44 {
45 }
46
47 void GL_Models_Init(void)
48 {
49         R_RegisterModule("GL_Models", gl_models_start, gl_models_shutdown, gl_models_newmap);
50 }
51
52 void R_Model_Alias_GetMeshVerts(const entity_render_t *ent, aliasmesh_t *mesh, float *vertices, float *normals, float *svectors, float *tvectors)
53 {
54         int i, vertcount;
55         float lerp1, lerp2, lerp3, lerp4;
56         const aliasvertex_t *verts1, *verts2, *verts3, *verts4;
57
58         if (vertices == NULL)
59                 Host_Error("R_Model_Alias_GetMeshVerts: vertices == NULL.\n");
60         if (svectors != NULL && (tvectors == NULL || normals == NULL))
61                 Host_Error("R_Model_Alias_GetMeshVerts: svectors requires tvectors and normals.\n");
62         if (tvectors != NULL && (svectors == NULL || normals == NULL))
63                 Host_Error("R_Model_Alias_GetMeshVerts: tvectors requires svectors and normals.\n");
64
65         vertcount = mesh->num_vertices;
66         verts1 = mesh->data_vertices + ent->frameblend[0].frame * vertcount;
67         lerp1 = ent->frameblend[0].lerp;
68         if (ent->frameblend[1].lerp)
69         {
70                 verts2 = mesh->data_vertices + ent->frameblend[1].frame * vertcount;
71                 lerp2 = ent->frameblend[1].lerp;
72                 if (ent->frameblend[2].lerp)
73                 {
74                         verts3 = mesh->data_vertices + ent->frameblend[2].frame * vertcount;
75                         lerp3 = ent->frameblend[2].lerp;
76                         if (ent->frameblend[3].lerp)
77                         {
78                                 verts4 = mesh->data_vertices + ent->frameblend[3].frame * vertcount;
79                                 lerp4 = ent->frameblend[3].lerp;
80                                 // generate vertices
81                                 if (svectors != NULL)
82                                 {
83                                         for (i = 0;i < vertcount;i++, vertices += 4, normals += 4, svectors += 4, tvectors += 4, verts1++, verts2++, verts3++, verts4++)
84                                         {
85                                                 VectorMAMAMAM(lerp1, verts1->origin, lerp2, verts2->origin, lerp3, verts3->origin, lerp4, verts4->origin, vertices);
86                                                 VectorMAMAMAM(lerp1, verts1->normal, lerp2, verts2->normal, lerp3, verts3->normal, lerp4, verts4->normal, normals);
87                                                 VectorMAMAMAM(lerp1, verts1->svector, lerp2, verts2->svector, lerp3, verts3->svector, lerp4, verts4->svector, svectors);
88                                                 CrossProduct(svectors, normals, tvectors);
89                                         }
90                                 }
91                                 else if (normals != NULL)
92                                 {
93                                         for (i = 0;i < vertcount;i++, vertices += 4, normals += 4, verts1++, verts2++, verts3++, verts4++)
94                                         {
95                                                 VectorMAMAMAM(lerp1, verts1->origin, lerp2, verts2->origin, lerp3, verts3->origin, lerp4, verts4->origin, vertices);
96                                                 VectorMAMAMAM(lerp1, verts1->normal, lerp2, verts2->normal, lerp3, verts3->normal, lerp4, verts4->normal, normals);
97                                         }
98                                 }
99                                 else
100                                         for (i = 0;i < vertcount;i++, vertices += 4, verts1++, verts2++, verts3++, verts4++)
101                                                 VectorMAMAMAM(lerp1, verts1->origin, lerp2, verts2->origin, lerp3, verts3->origin, lerp4, verts4->origin, vertices);
102                         }
103                         else
104                         {
105                                 // generate vertices
106                                 if (svectors != NULL)
107                                 {
108                                         for (i = 0;i < vertcount;i++, vertices += 4, normals += 4, svectors += 4, tvectors += 4, verts1++, verts2++, verts3++)
109                                         {
110                                                 VectorMAMAM(lerp1, verts1->origin, lerp2, verts2->origin, lerp3, verts3->origin, vertices);
111                                                 VectorMAMAM(lerp1, verts1->normal, lerp2, verts2->normal, lerp3, verts3->normal, normals);
112                                                 VectorMAMAM(lerp1, verts1->svector, lerp2, verts2->svector, lerp3, verts3->svector, svectors);
113                                                 CrossProduct(svectors, normals, tvectors);
114                                         }
115                                 }
116                                 else if (normals != NULL)
117                                 {
118                                         for (i = 0;i < vertcount;i++, vertices += 4, normals += 4, verts1++, verts2++, verts3++)
119                                         {
120                                                 VectorMAMAM(lerp1, verts1->origin, lerp2, verts2->origin, lerp3, verts3->origin, vertices);
121                                                 VectorMAMAM(lerp1, verts1->normal, lerp2, verts2->normal, lerp3, verts3->normal, normals);
122                                         }
123                                 }
124                                 else
125                                         for (i = 0;i < vertcount;i++, vertices += 4, verts1++, verts2++, verts3++)
126                                                 VectorMAMAM(lerp1, verts1->origin, lerp2, verts2->origin, lerp3, verts3->origin, vertices);
127                         }
128                 }
129                 else
130                 {
131                         // generate vertices
132                         if (svectors != NULL)
133                         {
134                                 for (i = 0;i < vertcount;i++, vertices += 4, normals += 4, svectors += 4, tvectors += 4, verts1++, verts2++)
135                                 {
136                                         VectorMAM(lerp1, verts1->origin, lerp2, verts2->origin, vertices);
137                                         VectorMAM(lerp1, verts1->normal, lerp2, verts2->normal, normals);
138                                         VectorMAM(lerp1, verts1->svector, lerp2, verts2->svector, svectors);
139                                         CrossProduct(svectors, normals, tvectors);
140                                 }
141                         }
142                         else if (normals != NULL)
143                         {
144                                 for (i = 0;i < vertcount;i++, vertices += 4, normals += 4, verts1++, verts2++)
145                                 {
146                                         VectorMAM(lerp1, verts1->origin, lerp2, verts2->origin, vertices);
147                                         VectorMAM(lerp1, verts1->normal, lerp2, verts2->normal, normals);
148                                 }
149                         }
150                         else
151                                 for (i = 0;i < vertcount;i++, vertices += 4, verts1++, verts2++)
152                                         VectorMAM(lerp1, verts1->origin, lerp2, verts2->origin, vertices);
153                 }
154         }
155         else
156         {
157                 // generate vertices
158                 if (svectors != NULL)
159                 {
160                         for (i = 0;i < vertcount;i++, vertices += 4, normals += 4, svectors += 4, tvectors += 4, verts1++)
161                         {
162                                 VectorM(lerp1, verts1->origin, vertices);
163                                 VectorM(lerp1, verts1->normal, normals);
164                                 VectorM(lerp1, verts1->svector, svectors);
165                                 CrossProduct(svectors, normals, tvectors);
166                         }
167                 }
168                 else if (normals != NULL)
169                 {
170                         for (i = 0;i < vertcount;i++, vertices += 4, normals += 4, verts1++)
171                         {
172                                 VectorM(lerp1, verts1->origin, vertices);
173                                 VectorM(lerp1, verts1->normal, normals);
174                         }
175                 }
176                 else
177                         for (i = 0;i < vertcount;i++, vertices += 4, verts1++)
178                                 VectorM(lerp1, verts1->origin, vertices);
179         }
180 }
181
182 skinframe_t *R_FetchSkinFrame(const entity_render_t *ent)
183 {
184         model_t *model = ent->model;
185         int s = ent->skinnum;
186         if ((unsigned int)s >= (unsigned int)model->numskins)
187                 s = 0;
188         if (model->skinscenes[s].framecount > 1)
189                 return &model->skinframes[model->skinscenes[s].firstframe + (int) (cl.time * 10) % model->skinscenes[s].framecount];
190         else
191                 return &model->skinframes[model->skinscenes[s].firstframe];
192 }
193
194 aliasskin_t *R_FetchAliasSkin(const entity_render_t *ent, const aliasmesh_t *mesh)
195 {
196         model_t *model = ent->model;
197         int s = ent->skinnum;
198         if ((unsigned int)s >= (unsigned int)model->numskins)
199                 s = 0;
200         if (model->skinscenes[s].framecount > 1)
201                 s = model->skinscenes[s].firstframe + (int) (cl.time * model->skinscenes[s].framerate) % model->skinscenes[s].framecount;
202         else
203                 s = model->skinscenes[s].firstframe;
204         if (s > mesh->num_skins)
205                 return mesh->data_skins;
206         return mesh->data_skins + s;
207 }
208
209 void R_DrawAliasModelCallback (const void *calldata1, int calldata2)
210 {
211         int c, fullbright, layernum;
212         float tint[3], fog, ifog, colorscale;
213         vec3_t diff;
214         qbyte *bcolor;
215         rmeshstate_t m;
216         const entity_render_t *ent = calldata1;
217         aliasmesh_t *mesh = ent->model->aliasdata_meshes + calldata2;
218         aliaslayer_t *layer;
219         aliasskin_t *skin;
220
221         R_Mesh_Matrix(&ent->matrix);
222
223         fog = 0;
224         if (fogenabled)
225         {
226                 VectorSubtract(ent->origin, r_origin, diff);
227                 fog = DotProduct(diff,diff);
228                 if (fog < 0.01f)
229                         fog = 0.01f;
230                 fog = exp(fogdensity/fog);
231                 if (fog > 1)
232                         fog = 1;
233                 if (fog < 0.01f)
234                         fog = 0;
235                 // fog method: darken, additive fog
236                 // 1. render model as normal, scaled by inverse of fog alpha (darkens it)
237                 // 2. render fog as additive
238         }
239         ifog = 1 - fog;
240
241         memset(&m, 0, sizeof(m));
242         skin = R_FetchAliasSkin(ent, mesh);
243         R_Mesh_ResizeCheck(mesh->num_vertices);
244         R_Model_Alias_GetMeshVerts(ent, mesh, varray_vertex, aliasvert_normals, NULL, NULL);
245         memcpy(varray_texcoord[0], mesh->data_texcoords, mesh->num_vertices * sizeof(float[4]));
246         for (layernum = 0, layer = skin->data_layers;layernum < skin->num_layers;layernum++, layer++)
247         {
248                 if (((layer->flags & ALIASLAYER_NODRAW_IF_NOTCOLORMAPPED) && ent->colormap < 0)
249                  || ((layer->flags & ALIASLAYER_NODRAW_IF_COLORMAPPED) && ent->colormap >= 0)
250                  ||  (layer->flags & ALIASLAYER_DRAW_PER_LIGHT))
251                         continue;
252                 if (layer->flags & ALIASLAYER_FOG)
253                 {
254                         m.blendfunc1 = GL_SRC_ALPHA;
255                         m.blendfunc2 = GL_ONE;
256                         colorscale = r_colorscale;
257                         m.texrgbscale[0] = 1;
258                         m.tex[0] = R_GetTexture(layer->texture);
259                         R_Mesh_State(&m);
260                         GL_Color(fogcolor[0] * fog * colorscale, fogcolor[1] * fog * colorscale, fogcolor[2] * fog * colorscale, ent->alpha);
261                         c_alias_polys += mesh->num_triangles;
262                         R_Mesh_Draw(mesh->num_vertices, mesh->num_triangles, mesh->data_elements);
263                         continue;
264                 }
265                 if ((layer->flags & ALIASLAYER_ADD) || ((layer->flags & ALIASLAYER_ALPHA) && (ent->effects & EF_ADDITIVE)))
266                 {
267                         m.blendfunc1 = GL_SRC_ALPHA;
268                         m.blendfunc2 = GL_ONE;
269                 }
270                 else if ((layer->flags & ALIASLAYER_ALPHA) || ent->alpha != 1.0)
271                 {
272                         m.blendfunc1 = GL_SRC_ALPHA;
273                         m.blendfunc2 = GL_ONE_MINUS_SRC_ALPHA;
274                 }
275                 else
276                 {
277                         m.blendfunc1 = GL_ONE;
278                         m.blendfunc2 = GL_ZERO;
279                 }
280                 colorscale = r_colorscale;
281                 m.texrgbscale[0] = 1;
282                 if (gl_combine.integer)
283                 {
284                         colorscale *= 0.25f;
285                         m.texrgbscale[0] = 4;
286                 }
287                 m.tex[0] = R_GetTexture(layer->texture);
288                 R_Mesh_State(&m);
289                 if (layer->flags & ALIASLAYER_COLORMAP_PANTS)
290                 {
291                         // 128-224 are backwards ranges
292                         c = (ent->colormap & 0xF) << 4;c += (c >= 128 && c < 224) ? 4 : 12;
293                         bcolor = (qbyte *) (&palette_complete[c]);
294                         fullbright = c >= 224;
295                         VectorScale(bcolor, (1.0f / 255.0f), tint);
296                 }
297                 else if (layer->flags & ALIASLAYER_COLORMAP_SHIRT)
298                 {
299                         // 128-224 are backwards ranges
300                         c = (ent->colormap & 0xF0);c += (c >= 128 && c < 224) ? 4 : 12;
301                         bcolor = (qbyte *) (&palette_complete[c]);
302                         fullbright = c >= 224;
303                         VectorScale(bcolor, (1.0f / 255.0f), tint);
304                 }
305                 else
306                 {
307                         tint[0] = tint[1] = tint[2] = 1;
308                         fullbright = false;
309                 }
310                 VectorScale(tint, ifog * colorscale, tint);
311                 if (!(layer->flags & ALIASLAYER_DIFFUSE))
312                         fullbright = true;
313                 if (ent->effects & EF_FULLBRIGHT)
314                         fullbright = true;
315                 if (fullbright)
316                         GL_Color(tint[0], tint[1], tint[2], ent->alpha);
317                 else
318                         R_LightModel(ent, mesh->num_vertices, varray_vertex, aliasvert_normals, varray_color, tint[0], tint[1], tint[2], false);
319                 c_alias_polys += mesh->num_triangles;
320                 R_Mesh_Draw(mesh->num_vertices, mesh->num_triangles, mesh->data_elements);
321         }
322 }
323
324 void R_Model_Alias_Draw(entity_render_t *ent)
325 {
326         int meshnum;
327         aliasmesh_t *mesh;
328         if (ent->alpha < (1.0f / 64.0f))
329                 return; // basically completely transparent
330
331         c_models++;
332
333         for (meshnum = 0, mesh = ent->model->aliasdata_meshes;meshnum < ent->model->aliasnum_meshes;meshnum++, mesh++)
334         {
335                 if (ent->effects & EF_ADDITIVE || ent->alpha != 1.0 || R_FetchAliasSkin(ent, mesh)->flags & ALIASSKIN_TRANSPARENT)
336                         R_MeshQueue_AddTransparent(ent->origin, R_DrawAliasModelCallback, ent, meshnum);
337                 else
338                         R_DrawAliasModelCallback(ent, meshnum);
339         }
340 }
341
342 void R_Model_Alias_DrawFakeShadow (entity_render_t *ent)
343 {
344         int i, meshnum;
345         aliasmesh_t *mesh;
346         aliasskin_t *skin;
347         rmeshstate_t m;
348         float *v, planenormal[3], planedist, dist, projection[3], floororigin[3], surfnormal[3], lightdirection[3], v2[3];
349
350         if ((ent->effects & EF_ADDITIVE) || ent->alpha < 1)
351                 return;
352
353         lightdirection[0] = 0.5;
354         lightdirection[1] = 0.2;
355         lightdirection[2] = -1;
356         VectorNormalizeFast(lightdirection);
357
358         VectorMA(ent->origin, 65536.0f, lightdirection, v2);
359         if (CL_TraceLine(ent->origin, v2, floororigin, surfnormal, 0, false, NULL) == 1)
360                 return;
361
362         R_Mesh_Matrix(&ent->matrix);
363
364         memset(&m, 0, sizeof(m));
365         m.blendfunc1 = GL_SRC_ALPHA;
366         m.blendfunc2 = GL_ONE_MINUS_SRC_ALPHA;
367         R_Mesh_State(&m);
368         GL_Color(0, 0, 0, 0.5);
369
370         // put a light direction in the entity's coordinate space
371         Matrix4x4_Transform3x3(&ent->inversematrix, lightdirection, projection);
372         VectorNormalizeFast(projection);
373
374         // put the plane's normal in the entity's coordinate space
375         Matrix4x4_Transform3x3(&ent->inversematrix, surfnormal, planenormal);
376         VectorNormalizeFast(planenormal);
377
378         // put the plane's distance in the entity's coordinate space
379         VectorSubtract(floororigin, ent->origin, floororigin);
380         planedist = DotProduct(floororigin, surfnormal) + 2;
381
382         dist = -1.0f / DotProduct(projection, planenormal);
383         VectorScale(projection, dist, projection);
384         for (meshnum = 0, mesh = ent->model->aliasdata_meshes;meshnum < ent->model->aliasnum_meshes;meshnum++)
385         {
386                 skin = R_FetchAliasSkin(ent, mesh);
387                 if (skin->flags & ALIASSKIN_TRANSPARENT)
388                         continue;
389                 R_Mesh_ResizeCheck(mesh->num_vertices);
390                 R_Model_Alias_GetMeshVerts(ent, mesh, varray_vertex, NULL, NULL, NULL);
391                 for (i = 0, v = varray_vertex;i < mesh->num_vertices;i++, v += 4)
392                 {
393                         dist = DotProduct(v, planenormal) - planedist;
394                         if (dist > 0)
395                                 VectorMA(v, dist, projection, v);
396                 }
397                 c_alias_polys += mesh->num_triangles;
398                 R_Mesh_Draw(mesh->num_vertices, mesh->num_triangles, mesh->data_elements);
399         }
400 }
401
402 void R_Model_Alias_DrawShadowVolume(entity_render_t *ent, vec3_t relativelightorigin, float lightradius)
403 {
404         int meshnum;
405         aliasmesh_t *mesh;
406         aliasskin_t *skin;
407         float projectdistance;
408         if (ent->effects & EF_ADDITIVE || ent->alpha < 1)
409                 return;
410         projectdistance = lightradius + ent->model->radius - sqrt(DotProduct(relativelightorigin, relativelightorigin));
411         if (projectdistance > 0.1)
412         {
413                 R_Mesh_Matrix(&ent->matrix);
414                 for (meshnum = 0, mesh = ent->model->aliasdata_meshes;meshnum < ent->model->aliasnum_meshes;meshnum++)
415                 {
416                         skin = R_FetchAliasSkin(ent, mesh);
417                         if (skin->flags & ALIASSKIN_TRANSPARENT)
418                                 continue;
419                         R_Mesh_ResizeCheck(mesh->num_vertices * 2);
420                         R_Model_Alias_GetMeshVerts(ent, mesh, varray_vertex, NULL, NULL, NULL);
421                         R_Shadow_Volume(mesh->num_vertices, mesh->num_triangles, mesh->data_elements, mesh->data_neighbors, relativelightorigin, lightradius, projectdistance);
422                 }
423         }
424 }
425
426 void R_Model_Alias_DrawLight(entity_render_t *ent, vec3_t relativelightorigin, vec3_t relativeeyeorigin, float lightradius, float *lightcolor)
427 {
428         int c, meshnum, layernum;
429         float fog, ifog, lightcolor2[3];
430         vec3_t diff;
431         qbyte *bcolor;
432         aliasmesh_t *mesh;
433         aliaslayer_t *layer;
434         aliasskin_t *skin;
435
436         if (ent->effects & (EF_ADDITIVE | EF_FULLBRIGHT) || ent->alpha < 1)
437                 return;
438
439         R_Mesh_Matrix(&ent->matrix);
440
441         fog = 0;
442         if (fogenabled)
443         {
444                 VectorSubtract(ent->origin, r_origin, diff);
445                 fog = DotProduct(diff,diff);
446                 if (fog < 0.01f)
447                         fog = 0.01f;
448                 fog = exp(fogdensity/fog);
449                 if (fog > 1)
450                         fog = 1;
451                 if (fog < 0.01f)
452                         fog = 0;
453                 // fog method: darken, additive fog
454                 // 1. render model as normal, scaled by inverse of fog alpha (darkens it)
455                 // 2. render fog as additive
456         }
457         ifog = 1 - fog;
458
459         for (meshnum = 0, mesh = ent->model->aliasdata_meshes;meshnum < ent->model->aliasnum_meshes;meshnum++, mesh++)
460         {
461                 skin = R_FetchAliasSkin(ent, mesh);
462                 if (skin->flags & ALIASSKIN_TRANSPARENT)
463                         continue;
464                 R_Mesh_ResizeCheck(mesh->num_vertices);
465                 R_Model_Alias_GetMeshVerts(ent, mesh, varray_vertex, aliasvert_normals, aliasvert_svectors, aliasvert_tvectors);
466                 for (layernum = 0, layer = skin->data_layers;layernum < skin->num_layers;layernum++, layer++)
467                 {
468                         if (!(layer->flags & ALIASLAYER_DRAW_PER_LIGHT)
469                          || ((layer->flags & ALIASLAYER_NODRAW_IF_NOTCOLORMAPPED) && ent->colormap < 0)
470                          || ((layer->flags & ALIASLAYER_NODRAW_IF_COLORMAPPED) && ent->colormap >= 0))
471                                 continue;
472                         lightcolor2[0] = lightcolor[0] * ifog;
473                         lightcolor2[1] = lightcolor[1] * ifog;
474                         lightcolor2[2] = lightcolor[2] * ifog;
475                         if (layer->flags & ALIASLAYER_SPECULAR)
476                         {
477                                 c_alias_polys += mesh->num_triangles;
478                                 R_Shadow_SpecularLighting(mesh->num_vertices, mesh->num_triangles, mesh->data_elements, aliasvert_svectors, aliasvert_tvectors, aliasvert_normals, mesh->data_texcoords, relativelightorigin, relativeeyeorigin, lightradius, lightcolor2, layer->texture, layer->nmap, NULL);
479                         }
480                         else if (layer->flags & ALIASLAYER_DIFFUSE)
481                         {
482                                 if (layer->flags & ALIASLAYER_COLORMAP_PANTS)
483                                 {
484                                         // 128-224 are backwards ranges
485                                         c = (ent->colormap & 0xF) << 4;c += (c >= 128 && c < 224) ? 4 : 12;
486                                         // fullbright passes were already taken care of, so skip them in realtime lighting passes
487                                         if (c >= 224)
488                                                 continue;
489                                         bcolor = (qbyte *) (&palette_complete[c]);
490                                         lightcolor2[0] *= bcolor[0] * (1.0f / 255.0f);
491                                         lightcolor2[1] *= bcolor[1] * (1.0f / 255.0f);
492                                         lightcolor2[2] *= bcolor[2] * (1.0f / 255.0f);
493                                 }
494                                 else if (layer->flags & ALIASLAYER_COLORMAP_SHIRT)
495                                 {
496                                         // 128-224 are backwards ranges
497                                         c = (ent->colormap & 0xF0);c += (c >= 128 && c < 224) ? 4 : 12;
498                                         // fullbright passes were already taken care of, so skip them in realtime lighting passes
499                                         if (c >= 224)
500                                                 continue;
501                                         bcolor = (qbyte *) (&palette_complete[c]);
502                                         lightcolor2[0] *= bcolor[0] * (1.0f / 255.0f);
503                                         lightcolor2[1] *= bcolor[1] * (1.0f / 255.0f);
504                                         lightcolor2[2] *= bcolor[2] * (1.0f / 255.0f);
505                                 }
506                                 c_alias_polys += mesh->num_triangles;
507                                 R_Shadow_DiffuseLighting(mesh->num_vertices, mesh->num_triangles, mesh->data_elements, aliasvert_svectors, aliasvert_tvectors, aliasvert_normals, mesh->data_texcoords, relativelightorigin, lightradius, lightcolor2, layer->texture, layer->nmap, NULL);
508                         }
509                 }
510         }
511 }
512
513 int ZymoticLerpBones(int count, const zymbonematrix *bonebase, const frameblend_t *blend, const zymbone_t *bone)
514 {
515         int i;
516         float lerp1, lerp2, lerp3, lerp4;
517         zymbonematrix *out, rootmatrix, m;
518         const zymbonematrix *bone1, *bone2, *bone3, *bone4;
519
520         rootmatrix.m[0][0] = 1;
521         rootmatrix.m[0][1] = 0;
522         rootmatrix.m[0][2] = 0;
523         rootmatrix.m[0][3] = 0;
524         rootmatrix.m[1][0] = 0;
525         rootmatrix.m[1][1] = 1;
526         rootmatrix.m[1][2] = 0;
527         rootmatrix.m[1][3] = 0;
528         rootmatrix.m[2][0] = 0;
529         rootmatrix.m[2][1] = 0;
530         rootmatrix.m[2][2] = 1;
531         rootmatrix.m[2][3] = 0;
532
533         bone1 = bonebase + blend[0].frame * count;
534         lerp1 = blend[0].lerp;
535         if (blend[1].lerp)
536         {
537                 bone2 = bonebase + blend[1].frame * count;
538                 lerp2 = blend[1].lerp;
539                 if (blend[2].lerp)
540                 {
541                         bone3 = bonebase + blend[2].frame * count;
542                         lerp3 = blend[2].lerp;
543                         if (blend[3].lerp)
544                         {
545                                 // 4 poses
546                                 bone4 = bonebase + blend[3].frame * count;
547                                 lerp4 = blend[3].lerp;
548                                 for (i = 0, out = zymbonepose;i < count;i++, out++)
549                                 {
550                                         // interpolate matrices
551                                         m.m[0][0] = bone1->m[0][0] * lerp1 + bone2->m[0][0] * lerp2 + bone3->m[0][0] * lerp3 + bone4->m[0][0] * lerp4;
552                                         m.m[0][1] = bone1->m[0][1] * lerp1 + bone2->m[0][1] * lerp2 + bone3->m[0][1] * lerp3 + bone4->m[0][1] * lerp4;
553                                         m.m[0][2] = bone1->m[0][2] * lerp1 + bone2->m[0][2] * lerp2 + bone3->m[0][2] * lerp3 + bone4->m[0][2] * lerp4;
554                                         m.m[0][3] = bone1->m[0][3] * lerp1 + bone2->m[0][3] * lerp2 + bone3->m[0][3] * lerp3 + bone4->m[0][3] * lerp4;
555                                         m.m[1][0] = bone1->m[1][0] * lerp1 + bone2->m[1][0] * lerp2 + bone3->m[1][0] * lerp3 + bone4->m[1][0] * lerp4;
556                                         m.m[1][1] = bone1->m[1][1] * lerp1 + bone2->m[1][1] * lerp2 + bone3->m[1][1] * lerp3 + bone4->m[1][1] * lerp4;
557                                         m.m[1][2] = bone1->m[1][2] * lerp1 + bone2->m[1][2] * lerp2 + bone3->m[1][2] * lerp3 + bone4->m[1][2] * lerp4;
558                                         m.m[1][3] = bone1->m[1][3] * lerp1 + bone2->m[1][3] * lerp2 + bone3->m[1][3] * lerp3 + bone4->m[1][3] * lerp4;
559                                         m.m[2][0] = bone1->m[2][0] * lerp1 + bone2->m[2][0] * lerp2 + bone3->m[2][0] * lerp3 + bone4->m[2][0] * lerp4;
560                                         m.m[2][1] = bone1->m[2][1] * lerp1 + bone2->m[2][1] * lerp2 + bone3->m[2][1] * lerp3 + bone4->m[2][1] * lerp4;
561                                         m.m[2][2] = bone1->m[2][2] * lerp1 + bone2->m[2][2] * lerp2 + bone3->m[2][2] * lerp3 + bone4->m[2][2] * lerp4;
562                                         m.m[2][3] = bone1->m[2][3] * lerp1 + bone2->m[2][3] * lerp2 + bone3->m[2][3] * lerp3 + bone4->m[2][3] * lerp4;
563                                         if (bone->parent >= 0)
564                                                 R_ConcatTransforms(&zymbonepose[bone->parent].m[0][0], &m.m[0][0], &out->m[0][0]);
565                                         else
566                                                 R_ConcatTransforms(&rootmatrix.m[0][0], &m.m[0][0], &out->m[0][0]);
567                                         bone1++;
568                                         bone2++;
569                                         bone3++;
570                                         bone4++;
571                                         bone++;
572                                 }
573                         }
574                         else
575                         {
576                                 // 3 poses
577                                 for (i = 0, out = zymbonepose;i < count;i++, out++)
578                                 {
579                                         // interpolate matrices
580                                         m.m[0][0] = bone1->m[0][0] * lerp1 + bone2->m[0][0] * lerp2 + bone3->m[0][0] * lerp3;
581                                         m.m[0][1] = bone1->m[0][1] * lerp1 + bone2->m[0][1] * lerp2 + bone3->m[0][1] * lerp3;
582                                         m.m[0][2] = bone1->m[0][2] * lerp1 + bone2->m[0][2] * lerp2 + bone3->m[0][2] * lerp3;
583                                         m.m[0][3] = bone1->m[0][3] * lerp1 + bone2->m[0][3] * lerp2 + bone3->m[0][3] * lerp3;
584                                         m.m[1][0] = bone1->m[1][0] * lerp1 + bone2->m[1][0] * lerp2 + bone3->m[1][0] * lerp3;
585                                         m.m[1][1] = bone1->m[1][1] * lerp1 + bone2->m[1][1] * lerp2 + bone3->m[1][1] * lerp3;
586                                         m.m[1][2] = bone1->m[1][2] * lerp1 + bone2->m[1][2] * lerp2 + bone3->m[1][2] * lerp3;
587                                         m.m[1][3] = bone1->m[1][3] * lerp1 + bone2->m[1][3] * lerp2 + bone3->m[1][3] * lerp3;
588                                         m.m[2][0] = bone1->m[2][0] * lerp1 + bone2->m[2][0] * lerp2 + bone3->m[2][0] * lerp3;
589                                         m.m[2][1] = bone1->m[2][1] * lerp1 + bone2->m[2][1] * lerp2 + bone3->m[2][1] * lerp3;
590                                         m.m[2][2] = bone1->m[2][2] * lerp1 + bone2->m[2][2] * lerp2 + bone3->m[2][2] * lerp3;
591                                         m.m[2][3] = bone1->m[2][3] * lerp1 + bone2->m[2][3] * lerp2 + bone3->m[2][3] * lerp3;
592                                         if (bone->parent >= 0)
593                                                 R_ConcatTransforms(&zymbonepose[bone->parent].m[0][0], &m.m[0][0], &out->m[0][0]);
594                                         else
595                                                 R_ConcatTransforms(&rootmatrix.m[0][0], &m.m[0][0], &out->m[0][0]);
596                                         bone1++;
597                                         bone2++;
598                                         bone3++;
599                                         bone++;
600                                 }
601                         }
602                 }
603                 else
604                 {
605                         // 2 poses
606                         for (i = 0, out = zymbonepose;i < count;i++, out++)
607                         {
608                                 // interpolate matrices
609                                 m.m[0][0] = bone1->m[0][0] * lerp1 + bone2->m[0][0] * lerp2;
610                                 m.m[0][1] = bone1->m[0][1] * lerp1 + bone2->m[0][1] * lerp2;
611                                 m.m[0][2] = bone1->m[0][2] * lerp1 + bone2->m[0][2] * lerp2;
612                                 m.m[0][3] = bone1->m[0][3] * lerp1 + bone2->m[0][3] * lerp2;
613                                 m.m[1][0] = bone1->m[1][0] * lerp1 + bone2->m[1][0] * lerp2;
614                                 m.m[1][1] = bone1->m[1][1] * lerp1 + bone2->m[1][1] * lerp2;
615                                 m.m[1][2] = bone1->m[1][2] * lerp1 + bone2->m[1][2] * lerp2;
616                                 m.m[1][3] = bone1->m[1][3] * lerp1 + bone2->m[1][3] * lerp2;
617                                 m.m[2][0] = bone1->m[2][0] * lerp1 + bone2->m[2][0] * lerp2;
618                                 m.m[2][1] = bone1->m[2][1] * lerp1 + bone2->m[2][1] * lerp2;
619                                 m.m[2][2] = bone1->m[2][2] * lerp1 + bone2->m[2][2] * lerp2;
620                                 m.m[2][3] = bone1->m[2][3] * lerp1 + bone2->m[2][3] * lerp2;
621                                 if (bone->parent >= 0)
622                                         R_ConcatTransforms(&zymbonepose[bone->parent].m[0][0], &m.m[0][0], &out->m[0][0]);
623                                 else
624                                         R_ConcatTransforms(&rootmatrix.m[0][0], &m.m[0][0], &out->m[0][0]);
625                                 bone1++;
626                                 bone2++;
627                                 bone++;
628                         }
629                 }
630         }
631         else
632         {
633                 // 1 pose
634                 if (lerp1 != 1)
635                 {
636                         // lerp != 1.0
637                         for (i = 0, out = zymbonepose;i < count;i++, out++)
638                         {
639                                 // interpolate matrices
640                                 m.m[0][0] = bone1->m[0][0] * lerp1;
641                                 m.m[0][1] = bone1->m[0][1] * lerp1;
642                                 m.m[0][2] = bone1->m[0][2] * lerp1;
643                                 m.m[0][3] = bone1->m[0][3] * lerp1;
644                                 m.m[1][0] = bone1->m[1][0] * lerp1;
645                                 m.m[1][1] = bone1->m[1][1] * lerp1;
646                                 m.m[1][2] = bone1->m[1][2] * lerp1;
647                                 m.m[1][3] = bone1->m[1][3] * lerp1;
648                                 m.m[2][0] = bone1->m[2][0] * lerp1;
649                                 m.m[2][1] = bone1->m[2][1] * lerp1;
650                                 m.m[2][2] = bone1->m[2][2] * lerp1;
651                                 m.m[2][3] = bone1->m[2][3] * lerp1;
652                                 if (bone->parent >= 0)
653                                         R_ConcatTransforms(&zymbonepose[bone->parent].m[0][0], &m.m[0][0], &out->m[0][0]);
654                                 else
655                                         R_ConcatTransforms(&rootmatrix.m[0][0], &m.m[0][0], &out->m[0][0]);
656                                 bone1++;
657                                 bone++;
658                         }
659                 }
660                 else
661                 {
662                         // lerp == 1.0
663                         for (i = 0, out = zymbonepose;i < count;i++, out++)
664                         {
665                                 if (bone->parent >= 0)
666                                         R_ConcatTransforms(&zymbonepose[bone->parent].m[0][0], &bone1->m[0][0], &out->m[0][0]);
667                                 else
668                                         R_ConcatTransforms(&rootmatrix.m[0][0], &bone1->m[0][0], &out->m[0][0]);
669                                 bone1++;
670                                 bone++;
671                         }
672                 }
673         }
674         return true;
675 }
676
677 void ZymoticTransformVerts(int vertcount, float *vertex, int *bonecounts, zymvertex_t *vert)
678 {
679         int c;
680         float *out = vertex;
681         zymbonematrix *matrix;
682         while(vertcount--)
683         {
684                 c = *bonecounts++;
685                 // FIXME: validate bonecounts at load time (must be >= 1)
686                 // FIXME: need 4th component in origin, for how much of the translate to blend in
687                 if (c == 1)
688                 {
689                         matrix = &zymbonepose[vert->bonenum];
690                         out[0] = vert->origin[0] * matrix->m[0][0] + vert->origin[1] * matrix->m[0][1] + vert->origin[2] * matrix->m[0][2] + matrix->m[0][3];
691                         out[1] = vert->origin[0] * matrix->m[1][0] + vert->origin[1] * matrix->m[1][1] + vert->origin[2] * matrix->m[1][2] + matrix->m[1][3];
692                         out[2] = vert->origin[0] * matrix->m[2][0] + vert->origin[1] * matrix->m[2][1] + vert->origin[2] * matrix->m[2][2] + matrix->m[2][3];
693                         vert++;
694                 }
695                 else
696                 {
697                         VectorClear(out);
698                         while(c--)
699                         {
700                                 matrix = &zymbonepose[vert->bonenum];
701                                 out[0] += vert->origin[0] * matrix->m[0][0] + vert->origin[1] * matrix->m[0][1] + vert->origin[2] * matrix->m[0][2] + matrix->m[0][3];
702                                 out[1] += vert->origin[0] * matrix->m[1][0] + vert->origin[1] * matrix->m[1][1] + vert->origin[2] * matrix->m[1][2] + matrix->m[1][3];
703                                 out[2] += vert->origin[0] * matrix->m[2][0] + vert->origin[1] * matrix->m[2][1] + vert->origin[2] * matrix->m[2][2] + matrix->m[2][3];
704                                 vert++;
705                         }
706                 }
707                 out += 4;
708         }
709 }
710
711 void ZymoticCalcNormals(int vertcount, float *vertex, float *normals, int shadercount, int *renderlist)
712 {
713         int a, b, c, d;
714         float *out, v1[3], v2[3], normal[3], s;
715         int *u;
716         // clear normals
717         memset(normals, 0, sizeof(float) * vertcount * 3);
718         memset(aliasvertusage, 0, sizeof(int) * vertcount);
719         // parse render list and accumulate surface normals
720         while(shadercount--)
721         {
722                 d = *renderlist++;
723                 while (d--)
724                 {
725                         a = renderlist[0]*4;
726                         b = renderlist[1]*4;
727                         c = renderlist[2]*4;
728                         v1[0] = vertex[a+0] - vertex[b+0];
729                         v1[1] = vertex[a+1] - vertex[b+1];
730                         v1[2] = vertex[a+2] - vertex[b+2];
731                         v2[0] = vertex[c+0] - vertex[b+0];
732                         v2[1] = vertex[c+1] - vertex[b+1];
733                         v2[2] = vertex[c+2] - vertex[b+2];
734                         CrossProduct(v1, v2, normal);
735                         VectorNormalizeFast(normal);
736                         // add surface normal to vertices
737                         a = renderlist[0] * 3;
738                         normals[a+0] += normal[0];
739                         normals[a+1] += normal[1];
740                         normals[a+2] += normal[2];
741                         aliasvertusage[renderlist[0]]++;
742                         a = renderlist[1] * 3;
743                         normals[a+0] += normal[0];
744                         normals[a+1] += normal[1];
745                         normals[a+2] += normal[2];
746                         aliasvertusage[renderlist[1]]++;
747                         a = renderlist[2] * 3;
748                         normals[a+0] += normal[0];
749                         normals[a+1] += normal[1];
750                         normals[a+2] += normal[2];
751                         aliasvertusage[renderlist[2]]++;
752                         renderlist += 3;
753                 }
754         }
755         // FIXME: precalc this
756         // average surface normals
757         out = normals;
758         u = aliasvertusage;
759         while(vertcount--)
760         {
761                 if (*u > 1)
762                 {
763                         s = ixtable[*u];
764                         out[0] *= s;
765                         out[1] *= s;
766                         out[2] *= s;
767                 }
768                 u++;
769                 out += 3;
770         }
771 }
772
773 void R_DrawZymoticModelMeshCallback (const void *calldata1, int calldata2)
774 {
775         float fog, ifog, colorscale;
776         vec3_t diff;
777         int i, *renderlist, *elements;
778         rtexture_t *texture;
779         rmeshstate_t mstate;
780         const entity_render_t *ent = calldata1;
781         int shadernum = calldata2;
782         int numverts, numtriangles;
783
784         R_Mesh_Matrix(&ent->matrix);
785
786         // find the vertex index list and texture
787         renderlist = ent->model->zymdata_renderlist;
788         for (i = 0;i < shadernum;i++)
789                 renderlist += renderlist[0] * 3 + 1;
790         texture = ent->model->zymdata_textures[shadernum];
791
792         numverts = ent->model->zymnum_verts;
793         numtriangles = *renderlist++;
794         elements = renderlist;
795         R_Mesh_ResizeCheck(numverts);
796
797         fog = 0;
798         if (fogenabled)
799         {
800                 VectorSubtract(ent->origin, r_origin, diff);
801                 fog = DotProduct(diff,diff);
802                 if (fog < 0.01f)
803                         fog = 0.01f;
804                 fog = exp(fogdensity/fog);
805                 if (fog > 1)
806                         fog = 1;
807                 if (fog < 0.01f)
808                         fog = 0;
809                 // fog method: darken, additive fog
810                 // 1. render model as normal, scaled by inverse of fog alpha (darkens it)
811                 // 2. render fog as additive
812         }
813         ifog = 1 - fog;
814
815         memset(&mstate, 0, sizeof(mstate));
816         if (ent->effects & EF_ADDITIVE)
817         {
818                 mstate.blendfunc1 = GL_SRC_ALPHA;
819                 mstate.blendfunc2 = GL_ONE;
820         }
821         else if (ent->alpha != 1.0 || R_TextureHasAlpha(texture))
822         {
823                 mstate.blendfunc1 = GL_SRC_ALPHA;
824                 mstate.blendfunc2 = GL_ONE_MINUS_SRC_ALPHA;
825         }
826         else
827         {
828                 mstate.blendfunc1 = GL_ONE;
829                 mstate.blendfunc2 = GL_ZERO;
830         }
831         colorscale = r_colorscale;
832         if (gl_combine.integer)
833         {
834                 mstate.texrgbscale[0] = 4;
835                 colorscale *= 0.25f;
836         }
837         mstate.tex[0] = R_GetTexture(texture);
838         R_Mesh_State(&mstate);
839         ZymoticLerpBones(ent->model->zymnum_bones, (zymbonematrix *) ent->model->zymdata_poses, ent->frameblend, ent->model->zymdata_bones);
840         ZymoticTransformVerts(numverts, varray_vertex, ent->model->zymdata_vertbonecounts, ent->model->zymdata_verts);
841         ZymoticCalcNormals(numverts, varray_vertex, aliasvert_normals, ent->model->zymnum_shaders, ent->model->zymdata_renderlist);
842         memcpy(varray_texcoord[0], ent->model->zymdata_texcoords, ent->model->zymnum_verts * sizeof(float[4]));
843         GL_UseColorArray();
844         R_LightModel(ent, numverts, varray_vertex, aliasvert_normals, varray_color, ifog * colorscale, ifog * colorscale, ifog * colorscale, false);
845         R_Mesh_Draw(numverts, numtriangles, elements);
846         c_alias_polys += numtriangles;
847
848         if (fog)
849         {
850                 memset(&mstate, 0, sizeof(mstate));
851                 mstate.blendfunc1 = GL_SRC_ALPHA;
852                 mstate.blendfunc2 = GL_ONE_MINUS_SRC_ALPHA;
853                 // FIXME: need alpha mask for fogging...
854                 //mstate.tex[0] = R_GetTexture(texture);
855                 R_Mesh_State(&mstate);
856                 GL_Color(fogcolor[0] * r_colorscale, fogcolor[1] * r_colorscale, fogcolor[2] * r_colorscale, ent->alpha * fog);
857                 R_Mesh_Draw(numverts, numtriangles, elements);
858                 c_alias_polys += numtriangles;
859         }
860 }
861
862 void R_Model_Zymotic_Draw(entity_render_t *ent)
863 {
864         int i;
865
866         if (ent->alpha < (1.0f / 64.0f))
867                 return; // basically completely transparent
868
869         c_models++;
870
871         for (i = 0;i < ent->model->zymnum_shaders;i++)
872         {
873                 if (ent->effects & EF_ADDITIVE || ent->alpha != 1.0 || R_TextureHasAlpha(ent->model->zymdata_textures[i]))
874                         R_MeshQueue_AddTransparent(ent->origin, R_DrawZymoticModelMeshCallback, ent, i);
875                 else
876                         R_DrawZymoticModelMeshCallback(ent, i);
877         }
878 }
879
880 void R_Model_Zymotic_DrawFakeShadow(entity_render_t *ent)
881 {
882         // FIXME
883 }
884
885 void R_Model_Zymotic_DrawLight(entity_render_t *ent, vec3_t relativelightorigin, float lightradius2, float lightdistbias, float lightsubtract, float *lightcolor)
886 {
887         // FIXME
888 }
889
890 void R_Model_Zymotic_DrawOntoLight(entity_render_t *ent)
891 {
892         // FIXME
893 }