misc macro --> bli math lib functions
[blender.git] / source / blender / blenkernel / intern / DerivedMesh.c
1 /*
2  * ***** BEGIN GPL LICENSE BLOCK *****
3  *
4  * This program is free software; you can redistribute it and/or
5  * modify it under the terms of the GNU General Public License
6  * as published by the Free Software Foundation; either version 2
7  * of the License, or (at your option) any later version.
8  *
9  * This program is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12  * GNU General Public License for more details.
13  *
14  * You should have received a copy of the GNU General Public License
15  * along with this program; if not, write to the Free Software Foundation,
16  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
17  *
18  * The Original Code is Copyright (C) 2005 Blender Foundation.
19  * All rights reserved.
20  *
21  * The Original Code is: all of this file.
22  *
23  * Contributor(s): none yet.
24  *
25  * ***** END GPL LICENSE BLOCK *****
26  */
27
28 /** \file blender/blenkernel/intern/DerivedMesh.c
29  *  \ingroup bke
30  */
31
32
33 #include <string.h>
34
35
36 #include "MEM_guardedalloc.h"
37
38 #include "DNA_cloth_types.h"
39 #include "DNA_key_types.h"
40 #include "DNA_meshdata_types.h"
41 #include "DNA_armature_types.h"
42 #include "DNA_object_types.h"
43 #include "DNA_scene_types.h" // N_T
44
45 #include "BLI_blenlib.h"
46 #include "BLI_editVert.h"
47 #include "BLI_math.h"
48 #include "BLI_memarena.h"
49 #include "BLI_pbvh.h"
50 #include "BLI_utildefines.h"
51 #include "BLI_linklist.h"
52
53 #include "BKE_cdderivedmesh.h"
54 #include "BKE_displist.h"
55 #include "BKE_key.h"
56 #include "BKE_modifier.h"
57 #include "BKE_mesh.h"
58 #include "BKE_object.h"
59 #include "BKE_paint.h"
60 #include "BKE_texture.h"
61 #include "BKE_multires.h"
62 #include "BKE_armature.h"
63
64 #ifdef WITH_GAMEENGINE
65 #include "BKE_navmesh_conversion.h"
66 static DerivedMesh *navmesh_dm_createNavMeshForVisualization(DerivedMesh *dm);
67 #endif
68
69 #include "BLO_sys_types.h" // for intptr_t support
70
71 #include "GL/glew.h"
72
73 #include "GPU_buffers.h"
74 #include "GPU_draw.h"
75 #include "GPU_extensions.h"
76 #include "GPU_material.h"
77
78 extern GLubyte stipple_quarttone[128]; /* glutil.c, bad level data */
79
80 ///////////////////////////////////
81 ///////////////////////////////////
82
83 static MVert *dm_getVertArray(DerivedMesh *dm)
84 {
85         MVert *mvert = CustomData_get_layer(&dm->vertData, CD_MVERT);
86
87         if (!mvert) {
88                 mvert = CustomData_add_layer(&dm->vertData, CD_MVERT, CD_CALLOC, NULL,
89                         dm->getNumVerts(dm));
90                 CustomData_set_layer_flag(&dm->vertData, CD_MVERT, CD_FLAG_TEMPORARY);
91                 dm->copyVertArray(dm, mvert);
92         }
93
94         return mvert;
95 }
96
97 static MEdge *dm_getEdgeArray(DerivedMesh *dm)
98 {
99         MEdge *medge = CustomData_get_layer(&dm->edgeData, CD_MEDGE);
100
101         if (!medge) {
102                 medge = CustomData_add_layer(&dm->edgeData, CD_MEDGE, CD_CALLOC, NULL,
103                         dm->getNumEdges(dm));
104                 CustomData_set_layer_flag(&dm->edgeData, CD_MEDGE, CD_FLAG_TEMPORARY);
105                 dm->copyEdgeArray(dm, medge);
106         }
107
108         return medge;
109 }
110
111 static MFace *dm_getFaceArray(DerivedMesh *dm)
112 {
113         MFace *mface = CustomData_get_layer(&dm->faceData, CD_MFACE);
114
115         if (!mface) {
116                 mface = CustomData_add_layer(&dm->faceData, CD_MFACE, CD_CALLOC, NULL,
117                         dm->getNumFaces(dm));
118                 CustomData_set_layer_flag(&dm->faceData, CD_MFACE, CD_FLAG_TEMPORARY);
119                 dm->copyFaceArray(dm, mface);
120         }
121
122         return mface;
123 }
124
125 static MVert *dm_dupVertArray(DerivedMesh *dm)
126 {
127         MVert *tmp = MEM_callocN(sizeof(*tmp) * dm->getNumVerts(dm),
128                                                          "dm_dupVertArray tmp");
129
130         if(tmp) dm->copyVertArray(dm, tmp);
131
132         return tmp;
133 }
134
135 static MEdge *dm_dupEdgeArray(DerivedMesh *dm)
136 {
137         MEdge *tmp = MEM_callocN(sizeof(*tmp) * dm->getNumEdges(dm),
138                                                          "dm_dupEdgeArray tmp");
139
140         if(tmp) dm->copyEdgeArray(dm, tmp);
141
142         return tmp;
143 }
144
145 static MFace *dm_dupFaceArray(DerivedMesh *dm)
146 {
147         MFace *tmp = MEM_callocN(sizeof(*tmp) * dm->getNumFaces(dm),
148                                                          "dm_dupFaceArray tmp");
149
150         if(tmp) dm->copyFaceArray(dm, tmp);
151
152         return tmp;
153 }
154
155 void DM_init_funcs(DerivedMesh *dm)
156 {
157         /* default function implementations */
158         dm->getVertArray = dm_getVertArray;
159         dm->getEdgeArray = dm_getEdgeArray;
160         dm->getFaceArray = dm_getFaceArray;
161         dm->dupVertArray = dm_dupVertArray;
162         dm->dupEdgeArray = dm_dupEdgeArray;
163         dm->dupFaceArray = dm_dupFaceArray;
164
165         dm->getVertData = DM_get_vert_data;
166         dm->getEdgeData = DM_get_edge_data;
167         dm->getFaceData = DM_get_face_data;
168         dm->getVertDataArray = DM_get_vert_data_layer;
169         dm->getEdgeDataArray = DM_get_edge_data_layer;
170         dm->getFaceDataArray = DM_get_face_data_layer;
171
172         bvhcache_init(&dm->bvhCache);
173 }
174
175 void DM_init(DerivedMesh *dm, DerivedMeshType type,
176                          int numVerts, int numEdges, int numFaces)
177 {
178         dm->type = type;
179         dm->numVertData = numVerts;
180         dm->numEdgeData = numEdges;
181         dm->numFaceData = numFaces;
182
183         DM_init_funcs(dm);
184         
185         dm->needsFree = 1;
186 }
187
188 void DM_from_template(DerivedMesh *dm, DerivedMesh *source, DerivedMeshType type,
189                                           int numVerts, int numEdges, int numFaces)
190 {
191         CustomData_copy(&source->vertData, &dm->vertData, CD_MASK_DERIVEDMESH,
192                                         CD_CALLOC, numVerts);
193         CustomData_copy(&source->edgeData, &dm->edgeData, CD_MASK_DERIVEDMESH,
194                                         CD_CALLOC, numEdges);
195         CustomData_copy(&source->faceData, &dm->faceData, CD_MASK_DERIVEDMESH,
196                                         CD_CALLOC, numFaces);
197
198         dm->type = type;
199         dm->numVertData = numVerts;
200         dm->numEdgeData = numEdges;
201         dm->numFaceData = numFaces;
202
203         DM_init_funcs(dm);
204
205         dm->needsFree = 1;
206 }
207
208 int DM_release(DerivedMesh *dm)
209 {
210         if (dm->needsFree) {
211                 bvhcache_free(&dm->bvhCache);
212                 GPU_drawobject_free( dm );
213                 CustomData_free(&dm->vertData, dm->numVertData);
214                 CustomData_free(&dm->edgeData, dm->numEdgeData);
215                 CustomData_free(&dm->faceData, dm->numFaceData);
216
217                 return 1;
218         }
219         else {
220                 CustomData_free_temporary(&dm->vertData, dm->numVertData);
221                 CustomData_free_temporary(&dm->edgeData, dm->numEdgeData);
222                 CustomData_free_temporary(&dm->faceData, dm->numFaceData);
223
224                 return 0;
225         }
226 }
227
228 void DM_to_mesh(DerivedMesh *dm, Mesh *me)
229 {
230         /* dm might depend on me, so we need to do everything with a local copy */
231         Mesh tmp = *me;
232         int totvert, totedge, totface;
233
234         memset(&tmp.vdata, 0, sizeof(tmp.vdata));
235         memset(&tmp.edata, 0, sizeof(tmp.edata));
236         memset(&tmp.fdata, 0, sizeof(tmp.fdata));
237
238         totvert = tmp.totvert = dm->getNumVerts(dm);
239         totedge = tmp.totedge = dm->getNumEdges(dm);
240         totface = tmp.totface = dm->getNumFaces(dm);
241
242         CustomData_copy(&dm->vertData, &tmp.vdata, CD_MASK_MESH, CD_DUPLICATE, totvert);
243         CustomData_copy(&dm->edgeData, &tmp.edata, CD_MASK_MESH, CD_DUPLICATE, totedge);
244         CustomData_copy(&dm->faceData, &tmp.fdata, CD_MASK_MESH, CD_DUPLICATE, totface);
245
246         /* not all DerivedMeshes store their verts/edges/faces in CustomData, so
247            we set them here in case they are missing */
248         if(!CustomData_has_layer(&tmp.vdata, CD_MVERT))
249                 CustomData_add_layer(&tmp.vdata, CD_MVERT, CD_ASSIGN, dm->dupVertArray(dm), totvert);
250         if(!CustomData_has_layer(&tmp.edata, CD_MEDGE))
251                 CustomData_add_layer(&tmp.edata, CD_MEDGE, CD_ASSIGN, dm->dupEdgeArray(dm), totedge);
252         if(!CustomData_has_layer(&tmp.fdata, CD_MFACE))
253                 CustomData_add_layer(&tmp.fdata, CD_MFACE, CD_ASSIGN, dm->dupFaceArray(dm), totface);
254
255         /* object had got displacement layer, should copy this layer to save sculpted data */
256         /* NOTE: maybe some other layers should be copied? nazgul */
257         if(CustomData_has_layer(&me->fdata, CD_MDISPS)) {
258                 if (totface == me->totface) {
259                         MDisps *mdisps = CustomData_get_layer(&me->fdata, CD_MDISPS);
260                         CustomData_add_layer(&tmp.fdata, CD_MDISPS, CD_DUPLICATE, mdisps, totface);
261                 }
262         }
263
264         mesh_update_customdata_pointers(&tmp);
265
266         CustomData_free(&me->vdata, me->totvert);
267         CustomData_free(&me->edata, me->totedge);
268         CustomData_free(&me->fdata, me->totface);
269
270         /* if the number of verts has changed, remove invalid data */
271         if(tmp.totvert != me->totvert) {
272                 if(tmp.key) tmp.key->id.us--;
273                 tmp.key = NULL;
274         }
275
276         *me = tmp;
277 }
278
279 void DM_to_meshkey(DerivedMesh *dm, Mesh *me, KeyBlock *kb)
280 {
281         int a, totvert = dm->getNumVerts(dm);
282         float *fp;
283         MVert *mvert;
284         
285         if(totvert==0 || me->totvert==0 || me->totvert!=totvert) return;
286         
287         if(kb->data) MEM_freeN(kb->data);
288         kb->data= MEM_callocN(me->key->elemsize*me->totvert, "kb->data");
289         kb->totelem= totvert;
290         
291         fp= kb->data;
292         mvert=dm->getVertDataArray(dm, CD_MVERT);
293         
294         for(a=0; a<kb->totelem; a++, fp+=3, mvert++) {
295                 VECCOPY(fp, mvert->co);
296         }
297 }
298
299 void DM_set_only_copy(DerivedMesh *dm, CustomDataMask mask)
300 {
301         CustomData_set_only_copy(&dm->vertData, mask);
302         CustomData_set_only_copy(&dm->edgeData, mask);
303         CustomData_set_only_copy(&dm->faceData, mask);
304 }
305
306 void DM_add_vert_layer(DerivedMesh *dm, int type, int alloctype, void *layer)
307 {
308         CustomData_add_layer(&dm->vertData, type, alloctype, layer, dm->numVertData);
309 }
310
311 void DM_add_edge_layer(DerivedMesh *dm, int type, int alloctype, void *layer)
312 {
313         CustomData_add_layer(&dm->edgeData, type, alloctype, layer, dm->numEdgeData);
314 }
315
316 void DM_add_face_layer(DerivedMesh *dm, int type, int alloctype, void *layer)
317 {
318         CustomData_add_layer(&dm->faceData, type, alloctype, layer, dm->numFaceData);
319 }
320
321 void *DM_get_vert_data(DerivedMesh *dm, int index, int type)
322 {
323         return CustomData_get(&dm->vertData, index, type);
324 }
325
326 void *DM_get_edge_data(DerivedMesh *dm, int index, int type)
327 {
328         return CustomData_get(&dm->edgeData, index, type);
329 }
330
331 void *DM_get_face_data(DerivedMesh *dm, int index, int type)
332 {
333         return CustomData_get(&dm->faceData, index, type);
334 }
335
336 void *DM_get_vert_data_layer(DerivedMesh *dm, int type)
337 {
338         if(type == CD_MVERT)
339                 return dm->getVertArray(dm);
340
341         return CustomData_get_layer(&dm->vertData, type);
342 }
343
344 void *DM_get_edge_data_layer(DerivedMesh *dm, int type)
345 {
346         if(type == CD_MEDGE)
347                 return dm->getEdgeArray(dm);
348
349         return CustomData_get_layer(&dm->edgeData, type);
350 }
351
352 void *DM_get_face_data_layer(DerivedMesh *dm, int type)
353 {
354         if(type == CD_MFACE)
355                 return dm->getFaceArray(dm);
356
357         return CustomData_get_layer(&dm->faceData, type);
358 }
359
360 void DM_set_vert_data(DerivedMesh *dm, int index, int type, void *data)
361 {
362         CustomData_set(&dm->vertData, index, type, data);
363 }
364
365 void DM_set_edge_data(DerivedMesh *dm, int index, int type, void *data)
366 {
367         CustomData_set(&dm->edgeData, index, type, data);
368 }
369
370 void DM_set_face_data(DerivedMesh *dm, int index, int type, void *data)
371 {
372         CustomData_set(&dm->faceData, index, type, data);
373 }
374
375 void DM_copy_vert_data(DerivedMesh *source, DerivedMesh *dest,
376                                            int source_index, int dest_index, int count)
377 {
378         CustomData_copy_data(&source->vertData, &dest->vertData,
379                                                  source_index, dest_index, count);
380 }
381
382 void DM_copy_edge_data(DerivedMesh *source, DerivedMesh *dest,
383                                            int source_index, int dest_index, int count)
384 {
385         CustomData_copy_data(&source->edgeData, &dest->edgeData,
386                                                  source_index, dest_index, count);
387 }
388
389 void DM_copy_face_data(DerivedMesh *source, DerivedMesh *dest,
390                                            int source_index, int dest_index, int count)
391 {
392         CustomData_copy_data(&source->faceData, &dest->faceData,
393                                                  source_index, dest_index, count);
394 }
395
396 void DM_free_vert_data(struct DerivedMesh *dm, int index, int count)
397 {
398         CustomData_free_elem(&dm->vertData, index, count);
399 }
400
401 void DM_free_edge_data(struct DerivedMesh *dm, int index, int count)
402 {
403         CustomData_free_elem(&dm->edgeData, index, count);
404 }
405
406 void DM_free_face_data(struct DerivedMesh *dm, int index, int count)
407 {
408         CustomData_free_elem(&dm->faceData, index, count);
409 }
410
411 void DM_interp_vert_data(DerivedMesh *source, DerivedMesh *dest,
412                                                  int *src_indices, float *weights,
413                                                  int count, int dest_index)
414 {
415         CustomData_interp(&source->vertData, &dest->vertData, src_indices,
416                                           weights, NULL, count, dest_index);
417 }
418
419 void DM_interp_edge_data(DerivedMesh *source, DerivedMesh *dest,
420                                                  int *src_indices,
421                                                  float *weights, EdgeVertWeight *vert_weights,
422                                                  int count, int dest_index)
423 {
424         CustomData_interp(&source->edgeData, &dest->edgeData, src_indices,
425                                           weights, (float*)vert_weights, count, dest_index);
426 }
427
428 void DM_interp_face_data(DerivedMesh *source, DerivedMesh *dest,
429                                                  int *src_indices,
430                                                  float *weights, FaceVertWeight *vert_weights,
431                                                  int count, int dest_index)
432 {
433         CustomData_interp(&source->faceData, &dest->faceData, src_indices,
434                                           weights, (float*)vert_weights, count, dest_index);
435 }
436
437 void DM_swap_face_data(DerivedMesh *dm, int index, const int *corner_indices)
438 {
439         CustomData_swap(&dm->faceData, index, corner_indices);
440 }
441
442 ///
443
444 DerivedMesh *mesh_create_derived(Mesh *me, Object *ob, float (*vertCos)[3])
445 {
446         DerivedMesh *dm = CDDM_from_mesh(me, ob);
447         
448         if(!dm)
449                 return NULL;
450         
451         if (vertCos)
452                 CDDM_apply_vert_coords(dm, vertCos);
453
454         CDDM_calc_normals(dm);
455
456         return dm;
457 }
458
459 ///
460
461 typedef struct {
462         DerivedMesh dm;
463
464         EditMesh *em;
465         float (*vertexCos)[3];
466         float (*vertexNos)[3];
467         float (*faceNos)[3];
468 } EditMeshDerivedMesh;
469
470 static void emDM_foreachMappedVert(DerivedMesh *dm, void (*func)(void *userData, int index, float *co, float *no_f, short *no_s), void *userData)
471 {
472         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
473         EditVert *eve;
474         int i;
475
476         for (i=0,eve= emdm->em->verts.first; eve; i++,eve=eve->next) {
477                 if (emdm->vertexCos) {
478                         func(userData, i, emdm->vertexCos[i], emdm->vertexNos[i], NULL);
479                 } else {
480                         func(userData, i, eve->co, eve->no, NULL);
481                 }
482         }
483 }
484 static void emDM_foreachMappedEdge(DerivedMesh *dm, void (*func)(void *userData, int index, float *v0co, float *v1co), void *userData)
485 {
486         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
487         EditEdge *eed;
488         int i;
489
490         if (emdm->vertexCos) {
491                 EditVert *eve;
492
493                 for (i=0,eve=emdm->em->verts.first; eve; eve= eve->next)
494                         eve->tmp.l = (intptr_t) i++;
495                 for(i=0,eed= emdm->em->edges.first; eed; i++,eed= eed->next)
496                         func(userData, i, emdm->vertexCos[(int) eed->v1->tmp.l], emdm->vertexCos[(int) eed->v2->tmp.l]);
497         } else {
498                 for(i=0,eed= emdm->em->edges.first; eed; i++,eed= eed->next)
499                         func(userData, i, eed->v1->co, eed->v2->co);
500         }
501 }
502 static void emDM_drawMappedEdges(DerivedMesh *dm, int (*setDrawOptions)(void *userData, int index), void *userData) 
503 {
504         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
505         EditEdge *eed;
506         int i;
507
508         if (emdm->vertexCos) {
509                 EditVert *eve;
510
511                 for (i=0,eve=emdm->em->verts.first; eve; eve= eve->next)
512                         eve->tmp.l = (intptr_t) i++;
513
514                 glBegin(GL_LINES);
515                 for(i=0,eed= emdm->em->edges.first; eed; i++,eed= eed->next) {
516                         if(!setDrawOptions || setDrawOptions(userData, i)) {
517                                 glVertex3fv(emdm->vertexCos[(int) eed->v1->tmp.l]);
518                                 glVertex3fv(emdm->vertexCos[(int) eed->v2->tmp.l]);
519                         }
520                 }
521                 glEnd();
522         } else {
523                 glBegin(GL_LINES);
524                 for(i=0,eed= emdm->em->edges.first; eed; i++,eed= eed->next) {
525                         if(!setDrawOptions || setDrawOptions(userData, i)) {
526                                 glVertex3fv(eed->v1->co);
527                                 glVertex3fv(eed->v2->co);
528                         }
529                 }
530                 glEnd();
531         }
532 }
533 static void emDM_drawEdges(DerivedMesh *dm, int UNUSED(drawLooseEdges), int UNUSED(drawAllEdges))
534 {
535         emDM_drawMappedEdges(dm, NULL, NULL);
536 }
537 static void emDM_drawMappedEdgesInterp(DerivedMesh *dm, int (*setDrawOptions)(void *userData, int index), void (*setDrawInterpOptions)(void *userData, int index, float t), void *userData) 
538 {
539         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
540         EditEdge *eed;
541         int i;
542
543         if (emdm->vertexCos) {
544                 EditVert *eve;
545
546                 for (i=0,eve=emdm->em->verts.first; eve; eve= eve->next)
547                         eve->tmp.l = (intptr_t) i++;
548
549                 glBegin(GL_LINES);
550                 for (i=0,eed= emdm->em->edges.first; eed; i++,eed= eed->next) {
551                         if(!setDrawOptions || setDrawOptions(userData, i)) {
552                                 setDrawInterpOptions(userData, i, 0.0);
553                                 glVertex3fv(emdm->vertexCos[(int) eed->v1->tmp.l]);
554                                 setDrawInterpOptions(userData, i, 1.0);
555                                 glVertex3fv(emdm->vertexCos[(int) eed->v2->tmp.l]);
556                         }
557                 }
558                 glEnd();
559         } else {
560                 glBegin(GL_LINES);
561                 for (i=0,eed= emdm->em->edges.first; eed; i++,eed= eed->next) {
562                         if(!setDrawOptions || setDrawOptions(userData, i)) {
563                                 setDrawInterpOptions(userData, i, 0.0);
564                                 glVertex3fv(eed->v1->co);
565                                 setDrawInterpOptions(userData, i, 1.0);
566                                 glVertex3fv(eed->v2->co);
567                         }
568                 }
569                 glEnd();
570         }
571 }
572
573 static void emDM_drawUVEdges(DerivedMesh *dm)
574 {
575         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
576         EditFace *efa;
577         MTFace *tf;
578
579         glBegin(GL_LINES);
580         for(efa= emdm->em->faces.first; efa; efa= efa->next) {
581                 tf = CustomData_em_get(&emdm->em->fdata, efa->data, CD_MTFACE);
582
583                 if(tf && !(efa->h)) {
584                         glVertex2fv(tf->uv[0]);
585                         glVertex2fv(tf->uv[1]);
586
587                         glVertex2fv(tf->uv[1]);
588                         glVertex2fv(tf->uv[2]);
589
590                         if (!efa->v4) {
591                                 glVertex2fv(tf->uv[2]);
592                                 glVertex2fv(tf->uv[0]);
593                         } else {
594                                 glVertex2fv(tf->uv[2]);
595                                 glVertex2fv(tf->uv[3]);
596                                 glVertex2fv(tf->uv[3]);
597                                 glVertex2fv(tf->uv[0]);
598                         }
599                 }
600         }
601         glEnd();
602 }
603
604 static void emDM__calcFaceCent(EditFace *efa, float cent[3], float (*vertexCos)[3])
605 {
606         if (vertexCos) {
607                 VECCOPY(cent, vertexCos[(int) efa->v1->tmp.l]);
608                 add_v3_v3(cent, vertexCos[(int) efa->v2->tmp.l]);
609                 add_v3_v3(cent, vertexCos[(int) efa->v3->tmp.l]);
610                 if (efa->v4) add_v3_v3(cent, vertexCos[(int) efa->v4->tmp.l]);
611         } else {
612                 VECCOPY(cent, efa->v1->co);
613                 add_v3_v3(cent, efa->v2->co);
614                 add_v3_v3(cent, efa->v3->co);
615                 if (efa->v4) add_v3_v3(cent, efa->v4->co);
616         }
617
618         if (efa->v4) {
619                 mul_v3_fl(cent, 0.25f);
620         } else {
621                 mul_v3_fl(cent, 0.33333333333f);
622         }
623 }
624 static void emDM_foreachMappedFaceCenter(DerivedMesh *dm, void (*func)(void *userData, int index, float *co, float *no), void *userData)
625 {
626         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
627         EditVert *eve;
628         EditFace *efa;
629         float cent[3];
630         int i;
631
632         if (emdm->vertexCos) {
633                 for (i=0,eve=emdm->em->verts.first; eve; eve= eve->next)
634                         eve->tmp.l = (intptr_t) i++;
635         }
636
637         for(i=0,efa= emdm->em->faces.first; efa; i++,efa= efa->next) {
638                 emDM__calcFaceCent(efa, cent, emdm->vertexCos);
639                 func(userData, i, cent, emdm->vertexCos?emdm->faceNos[i]:efa->n);
640         }
641 }
642
643 /* note, material function is ignored for now. */
644 static void emDM_drawMappedFaces(DerivedMesh *dm, int (*setDrawOptions)(void *userData, int index, int *drawSmooth_r), void *userData, int UNUSED(useColors), int (*setMaterial)(int, void *attribs),
645                         int (*compareDrawOptions)(void *userData, int cur_index, int next_index))
646 {
647         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
648         EditFace *efa;
649         int i, draw;
650         const int skip_normals= !glIsEnabled(GL_LIGHTING); /* could be passed as an arg */
651
652         /* GL_ZERO is used to detect if drawing has started or not */
653         GLenum poly_prev= GL_ZERO;
654         GLenum shade_prev= GL_ZERO;
655
656         (void)setMaterial; /* unused */
657
658         /* currently unused -- each original face is handled separately */
659         (void)compareDrawOptions;
660
661         if (emdm->vertexCos) {
662                 /* add direct access */
663                 float (*vertexCos)[3]= emdm->vertexCos;
664                 float (*vertexNos)[3]= emdm->vertexNos;
665                 float (*faceNos)[3]=   emdm->faceNos;
666                 EditVert *eve;
667
668                 for (i=0,eve=emdm->em->verts.first; eve; eve= eve->next)
669                         eve->tmp.l = (intptr_t) i++;
670
671                 for (i=0,efa= emdm->em->faces.first; efa; i++,efa= efa->next) {
672                         int drawSmooth = (efa->flag & ME_SMOOTH);
673                         draw = setDrawOptions==NULL ? 1 : setDrawOptions(userData, i, &drawSmooth);
674                         if(draw) {
675                                 const GLenum poly_type= efa->v4 ? GL_QUADS:GL_TRIANGLES;
676                                 if (draw==2) { /* enabled with stipple */
677
678                                         if(poly_prev != GL_ZERO) glEnd();
679                                         poly_prev= GL_ZERO; /* force glBegin */
680
681                                         glEnable(GL_POLYGON_STIPPLE);
682                                         glPolygonStipple(stipple_quarttone);
683                                 }
684                                 
685                                 if(skip_normals) {
686                                         if(poly_type != poly_prev) {
687                                                 if(poly_prev != GL_ZERO) glEnd();
688                                                 glBegin((poly_prev= poly_type));
689                                         }
690                                         glVertex3fv(vertexCos[(int) efa->v1->tmp.l]);
691                                         glVertex3fv(vertexCos[(int) efa->v2->tmp.l]);
692                                         glVertex3fv(vertexCos[(int) efa->v3->tmp.l]);
693                                         if(poly_type == GL_QUADS) glVertex3fv(vertexCos[(int) efa->v4->tmp.l]);
694                                 }
695                                 else {
696                                         const GLenum shade_type= drawSmooth ? GL_SMOOTH : GL_FLAT;
697                                         if (shade_type != shade_prev) {
698                                                 if(poly_prev != GL_ZERO) glEnd();
699                                                 glShadeModel((shade_prev= shade_type)); /* same as below but switch shading */
700                                                 glBegin((poly_prev= poly_type));
701                                         }
702                                         else if(poly_type != poly_prev) {
703                                                 if(poly_prev != GL_ZERO) glEnd();
704                                                 glBegin((poly_prev= poly_type));
705                                         }
706
707                                         if (!drawSmooth) {
708                                                 glNormal3fv(faceNos[i]);
709                                                 glVertex3fv(vertexCos[(int) efa->v1->tmp.l]);
710                                                 glVertex3fv(vertexCos[(int) efa->v2->tmp.l]);
711                                                 glVertex3fv(vertexCos[(int) efa->v3->tmp.l]);
712                                                 if(poly_type == GL_QUADS) glVertex3fv(vertexCos[(int) efa->v4->tmp.l]);
713                                         } else {
714                                                 glNormal3fv(vertexNos[(int) efa->v1->tmp.l]);
715                                                 glVertex3fv(vertexCos[(int) efa->v1->tmp.l]);
716                                                 glNormal3fv(vertexNos[(int) efa->v2->tmp.l]);
717                                                 glVertex3fv(vertexCos[(int) efa->v2->tmp.l]);
718                                                 glNormal3fv(vertexNos[(int) efa->v3->tmp.l]);
719                                                 glVertex3fv(vertexCos[(int) efa->v3->tmp.l]);
720                                                 if(poly_type == GL_QUADS) {
721                                                         glNormal3fv(vertexNos[(int) efa->v4->tmp.l]);
722                                                         glVertex3fv(vertexCos[(int) efa->v4->tmp.l]);
723                                                 }
724                                         }
725                                 }
726
727                                 
728                                 if (draw==2) {
729                                         glEnd();
730                                         poly_prev= GL_ZERO; /* force glBegin */
731
732                                         glDisable(GL_POLYGON_STIPPLE);
733                                 }
734                         }
735                 }
736         }
737         else {
738                 for (i=0,efa= emdm->em->faces.first; efa; i++,efa= efa->next) {
739                         int drawSmooth = (efa->flag & ME_SMOOTH);
740                         draw = setDrawOptions==NULL ? 1 : setDrawOptions(userData, i, &drawSmooth);
741                         if(draw) {
742                                 const GLenum poly_type= efa->v4 ? GL_QUADS:GL_TRIANGLES;
743                                 if (draw==2) { /* enabled with stipple */
744
745                                         if(poly_prev != GL_ZERO) glEnd();
746                                         poly_prev= GL_ZERO; /* force glBegin */
747
748                                         glEnable(GL_POLYGON_STIPPLE);
749                                         glPolygonStipple(stipple_quarttone);
750                                 }
751
752                                 if(skip_normals) {
753                                         if(poly_type != poly_prev) {
754                                                 if(poly_prev != GL_ZERO) glEnd();
755                                                 glBegin((poly_prev= poly_type));
756                                         }
757                                         glVertex3fv(efa->v1->co);
758                                         glVertex3fv(efa->v2->co);
759                                         glVertex3fv(efa->v3->co);
760                                         if(poly_type == GL_QUADS) glVertex3fv(efa->v4->co);
761                                 }
762                                 else {
763                                         const GLenum shade_type= drawSmooth ? GL_SMOOTH : GL_FLAT;
764                                         if (shade_type != shade_prev) {
765                                                 if(poly_prev != GL_ZERO) glEnd();
766                                                 glShadeModel((shade_prev= shade_type)); /* same as below but switch shading */
767                                                 glBegin((poly_prev= poly_type));
768                                         }
769                                         else if(poly_type != poly_prev) {
770                                                 if(poly_prev != GL_ZERO) glEnd();
771                                                 glBegin((poly_prev= poly_type));
772                                         }
773
774                                         if (!drawSmooth) {
775                                                 glNormal3fv(efa->n);
776                                                 glVertex3fv(efa->v1->co);
777                                                 glVertex3fv(efa->v2->co);
778                                                 glVertex3fv(efa->v3->co);
779                                                 if(poly_type == GL_QUADS) glVertex3fv(efa->v4->co);
780                                         } else {
781                                                 glNormal3fv(efa->v1->no);
782                                                 glVertex3fv(efa->v1->co);
783                                                 glNormal3fv(efa->v2->no);
784                                                 glVertex3fv(efa->v2->co);
785                                                 glNormal3fv(efa->v3->no);
786                                                 glVertex3fv(efa->v3->co);
787                                                 if(poly_type == GL_QUADS) {
788                                                         glNormal3fv(efa->v4->no);
789                                                         glVertex3fv(efa->v4->co);
790                                                 }
791                                         }
792                                 }
793
794                                 
795                                 if (draw==2) {
796                                         glEnd();
797                                         poly_prev= GL_ZERO;
798
799                                         glDisable(GL_POLYGON_STIPPLE);
800                                 }
801                         }
802                 }
803         }
804
805         /* if non zero we know a face was rendered */
806         if(poly_prev != GL_ZERO) glEnd();
807 }
808
809 static void emDM_drawFacesTex_common(DerivedMesh *dm,
810                            int (*drawParams)(MTFace *tface, int has_mcol, int matnr),
811                            int (*drawParamsMapped)(void *userData, int index),
812                            void *userData) 
813 {
814         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
815         EditMesh *em= emdm->em;
816         float (*vertexCos)[3]= emdm->vertexCos;
817         float (*vertexNos)[3]= emdm->vertexNos;
818         EditFace *efa;
819         int i;
820
821         /* always use smooth shading even for flat faces, else vertex colors wont interpolate */
822         glShadeModel(GL_SMOOTH);
823         
824         if (vertexCos) {
825                 EditVert *eve;
826
827                 for (i=0,eve=em->verts.first; eve; eve= eve->next)
828                         eve->tmp.l = (intptr_t) i++;
829
830                 for (i=0,efa= em->faces.first; efa; i++,efa= efa->next) {
831                         MTFace *tf= CustomData_em_get(&em->fdata, efa->data, CD_MTFACE);
832                         MCol *mcol= CustomData_em_get(&em->fdata, efa->data, CD_MCOL);
833                         unsigned char *cp= NULL;
834                         int drawSmooth= (efa->flag & ME_SMOOTH);
835                         int flag;
836
837                         if(drawParams)
838                                 flag= drawParams(tf, (mcol != NULL), efa->mat_nr);
839                         else if(drawParamsMapped)
840                                 flag= drawParamsMapped(userData, i);
841                         else
842                                 flag= 1;
843
844                         if(flag != 0) { /* flag 0 == the face is hidden or invisible */
845                                 
846                                 /* we always want smooth here since otherwise vertex colors dont interpolate */
847                                 if (mcol) {
848                                         if (flag==1) {
849                                                 cp= (unsigned char*)mcol;
850                                         }
851                                 } else {
852                                         glShadeModel(drawSmooth?GL_SMOOTH:GL_FLAT);
853                                 } 
854                                 
855                                 glBegin(efa->v4?GL_QUADS:GL_TRIANGLES);
856                                 if (!drawSmooth) {
857                                         glNormal3fv(emdm->faceNos[i]);
858
859                                         if(tf) glTexCoord2fv(tf->uv[0]);
860                                         if(cp) glColor3ub(cp[3], cp[2], cp[1]);
861                                         glVertex3fv(vertexCos[(int) efa->v1->tmp.l]);
862
863                                         if(tf) glTexCoord2fv(tf->uv[1]);
864                                         if(cp) glColor3ub(cp[7], cp[6], cp[5]);
865                                         glVertex3fv(vertexCos[(int) efa->v2->tmp.l]);
866
867                                         if(tf) glTexCoord2fv(tf->uv[2]);
868                                         if(cp) glColor3ub(cp[11], cp[10], cp[9]);
869                                         glVertex3fv(vertexCos[(int) efa->v3->tmp.l]);
870
871                                         if(efa->v4) {
872                                                 if(tf) glTexCoord2fv(tf->uv[3]);
873                                                 if(cp) glColor3ub(cp[15], cp[14], cp[13]);
874                                                 glVertex3fv(vertexCos[(int) efa->v4->tmp.l]);
875                                         }
876                                 } else {
877                                         if(tf) glTexCoord2fv(tf->uv[0]);
878                                         if(cp) glColor3ub(cp[3], cp[2], cp[1]);
879                                         glNormal3fv(vertexNos[(int) efa->v1->tmp.l]);
880                                         glVertex3fv(vertexCos[(int) efa->v1->tmp.l]);
881
882                                         if(tf) glTexCoord2fv(tf->uv[1]);
883                                         if(cp) glColor3ub(cp[7], cp[6], cp[5]);
884                                         glNormal3fv(vertexNos[(int) efa->v2->tmp.l]);
885                                         glVertex3fv(vertexCos[(int) efa->v2->tmp.l]);
886
887                                         if(tf) glTexCoord2fv(tf->uv[2]);
888                                         if(cp) glColor3ub(cp[11], cp[10], cp[9]);
889                                         glNormal3fv(vertexNos[(int) efa->v3->tmp.l]);
890                                         glVertex3fv(vertexCos[(int) efa->v3->tmp.l]);
891
892                                         if(efa->v4) {
893                                                 if(tf) glTexCoord2fv(tf->uv[3]);
894                                                 if(cp) glColor3ub(cp[15], cp[14], cp[13]);
895                                                 glNormal3fv(vertexNos[(int) efa->v4->tmp.l]);
896                                                 glVertex3fv(vertexCos[(int) efa->v4->tmp.l]);
897                                         }
898                                 }
899                                 glEnd();
900                         }
901                 }
902         } else {
903                 for (i=0,efa= em->faces.first; efa; i++,efa= efa->next) {
904                         MTFace *tf= CustomData_em_get(&em->fdata, efa->data, CD_MTFACE);
905                         MCol *mcol= CustomData_em_get(&em->fdata, efa->data, CD_MCOL);
906                         unsigned char *cp= NULL;
907                         int drawSmooth= (efa->flag & ME_SMOOTH);
908                         int flag;
909
910                         if(drawParams)
911                                 flag= drawParams(tf, (mcol != NULL), efa->mat_nr);
912                         else if(drawParamsMapped)
913                                 flag= drawParamsMapped(userData, i);
914                         else
915                                 flag= 1;
916
917                         if(flag != 0) { /* flag 0 == the face is hidden or invisible */
918                                 /* we always want smooth here since otherwise vertex colors dont interpolate */
919                                 if (mcol) {
920                                         if (flag==1) {
921                                                 cp= (unsigned char*)mcol;
922                                         }
923                                 } else {
924                                         glShadeModel(drawSmooth?GL_SMOOTH:GL_FLAT);
925                                 } 
926
927                                 glBegin(efa->v4?GL_QUADS:GL_TRIANGLES);
928                                 if (!drawSmooth) {
929                                         glNormal3fv(efa->n);
930
931                                         if(tf) glTexCoord2fv(tf->uv[0]);
932                                         if(cp) glColor3ub(cp[3], cp[2], cp[1]);
933                                         glVertex3fv(efa->v1->co);
934
935                                         if(tf) glTexCoord2fv(tf->uv[1]);
936                                         if(cp) glColor3ub(cp[7], cp[6], cp[5]);
937                                         glVertex3fv(efa->v2->co);
938
939                                         if(tf) glTexCoord2fv(tf->uv[2]);
940                                         if(cp) glColor3ub(cp[11], cp[10], cp[9]);
941                                         glVertex3fv(efa->v3->co);
942
943                                         if(efa->v4) {
944                                                 if(tf) glTexCoord2fv(tf->uv[3]);
945                                                 if(cp) glColor3ub(cp[15], cp[14], cp[13]);
946                                                 glVertex3fv(efa->v4->co);
947                                         }
948                                 } else {
949                                         if(tf) glTexCoord2fv(tf->uv[0]);
950                                         if(cp) glColor3ub(cp[3], cp[2], cp[1]);
951                                         glNormal3fv(efa->v1->no);
952                                         glVertex3fv(efa->v1->co);
953
954                                         if(tf) glTexCoord2fv(tf->uv[1]);
955                                         if(cp) glColor3ub(cp[7], cp[6], cp[5]);
956                                         glNormal3fv(efa->v2->no);
957                                         glVertex3fv(efa->v2->co);
958
959                                         if(tf) glTexCoord2fv(tf->uv[2]);
960                                         if(cp) glColor3ub(cp[11], cp[10], cp[9]);
961                                         glNormal3fv(efa->v3->no);
962                                         glVertex3fv(efa->v3->co);
963
964                                         if(efa->v4) {
965                                                 if(tf) glTexCoord2fv(tf->uv[3]);
966                                                 if(cp) glColor3ub(cp[15], cp[14], cp[13]);
967                                                 glNormal3fv(efa->v4->no);
968                                                 glVertex3fv(efa->v4->co);
969                                         }
970                                 }
971                                 glEnd();
972                         }
973                 }
974         }
975 }
976
977 static void emDM_drawFacesTex(DerivedMesh *dm, int (*setDrawOptions)(MTFace *tface, int has_mcol, int matnr))
978 {
979         emDM_drawFacesTex_common(dm, setDrawOptions, NULL, NULL);
980 }
981
982 static void emDM_drawMappedFacesTex(DerivedMesh *dm, int (*setDrawOptions)(void *userData, int index), void *userData)
983 {
984         emDM_drawFacesTex_common(dm, NULL, setDrawOptions, userData);
985 }
986
987 static void emDM_drawMappedFacesGLSL(DerivedMesh *dm,
988                            int (*setMaterial)(int, void *attribs),
989                            int (*setDrawOptions)(void *userData, int index), void *userData) 
990 {
991         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
992         EditMesh *em= emdm->em;
993         float (*vertexCos)[3]= emdm->vertexCos;
994         float (*vertexNos)[3]= emdm->vertexNos;
995         EditVert *eve;
996         EditFace *efa;
997         DMVertexAttribs attribs= {{{0}}};
998         GPUVertexAttribs gattribs;
999         /* int tfoffset; */ /* UNUSED */
1000         int i, b, matnr, new_matnr, dodraw /* , layer */ /* UNUSED */;
1001
1002         dodraw = 0;
1003         matnr = -1;
1004
1005         /* layer = CustomData_get_layer_index(&em->fdata, CD_MTFACE); */ /* UNUSED */
1006         /* tfoffset = (layer == -1)? -1: em->fdata.layers[layer].offset; */ /* UNUSED */
1007
1008         /* always use smooth shading even for flat faces, else vertex colors wont interpolate */
1009         glShadeModel(GL_SMOOTH);
1010
1011         for (i=0,eve=em->verts.first; eve; eve= eve->next)
1012                 eve->tmp.l = (intptr_t) i++;
1013
1014 #define PASSATTRIB(efa, eve, vert) {                                                                                    \
1015         if(attribs.totorco) {                                                                                                           \
1016                 float *orco = attribs.orco.array[eve->tmp.l];                                                   \
1017                 glVertexAttrib3fvARB(attribs.orco.glIndex, orco);                                               \
1018         }                                                                                                                                                       \
1019         for(b = 0; b < attribs.tottface; b++) {                                                                         \
1020                 MTFace *_tf = (MTFace*)((char*)efa->data + attribs.tface[b].emOffset);  \
1021                 glVertexAttrib2fvARB(attribs.tface[b].glIndex, _tf->uv[vert]);                  \
1022         }                                                                                                                                                       \
1023         for(b = 0; b < attribs.totmcol; b++) {                                                                          \
1024                 MCol *cp = (MCol*)((char*)efa->data + attribs.mcol[b].emOffset);                \
1025                 GLubyte col[4];                                                                                                                 \
1026                 col[0]= cp->b; col[1]= cp->g; col[2]= cp->r; col[3]= cp->a;                             \
1027                 glVertexAttrib4ubvARB(attribs.mcol[b].glIndex, col);                                    \
1028         }                                                                                                                                                       \
1029         if(attribs.tottang) {                                                                                                           \
1030                 float *tang = attribs.tang.array[i*4 + vert];                                                   \
1031                 glVertexAttrib4fvARB(attribs.tang.glIndex, tang);                                               \
1032         }                                                                                                                                                       \
1033 }
1034
1035         for (i=0,efa= em->faces.first; efa; i++,efa= efa->next) {
1036                 int drawSmooth= (efa->flag & ME_SMOOTH);
1037
1038                 if(setDrawOptions && !setDrawOptions(userData, i))
1039                         continue;
1040
1041                 new_matnr = efa->mat_nr + 1;
1042                 if(new_matnr != matnr) {
1043                         dodraw = setMaterial(matnr = new_matnr, &gattribs);
1044                         if(dodraw)
1045                                 DM_vertex_attributes_from_gpu(dm, &gattribs, &attribs);
1046                 }
1047
1048                 if(dodraw) {
1049                         glBegin(efa->v4?GL_QUADS:GL_TRIANGLES);
1050                         if (!drawSmooth) {
1051                                 if(vertexCos) glNormal3fv(emdm->faceNos[i]);
1052                                 else glNormal3fv(efa->n);
1053
1054                                 PASSATTRIB(efa, efa->v1, 0);
1055                                 if(vertexCos) glVertex3fv(vertexCos[(int) efa->v1->tmp.l]);
1056                                 else glVertex3fv(efa->v1->co);
1057
1058                                 PASSATTRIB(efa, efa->v2, 1);
1059                                 if(vertexCos) glVertex3fv(vertexCos[(int) efa->v2->tmp.l]);
1060                                 else glVertex3fv(efa->v2->co);
1061
1062                                 PASSATTRIB(efa, efa->v3, 2);
1063                                 if(vertexCos) glVertex3fv(vertexCos[(int) efa->v3->tmp.l]);
1064                                 else glVertex3fv(efa->v3->co);
1065
1066                                 if(efa->v4) {
1067                                         PASSATTRIB(efa, efa->v4, 3);
1068                                         if(vertexCos) glVertex3fv(vertexCos[(int) efa->v4->tmp.l]);
1069                                         else glVertex3fv(efa->v4->co);
1070                                 }
1071                         } else {
1072                                 PASSATTRIB(efa, efa->v1, 0);
1073                                 if(vertexCos) {
1074                                         glNormal3fv(vertexNos[(int) efa->v1->tmp.l]);
1075                                         glVertex3fv(vertexCos[(int) efa->v1->tmp.l]);
1076                                 }
1077                                 else {
1078                                         glNormal3fv(efa->v1->no);
1079                                         glVertex3fv(efa->v1->co);
1080                                 }
1081
1082                                 PASSATTRIB(efa, efa->v2, 1);
1083                                 if(vertexCos) {
1084                                         glNormal3fv(vertexNos[(int) efa->v2->tmp.l]);
1085                                         glVertex3fv(vertexCos[(int) efa->v2->tmp.l]);
1086                                 }
1087                                 else {
1088                                         glNormal3fv(efa->v2->no);
1089                                         glVertex3fv(efa->v2->co);
1090                                 }
1091
1092                                 PASSATTRIB(efa, efa->v3, 2);
1093                                 if(vertexCos) {
1094                                         glNormal3fv(vertexNos[(int) efa->v3->tmp.l]);
1095                                         glVertex3fv(vertexCos[(int) efa->v3->tmp.l]);
1096                                 }
1097                                 else {
1098                                         glNormal3fv(efa->v3->no);
1099                                         glVertex3fv(efa->v3->co);
1100                                 }
1101
1102                                 if(efa->v4) {
1103                                         PASSATTRIB(efa, efa->v4, 3);
1104                                         if(vertexCos) {
1105                                                 glNormal3fv(vertexNos[(int) efa->v4->tmp.l]);
1106                                                 glVertex3fv(vertexCos[(int) efa->v4->tmp.l]);
1107                                         }
1108                                         else {
1109                                                 glNormal3fv(efa->v4->no);
1110                                                 glVertex3fv(efa->v4->co);
1111                                         }
1112                                 }
1113                         }
1114                         glEnd();
1115                 }
1116         }
1117 #undef PASSATTRIB
1118 }
1119
1120 static void emDM_drawFacesGLSL(DerivedMesh *dm,
1121                            int (*setMaterial)(int, void *attribs))
1122 {
1123         dm->drawMappedFacesGLSL(dm, setMaterial, NULL, NULL);
1124 }
1125
1126 static void emDM_getMinMax(DerivedMesh *dm, float min_r[3], float max_r[3])
1127 {
1128         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
1129         EditVert *eve;
1130         int i;
1131
1132         if (emdm->em->verts.first) {
1133                 for (i=0,eve= emdm->em->verts.first; eve; i++,eve= eve->next) {
1134                         if (emdm->vertexCos) {
1135                                 DO_MINMAX(emdm->vertexCos[i], min_r, max_r);
1136                         } else {
1137                                 DO_MINMAX(eve->co, min_r, max_r);
1138                         }
1139                 }
1140         } else {
1141                 min_r[0] = min_r[1] = min_r[2] = max_r[0] = max_r[1] = max_r[2] = 0.0;
1142         }
1143 }
1144 static int emDM_getNumVerts(DerivedMesh *dm)
1145 {
1146         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
1147
1148         return BLI_countlist(&emdm->em->verts);
1149 }
1150
1151 static int emDM_getNumEdges(DerivedMesh *dm)
1152 {
1153         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
1154
1155         return BLI_countlist(&emdm->em->edges);
1156 }
1157
1158 static int emDM_getNumFaces(DerivedMesh *dm)
1159 {
1160         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
1161
1162         return BLI_countlist(&emdm->em->faces);
1163 }
1164
1165 static void emDM_getVertCos(DerivedMesh *dm, float (*cos_r)[3])
1166 {
1167         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
1168         EditVert *eve;
1169         int i;
1170
1171         for (i=0,eve= emdm->em->verts.first; eve; i++,eve=eve->next) {
1172                 if (emdm->vertexCos) {
1173                         copy_v3_v3(cos_r[i], emdm->vertexCos[i]);
1174                 } else {
1175                         copy_v3_v3(cos_r[i], eve->co);
1176                 }
1177         }
1178 }
1179
1180 static void emDM_getVert(DerivedMesh *dm, int index, MVert *vert_r)
1181 {
1182         EditVert *ev = ((EditMeshDerivedMesh *)dm)->em->verts.first;
1183         int i;
1184
1185         for(i = 0; i < index; ++i) ev = ev->next;
1186
1187         VECCOPY(vert_r->co, ev->co);
1188
1189         normal_float_to_short_v3(vert_r->no, ev->no);
1190
1191         /* TODO what to do with vert_r->flag? */
1192         vert_r->bweight = (unsigned char) (ev->bweight*255.0f);
1193 }
1194
1195 static void emDM_getEdge(DerivedMesh *dm, int index, MEdge *edge_r)
1196 {
1197         EditMesh *em = ((EditMeshDerivedMesh *)dm)->em;
1198         EditEdge *ee = em->edges.first;
1199         EditVert *ev, *v1, *v2;
1200         int i;
1201
1202         for(i = 0; i < index; ++i) ee = ee->next;
1203
1204         edge_r->crease = (unsigned char) (ee->crease*255.0f);
1205         edge_r->bweight = (unsigned char) (ee->bweight*255.0f);
1206         /* TODO what to do with edge_r->flag? */
1207         edge_r->flag = ME_EDGEDRAW|ME_EDGERENDER;
1208         if (ee->seam) edge_r->flag |= ME_SEAM;
1209         if (ee->sharp) edge_r->flag |= ME_SHARP;
1210 #if 0
1211         /* this needs setup of f2 field */
1212         if (!ee->f2) edge_r->flag |= ME_LOOSEEDGE;
1213 #endif
1214
1215         /* goddamn, we have to search all verts to find indices */
1216         v1 = ee->v1;
1217         v2 = ee->v2;
1218         for(i = 0, ev = em->verts.first; v1 || v2; i++, ev = ev->next) {
1219                 if(ev == v1) {
1220                         edge_r->v1 = i;
1221                         v1 = NULL;
1222                 }
1223                 if(ev == v2) {
1224                         edge_r->v2 = i;
1225                         v2 = NULL;
1226                 }
1227         }
1228 }
1229
1230 static void emDM_getFace(DerivedMesh *dm, int index, MFace *face_r)
1231 {
1232         EditMesh *em = ((EditMeshDerivedMesh *)dm)->em;
1233         EditFace *ef = em->faces.first;
1234         EditVert *ev, *v1, *v2, *v3, *v4;
1235         int i;
1236
1237         for(i = 0; i < index; ++i) ef = ef->next;
1238
1239         face_r->mat_nr = ef->mat_nr;
1240         face_r->flag = ef->flag;
1241
1242         /* goddamn, we have to search all verts to find indices */
1243         v1 = ef->v1;
1244         v2 = ef->v2;
1245         v3 = ef->v3;
1246         v4 = ef->v4;
1247         if(!v4) face_r->v4 = 0;
1248
1249         for(i = 0, ev = em->verts.first; v1 || v2 || v3 || v4;
1250                 i++, ev = ev->next) {
1251                 if(ev == v1) {
1252                         face_r->v1 = i;
1253                         v1 = NULL;
1254                 }
1255                 if(ev == v2) {
1256                         face_r->v2 = i;
1257                         v2 = NULL;
1258                 }
1259                 if(ev == v3) {
1260                         face_r->v3 = i;
1261                         v3 = NULL;
1262                 }
1263                 if(ev == v4) {
1264                         face_r->v4 = i;
1265                         v4 = NULL;
1266                 }
1267         }
1268
1269         test_index_face(face_r, NULL, 0, ef->v4?4:3);
1270 }
1271
1272 static void emDM_copyVertArray(DerivedMesh *dm, MVert *vert_r)
1273 {
1274         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
1275         EditVert *ev = emdm->em->verts.first;
1276         int i;
1277
1278         for(i=0; ev; ev = ev->next, ++vert_r, ++i) {
1279                 if(emdm->vertexCos)
1280                         copy_v3_v3(vert_r->co, emdm->vertexCos[i]);
1281                 else
1282                         copy_v3_v3(vert_r->co, ev->co);
1283
1284                 normal_float_to_short_v3(vert_r->no, ev->no);
1285
1286                 /* TODO what to do with vert_r->flag? */
1287                 vert_r->flag = 0;
1288                 vert_r->bweight = (unsigned char) (ev->bweight*255.0f);
1289         }
1290 }
1291
1292 static void emDM_copyEdgeArray(DerivedMesh *dm, MEdge *edge_r)
1293 {
1294         EditMesh *em = ((EditMeshDerivedMesh *)dm)->em;
1295         EditEdge *ee = em->edges.first;
1296         EditVert *ev;
1297         int i;
1298
1299         /* store vertex indices in tmp union */
1300         for(ev = em->verts.first, i = 0; ev; ev = ev->next, ++i)
1301                 ev->tmp.l = (intptr_t) i;
1302
1303         for( ; ee; ee = ee->next, ++edge_r) {
1304                 edge_r->crease = (unsigned char) (ee->crease*255.0f);
1305                 edge_r->bweight = (unsigned char) (ee->bweight*255.0f);
1306                 /* TODO what to do with edge_r->flag? */
1307                 edge_r->flag = ME_EDGEDRAW|ME_EDGERENDER;
1308                 if (ee->seam) edge_r->flag |= ME_SEAM;
1309                 if (ee->sharp) edge_r->flag |= ME_SHARP;
1310 #if 0
1311                 /* this needs setup of f2 field */
1312                 if (!ee->f2) edge_r->flag |= ME_LOOSEEDGE;
1313 #endif
1314
1315                 edge_r->v1 = (int)ee->v1->tmp.l;
1316                 edge_r->v2 = (int)ee->v2->tmp.l;
1317         }
1318 }
1319
1320 static void emDM_copyFaceArray(DerivedMesh *dm, MFace *face_r)
1321 {
1322         EditMesh *em = ((EditMeshDerivedMesh *)dm)->em;
1323         EditFace *ef = em->faces.first;
1324         EditVert *ev;
1325         int i;
1326
1327         /* store vertexes indices in tmp union */
1328         for(ev = em->verts.first, i = 0; ev; ev = ev->next, ++i)
1329                 ev->tmp.l = (intptr_t) i;
1330
1331         for( ; ef; ef = ef->next, ++face_r) {
1332                 face_r->mat_nr = ef->mat_nr;
1333                 face_r->flag = ef->flag;
1334
1335                 face_r->v1 = (int)ef->v1->tmp.l;
1336                 face_r->v2 = (int)ef->v2->tmp.l;
1337                 face_r->v3 = (int)ef->v3->tmp.l;
1338                 if(ef->v4) face_r->v4 = (int)ef->v4->tmp.l;
1339                 else face_r->v4 = 0;
1340
1341                 test_index_face(face_r, NULL, 0, ef->v4?4:3);
1342         }
1343 }
1344
1345 static void *emDM_getFaceDataArray(DerivedMesh *dm, int type)
1346 {
1347         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
1348         EditMesh *em= emdm->em;
1349         EditFace *efa;
1350         char *data, *emdata;
1351         void *datalayer;
1352         int index, size;
1353
1354         datalayer = DM_get_face_data_layer(dm, type);
1355         if(datalayer)
1356                 return datalayer;
1357
1358         /* layers are store per face for editmesh, we convert to a temporary
1359          * data layer array in the derivedmesh when these are requested */
1360         if(type == CD_MTFACE || type == CD_MCOL) {
1361                 index = CustomData_get_layer_index(&em->fdata, type);
1362
1363                 if(index != -1) {
1364                         /* int offset = em->fdata.layers[index].offset; */ /* UNUSED */
1365                         size = CustomData_sizeof(type);
1366
1367                         DM_add_face_layer(dm, type, CD_CALLOC, NULL);
1368                         index = CustomData_get_layer_index(&dm->faceData, type);
1369                         dm->faceData.layers[index].flag |= CD_FLAG_TEMPORARY;
1370
1371                         data = datalayer = DM_get_face_data_layer(dm, type);
1372                         for(efa=em->faces.first; efa; efa=efa->next, data+=size) {
1373                                 emdata = CustomData_em_get(&em->fdata, efa->data, type);
1374                                 memcpy(data, emdata, size);
1375                         }
1376                 }
1377         }
1378
1379         return datalayer;
1380 }
1381
1382 static void emDM_release(DerivedMesh *dm)
1383 {
1384         EditMeshDerivedMesh *emdm= (EditMeshDerivedMesh*) dm;
1385
1386         if (DM_release(dm)) {
1387                 if (emdm->vertexCos) {
1388                         MEM_freeN(emdm->vertexCos);
1389                         MEM_freeN(emdm->vertexNos);
1390                         MEM_freeN(emdm->faceNos);
1391                 }
1392
1393                 MEM_freeN(emdm);
1394         }
1395 }
1396
1397 DerivedMesh *editmesh_get_derived(EditMesh *em, float (*vertexCos)[3])
1398 {
1399         EditMeshDerivedMesh *emdm = MEM_callocN(sizeof(*emdm), "emdm");
1400
1401         DM_init(&emdm->dm, DM_TYPE_EDITMESH, BLI_countlist(&em->verts),
1402                                          BLI_countlist(&em->edges), BLI_countlist(&em->faces));
1403
1404         emdm->dm.getMinMax = emDM_getMinMax;
1405
1406         emdm->dm.getNumVerts = emDM_getNumVerts;
1407         emdm->dm.getNumEdges = emDM_getNumEdges;
1408         emdm->dm.getNumFaces = emDM_getNumFaces;
1409
1410         emdm->dm.getVertCos = emDM_getVertCos;
1411
1412         emdm->dm.getVert = emDM_getVert;
1413         emdm->dm.getEdge = emDM_getEdge;
1414         emdm->dm.getFace = emDM_getFace;
1415         emdm->dm.copyVertArray = emDM_copyVertArray;
1416         emdm->dm.copyEdgeArray = emDM_copyEdgeArray;
1417         emdm->dm.copyFaceArray = emDM_copyFaceArray;
1418         emdm->dm.getFaceDataArray = emDM_getFaceDataArray;
1419
1420         emdm->dm.foreachMappedVert = emDM_foreachMappedVert;
1421         emdm->dm.foreachMappedEdge = emDM_foreachMappedEdge;
1422         emdm->dm.foreachMappedFaceCenter = emDM_foreachMappedFaceCenter;
1423
1424         emdm->dm.drawEdges = emDM_drawEdges;
1425         emdm->dm.drawMappedEdges = emDM_drawMappedEdges;
1426         emdm->dm.drawMappedEdgesInterp = emDM_drawMappedEdgesInterp;
1427         emdm->dm.drawMappedFaces = emDM_drawMappedFaces;
1428         emdm->dm.drawMappedFacesTex = emDM_drawMappedFacesTex;
1429         emdm->dm.drawMappedFacesGLSL = emDM_drawMappedFacesGLSL;
1430         emdm->dm.drawFacesTex = emDM_drawFacesTex;
1431         emdm->dm.drawFacesGLSL = emDM_drawFacesGLSL;
1432         emdm->dm.drawUVEdges = emDM_drawUVEdges;
1433
1434         emdm->dm.release = emDM_release;
1435         
1436         emdm->em = em;
1437         emdm->vertexCos = vertexCos;
1438
1439         if(CustomData_has_layer(&em->vdata, CD_MDEFORMVERT)) {
1440                 EditVert *eve;
1441                 int i;
1442
1443                 DM_add_vert_layer(&emdm->dm, CD_MDEFORMVERT, CD_CALLOC, NULL);
1444
1445                 for(eve = em->verts.first, i = 0; eve; eve = eve->next, ++i)
1446                         DM_set_vert_data(&emdm->dm, i, CD_MDEFORMVERT,
1447                                                          CustomData_em_get(&em->vdata, eve->data, CD_MDEFORMVERT));
1448         }
1449
1450         if(vertexCos) {
1451                 EditVert *eve;
1452                 EditFace *efa;
1453                 int totface = BLI_countlist(&em->faces);
1454                 int i;
1455
1456                 for (i=0,eve=em->verts.first; eve; eve= eve->next)
1457                         eve->tmp.l = (intptr_t) i++;
1458
1459                 emdm->vertexNos = MEM_callocN(sizeof(*emdm->vertexNos)*i, "emdm_vno");
1460                 emdm->faceNos = MEM_mallocN(sizeof(*emdm->faceNos)*totface, "emdm_vno");
1461
1462                 for(i=0, efa= em->faces.first; efa; i++, efa=efa->next) {
1463                         float *v1 = vertexCos[(int) efa->v1->tmp.l];
1464                         float *v2 = vertexCos[(int) efa->v2->tmp.l];
1465                         float *v3 = vertexCos[(int) efa->v3->tmp.l];
1466                         float *no = emdm->faceNos[i];
1467                         
1468                         if(efa->v4) {
1469                                 float *v4 = vertexCos[(int) efa->v4->tmp.l];
1470
1471                                 normal_quad_v3( no,v1, v2, v3, v4);
1472                                 add_v3_v3(emdm->vertexNos[(int) efa->v4->tmp.l], no);
1473                         }
1474                         else {
1475                                 normal_tri_v3( no,v1, v2, v3);
1476                         }
1477
1478                         add_v3_v3(emdm->vertexNos[(int) efa->v1->tmp.l], no);
1479                         add_v3_v3(emdm->vertexNos[(int) efa->v2->tmp.l], no);
1480                         add_v3_v3(emdm->vertexNos[(int) efa->v3->tmp.l], no);
1481                 }
1482
1483                 for(i=0, eve= em->verts.first; eve; i++, eve=eve->next) {
1484                         float *no = emdm->vertexNos[i];
1485                         /* following Mesh convention; we use vertex coordinate itself
1486                          * for normal in this case */
1487                         if (normalize_v3(no) == 0.0f) {
1488                                 normalize_v3_v3(no, vertexCos[i]);
1489                         }
1490                 }
1491         }
1492
1493         return (DerivedMesh*) emdm;
1494 }
1495
1496 /***/
1497
1498 DerivedMesh *mesh_create_derived_for_modifier(Scene *scene, Object *ob, ModifierData *md)
1499 {
1500         Mesh *me = ob->data;
1501         ModifierTypeInfo *mti = modifierType_getInfo(md->type);
1502         DerivedMesh *dm;
1503
1504         md->scene= scene;
1505         
1506         if (!(md->mode&eModifierMode_Realtime)) return NULL;
1507         if (mti->isDisabled && mti->isDisabled(md, 0)) return NULL;
1508
1509         if (mti->type==eModifierTypeType_OnlyDeform) {
1510                 int numVerts;
1511                 float (*deformedVerts)[3] = mesh_getVertexCos(me, &numVerts);
1512
1513                 mti->deformVerts(md, ob, NULL, deformedVerts, numVerts, 0, 0);
1514                 dm = mesh_create_derived(me, ob, deformedVerts);
1515
1516                 MEM_freeN(deformedVerts);
1517         } else {
1518                 DerivedMesh *tdm = mesh_create_derived(me, ob, NULL);
1519                 dm = mti->applyModifier(md, ob, tdm, 0, 0);
1520
1521                 if(tdm != dm) tdm->release(tdm);
1522         }
1523
1524         return dm;
1525 }
1526
1527 static float *get_editmesh_orco_verts(EditMesh *em)
1528 {
1529         EditVert *eve;
1530         float *orco;
1531         int a, totvert;
1532
1533         /* these may not really be the orco's, but it's only for preview.
1534          * could be solver better once, but isn't simple */
1535
1536         totvert= 0;
1537         for(eve=em->verts.first; eve; eve=eve->next)
1538                 totvert++;
1539         
1540         orco = MEM_mallocN(sizeof(float)*3*totvert, "EditMesh Orco");
1541
1542         for(a=0, eve=em->verts.first; eve; eve=eve->next, a+=3)
1543                 VECCOPY(orco+a, eve->co);
1544         
1545         return orco;
1546 }
1547
1548 /* orco custom data layer */
1549
1550 static void *get_orco_coords_dm(Object *ob, EditMesh *em, int layer, int *free)
1551 {
1552         *free= 0;
1553
1554         if(layer == CD_ORCO) {
1555                 /* get original coordinates */
1556                 *free= 1;
1557
1558                 if(em)
1559                         return (float(*)[3])get_editmesh_orco_verts(em);
1560                 else
1561                         return (float(*)[3])get_mesh_orco_verts(ob);
1562         }
1563         else if(layer == CD_CLOTH_ORCO) {
1564                 /* apply shape key for cloth, this should really be solved
1565                    by a more flexible customdata system, but not simple */
1566                 if(!em) {
1567                         ClothModifierData *clmd = (ClothModifierData *)modifiers_findByType(ob, eModifierType_Cloth);
1568                         KeyBlock *kb= key_get_keyblock(ob_get_key(ob), clmd->sim_parms->shapekey_rest);
1569
1570                         if(kb->data)
1571                                 return kb->data;
1572                 }
1573
1574                 return NULL;
1575         }
1576
1577         return NULL;
1578 }
1579
1580 static DerivedMesh *create_orco_dm(Object *ob, Mesh *me, EditMesh *em, int layer)
1581 {
1582         DerivedMesh *dm;
1583         float (*orco)[3];
1584         int free;
1585
1586         if(em) dm= CDDM_from_editmesh(em, me);
1587         else dm= CDDM_from_mesh(me, ob);
1588
1589         orco= get_orco_coords_dm(ob, em, layer, &free);
1590
1591         if(orco) {
1592                 CDDM_apply_vert_coords(dm, orco);
1593                 if(free) MEM_freeN(orco);
1594         }
1595
1596         CDDM_calc_normals(dm);
1597
1598         return dm;
1599 }
1600
1601 static void add_orco_dm(Object *ob, EditMesh *em, DerivedMesh *dm, DerivedMesh *orcodm, int layer)
1602 {
1603         float (*orco)[3], (*layerorco)[3];
1604         int totvert, free;
1605
1606         totvert= dm->getNumVerts(dm);
1607
1608         if(orcodm) {
1609                 orco= MEM_callocN(sizeof(float)*3*totvert, "dm orco");
1610                 free= 1;
1611
1612                 if(orcodm->getNumVerts(orcodm) == totvert)
1613                         orcodm->getVertCos(orcodm, orco);
1614                 else
1615                         dm->getVertCos(dm, orco);
1616         }
1617         else
1618                 orco= get_orco_coords_dm(ob, em, layer, &free);
1619
1620         if(orco) {
1621                 if(layer == CD_ORCO)
1622                         transform_mesh_orco_verts(ob->data, orco, totvert, 0);
1623
1624                 if(!(layerorco = DM_get_vert_data_layer(dm, layer))) {
1625                         DM_add_vert_layer(dm, layer, CD_CALLOC, NULL);
1626                         layerorco = DM_get_vert_data_layer(dm, layer);
1627                 }
1628
1629                 memcpy(layerorco, orco, sizeof(float)*3*totvert);
1630                 if(free) MEM_freeN(orco);
1631         }
1632 }
1633
1634 /* weight paint colors */
1635
1636 /* Something of a hack, at the moment deal with weightpaint
1637  * by tucking into colors during modifier eval, only in
1638  * wpaint mode. Works ok but need to make sure recalc
1639  * happens on enter/exit wpaint.
1640  */
1641
1642 void weight_to_rgb(float input, float *fr, float *fg, float *fb)
1643 {
1644         float blend;
1645         
1646         blend= ((input/2.0f)+0.5f);
1647         
1648         if (input<=0.25f){      // blue->cyan
1649                 *fr= 0.0f;
1650                 *fg= blend*input*4.0f;
1651                 *fb= blend;
1652         }
1653         else if (input<=0.50f){ // cyan->green
1654                 *fr= 0.0f;
1655                 *fg= blend;
1656                 *fb= blend*(1.0f-((input-0.25f)*4.0f)); 
1657         }
1658         else if (input <= 0.75f){       // green->yellow
1659                 *fr= blend * ((input-0.50f)*4.0f);
1660                 *fg= blend;
1661                 *fb= 0.0f;
1662         }
1663         else if (input <= 1.0f){ // yellow->red
1664                 *fr= blend;
1665                 *fg= blend * (1.0f-((input-0.75f)*4.0f)); 
1666                 *fb= 0.0f;
1667         }
1668 }
1669
1670 /* draw_flag's for calc_weightpaint_vert_color */
1671 enum {
1672         CALC_WP_MULTIPAINT= (1<<0),
1673         CALC_WP_AUTO_NORMALIZE= (1<<1),
1674 };
1675
1676 static void calc_weightpaint_vert_color(Object *ob, ColorBand *coba, int vert, unsigned char *col, char *dg_flags, int selected, int UNUSED(unselected), const int draw_flag)
1677 {
1678         Mesh *me = ob->data;
1679         float colf[4], input = 0.0f;
1680         int i;
1681
1682         
1683         int make_black= FALSE;
1684
1685         if (me->dvert) {
1686                 if ((selected > 1) && (draw_flag & CALC_WP_MULTIPAINT)) {
1687                         
1688                         int was_a_nonzero= FALSE;
1689                         for (i=0; i<me->dvert[vert].totweight; i++) {
1690                                 /* in multipaint, get the average if auto normalize is inactive
1691                                  * get the sum if it is active */
1692                                 if(dg_flags[me->dvert[vert].dw[i].def_nr]) {
1693                                         if(me->dvert[vert].dw[i].weight) {
1694                                                 input+= me->dvert[vert].dw[i].weight;
1695                                                 was_a_nonzero= TRUE;
1696                                         }
1697                                 }
1698                         }
1699
1700                         /* make it black if the selected groups have no weight on a vertex */
1701                         if(was_a_nonzero == FALSE) {
1702                                 make_black = TRUE;
1703                         }
1704                         else if ((draw_flag & CALC_WP_AUTO_NORMALIZE) == FALSE) {
1705                                 input /= selected; /* get the average */
1706                         }
1707                 }
1708                 else {
1709                         /* default, non tricky behavior */
1710                         for (i=0; i<me->dvert[vert].totweight; i++) {
1711                                 if (me->dvert[vert].dw[i].def_nr==ob->actdef-1) {
1712                                         input+=me->dvert[vert].dw[i].weight;
1713                                 }
1714                         }
1715                 }
1716         }
1717         
1718         if (make_black) {
1719                 col[3] = 0;
1720                 col[2] = 0;
1721                 col[1] = 0;
1722                 col[0] = 255;
1723                 return;
1724         }
1725
1726         CLAMP(input, 0.0f, 1.0f);       
1727
1728         if(coba)
1729                 do_colorband(coba, input, colf);
1730         else
1731                 weight_to_rgb(input, colf, colf+1, colf+2);
1732         
1733         col[3] = (unsigned char)(colf[0] * 255.0f);
1734         col[2] = (unsigned char)(colf[1] * 255.0f);
1735         col[1] = (unsigned char)(colf[2] * 255.0f);
1736         col[0] = 255;
1737 }
1738
1739 static ColorBand *stored_cb= NULL;
1740
1741 void vDM_ColorBand_store(ColorBand *coba)
1742 {
1743         stored_cb= coba;
1744 }
1745
1746 static void add_weight_mcol_dm(Object *ob, DerivedMesh *dm, int const draw_flag)
1747 {
1748         Mesh *me = ob->data;
1749         MFace *mf = me->mface;
1750         ColorBand *coba= stored_cb;     /* warning, not a local var */
1751         unsigned char *wtcol;
1752         int i;
1753         
1754         int defbase_len = BLI_countlist(&ob->defbase);
1755         char *defbase_sel = MEM_mallocN(defbase_len * sizeof(char), __func__);
1756         int selected = get_selected_defgroups(ob, defbase_sel, defbase_len);
1757         int unselected = defbase_len - selected;
1758
1759         wtcol = MEM_callocN (sizeof (unsigned char) * me->totface*4*4, "weightmap");
1760         
1761         memset(wtcol, 0x55, sizeof (unsigned char) * me->totface*4*4);
1762         for (i=0; i<me->totface; i++, mf++) {
1763                 calc_weightpaint_vert_color(ob, coba, mf->v1, &wtcol[(i*4 + 0)*4], defbase_sel, selected, unselected, draw_flag);
1764                 calc_weightpaint_vert_color(ob, coba, mf->v2, &wtcol[(i*4 + 1)*4], defbase_sel, selected, unselected, draw_flag);
1765                 calc_weightpaint_vert_color(ob, coba, mf->v3, &wtcol[(i*4 + 2)*4], defbase_sel, selected, unselected, draw_flag);
1766                 if (mf->v4)
1767                         calc_weightpaint_vert_color(ob, coba, mf->v4, &wtcol[(i*4 + 3)*4], defbase_sel, selected, unselected, draw_flag);
1768         }
1769         
1770         MEM_freeN(defbase_sel);
1771
1772         CustomData_add_layer(&dm->faceData, CD_WEIGHT_MCOL, CD_ASSIGN, wtcol, dm->numFaceData);
1773 }
1774
1775 /* new value for useDeform -1  (hack for the gameengine):
1776  * - apply only the modifier stack of the object, skipping the virtual modifiers,
1777  * - don't apply the key
1778  * - apply deform modifiers and input vertexco
1779  */
1780 static void mesh_calc_modifiers(Scene *scene, Object *ob, float (*inputVertexCos)[3],
1781                                                                 DerivedMesh **deform_r, DerivedMesh **final_r,
1782                                                                 int useRenderParams, int useDeform,
1783                                                                 int needMapping, CustomDataMask dataMask, int index, int useCache)
1784 {
1785         Mesh *me = ob->data;
1786         ModifierData *firstmd, *md;
1787         LinkNode *datamasks, *curr;
1788         CustomDataMask mask, nextmask;
1789         float (*deformedVerts)[3] = NULL;
1790         DerivedMesh *dm, *orcodm, *clothorcodm, *finaldm;
1791         int numVerts = me->totvert;
1792         int required_mode;
1793         int isPrevDeform= FALSE;
1794         int skipVirtualArmature = (useDeform < 0);
1795         MultiresModifierData *mmd= get_multires_modifier(scene, ob, 0);
1796         int has_multires = mmd != NULL, multires_applied = 0;
1797         int sculpt_mode = ob->mode & OB_MODE_SCULPT && ob->sculpt;
1798
1799         int draw_flag= ((scene->toolsettings->multipaint ? CALC_WP_MULTIPAINT : 0) |
1800                         (scene->toolsettings->auto_normalize ? CALC_WP_AUTO_NORMALIZE : 0));
1801
1802         if(mmd && !mmd->sculptlvl)
1803                 has_multires = 0;
1804
1805         if(!skipVirtualArmature) {
1806                 firstmd = modifiers_getVirtualModifierList(ob);
1807         }
1808         else {
1809                 /* game engine exception */
1810                 firstmd = ob->modifiers.first;
1811                 if(firstmd && firstmd->type == eModifierType_Armature)
1812                         firstmd = firstmd->next;
1813         }
1814
1815         md = firstmd;
1816
1817         modifiers_clearErrors(ob);
1818
1819         if(useRenderParams) required_mode = eModifierMode_Render;
1820         else required_mode = eModifierMode_Realtime;
1821
1822         datamasks = modifiers_calcDataMasks(scene, ob, md, dataMask, required_mode);
1823         curr = datamasks;
1824
1825         if(deform_r) *deform_r = NULL;
1826         *final_r = NULL;
1827
1828         if(useDeform) {
1829                 if(inputVertexCos)
1830                         deformedVerts = inputVertexCos;
1831                 
1832                 /* Apply all leading deforming modifiers */
1833                 for(;md; md = md->next, curr = curr->next) {
1834                         ModifierTypeInfo *mti = modifierType_getInfo(md->type);
1835
1836                         md->scene= scene;
1837                         
1838                         if(!modifier_isEnabled(scene, md, required_mode)) continue;
1839                         if(useDeform < 0 && mti->dependsOnTime && mti->dependsOnTime(md)) continue;
1840
1841                         if(mti->type == eModifierTypeType_OnlyDeform) {
1842                                 if(!deformedVerts)
1843                                         deformedVerts = mesh_getVertexCos(me, &numVerts);
1844
1845                                 mti->deformVerts(md, ob, NULL, deformedVerts, numVerts, useRenderParams, useDeform);
1846                         } else {
1847                                 break;
1848                         }
1849                         
1850                         /* grab modifiers until index i */
1851                         if((index >= 0) && (modifiers_indexInObject(ob, md) >= index))
1852                                 break;
1853                 }
1854
1855                 /* Result of all leading deforming modifiers is cached for
1856                  * places that wish to use the original mesh but with deformed
1857                  * coordinates (vpaint, etc.)
1858                  */
1859                 if (deform_r) {
1860                         *deform_r = CDDM_from_mesh(me, ob);
1861
1862                         if(deformedVerts) {
1863                                 CDDM_apply_vert_coords(*deform_r, deformedVerts);
1864                                 CDDM_calc_normals(*deform_r);
1865                         }
1866                 }
1867         } else {
1868                 /* default behaviour for meshes */
1869                 if(inputVertexCos)
1870                         deformedVerts = inputVertexCos;
1871                 else
1872                         deformedVerts = mesh_getVertexCos(me, &numVerts);
1873         }
1874
1875
1876         /* Now apply all remaining modifiers. If useDeform is off then skip
1877          * OnlyDeform ones. 
1878          */
1879         dm = NULL;
1880         orcodm = NULL;
1881         clothorcodm = NULL;
1882
1883         for(;md; md = md->next, curr = curr->next) {
1884                 ModifierTypeInfo *mti = modifierType_getInfo(md->type);
1885
1886                 md->scene= scene;
1887
1888                 if(!modifier_isEnabled(scene, md, required_mode)) continue;
1889                 if(mti->type == eModifierTypeType_OnlyDeform && !useDeform) continue;
1890                 if((mti->flags & eModifierTypeFlag_RequiresOriginalData) && dm) {
1891                         modifier_setError(md, "Modifier requires original data, bad stack position.");
1892                         continue;
1893                 }
1894                 if(sculpt_mode && (!has_multires || multires_applied)) {
1895                         int unsupported= 0;
1896
1897                         if(scene->toolsettings->sculpt->flags & SCULPT_ONLY_DEFORM)
1898                                 unsupported|= mti->type != eModifierTypeType_OnlyDeform;
1899
1900                         unsupported|= md->type == eModifierType_Multires && ((MultiresModifierData*)md)->sculptlvl==0;
1901                         unsupported|= multires_applied;
1902
1903                         if(unsupported) {
1904                                 modifier_setError(md, "Not supported in sculpt mode.");
1905                                 continue;
1906                         }
1907                 }
1908                 if(needMapping && !modifier_supportsMapping(md)) continue;
1909                 if(useDeform < 0 && mti->dependsOnTime && mti->dependsOnTime(md)) continue;
1910
1911                 /* add an orco layer if needed by this modifier */
1912                 if(mti->requiredDataMask)
1913                         mask = mti->requiredDataMask(ob, md);
1914                 else
1915                         mask = 0;
1916
1917                 if(dm && (mask & CD_MASK_ORCO))
1918                         add_orco_dm(ob, NULL, dm, orcodm, CD_ORCO);
1919
1920                 /* How to apply modifier depends on (a) what we already have as
1921                  * a result of previous modifiers (could be a DerivedMesh or just
1922                  * deformed vertices) and (b) what type the modifier is.
1923                  */
1924
1925                 if(mti->type == eModifierTypeType_OnlyDeform) {
1926                         /* No existing verts to deform, need to build them. */
1927                         if(!deformedVerts) {
1928                                 if(dm) {
1929                                         /* Deforming a derived mesh, read the vertex locations
1930                                          * out of the mesh and deform them. Once done with this
1931                                          * run of deformers verts will be written back.
1932                                          */
1933                                         numVerts = dm->getNumVerts(dm);
1934                                         deformedVerts =
1935                                                 MEM_mallocN(sizeof(*deformedVerts) * numVerts, "dfmv");
1936                                         dm->getVertCos(dm, deformedVerts);
1937                                 } else {
1938                                         deformedVerts = mesh_getVertexCos(me, &numVerts);
1939                                 }
1940                         }
1941
1942                         /* if this is not the last modifier in the stack then recalculate the normals
1943                          * to avoid giving bogus normals to the next modifier see: [#23673] */
1944                         if(isPrevDeform &&  mti->dependsOnNormals && mti->dependsOnNormals(md)) {
1945                                 /* XXX, this covers bug #23673, but we may need normal calc for other types */
1946                                 if(dm && dm->type == DM_TYPE_CDDM) {
1947                                         CDDM_apply_vert_coords(dm, deformedVerts);
1948                                         CDDM_calc_normals(dm);
1949                                 }
1950                         }
1951
1952                         mti->deformVerts(md, ob, dm, deformedVerts, numVerts, useRenderParams, useDeform);
1953                 } else {
1954                         DerivedMesh *ndm;
1955
1956                         /* determine which data layers are needed by following modifiers */
1957                         if(curr->next)
1958                                 nextmask= (CustomDataMask)GET_INT_FROM_POINTER(curr->next->link);
1959                         else
1960                                 nextmask= dataMask;
1961
1962                         /* apply vertex coordinates or build a DerivedMesh as necessary */
1963                         if(dm) {
1964                                 if(deformedVerts) {
1965                                         DerivedMesh *tdm = CDDM_copy(dm);
1966                                         dm->release(dm);
1967                                         dm = tdm;
1968
1969                                         CDDM_apply_vert_coords(dm, deformedVerts);
1970                                         CDDM_calc_normals(dm);
1971                                 }
1972                         } else {
1973                                 dm = CDDM_from_mesh(me, ob);
1974
1975                                 if(deformedVerts) {
1976                                         CDDM_apply_vert_coords(dm, deformedVerts);
1977                                         CDDM_calc_normals(dm);
1978                                 }
1979
1980                                 if((dataMask & CD_MASK_WEIGHT_MCOL) && (ob->mode & OB_MODE_WEIGHT_PAINT))
1981                                         add_weight_mcol_dm(ob, dm, draw_flag); 
1982
1983                                 /* Constructive modifiers need to have an origindex
1984                                  * otherwise they wont have anywhere to copy the data from.
1985                                  *
1986                                  * Also create ORIGINDEX data if any of the following modifiers
1987                                  * requests it, this way Mirror, Solidify etc will keep ORIGINDEX
1988                                  * data by using generic DM_copy_vert_data() functions.
1989                                  */
1990                                 if(needMapping || (nextmask & CD_MASK_ORIGINDEX)) {
1991                                         /* calc */
1992                                         DM_add_vert_layer(dm, CD_ORIGINDEX, CD_CALLOC, NULL);
1993                                         DM_add_edge_layer(dm, CD_ORIGINDEX, CD_CALLOC, NULL);
1994                                         DM_add_face_layer(dm, CD_ORIGINDEX, CD_CALLOC, NULL);
1995
1996                                         range_vni(DM_get_vert_data_layer(dm, CD_ORIGINDEX), dm->numVertData, 0);
1997                                         range_vni(DM_get_edge_data_layer(dm, CD_ORIGINDEX), dm->numEdgeData, 0);
1998                                         range_vni(DM_get_face_data_layer(dm, CD_ORIGINDEX), dm->numFaceData, 0);
1999                                 }
2000                         }
2001
2002                         
2003                         /* set the DerivedMesh to only copy needed data */
2004                         mask= (CustomDataMask)GET_INT_FROM_POINTER(curr->link);
2005                         /* needMapping check here fixes bug [#28112], otherwise its
2006                          * possible that it wont be copied */
2007                         DM_set_only_copy(dm, mask | (needMapping ? CD_MASK_ORIGINDEX : 0));
2008                         
2009                         /* add cloth rest shape key if need */
2010                         if(mask & CD_MASK_CLOTH_ORCO)
2011                                 add_orco_dm(ob, NULL, dm, clothorcodm, CD_CLOTH_ORCO);
2012
2013                         /* add an origspace layer if needed */
2014                         if(((CustomDataMask)GET_INT_FROM_POINTER(curr->link)) & CD_MASK_ORIGSPACE)
2015                                 if(!CustomData_has_layer(&dm->faceData, CD_ORIGSPACE))
2016                                         DM_add_face_layer(dm, CD_ORIGSPACE, CD_DEFAULT, NULL);
2017
2018                         ndm = mti->applyModifier(md, ob, dm, useRenderParams, useCache);
2019
2020                         if(ndm) {
2021                                 /* if the modifier returned a new dm, release the old one */
2022                                 if(dm && dm != ndm) dm->release(dm);
2023
2024                                 dm = ndm;
2025
2026                                 if(deformedVerts) {
2027                                         if(deformedVerts != inputVertexCos)
2028                                                 MEM_freeN(deformedVerts);
2029
2030                                         deformedVerts = NULL;
2031                                 }
2032                         } 
2033
2034                         /* create an orco derivedmesh in parallel */
2035                         if(nextmask & CD_MASK_ORCO) {
2036                                 if(!orcodm)
2037                                         orcodm= create_orco_dm(ob, me, NULL, CD_ORCO);
2038
2039                                 nextmask &= ~CD_MASK_ORCO;
2040                                 DM_set_only_copy(orcodm, nextmask | CD_MASK_ORIGINDEX);
2041                                 ndm = mti->applyModifier(md, ob, orcodm, useRenderParams, 0);
2042
2043                                 if(ndm) {
2044                                         /* if the modifier returned a new dm, release the old one */
2045                                         if(orcodm && orcodm != ndm) orcodm->release(orcodm);
2046                                         orcodm = ndm;
2047                                 }
2048                         }
2049
2050                         /* create cloth orco derivedmesh in parallel */
2051                         if(nextmask & CD_MASK_CLOTH_ORCO) {
2052                                 if(!clothorcodm)
2053                                         clothorcodm= create_orco_dm(ob, me, NULL, CD_CLOTH_ORCO);
2054
2055                                 nextmask &= ~CD_MASK_CLOTH_ORCO;
2056                                 DM_set_only_copy(clothorcodm, nextmask | CD_MASK_ORIGINDEX);
2057                                 ndm = mti->applyModifier(md, ob, clothorcodm, useRenderParams, 0);
2058
2059                                 if(ndm) {
2060                                         /* if the modifier returned a new dm, release the old one */
2061                                         if(clothorcodm && clothorcodm != ndm) clothorcodm->release(clothorcodm);
2062                                         clothorcodm = ndm;
2063                                 }
2064                         }
2065                 }
2066
2067                 isPrevDeform= (mti->type == eModifierTypeType_OnlyDeform);
2068
2069                 /* grab modifiers until index i */
2070                 if((index >= 0) && (modifiers_indexInObject(ob, md) >= index))
2071                         break;
2072
2073                 if(sculpt_mode && md->type == eModifierType_Multires)
2074                         multires_applied = 1;
2075         }
2076
2077         for(md=firstmd; md; md=md->next)
2078                 modifier_freeTemporaryData(md);
2079
2080         /* Yay, we are done. If we have a DerivedMesh and deformed vertices
2081          * need to apply these back onto the DerivedMesh. If we have no
2082          * DerivedMesh then we need to build one.
2083          */
2084         if(dm && deformedVerts) {
2085                 finaldm = CDDM_copy(dm);
2086
2087                 dm->release(dm);
2088
2089                 CDDM_apply_vert_coords(finaldm, deformedVerts);
2090                 CDDM_calc_normals(finaldm);
2091
2092                 if((dataMask & CD_MASK_WEIGHT_MCOL) && (ob->mode & OB_MODE_WEIGHT_PAINT))
2093                         add_weight_mcol_dm(ob, finaldm, draw_flag);
2094         } else if(dm) {
2095                 finaldm = dm;
2096         } else {
2097                 finaldm = CDDM_from_mesh(me, ob);
2098
2099                 if(deformedVerts) {
2100                         CDDM_apply_vert_coords(finaldm, deformedVerts);
2101                         CDDM_calc_normals(finaldm);
2102                 }
2103
2104                 if((dataMask & CD_MASK_WEIGHT_MCOL) && (ob->mode & OB_MODE_WEIGHT_PAINT))
2105                         add_weight_mcol_dm(ob, finaldm, draw_flag);
2106         }
2107
2108         /* add an orco layer if needed */
2109         if(dataMask & CD_MASK_ORCO) {
2110                 add_orco_dm(ob, NULL, finaldm, orcodm, CD_ORCO);
2111
2112                 if(deform_r && *deform_r)
2113                         add_orco_dm(ob, NULL, *deform_r, NULL, CD_ORCO);
2114         }
2115
2116 #ifdef WITH_GAMEENGINE
2117         /* NavMesh - this is a hack but saves having a NavMesh modifier */
2118         if ((ob->gameflag & OB_NAVMESH) && (finaldm->type == DM_TYPE_CDDM)) {
2119                 DerivedMesh *tdm;
2120                 tdm= navmesh_dm_createNavMeshForVisualization(finaldm);
2121                 if (finaldm != tdm) {
2122                         finaldm->release(finaldm);
2123                         finaldm= tdm;
2124                 }
2125         }
2126 #endif /* WITH_GAMEENGINE */
2127
2128         *final_r = finaldm;
2129
2130         if(orcodm)
2131                 orcodm->release(orcodm);
2132         if(clothorcodm)
2133                 clothorcodm->release(clothorcodm);
2134
2135         if(deformedVerts && deformedVerts != inputVertexCos)
2136                 MEM_freeN(deformedVerts);
2137
2138         BLI_linklist_free(datamasks, NULL);
2139 }
2140
2141 float (*editmesh_get_vertex_cos(EditMesh *em, int *numVerts_r))[3]
2142 {
2143         int i, numVerts = *numVerts_r = BLI_countlist(&em->verts);
2144         float (*cos)[3];
2145         EditVert *eve;
2146
2147         cos = MEM_mallocN(sizeof(*cos)*numVerts, "vertexcos");
2148         for (i=0,eve=em->verts.first; i<numVerts; i++,eve=eve->next) {
2149                 VECCOPY(cos[i], eve->co);
2150         }
2151
2152         return cos;
2153 }
2154
2155 int editmesh_modifier_is_enabled(Scene *scene, ModifierData *md, DerivedMesh *dm)
2156 {
2157         ModifierTypeInfo *mti = modifierType_getInfo(md->type);
2158         int required_mode = eModifierMode_Realtime | eModifierMode_Editmode;
2159
2160         if(!modifier_isEnabled(scene, md, required_mode)) return 0;
2161         if((mti->flags & eModifierTypeFlag_RequiresOriginalData) && dm) {
2162                 modifier_setError(md, "Modifier requires original data, bad stack position.");
2163                 return 0;
2164         }
2165         
2166         return 1;
2167 }
2168
2169 static void editmesh_calc_modifiers(Scene *scene, Object *ob, EditMesh *em, DerivedMesh **cage_r,
2170                                                                         DerivedMesh **final_r,
2171                                                                         CustomDataMask dataMask)
2172 {
2173         ModifierData *md;
2174         float (*deformedVerts)[3] = NULL;
2175         CustomDataMask mask;
2176         DerivedMesh *dm, *orcodm = NULL;
2177         int i, numVerts = 0, cageIndex = modifiers_getCageIndex(scene, ob, NULL, 1);
2178         LinkNode *datamasks, *curr;
2179         int required_mode = eModifierMode_Realtime | eModifierMode_Editmode;
2180
2181         modifiers_clearErrors(ob);
2182
2183         if(cage_r && cageIndex == -1) {
2184                 *cage_r = editmesh_get_derived(em, NULL);
2185         }
2186
2187         dm = NULL;
2188         md = modifiers_getVirtualModifierList(ob);
2189
2190         datamasks = modifiers_calcDataMasks(scene, ob, md, dataMask, required_mode);
2191
2192         curr = datamasks;
2193         for(i = 0; md; i++, md = md->next, curr = curr->next) {
2194                 ModifierTypeInfo *mti = modifierType_getInfo(md->type);
2195
2196                 md->scene= scene;
2197                 
2198                 if(!editmesh_modifier_is_enabled(scene, md, dm))
2199                         continue;
2200
2201                 /* add an orco layer if needed by this modifier */
2202                 if(dm && mti->requiredDataMask) {
2203                         mask = mti->requiredDataMask(ob, md);
2204                         if(mask & CD_MASK_ORCO)
2205                                 add_orco_dm(ob, em, dm, orcodm, CD_ORCO);
2206                 }
2207
2208                 /* How to apply modifier depends on (a) what we already have as
2209                  * a result of previous modifiers (could be a DerivedMesh or just
2210                  * deformed vertices) and (b) what type the modifier is.
2211                  */
2212
2213                 if(mti->type == eModifierTypeType_OnlyDeform) {
2214                         /* No existing verts to deform, need to build them. */
2215                         if(!deformedVerts) {
2216                                 if(dm) {
2217                                         /* Deforming a derived mesh, read the vertex locations
2218                                          * out of the mesh and deform them. Once done with this
2219                                          * run of deformers verts will be written back.
2220                                          */
2221                                         numVerts = dm->getNumVerts(dm);
2222                                         deformedVerts =
2223                                                 MEM_mallocN(sizeof(*deformedVerts) * numVerts, "dfmv");
2224                                         dm->getVertCos(dm, deformedVerts);
2225                                 } else {
2226                                         deformedVerts = editmesh_get_vertex_cos(em, &numVerts);
2227                                 }
2228                         }
2229
2230                         if (mti->deformVertsEM)
2231                                 mti->deformVertsEM(md, ob, em, dm, deformedVerts, numVerts);
2232                         else mti->deformVerts(md, ob, dm, deformedVerts, numVerts, 0, 0);
2233                 } else {
2234                         DerivedMesh *ndm;
2235
2236                         /* apply vertex coordinates or build a DerivedMesh as necessary */
2237                         if(dm) {
2238                                 if(deformedVerts) {
2239                                         DerivedMesh *tdm = CDDM_copy(dm);
2240                                         if(!(cage_r && dm == *cage_r)) dm->release(dm);
2241                                         dm = tdm;
2242
2243                                         CDDM_apply_vert_coords(dm, deformedVerts);
2244                                         CDDM_calc_normals(dm);
2245                                 } else if(cage_r && dm == *cage_r) {
2246                                         /* dm may be changed by this modifier, so we need to copy it
2247                                          */
2248                                         dm = CDDM_copy(dm);
2249                                 }
2250
2251                         } else {
2252                                 dm = CDDM_from_editmesh(em, ob->data);
2253
2254                                 if(deformedVerts) {
2255                                         CDDM_apply_vert_coords(dm, deformedVerts);
2256                                         CDDM_calc_normals(dm);
2257                                 }
2258                         }
2259
2260                         /* create an orco derivedmesh in parallel */
2261                         mask= (CustomDataMask)GET_INT_FROM_POINTER(curr->link);
2262                         if(mask & CD_MASK_ORCO) {
2263                                 if(!orcodm)
2264                                         orcodm= create_orco_dm(ob, ob->data, em, CD_ORCO);
2265
2266                                 mask &= ~CD_MASK_ORCO;
2267                                 DM_set_only_copy(orcodm, mask | CD_MASK_ORIGINDEX);
2268
2269                                 if (mti->applyModifierEM)
2270                                         ndm = mti->applyModifierEM(md, ob, em, orcodm);
2271                                 else
2272                                         ndm = mti->applyModifier(md, ob, orcodm, 0, 0);
2273
2274                                 if(ndm) {
2275                                         /* if the modifier returned a new dm, release the old one */
2276                                         if(orcodm && orcodm != ndm) orcodm->release(orcodm);
2277                                         orcodm = ndm;
2278                                 }
2279                         }
2280
2281                         /* set the DerivedMesh to only copy needed data */
2282                         mask= (CustomDataMask)GET_INT_FROM_POINTER(curr->link); /* CD_MASK_ORCO may have been cleared above */
2283
2284                         DM_set_only_copy(dm, mask | CD_MASK_ORIGINDEX);
2285
2286                         if(mask & CD_MASK_ORIGSPACE)
2287                                 if(!CustomData_has_layer(&dm->faceData, CD_ORIGSPACE))
2288                                         DM_add_face_layer(dm, CD_ORIGSPACE, CD_DEFAULT, NULL);
2289                         
2290                         if (mti->applyModifierEM)
2291                                 ndm = mti->applyModifierEM(md, ob, em, dm);
2292                         else
2293                                 ndm = mti->applyModifier(md, ob, dm, 0, 0);
2294
2295                         if (ndm) {
2296                                 if(dm && dm != ndm)
2297                                         dm->release(dm);
2298
2299                                 dm = ndm;
2300
2301                                 if (deformedVerts) {
2302                                         MEM_freeN(deformedVerts);
2303                                         deformedVerts = NULL;
2304                                 }
2305                         }
2306                 }
2307
2308                 if(cage_r && i == cageIndex) {
2309                         if(dm && deformedVerts) {
2310                                 *cage_r = CDDM_copy(dm);
2311                                 CDDM_apply_vert_coords(*cage_r, deformedVerts);
2312                         } else if(dm) {
2313                                 *cage_r = dm;
2314                         } else {
2315                                 *cage_r =
2316                                         editmesh_get_derived(em,
2317                                                 deformedVerts ? MEM_dupallocN(deformedVerts) : NULL);
2318                         }
2319                 }
2320         }
2321
2322         BLI_linklist_free(datamasks, NULL);
2323
2324         /* Yay, we are done. If we have a DerivedMesh and deformed vertices need
2325          * to apply these back onto the DerivedMesh. If we have no DerivedMesh
2326          * then we need to build one.
2327          */
2328         if(dm && deformedVerts) {
2329                 *final_r = CDDM_copy(dm);
2330
2331                 if(!(cage_r && dm == *cage_r)) dm->release(dm);
2332
2333                 CDDM_apply_vert_coords(*final_r, deformedVerts);
2334                 CDDM_calc_normals(*final_r);
2335         } else if (dm) {
2336                 *final_r = dm;
2337         } else if (!deformedVerts && cage_r && *cage_r) {
2338                 *final_r = *cage_r;
2339         } else {
2340                 *final_r = editmesh_get_derived(em, deformedVerts);
2341                 deformedVerts = NULL;
2342         }
2343
2344         /* add an orco layer if needed */
2345         if(dataMask & CD_MASK_ORCO)
2346                 add_orco_dm(ob, em, *final_r, orcodm, CD_ORCO);
2347
2348         if(orcodm)
2349                 orcodm->release(orcodm);
2350
2351         if(deformedVerts)
2352                 MEM_freeN(deformedVerts);
2353 }
2354
2355 static void clear_mesh_caches(Object *ob)
2356 {
2357         Mesh *me= ob->data;
2358
2359                 /* also serves as signal to remake texspace */
2360         if (ob->bb) {
2361                 MEM_freeN(ob->bb);
2362                 ob->bb = NULL;
2363         }
2364         if (me->bb) {
2365                 MEM_freeN(me->bb);
2366                 me->bb = NULL;
2367         }
2368
2369         freedisplist(&ob->disp);
2370
2371         if (ob->derivedFinal) {
2372                 ob->derivedFinal->needsFree = 1;
2373                 ob->derivedFinal->release(ob->derivedFinal);
2374                 ob->derivedFinal= NULL;
2375         }
2376         if (ob->derivedDeform) {
2377                 ob->derivedDeform->needsFree = 1;
2378                 ob->derivedDeform->release(ob->derivedDeform);
2379                 ob->derivedDeform= NULL;
2380         }
2381
2382         if(ob->sculpt) {
2383                 object_sculpt_modifiers_changed(ob);
2384         }
2385 }
2386
2387 static void mesh_build_data(Scene *scene, Object *ob, CustomDataMask dataMask)
2388 {
2389         Object *obact = scene->basact?scene->basact->object:NULL;
2390         int editing = paint_facesel_test(ob) || paint_vertsel_test(ob);/* paint_vertsel_test */
2391         /* weight paint and face select need original indices because of selection buffer drawing */
2392         int needMapping = (ob==obact) && (editing || (ob->mode & (OB_MODE_WEIGHT_PAINT|OB_MODE_VERTEX_PAINT)));
2393
2394         clear_mesh_caches(ob);
2395
2396         mesh_calc_modifiers(scene, ob, NULL, &ob->derivedDeform,
2397                                                 &ob->derivedFinal, 0, 1,
2398                                                 needMapping, dataMask, -1, 1);
2399
2400         DM_set_object_boundbox (ob, ob->derivedFinal);
2401
2402         ob->derivedFinal->needsFree = 0;
2403         ob->derivedDeform->needsFree = 0;
2404         ob->lastDataMask = dataMask;
2405 }
2406
2407 static void editmesh_build_data(Scene *scene, Object *obedit, EditMesh *em, CustomDataMask dataMask)
2408 {
2409         clear_mesh_caches(obedit);
2410
2411         if (em->derivedFinal) {
2412                 if (em->derivedFinal!=em->derivedCage) {
2413                         em->derivedFinal->needsFree = 1;
2414                         em->derivedFinal->release(em->derivedFinal);
2415                 }
2416                 em->derivedFinal = NULL;
2417         }
2418         if (em->derivedCage) {
2419                 em->derivedCage->needsFree = 1;
2420                 em->derivedCage->release(em->derivedCage);
2421                 em->derivedCage = NULL;
2422         }
2423
2424         editmesh_calc_modifiers(scene, obedit, em, &em->derivedCage, &em->derivedFinal, dataMask);
2425         DM_set_object_boundbox (obedit, em->derivedFinal);
2426
2427         em->lastDataMask = dataMask;
2428         em->derivedFinal->needsFree = 0;
2429         em->derivedCage->needsFree = 0;
2430 }
2431
2432 void makeDerivedMesh(Scene *scene, Object *ob, EditMesh *em, CustomDataMask dataMask)
2433 {
2434         if (em) {
2435                 editmesh_build_data(scene, ob, em, dataMask);
2436         } else {
2437                 mesh_build_data(scene, ob, dataMask);
2438         }
2439 }
2440
2441 /***/
2442
2443 DerivedMesh *mesh_get_derived_final(Scene *scene, Object *ob, CustomDataMask dataMask)
2444 {
2445         /* if there's no derived mesh or the last data mask used doesn't include
2446          * the data we need, rebuild the derived mesh
2447          */
2448         if(!ob->derivedFinal || (dataMask & ob->lastDataMask) != dataMask)
2449                 mesh_build_data(scene, ob, dataMask);
2450
2451         return ob->derivedFinal;
2452 }
2453
2454 DerivedMesh *mesh_get_derived_deform(Scene *scene, Object *ob, CustomDataMask dataMask)
2455 {
2456         /* if there's no derived mesh or the last data mask used doesn't include
2457          * the data we need, rebuild the derived mesh
2458          */
2459         if(!ob->derivedDeform || (dataMask & ob->lastDataMask) != dataMask)
2460                 mesh_build_data(scene, ob, dataMask);
2461
2462         return ob->derivedDeform;
2463 }
2464
2465 DerivedMesh *mesh_create_derived_render(Scene *scene, Object *ob, CustomDataMask dataMask)
2466 {
2467         DerivedMesh *final;
2468         
2469         mesh_calc_modifiers(scene, ob, NULL, NULL, &final, 1, 1, 0, dataMask, -1, 0);
2470
2471         return final;
2472 }
2473
2474 DerivedMesh *mesh_create_derived_index_render(Scene *scene, Object *ob, CustomDataMask dataMask, int index)
2475 {
2476         DerivedMesh *final;
2477         
2478         mesh_calc_modifiers(scene, ob, NULL, NULL, &final, 1, 1, 0, dataMask, index, 0);
2479
2480         return final;
2481 }
2482
2483 DerivedMesh *mesh_create_derived_view(Scene *scene, Object *ob, CustomDataMask dataMask)
2484 {
2485         DerivedMesh *final;
2486
2487         mesh_calc_modifiers(scene, ob, NULL, NULL, &final, 0, 1, 0, dataMask, -1, 0);
2488
2489         return final;
2490 }
2491
2492 DerivedMesh *mesh_create_derived_no_deform(Scene *scene, Object *ob, float (*vertCos)[3],
2493                                                                                    CustomDataMask dataMask)
2494 {
2495         DerivedMesh *final;
2496         
2497         mesh_calc_modifiers(scene, ob, vertCos, NULL, &final, 0, 0, 0, dataMask, -1, 0);
2498
2499         return final;
2500 }
2501
2502 DerivedMesh *mesh_create_derived_no_virtual(Scene *scene, Object *ob, float (*vertCos)[3],
2503                                                                                         CustomDataMask dataMask)
2504 {
2505         DerivedMesh *final;
2506         
2507         mesh_calc_modifiers(scene, ob, vertCos, NULL, &final, 0, -1, 0, dataMask, -1, 0);
2508
2509         return final;
2510 }
2511
2512 DerivedMesh *mesh_create_derived_physics(Scene *scene, Object *ob, float (*vertCos)[3],
2513                                                                                         CustomDataMask dataMask)
2514 {
2515         DerivedMesh *final;
2516         
2517         mesh_calc_modifiers(scene, ob, vertCos, NULL, &final, 0, -1, 1, dataMask, -1, 0);
2518
2519         return final;
2520 }
2521
2522 DerivedMesh *mesh_create_derived_no_deform_render(Scene *scene, Object *ob,
2523                                                                                                   float (*vertCos)[3],
2524                                                                                                   CustomDataMask dataMask)
2525 {
2526         DerivedMesh *final;
2527
2528         mesh_calc_modifiers(scene, ob, vertCos, NULL, &final, 1, 0, 0, dataMask, -1, 0);
2529
2530         return final;
2531 }
2532
2533 /***/
2534
2535 DerivedMesh *editmesh_get_derived_cage_and_final(Scene *scene, Object *obedit, EditMesh *em, DerivedMesh **final_r,
2536                                                                                                  CustomDataMask dataMask)
2537 {
2538         /* if there's no derived mesh or the last data mask used doesn't include
2539          * the data we need, rebuild the derived mesh
2540          */
2541         if(!em->derivedCage ||
2542            (em->lastDataMask & dataMask) != dataMask)
2543                 editmesh_build_data(scene, obedit, em, dataMask);
2544
2545         *final_r = em->derivedFinal;
2546         return em->derivedCage;
2547 }
2548
2549 DerivedMesh *editmesh_get_derived_cage(Scene *scene, Object *obedit, EditMesh *em, CustomDataMask dataMask)
2550 {
2551         /* if there's no derived mesh or the last data mask used doesn't include
2552          * the data we need, rebuild the derived mesh
2553          */
2554         if(!em->derivedCage ||
2555            (em->lastDataMask & dataMask) != dataMask)
2556                 editmesh_build_data(scene, obedit, em, dataMask);
2557
2558         return em->derivedCage;
2559 }
2560
2561 DerivedMesh *editmesh_get_derived_base(Object *UNUSED(obedit), EditMesh *em)
2562 {
2563         return editmesh_get_derived(em, NULL);
2564 }
2565
2566
2567 /* ********* For those who don't grasp derived stuff! (ton) :) *************** */
2568
2569 static void make_vertexcosnos__mapFunc(void *userData, int index, float *co, float *no_f, short *no_s)
2570 {
2571         float *vec = userData;
2572         
2573         vec+= 6*index;
2574
2575         /* check if we've been here before (normal should not be 0) */
2576         if(vec[3] || vec[4] || vec[5]) return;
2577
2578         copy_v3_v3(vec, co);
2579         vec+= 3;
2580         if(no_f) {
2581                 copy_v3_v3(vec, no_f);
2582         }
2583         else {
2584                 normal_short_to_float_v3(vec, no_s);
2585         }
2586 }
2587
2588 /* always returns original amount me->totvert of vertices and normals, but fully deformed and subsurfered */
2589 /* this is needed for all code using vertexgroups (no subsurf support) */
2590 /* it stores the normals as floats, but they can still be scaled as shorts (32767 = unit) */
2591 /* in use now by vertex/weight paint and particle generating */
2592
2593 float *mesh_get_mapped_verts_nors(Scene *scene, Object *ob)
2594 {
2595         Mesh *me= ob->data;
2596         DerivedMesh *dm;
2597         float *vertexcosnos;
2598         
2599         /* lets prevent crashing... */
2600         if(ob->type!=OB_MESH || me->totvert==0)
2601                 return NULL;
2602         
2603         dm= mesh_get_derived_final(scene, ob, CD_MASK_BAREMESH);
2604         vertexcosnos= MEM_callocN(6*sizeof(float)*me->totvert, "vertexcosnos map");
2605         
2606         if(dm->foreachMappedVert) {
2607                 dm->foreachMappedVert(dm, make_vertexcosnos__mapFunc, vertexcosnos);
2608         }
2609         else {
2610                 float *fp= vertexcosnos;
2611                 int a;
2612                 
2613                 for(a=0; a< me->totvert; a++, fp+=6) {
2614                         dm->getVertCo(dm, a, fp);
2615                         dm->getVertNo(dm, a, fp+3);
2616                 }
2617         }
2618         
2619         dm->release(dm);
2620         return vertexcosnos;
2621 }
2622
2623 /* ******************* GLSL ******************** */
2624
2625 typedef struct
2626 {
2627         float * precomputedFaceNormals;
2628         MTFace * mtface;        // texture coordinates
2629         MFace * mface;          // indices
2630         MVert * mvert;          // vertices & normals
2631         float (*orco)[3];
2632         float (*tangent)[4];    // destination
2633         int numFaces;
2634
2635 } SGLSLMeshToTangent;
2636
2637 // interface
2638 #include "mikktspace.h"
2639
2640 static int GetNumFaces(const SMikkTSpaceContext * pContext)
2641 {
2642         SGLSLMeshToTangent * pMesh = (SGLSLMeshToTangent *) pContext->m_pUserData;
2643         return pMesh->numFaces;
2644 }
2645
2646 static int GetNumVertsOfFace(const SMikkTSpaceContext * pContext, const int face_num)
2647 {
2648         SGLSLMeshToTangent * pMesh = (SGLSLMeshToTangent *) pContext->m_pUserData;
2649         return pMesh->mface[face_num].v4!=0 ? 4 : 3;
2650 }
2651
2652 static void GetPosition(const SMikkTSpaceContext * pContext, float fPos[], const int face_num, const int vert_index)
2653 {
2654         //assert(vert_index>=0 && vert_index<4);
2655         SGLSLMeshToTangent * pMesh = (SGLSLMeshToTangent *) pContext->m_pUserData;
2656         const float *co= pMesh->mvert[(&pMesh->mface[face_num].v1)[vert_index]].co;
2657         VECCOPY(fPos, co);
2658 }
2659
2660 static void GetTextureCoordinate(const SMikkTSpaceContext * pContext, float fUV[], const int face_num, const int vert_index)
2661 {
2662         //assert(vert_index>=0 && vert_index<4);
2663         SGLSLMeshToTangent * pMesh = (SGLSLMeshToTangent *) pContext->m_pUserData;
2664
2665         if(pMesh->mtface!=NULL) {
2666                 float * uv = pMesh->mtface[face_num].uv[vert_index];
2667                 fUV[0]=uv[0]; fUV[1]=uv[1];
2668         }
2669         else {
2670                 const float *orco= pMesh->orco[(&pMesh->mface[face_num].v1)[vert_index]];
2671                 map_to_sphere( &fUV[0], &fUV[1], orco[0], orco[1], orco[2]);
2672         }
2673 }
2674
2675 static void GetNormal(const SMikkTSpaceContext * pContext, float fNorm[], const int face_num, const int vert_index)
2676 {
2677         //assert(vert_index>=0 && vert_index<4);
2678         SGLSLMeshToTangent * pMesh = (SGLSLMeshToTangent *) pContext->m_pUserData;
2679
2680         const int smoothnormal = (pMesh->mface[face_num].flag & ME_SMOOTH);
2681         if(!smoothnormal) {     // flat
2682                 if(pMesh->precomputedFaceNormals) {
2683                         VECCOPY(fNorm, &pMesh->precomputedFaceNormals[3*face_num]);
2684                 }
2685                 else {
2686                         MFace *mf= &pMesh->mface[face_num];
2687                         float *p0= pMesh->mvert[mf->v1].co;
2688                         float *p1= pMesh->mvert[mf->v2].co;
2689                         float *p2= pMesh->mvert[mf->v3].co;
2690
2691                         if(mf->v4) {
2692                                 float *p3 = pMesh->mvert[mf->v4].co;
2693                                 normal_quad_v3(fNorm, p0, p1, p2, p3);
2694                         }
2695                         else {
2696                                 normal_tri_v3(fNorm, p0, p1, p2);
2697                         }
2698                 }
2699         }
2700         else {
2701                 const short *no= pMesh->mvert[(&pMesh->mface[face_num].v1)[vert_index]].no;
2702                 normal_short_to_float_v3(fNorm, no);
2703         }
2704 }
2705 static void SetTSpace(const SMikkTSpaceContext * pContext, const float fvTangent[], const float fSign, const int face_num, const int iVert)
2706 {
2707         //assert(vert_index>=0 && vert_index<4);
2708         SGLSLMeshToTangent * pMesh = (SGLSLMeshToTangent *) pContext->m_pUserData;
2709         float * pRes = pMesh->tangent[4*face_num+iVert];
2710         VECCOPY(pRes, fvTangent);
2711         pRes[3]=fSign;
2712 }
2713
2714
2715 void DM_add_tangent_layer(DerivedMesh *dm)
2716 {
2717         /* mesh vars */
2718         MTFace *mtface, *tf;
2719         MFace *mface, *mf;
2720         MVert *mvert, *v1, *v2, *v3, *v4;
2721         MemArena *arena= NULL;
2722         VertexTangent **vtangents= NULL;
2723         float (*orco)[3]= NULL, (*tangent)[4];
2724         float *uv1, *uv2, *uv3, *uv4, *vtang;
2725         float fno[3], tang[3], uv[4][2];
2726         int i, j, len, mf_vi[4], totvert, totface, iCalcNewMethod;
2727         float *nors;
2728
2729         if(CustomData_get_layer_index(&dm->faceData, CD_TANGENT) != -1)
2730                 return;
2731
2732         nors = dm->getFaceDataArray(dm, CD_NORMAL);
2733
2734         /* check we have all the needed layers */
2735         totvert= dm->getNumVerts(dm);
2736         totface= dm->getNumFaces(dm);
2737
2738         mvert= dm->getVertArray(dm);
2739         mface= dm->getFaceArray(dm);
2740         mtface= dm->getFaceDataArray(dm, CD_MTFACE);
2741
2742         if(!mtface) {
2743                 orco= dm->getVertDataArray(dm, CD_ORCO);
2744                 if(!orco)
2745                         return;
2746         }
2747         
2748         /* create tangent layer */
2749         DM_add_face_layer(dm, CD_TANGENT, CD_CALLOC, NULL);
2750         tangent= DM_get_face_data_layer(dm, CD_TANGENT);
2751         
2752         /* allocate some space */
2753         arena= BLI_memarena_new(BLI_MEMARENA_STD_BUFSIZE, "tangent layer arena");
2754         BLI_memarena_use_calloc(arena);
2755         vtangents= MEM_callocN(sizeof(VertexTangent*)*totvert, "VertexTangent");
2756
2757         // new computation method
2758         iCalcNewMethod = 1;
2759         if(iCalcNewMethod != 0) {
2760                 SGLSLMeshToTangent mesh2tangent= {0};
2761                 SMikkTSpaceContext sContext= {0};
2762                 SMikkTSpaceInterface sInterface= {0};
2763
2764                 mesh2tangent.precomputedFaceNormals = nors;
2765                 mesh2tangent.mtface = mtface;
2766                 mesh2tangent.mface = mface;
2767                 mesh2tangent.mvert = mvert;
2768                 mesh2tangent.orco = orco;
2769                 mesh2tangent.tangent = tangent;
2770                 mesh2tangent.numFaces = totface;
2771
2772                 sContext.m_pUserData = &mesh2tangent;
2773                 sContext.m_pInterface = &sInterface;
2774                 sInterface.m_getNumFaces = GetNumFaces;
2775                 sInterface.m_getNumVerticesOfFace = GetNumVertsOfFace;
2776                 sInterface.m_getPosition = GetPosition;
2777                 sInterface.m_getTexCoord = GetTextureCoordinate;
2778                 sInterface.m_getNormal = GetNormal;
2779                 sInterface.m_setTSpaceBasic = SetTSpace;
2780
2781                 // 0 if failed
2782                 iCalcNewMethod = genTangSpaceDefault(&sContext);
2783         }
2784
2785         if(!iCalcNewMethod) {
2786                 /* sum tangents at connected vertices */
2787                 for(i=0, tf=mtface, mf=mface; i < totface; mf++, tf++, i++) {
2788                         v1= &mvert[mf->v1];
2789                         v2= &mvert[mf->v2];
2790                         v3= &mvert[mf->v3];
2791
2792                         if (mf->v4) {
2793                                 v4= &mvert[mf->v4];
2794                                 normal_quad_v3( fno,v4->co, v3->co, v2->co, v1->co);
2795                         }
2796                         else {
2797                                 v4= NULL;
2798                                 normal_tri_v3( fno,v3->co, v2->co, v1->co);
2799                         }
2800                 
2801                         if(mtface) {
2802                                 uv1= tf->uv[0];
2803                                 uv2= tf->uv[1];
2804                                 uv3= tf->uv[2];
2805                                 uv4= tf->uv[3];
2806                         }
2807                         else {
2808                                 uv1= uv[0]; uv2= uv[1]; uv3= uv[2]; uv4= uv[3];
2809                                 map_to_sphere( &uv[0][0], &uv[0][1],orco[mf->v1][0], orco[mf->v1][1], orco[mf->v1][2]);
2810                                 map_to_sphere( &uv[1][0], &uv[1][1],orco[mf->v2][0], orco[mf->v2][1], orco[mf->v2][2]);
2811                                 map_to_sphere( &uv[2][0], &uv[2][1],orco[mf->v3][0], orco[mf->v3][1], orco[mf->v3][2]);
2812                                 if(v4)
2813                                         map_to_sphere( &uv[3][0], &uv[3][1],orco[mf->v4][0], orco[mf->v4][1], orco[mf->v4][2]);
2814                         }
2815                 
2816                         tangent_from_uv(uv1, uv2, uv3, v1->co, v2->co, v3->co, fno, tang);
2817                         sum_or_add_vertex_tangent(arena, &vtangents[mf->v1], tang, uv1);
2818                         sum_or_add_vertex_tangent(arena, &vtangents[mf->v2], tang, uv2);
2819                         sum_or_add_vertex_tangent(arena, &vtangents[mf->v3], tang, uv3);
2820                 
2821                         if(mf->v4) {
2822                                 v4= &mvert[mf->v4];
2823                         
2824                                 tangent_from_uv(uv1, uv3, uv4, v1->co, v3->co, v4->co, fno, tang);
2825                                 sum_or_add_vertex_tangent(arena, &vtangents[mf->v1], tang, uv1);
2826                                 sum_or_add_vertex_tangent(arena, &vtangents[mf->v3], tang, uv3);
2827                                 sum_or_add_vertex_tangent(arena, &vtangents[mf->v4], tang, uv4);
2828                         }
2829                 }
2830         
2831                 /* write tangent to layer */
2832                 for(i=0, tf=mtface, mf=mface; i < totface; mf++, tf++, i++, tangent+=4) {
2833                         len= (mf->v4)? 4 : 3; 
2834
2835                         if(mtface == NULL) {
2836                                 map_to_sphere( &uv[0][0], &uv[0][1],orco[mf->v1][0], orco[mf->v1][1], orco[mf->v1][2]);
2837                                 map_to_sphere( &uv[1][0], &uv[1][1],orco[mf->v2][0], orco[mf->v2][1], orco[mf->v2][2]);
2838                                 map_to_sphere( &uv[2][0], &uv[2][1],orco[mf->v3][0], orco[mf->v3][1], orco[mf->v3][2]);
2839                                 if(len==4)
2840                                         map_to_sphere( &uv[3][0], &uv[3][1],orco[mf->v4][0], orco[mf->v4][1], orco[mf->v4][2]);
2841                         }
2842                 
2843                         mf_vi[0]= mf->v1;
2844                         mf_vi[1]= mf->v2;
2845                         mf_vi[2]= mf->v3;
2846                         mf_vi[3]= mf->v4;
2847                 
2848                         for(j=0; j<len; j++) {
2849                                 vtang= find_vertex_tangent(vtangents[mf_vi[j]], mtface ? tf->uv[j] : uv[j]);
2850                                 normalize_v3_v3(tangent[j], vtang);
2851                                 ((float *) tangent[j])[3]=1.0f;
2852                         }
2853                 }
2854         }
2855         
2856         BLI_memarena_free(arena);
2857         MEM_freeN(vtangents);
2858 }
2859
2860 void DM_vertex_attributes_from_gpu(DerivedMesh *dm, GPUVertexAttribs *gattribs, DMVertexAttribs *attribs)
2861 {
2862         CustomData *vdata, *fdata, *tfdata = NULL;
2863         int a, b, layer;
2864
2865         /* From the layers requested by the GLSL shader, figure out which ones are
2866          * actually available for this derivedmesh, and retrieve the pointers */
2867
2868         memset(attribs, 0, sizeof(DMVertexAttribs));
2869
2870         vdata = &dm->vertData;
2871         fdata = &dm->faceData;
2872
2873         /* ugly hack, editmesh derivedmesh doesn't copy face data, this way we
2874          * can use offsets instead */
2875         if(dm->release == emDM_release)
2876                 tfdata = &((EditMeshDerivedMesh*)dm)->em->fdata;
2877         else
2878                 tfdata = fdata;
2879
2880         /* add a tangent layer if necessary */
2881         for(b = 0; b < gattribs->totlayer; b++)
2882                 if(gattribs->layer[b].type == CD_TANGENT)
2883                         if(CustomData_get_layer_index(fdata, CD_TANGENT) == -1)
2884                                 DM_add_tangent_layer(dm);
2885
2886         for(b = 0; b < gattribs->totlayer; b++) {
2887                 if(gattribs->layer[b].type == CD_MTFACE) {
2888                         /* uv coordinates */
2889                         if(gattribs->layer[b].name[0])
2890                                 layer = CustomData_get_named_layer_index(tfdata, CD_MTFACE,
2891                                         gattribs->layer[b].name);
2892                         else
2893                                 layer = CustomData_get_active_layer_index(tfdata, CD_MTFACE);
2894
2895                         if(layer != -1) {
2896                                 a = attribs->tottface++;
2897
2898                                 attribs->tface[a].array = tfdata->layers[layer].data;
2899                                 attribs->tface[a].emOffset = tfdata->layers[layer].offset;
2900                                 attribs->tface[a].glIndex = gattribs->layer[b].glindex;
2901                                 attribs->tface[a].glTexco = gattribs->layer[b].gltexco;
2902                         }
2903                 }
2904                 else if(gattribs->layer[b].type == CD_MCOL) {
2905                         /* vertex colors */
2906                         if(gattribs->layer[b].name[0])
2907                                 layer = CustomData_get_named_layer_index(tfdata, CD_MCOL,
2908                                         gattribs->layer[b].name);
2909                         else
2910                                 layer = CustomData_get_active_layer_index(tfdata, CD_MCOL);
2911
2912                         if(layer != -1) {
2913                                 a = attribs->totmcol++;
2914
2915                                 attribs->mcol[a].array = tfdata->layers[layer].data;
2916                                 attribs->mcol[a].emOffset = tfdata->layers[layer].offset;
2917                                 attribs->mcol[a].glIndex = gattribs->layer[b].glindex;
2918                         }
2919                 }
2920                 else if(gattribs->layer[b].type == CD_TANGENT) {
2921                         /* tangents */
2922                         layer = CustomData_get_layer_index(fdata, CD_TANGENT);
2923
2924                         if(layer != -1) {
2925                                 attribs->tottang = 1;
2926
2927                                 attribs->tang.array = fdata->layers[layer].data;
2928                                 attribs->tang.emOffset = fdata->layers[layer].offset;
2929                                 attribs->tang.glIndex = gattribs->layer[b].glindex;
2930                         }
2931                 }
2932                 else if(gattribs->layer[b].type == CD_ORCO) {
2933                         /* original coordinates */
2934                         layer = CustomData_get_layer_index(vdata, CD_ORCO);
2935
2936                         if(layer != -1) {
2937                                 attribs->totorco = 1;
2938
2939                                 attribs->orco.array = vdata->layers[layer].data;
2940                                 attribs->orco.emOffset = vdata->layers[layer].offset;
2941                                 attribs->orco.glIndex = gattribs->layer[b].glindex;
2942                                 attribs->orco.glTexco = gattribs->layer[b].gltexco;
2943                         }
2944                 }
2945         }
2946 }
2947
2948 /* Set object's bounding box based on DerivedMesh min/max data */
2949 void DM_set_object_boundbox(Object *ob, DerivedMesh *dm)
2950 {
2951         float min[3], max[3];
2952
2953         INIT_MINMAX(min, max);
2954
2955         dm->getMinMax(dm, min, max);
2956
2957         if(!ob->bb)
2958                 ob->bb= MEM_callocN(sizeof(BoundBox), "DM-BoundBox");
2959
2960         boundbox_set_from_min_max(ob->bb, min, max);
2961 }
2962
2963 /* --- NAVMESH (begin) --- */
2964 #ifdef WITH_GAMEENGINE
2965
2966 BM_INLINE int navmesh_bit(int a, int b)
2967 {
2968         return (a & (1 << b)) >> b;
2969 }
2970
2971 BM_INLINE void navmesh_intToCol(int i, float col[3])
2972 {
2973         int     r = navmesh_bit(i, 0) + navmesh_bit(i, 3) * 2 + 1;
2974         int     g = navmesh_bit(i, 1) + navmesh_bit(i, 4) * 2 + 1;
2975         int     b = navmesh_bit(i, 2) + navmesh_bit(i, 5) * 2 + 1;
2976         col[0] = 1 - r*63.0f/255.0f;
2977         col[1] = 1 - g*63.0f/255.0f;
2978         col[2] = 1 - b*63.0f/255.0f;
2979 }
2980
2981 static void navmesh_drawColored(DerivedMesh *dm)
2982 {
2983         int a, glmode;
2984         MVert *mvert = (MVert *)CustomData_get_layer(&dm->vertData, CD_MVERT);
2985         MFace *mface = (MFace *)CustomData_get_layer(&dm->faceData, CD_MFACE);
2986         int *polygonIdx = (int *)CustomData_get_layer(&dm->faceData, CD_RECAST);
2987         float col[3];
2988
2989         if (!polygonIdx)
2990                 return;
2991
2992         /*
2993         //UI_ThemeColor(TH_WIRE);
2994         glDisable(GL_LIGHTING);
2995         glLineWidth(2.0);
2996         dm->drawEdges(dm, 0, 1);
2997         glLineWidth(1.0);
2998         glEnable(GL_LIGHTING);*/
2999
3000         glDisable(GL_LIGHTING);
3001         /*  if(GPU_buffer_legacy(dm) ) */ { /* TODO - VBO draw code, not high priority - campbell */
3002                 DEBUG_VBO( "Using legacy code. drawNavMeshColored\n" );
3003                 //glShadeModel(GL_SMOOTH);
3004                 glBegin(glmode = GL_QUADS);
3005                 for(a = 0; a < dm->numFaceData; a++, mface++) {
3006                         int new_glmode = mface->v4?GL_QUADS:GL_TRIANGLES;
3007                         int pi = polygonIdx[a];
3008                         if (pi <= 0) {
3009                                 zero_v3(col);
3010                         }
3011                         else {
3012                                 navmesh_intToCol(pi, col);
3013            &nbs