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