svn merge ^/trunk/blender -r48095:48592
[blender.git] / source / blender / editors / mesh / meshtools.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) 2004 by 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/editors/mesh/meshtools.c
29  *  \ingroup edmesh
30  */
31
32
33 /*
34  * meshtools.c: no editmode (violated already :), tools operating on meshes
35  */
36
37 #include <stddef.h>
38 #include <stdlib.h>
39 #include <math.h>
40 #include <float.h>
41
42 #include "MEM_guardedalloc.h"
43
44 #include "DNA_mesh_types.h"
45 #include "DNA_key_types.h"
46 #include "DNA_material_types.h"
47 #include "DNA_meshdata_types.h"
48 #include "DNA_modifier_types.h"
49 #include "DNA_object_types.h"
50 #include "DNA_scene_types.h"
51
52 #include "BLI_math.h"
53 #include "BLI_blenlib.h"
54 #include "BLI_utildefines.h"
55 #include "BLI_ghash.h"
56 #include "BLI_rand.h" /* for randome face sorting */
57 #include "BLI_threads.h"
58
59
60 #include "BKE_context.h"
61 #include "BKE_depsgraph.h"
62 #include "BKE_deform.h"
63 #include "BKE_DerivedMesh.h"
64 #include "BKE_key.h"
65 #include "BKE_library.h"
66 #include "BKE_main.h"
67 #include "BKE_mesh.h"
68 #include "BKE_material.h"
69 #include "BKE_report.h"
70 #include "BKE_tessmesh.h"
71 #include "BKE_multires.h"
72
73 #include "BLO_sys_types.h" // for intptr_t support
74
75 #include "ED_mesh.h"
76 #include "ED_object.h"
77 #include "ED_view3d.h"
78
79 #include "WM_api.h"
80 #include "WM_types.h"
81
82 /* own include */
83 #include "mesh_intern.h"
84 #include "uvedit_intern.h"
85
86 /* * ********************** no editmode!!! *********** */
87
88 /*********************** JOIN ***************************/
89
90 /* join selected meshes into the active mesh, context sensitive
91  * return 0 if no join is made (error) and 1 if the join is done */
92
93 int join_mesh_exec(bContext *C, wmOperator *op)
94 {
95         Main *bmain = CTX_data_main(C);
96         Scene *scene = CTX_data_scene(C);
97         Object *ob = CTX_data_active_object(C);
98         Material **matar, *ma;
99         Mesh *me;
100         MVert *mvert, *mv;
101         MEdge *medge = NULL;
102         MPoly *mpoly = NULL;
103         MLoop *mloop = NULL;
104         Key *key, *nkey = NULL;
105         KeyBlock *kb, *okb, *kbn;
106         float imat[4][4], cmat[4][4], *fp1, *fp2, curpos;
107         int a, b, totcol, totmat = 0, totedge = 0, totvert = 0, ok = 0;
108         int totloop = 0, totpoly = 0, vertofs, *matmap = NULL;
109         int i, j, index, haskey = 0, edgeofs, loopofs, polyofs;
110         bDeformGroup *dg, *odg;
111         MDeformVert *dvert;
112         CustomData vdata, edata, fdata, ldata, pdata;
113
114         if (scene->obedit) {
115                 BKE_report(op->reports, RPT_WARNING, "Cant join while in editmode");
116                 return OPERATOR_CANCELLED;
117         }
118         
119         /* ob is the object we are adding geometry to */
120         if (!ob || ob->type != OB_MESH) {
121                 BKE_report(op->reports, RPT_WARNING, "Active object is not a mesh");
122                 return OPERATOR_CANCELLED;
123         }
124         
125         /* count & check */
126         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
127         {
128                 if (base->object->type == OB_MESH) {
129                         me = base->object->data;
130
131                         totvert += me->totvert;
132                         totedge += me->totedge;
133                         totloop += me->totloop;
134                         totpoly += me->totpoly;
135                         totmat += base->object->totcol;
136                         
137                         if (base->object == ob)
138                                 ok = 1;
139                         
140                         /* check for shapekeys */
141                         if (me->key)
142                                 haskey++;
143                 }
144         }
145         CTX_DATA_END;
146         
147         /* that way the active object is always selected */ 
148         if (ok == 0) {
149                 BKE_report(op->reports, RPT_WARNING, "Active object is not a selected mesh");
150                 return OPERATOR_CANCELLED;
151         }
152         
153         /* only join meshes if there are verts to join, there aren't too many, and we only had one mesh selected */
154         me = (Mesh *)ob->data;
155         key = me->key;
156
157         if (totvert == 0 || totvert == me->totvert) {
158                 BKE_report(op->reports, RPT_WARNING, "No mesh data to join");
159                 return OPERATOR_CANCELLED;
160         }
161         
162         if (totvert > MESH_MAX_VERTS) {
163                 BKE_reportf(op->reports, RPT_WARNING, "Joining results in %d vertices, limit is " STRINGIFY(MESH_MAX_VERTS), totvert);
164                 return OPERATOR_CANCELLED;              
165         }
166
167         /* new material indices and material array */
168         matar = MEM_callocN(sizeof(void *) * totmat, "join_mesh matar");
169         if (totmat) matmap = MEM_callocN(sizeof(int) * totmat, "join_mesh matmap");
170         totcol = ob->totcol;
171         
172         /* obact materials in new main array, is nicer start! */
173         for (a = 0; a < ob->totcol; a++) {
174                 matar[a] = give_current_material(ob, a + 1);
175                 id_us_plus((ID *)matar[a]);
176                 /* increase id->us : will be lowered later */
177         }
178         
179         /* - if destination mesh had shapekeys, move them somewhere safe, and set up placeholders
180          *  with arrays that are large enough to hold shapekey data for all meshes
181          * -    if destination mesh didn't have shapekeys, but we encountered some in the meshes we're 
182          *      joining, set up a new keyblock and assign to the mesh
183          */
184         if (key) {
185                 /* make a duplicate copy that will only be used here... (must remember to free it!) */
186                 nkey = BKE_key_copy(key);
187                 
188                 /* for all keys in old block, clear data-arrays */
189                 for (kb = key->block.first; kb; kb = kb->next) {
190                         if (kb->data) MEM_freeN(kb->data);
191                         kb->data = MEM_callocN(sizeof(float) * 3 * totvert, "join_shapekey");
192                         kb->totelem = totvert;
193                         kb->weights = NULL;
194                 }
195         }
196         else if (haskey) {
197                 /* add a new key-block and add to the mesh */
198                 key = me->key = add_key((ID *)me);
199                 key->type = KEY_RELATIVE;
200         }
201         
202         /* first pass over objects - copying materials and vertexgroups across */
203         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
204         {
205                 /* only act if a mesh, and not the one we're joining to */
206                 if ((ob != base->object) && (base->object->type == OB_MESH)) {
207                         me = base->object->data;
208                         
209                         /* Join this object's vertex groups to the base one's */
210                         for (dg = base->object->defbase.first; dg; dg = dg->next) {
211                                 /* See if this group exists in the object (if it doesn't, add it to the end) */
212                                 if (!defgroup_find_name(ob, dg->name)) {
213                                         odg = MEM_callocN(sizeof(bDeformGroup), "join deformGroup");
214                                         memcpy(odg, dg, sizeof(bDeformGroup));
215                                         BLI_addtail(&ob->defbase, odg);
216                                 }
217                         }
218                         if (ob->defbase.first && ob->actdef == 0)
219                                 ob->actdef = 1;
220                         
221                         
222                         if (me->totvert) {
223                                 /* Add this object's materials to the base one's if they don't exist already (but only if limits not exceeded yet) */
224                                 if (totcol < MAXMAT) {
225                                         for (a = 1; a <= base->object->totcol; a++) {
226                                                 ma = give_current_material(base->object, a);
227
228                                                 for (b = 0; b < totcol; b++) {
229                                                         if (ma == matar[b]) break;
230                                                 }
231                                                 if (b == totcol) {
232                                                         matar[b] = ma;
233                                                         if (ma) {
234                                                                 id_us_plus(&ma->id);
235                                                         }
236                                                         totcol++;
237                                                 }
238                                                 if (totcol >= MAXMAT)
239                                                         break;
240                                         }
241                                 }
242                                 
243                                 /* if this mesh has shapekeys, check if destination mesh already has matching entries too */
244                                 if (me->key && key) {
245                                         for (kb = me->key->block.first; kb; kb = kb->next) {
246                                                 /* if key doesn't exist in destination mesh, add it */
247                                                 if (key_get_named_keyblock(key, kb->name) == NULL) {
248                                                         /* copy this existing one over to the new shapekey block */
249                                                         kbn = MEM_dupallocN(kb);
250                                                         kbn->prev = kbn->next = NULL;
251                                                         
252                                                         /* adjust settings to fit (allocate a new data-array) */
253                                                         kbn->data = MEM_callocN(sizeof(float) * 3 * totvert, "joined_shapekey");
254                                                         kbn->totelem = totvert;
255                                                         kbn->weights = NULL;
256                                                         
257                                                         okb = key->block.last;
258                                                         curpos = (okb) ? okb->pos : -0.1f;
259                                                         if (key->type == KEY_RELATIVE)
260                                                                 kbn->pos = curpos + 0.1f;
261                                                         else
262                                                                 kbn->pos = curpos;
263                                                         
264                                                         BLI_addtail(&key->block, kbn);
265                                                         key->totkey++;
266                                                         if (key->totkey == 1) key->refkey = kbn;
267                                                         
268                                                         // XXX 2.5 Animato
269 #if 0
270                                                         /* also, copy corresponding ipo-curve to ipo-block if applicable */
271                                                         if (me->key->ipo && key->ipo) {
272                                                                 // FIXME... this is a luxury item!
273                                                                 puts("FIXME: ignoring IPO's when joining shapekeys on Meshes for now...");
274                                                         }
275 #endif
276                                                 }
277                                         }
278                                 }
279                         }
280                 }
281         }
282         CTX_DATA_END;
283         
284         /* setup new data for destination mesh */
285         memset(&vdata, 0, sizeof(vdata));
286         memset(&edata, 0, sizeof(edata));
287         memset(&fdata, 0, sizeof(fdata));
288         memset(&ldata, 0, sizeof(ldata));
289         memset(&pdata, 0, sizeof(pdata));
290         
291         mvert = CustomData_add_layer(&vdata, CD_MVERT, CD_CALLOC, NULL, totvert);
292         medge = CustomData_add_layer(&edata, CD_MEDGE, CD_CALLOC, NULL, totedge);
293         mloop = CustomData_add_layer(&ldata, CD_MLOOP, CD_CALLOC, NULL, totloop);
294         mpoly = CustomData_add_layer(&pdata, CD_MPOLY, CD_CALLOC, NULL, totpoly);
295
296         vertofs = 0;
297         edgeofs = 0;
298         loopofs = 0;
299         polyofs = 0;
300         
301         /* inverse transform for all selected meshes in this object */
302         invert_m4_m4(imat, ob->obmat);
303         
304         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
305         {
306                 /* only join if this is a mesh */
307                 if (base->object->type == OB_MESH) {
308                         me = base->object->data;
309                         
310                         if (me->totvert) {
311                                 /* standard data */
312                                 CustomData_merge(&me->vdata, &vdata, CD_MASK_MESH, CD_DEFAULT, totvert);
313                                 CustomData_copy_data(&me->vdata, &vdata, 0, vertofs, me->totvert);
314                                 
315                                 /* vertex groups */
316                                 dvert = CustomData_get(&vdata, vertofs, CD_MDEFORMVERT);
317                                 
318                                 /* NB: vertex groups here are new version */
319                                 if (dvert) {
320                                         for (i = 0; i < me->totvert; i++) {
321                                                 for (j = 0; j < dvert[i].totweight; j++) {
322                                                         /*      Find the old vertex group */
323                                                         odg = BLI_findlink(&base->object->defbase, dvert[i].dw[j].def_nr);
324                                                         if (odg) {
325                                                                 /*      Search for a match in the new object, and set new index */
326                                                                 for (dg = ob->defbase.first, index = 0; dg; dg = dg->next, index++) {
327                                                                         if (!strcmp(dg->name, odg->name)) {
328                                                                                 dvert[i].dw[j].def_nr = index;
329                                                                                 break;
330                                                                         }
331                                                                 }
332                                                         }
333                                                 }
334                                         }
335                                 }
336                                 
337                                 /* if this is the object we're merging into, no need to do anything */
338                                 if (base->object != ob) {
339                                         /* watch this: switch matmul order really goes wrong */
340                                         mult_m4_m4m4(cmat, imat, base->object->obmat);
341                                         
342                                         /* transform vertex coordinates into new space */
343                                         for (a = 0, mv = mvert; a < me->totvert; a++, mv++) {
344                                                 mul_m4_v3(cmat, mv->co);
345                                         }
346                                         
347                                         /* for each shapekey in destination mesh:
348                                          *      - if there's a matching one, copy it across (will need to transform vertices into new space...)
349                                          *      - otherwise, just copy own coordinates of mesh (no need to transform vertex coordinates into new space)
350                                          */
351                                         if (key) {
352                                                 /* if this mesh has any shapekeys, check first, otherwise just copy coordinates */
353                                                 for (kb = key->block.first; kb; kb = kb->next) {
354                                                         /* get pointer to where to write data for this mesh in shapekey's data array */
355                                                         fp1 = ((float *)kb->data) + (vertofs * 3);
356                                                         
357                                                         /* check if this mesh has such a shapekey */
358                                                         okb = key_get_named_keyblock(me->key, kb->name);
359                                                         if (okb) {
360                                                                 /* copy this mesh's shapekey to the destination shapekey (need to transform first) */
361                                                                 fp2 = ((float *)(okb->data));
362                                                                 for (a = 0; a < me->totvert; a++, fp1 += 3, fp2 += 3) {
363                                                                         copy_v3_v3(fp1, fp2);
364                                                                         mul_m4_v3(cmat, fp1);
365                                                                 }
366                                                         }
367                                                         else {
368                                                                 /* copy this mesh's vertex coordinates to the destination shapekey */
369                                                                 mv = mvert;
370                                                                 for (a = 0; a < me->totvert; a++, fp1 += 3, mv++) {
371                                                                         copy_v3_v3(fp1, mv->co);
372                                                                 }
373                                                         }
374                                                 }
375                                         }
376                                 }
377                                 else {
378                                         /* for each shapekey in destination mesh:
379                                          *      - if it was an 'original', copy the appropriate data from nkey
380                                          *      - otherwise, copy across plain coordinates (no need to transform coordinates)
381                                          */
382                                         if (key) {
383                                                 for (kb = key->block.first; kb; kb = kb->next) {
384                                                         /* get pointer to where to write data for this mesh in shapekey's data array */
385                                                         fp1 = ((float *)kb->data) + (vertofs * 3);
386                                                         
387                                                         /* check if this was one of the original shapekeys */
388                                                         okb = key_get_named_keyblock(nkey, kb->name);
389                                                         if (okb) {
390                                                                 /* copy this mesh's shapekey to the destination shapekey */
391                                                                 fp2 = ((float *)(okb->data));
392                                                                 for (a = 0; a < me->totvert; a++, fp1 += 3, fp2 += 3) {
393                                                                         copy_v3_v3(fp1, fp2);
394                                                                 }
395                                                         }
396                                                         else {
397                                                                 /* copy base-coordinates to the destination shapekey */
398                                                                 mv = mvert;
399                                                                 for (a = 0; a < me->totvert; a++, fp1 += 3, mv++) {
400                                                                         copy_v3_v3(fp1, mv->co);
401                                                                 }
402                                                         }
403                                                 }
404                                         }
405                                 }
406                                 
407                                 /* advance mvert pointer to end of base mesh's data */
408                                 mvert += me->totvert;
409                         }
410                         
411                         if (me->totedge) {
412                                 CustomData_merge(&me->edata, &edata, CD_MASK_MESH, CD_DEFAULT, totedge);
413                                 CustomData_copy_data(&me->edata, &edata, 0, edgeofs, me->totedge);
414                                 
415                                 for (a = 0; a < me->totedge; a++, medge++) {
416                                         medge->v1 += vertofs;
417                                         medge->v2 += vertofs;
418                                 }
419                         }
420
421                         if (me->totloop) {
422                                 if (base->object != ob) {
423                                         MultiresModifierData *mmd;
424
425                                         multiresModifier_prepare_join(scene, base->object, ob);
426
427                                         if ((mmd = get_multires_modifier(scene, base->object, TRUE))) {
428                                                 ED_object_iter_other(bmain, base->object, TRUE,
429                                                                                          ED_object_multires_update_totlevels_cb,
430                                                                                          &mmd->totlvl);
431                                         }
432                                 }
433                                 
434                                 CustomData_merge(&me->ldata, &ldata, CD_MASK_MESH, CD_DEFAULT, totloop);
435                                 CustomData_copy_data(&me->ldata, &ldata, 0, loopofs, me->totloop);
436                                 
437                                 for (a = 0; a < me->totloop; a++, mloop++) {
438                                         mloop->v += vertofs;
439                                         mloop->e += edgeofs;
440                                 }
441                         }
442                         
443                         if (me->totpoly) {
444                                 /* make mapping for materials */
445                                 for (a = 1; a <= base->object->totcol; a++) {
446                                         ma = give_current_material(base->object, a);
447
448                                         for (b = 0; b < totcol; b++) {
449                                                 if (ma == matar[b]) {
450                                                         matmap[a - 1] = b;
451                                                         break;
452                                                 }
453                                         }
454                                 }
455                                 
456                                 CustomData_merge(&me->pdata, &pdata, CD_MASK_MESH, CD_DEFAULT, totpoly);
457                                 CustomData_copy_data(&me->pdata, &pdata, 0, polyofs, me->totpoly);
458                                 
459                                 for (a = 0; a < me->totpoly; a++, mpoly++) {
460                                         mpoly->loopstart += loopofs;
461                                         mpoly->mat_nr = matmap ? matmap[(int)mpoly->mat_nr] : 0;
462                                 }
463                                 
464                                 polyofs += me->totpoly;
465                         }
466
467                         /* these are used for relinking (cannot be set earlier, 
468                          * or else reattaching goes wrong)
469                          */
470                         vertofs += me->totvert;
471                         edgeofs += me->totedge;
472                         loopofs += me->totloop;
473                         
474                         /* free base, now that data is merged */
475                         if (base->object != ob)
476                                 ED_base_object_free_and_unlink(bmain, scene, base);
477                 }
478         }
479         CTX_DATA_END;
480         
481         /* return to mesh we're merging to */
482         me = ob->data;
483         
484         CustomData_free(&me->vdata, me->totvert);
485         CustomData_free(&me->edata, me->totedge);
486         CustomData_free(&me->ldata, me->totloop);
487         CustomData_free(&me->pdata, me->totpoly);
488
489         me->totvert = totvert;
490         me->totedge = totedge;
491         me->totloop = totloop;
492         me->totpoly = totpoly;
493
494         me->vdata = vdata;
495         me->edata = edata;
496         me->ldata = ldata;
497         me->pdata = pdata;
498
499         mesh_update_customdata_pointers(me, TRUE); /* BMESH_TODO, check if this arg can be failse, non urgent - campbell */
500         
501         /* old material array */
502         for (a = 1; a <= ob->totcol; a++) {
503                 ma = ob->mat[a - 1];
504                 if (ma) ma->id.us--;
505         }
506         for (a = 1; a <= me->totcol; a++) {
507                 ma = me->mat[a - 1];
508                 if (ma) ma->id.us--;
509         }
510         if (ob->mat) MEM_freeN(ob->mat);
511         if (ob->matbits) MEM_freeN(ob->matbits);
512         if (me->mat) MEM_freeN(me->mat);
513         ob->mat = me->mat = NULL;
514         ob->matbits = NULL;
515         
516         if (totcol) {
517                 me->mat = matar;
518                 ob->mat = MEM_callocN(sizeof(void *) * totcol, "join obmatar");
519                 ob->matbits = MEM_callocN(sizeof(char) * totcol, "join obmatbits");
520         }
521         else
522                 MEM_freeN(matar);
523         
524         ob->totcol = me->totcol = totcol;
525
526         if (matmap) MEM_freeN(matmap);
527         
528         /* other mesh users */
529         test_object_materials((ID *)me);
530         
531         /* free temp copy of destination shapekeys (if applicable) */
532         if (nkey) {
533                 // XXX 2.5 Animato
534 #if 0
535                 /* free it's ipo too - both are not actually freed from memory yet as ID-blocks */
536                 if (nkey->ipo) {
537                         BKE_ipo_free(nkey->ipo);
538                         BLI_remlink(&bmain->ipo, nkey->ipo);
539                         MEM_freeN(nkey->ipo);
540                 }
541 #endif
542                 
543                 BKE_key_free(nkey);
544                 BLI_remlink(&bmain->key, nkey);
545                 MEM_freeN(nkey);
546         }
547         
548         DAG_scene_sort(bmain, scene);   // removed objects, need to rebuild dag before editmode call
549
550 #if 0
551         ED_object_enter_editmode(C, EM_WAITCURSOR);
552         ED_object_exit_editmode(C, EM_FREEDATA | EM_WAITCURSOR | EM_DO_UNDO);
553 #else
554         /* toggle editmode using lower level functions so this can be called from python */
555         EDBM_mesh_make(scene->toolsettings, scene, ob);
556         EDBM_mesh_load(ob);
557         EDBM_mesh_free(me->edit_btmesh);
558         MEM_freeN(me->edit_btmesh);
559         me->edit_btmesh = NULL;
560         DAG_id_tag_update(&ob->id, OB_RECALC_OB | OB_RECALC_DATA);
561 #endif
562         WM_event_add_notifier(C, NC_SCENE | ND_OB_ACTIVE, scene);
563
564         return OPERATOR_FINISHED;
565 }
566
567 /*********************** JOIN AS SHAPES ***************************/
568
569 /* Append selected meshes vertex locations as shapes of the active mesh, 
570  * return 0 if no join is made (error) and 1 of the join is done */
571
572 int join_mesh_shapes_exec(bContext *C, wmOperator *op)
573 {
574         Scene *scene = CTX_data_scene(C);
575         Object *ob = CTX_data_active_object(C);
576         Mesh *me = (Mesh *)ob->data;
577         Mesh *selme = NULL;
578         DerivedMesh *dm = NULL;
579         Key *key = me->key;
580         KeyBlock *kb;
581         int ok = 0, nonequal_verts = 0;
582         
583         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
584         {
585                 if (base->object == ob) continue;
586                 
587                 if (base->object->type == OB_MESH) {
588                         selme = (Mesh *)base->object->data;
589                         
590                         if (selme->totvert == me->totvert)
591                                 ok++;
592                         else
593                                 nonequal_verts = 1;
594                 }
595         }
596         CTX_DATA_END;
597         
598         if (!ok) {
599                 if (nonequal_verts)
600                         BKE_report(op->reports, RPT_WARNING, "Selected meshes must have equal numbers of vertices");
601                 else
602                         BKE_report(op->reports, RPT_WARNING, "No additional selected meshes with equal vertex count to join");
603                 return OPERATOR_CANCELLED;
604         }
605         
606         if (key == NULL) {
607                 key = me->key = add_key((ID *)me);
608                 key->type = KEY_RELATIVE;
609
610                 /* first key added, so it was the basis. initialize it with the existing mesh */
611                 kb = add_keyblock(key, NULL);
612                 mesh_to_key(me, kb);
613         }
614         
615         /* now ready to add new keys from selected meshes */
616         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
617         {
618                 if (base->object == ob) continue;
619                 
620                 if (base->object->type == OB_MESH) {
621                         selme = (Mesh *)base->object->data;
622                         
623                         if (selme->totvert == me->totvert) {
624                                 dm = mesh_get_derived_deform(scene, base->object, CD_MASK_BAREMESH);
625                                 
626                                 if (!dm) continue;
627                                         
628                                 kb = add_keyblock(key, base->object->id.name + 2);
629                                 
630                                 DM_to_meshkey(dm, me, kb);
631                                 
632                                 dm->release(dm);
633                         }
634                 }
635         }
636         CTX_DATA_END;
637         
638         WM_event_add_notifier(C, NC_SCENE | ND_OB_ACTIVE, scene);
639         
640         return OPERATOR_FINISHED;
641 }
642
643 /* ********************* MESH VERTEX OCTREE LOOKUP ************* */
644
645 /* important note; this is unfinished, needs better API for editmode, and custom threshold */
646
647 #define MOC_RES         8
648 #define MOC_NODE_RES    8
649 #define MOC_THRESH      0.00002f
650
651 typedef struct MocNode {
652         struct MocNode *next;
653         intptr_t index[MOC_NODE_RES];
654 } MocNode;
655
656 static int mesh_octree_get_base_offs(const float co[3], const float offs[3], const float div[3])
657 {
658         int vx, vy, vz;
659         
660         vx = floor((co[0] - offs[0]) / div[0]);
661         vy = floor((co[1] - offs[1]) / div[1]);
662         vz = floor((co[2] - offs[2]) / div[2]);
663
664         CLAMP(vx, 0, MOC_RES - 1);
665         CLAMP(vy, 0, MOC_RES - 1);
666         CLAMP(vz, 0, MOC_RES - 1);
667
668         return (vx * MOC_RES * MOC_RES) + vy * MOC_RES + vz;
669 }
670
671 static void mesh_octree_add_node(MocNode **bt, intptr_t index)
672 {
673         if (*bt == NULL) {
674                 *bt = MEM_callocN(sizeof(MocNode), "MocNode");
675                 (*bt)->index[0] = index;
676         }
677         else {
678                 int a;
679                 for (a = 0; a < MOC_NODE_RES; a++) {
680                         if ((*bt)->index[a] == index)
681                                 return;
682                         else if ((*bt)->index[a] == 0) {
683                                 (*bt)->index[a] = index;
684                                 return;
685                         }
686                 }
687                 mesh_octree_add_node(&(*bt)->next, index);
688         }
689 }
690
691 static void mesh_octree_free_node(MocNode **bt)
692 {
693         if ( (*bt)->next) {
694                 mesh_octree_free_node(&(*bt)->next);
695         }
696         MEM_freeN(*bt);
697 }
698
699
700 /* temporal define, just to make nicer code below */
701 #define MOC_INDEX(vx, vy, vz)  (((vx) * MOC_RES * MOC_RES) + (vy) * MOC_RES + (vz))
702
703 static void mesh_octree_add_nodes(MocNode **basetable, const float co[3], const float offs[3],
704                                   const float div[3], intptr_t index)
705 {
706         float fx, fy, fz;
707         int vx, vy, vz;
708         
709         if (!finite(co[0]) ||
710             !finite(co[1]) ||
711             !finite(co[2]))
712         {
713                 return;
714         }
715         
716         fx = (co[0] - offs[0]) / div[0];
717         fy = (co[1] - offs[1]) / div[1];
718         fz = (co[2] - offs[2]) / div[2];
719         CLAMP(fx, 0.0f, MOC_RES - MOC_THRESH);
720         CLAMP(fy, 0.0f, MOC_RES - MOC_THRESH);
721         CLAMP(fz, 0.0f, MOC_RES - MOC_THRESH);
722
723         vx = (int)floorf(fx);
724         vy = (int)floorf(fy);
725         vz = (int)floorf(fz);
726
727         mesh_octree_add_node(basetable + MOC_INDEX(vx, vy, vz), index);
728
729         if (vx > 0)
730                 if (fx - ((float)vx) - MOC_THRESH < 0.0f)
731                         mesh_octree_add_node(basetable + MOC_INDEX(vx - 1, vy, vz), index);
732         if (vx < MOC_RES - 2)
733                 if (fx - ((float)vx) + MOC_THRESH > 1.0f)
734                         mesh_octree_add_node(basetable + MOC_INDEX(vx + 1, vy, vz), index);
735
736         if (vy > 0)
737                 if (fy - ((float)vy) - MOC_THRESH < 0.0f)
738                         mesh_octree_add_node(basetable + MOC_INDEX(vx, vy - 1, vz), index);
739         if (vy < MOC_RES - 2)
740                 if (fy - ((float)vy) + MOC_THRESH > 1.0f)
741                         mesh_octree_add_node(basetable + MOC_INDEX(vx, vy + 1, vz), index);
742
743         if (vz > 0)
744                 if (fz - ((float)vz) - MOC_THRESH < 0.0f)
745                         mesh_octree_add_node(basetable + MOC_INDEX(vx, vy, vz - 1), index);
746         if (vz < MOC_RES - 2)
747                 if (fz - ((float)vz) + MOC_THRESH > 1.0f)
748                         mesh_octree_add_node(basetable + MOC_INDEX(vx, vy, vz + 1), index);
749
750 }
751
752 static intptr_t mesh_octree_find_index(MocNode **bt, MVert *mvert, const float co[3])
753 {
754         float *vec;
755         int a;
756         
757         if (*bt == NULL)
758                 return -1;
759         
760         for (a = 0; a < MOC_NODE_RES; a++) {
761                 if ((*bt)->index[a]) {
762                         /* does mesh verts and editmode, code looks potential dangerous, octree should really be filled OK! */
763                         if (mvert) {
764                                 vec = (mvert + (*bt)->index[a] - 1)->co;
765                                 if (compare_v3v3(vec, co, MOC_THRESH))
766                                         return (*bt)->index[a] - 1;
767                         }
768                         else {
769                                 BMVert *eve = (BMVert *)((*bt)->index[a]);
770                                 if (compare_v3v3(eve->co, co, MOC_THRESH))
771                                         return (*bt)->index[a];
772                         }
773                 }
774                 else return -1;
775         }
776         if ( (*bt)->next)
777                 return mesh_octree_find_index(&(*bt)->next, mvert, co);
778         
779         return -1;
780 }
781
782 static struct {
783         MocNode **table;
784         float offs[3], div[3];
785 } MeshOctree = {NULL, {0, 0, 0}, {0, 0, 0}};
786
787 /* mode is 's' start, or 'e' end, or 'u' use */
788 /* if end, ob can be NULL */
789 intptr_t mesh_octree_table(Object *ob, BMEditMesh *em, const float co[3], char mode)
790 {
791         MocNode **bt;
792         
793         if (mode == 'u') {        /* use table */
794                 if (MeshOctree.table == NULL)
795                         mesh_octree_table(ob, em, NULL, 's');
796
797                 if (MeshOctree.table) {
798                         Mesh *me = ob->data;
799                         bt = MeshOctree.table + mesh_octree_get_base_offs(co, MeshOctree.offs, MeshOctree.div);
800                         if (em)
801                                 return mesh_octree_find_index(bt, NULL, co);
802                         else
803                                 return mesh_octree_find_index(bt, me->mvert, co);
804                 }
805                 return -1;
806         }
807         else if (mode == 's') {   /* start table */
808                 Mesh *me = ob->data;
809                 float min[3], max[3];
810
811                 /* we compute own bounding box and don't reuse ob->bb because
812                  * we are using the undeformed coordinates*/
813                 INIT_MINMAX(min, max);
814
815                 if (em && me->edit_btmesh == em) {
816                         BMIter iter;
817                         BMVert *eve;
818                         
819                         BM_ITER_MESH (eve, &iter, em->bm, BM_VERTS_OF_MESH) {
820                                 minmax_v3v3_v3(min, max, eve->co);
821                         }
822                 }
823                 else {          
824                         MVert *mvert;
825                         int a;
826                         
827                         for (a = 0, mvert = me->mvert; a < me->totvert; a++, mvert++)
828                                 minmax_v3v3_v3(min, max, mvert->co);
829                 }
830                 
831                 /* for quick unit coordinate calculus */
832                 copy_v3_v3(MeshOctree.offs, min);
833                 /* we offset it 1 threshold unit extra */
834                 add_v3_fl(MeshOctree.offs, -MOC_THRESH);
835                 
836                 sub_v3_v3v3(MeshOctree.div, max, min);
837                 /* and divide with 2 threshold unit more extra (try 8x8 unit grid on paint) */
838                 add_v3_fl(MeshOctree.div, 2.0f * MOC_THRESH);
839
840                 mul_v3_fl(MeshOctree.div, 1.0f / MOC_RES);
841                 if (MeshOctree.div[0] == 0.0f) MeshOctree.div[0] = 1.0f;
842                 if (MeshOctree.div[1] == 0.0f) MeshOctree.div[1] = 1.0f;
843                 if (MeshOctree.div[2] == 0.0f) MeshOctree.div[2] = 1.0f;
844                         
845                 if (MeshOctree.table) /* happens when entering this call without ending it */
846                         mesh_octree_table(ob, em, co, 'e');
847                 
848                 MeshOctree.table = MEM_callocN(MOC_RES * MOC_RES * MOC_RES * sizeof(void *), "sym table");
849                 
850                 if (em && me->edit_btmesh == em) {
851                         BMVert *eve;
852                         BMIter iter;
853
854                         BM_ITER_MESH (eve, &iter, em->bm, BM_VERTS_OF_MESH) {
855                                 mesh_octree_add_nodes(MeshOctree.table, eve->co, MeshOctree.offs, MeshOctree.div, (intptr_t)(eve));
856                         }
857                 }
858                 else {          
859                         MVert *mvert;
860                         int a;
861                         
862                         for (a = 0, mvert = me->mvert; a < me->totvert; a++, mvert++)
863                                 mesh_octree_add_nodes(MeshOctree.table, mvert->co, MeshOctree.offs, MeshOctree.div, a + 1);
864                 }
865         }
866         else if (mode == 'e') { /* end table */
867                 if (MeshOctree.table) {
868                         int a;
869                         
870                         for (a = 0, bt = MeshOctree.table; a < MOC_RES * MOC_RES * MOC_RES; a++, bt++) {
871                                 if (*bt) mesh_octree_free_node(bt);
872                         }
873                         MEM_freeN(MeshOctree.table);
874                         MeshOctree.table = NULL;
875                 }
876         }
877         return 0;
878 }
879
880 MirrTopoStore_t mesh_topo_store = {NULL, -1. - 1, -1};
881
882 /* mode is 's' start, or 'e' end, or 'u' use */
883 /* if end, ob can be NULL */
884 /* note, is supposed return -1 on error, which callers are currently checking for, but is not used so far */
885 int mesh_mirrtopo_table(Object *ob, char mode)
886 {
887         if (mode == 'u') {        /* use table */
888                 if (ED_mesh_mirrtopo_recalc_check(ob->data, ob->mode, &mesh_topo_store)) {
889                         mesh_mirrtopo_table(ob, 's');
890                 }
891         }
892         else if (mode == 's') { /* start table */
893                 ED_mesh_mirrtopo_init(ob->data, ob->mode, &mesh_topo_store, FALSE);
894         }
895         else if (mode == 'e') { /* end table */
896                 ED_mesh_mirrtopo_free(&mesh_topo_store);
897         }
898         return 0;
899 }
900
901 static int mesh_get_x_mirror_vert_spatial(Object *ob, int index)
902 {
903         Mesh *me = ob->data;
904         MVert *mvert;
905         float vec[3];
906         
907         mvert = me->mvert + index;
908         vec[0] = -mvert->co[0];
909         vec[1] = mvert->co[1];
910         vec[2] = mvert->co[2];
911         
912         return mesh_octree_table(ob, NULL, vec, 'u');
913 }
914
915 static int mesh_get_x_mirror_vert_topo(Object *ob, int index)
916 {
917         if (mesh_mirrtopo_table(ob, 'u') == -1)
918                 return -1;
919
920         return mesh_topo_store.index_lookup[index];
921 }
922
923 int mesh_get_x_mirror_vert(Object *ob, int index)
924 {
925         if (((Mesh *)ob->data)->editflag & ME_EDIT_MIRROR_TOPO) {
926                 return mesh_get_x_mirror_vert_topo(ob, index);
927         }
928         else {
929                 return mesh_get_x_mirror_vert_spatial(ob, index);
930         }
931         return 0;
932 }
933
934 static BMVert *editbmesh_get_x_mirror_vert_spatial(Object *ob, BMEditMesh *em, const float co[3])
935 {
936         float vec[3];
937         intptr_t poinval;
938         
939         /* ignore nan verts */
940         if (!finite(co[0]) ||
941             !finite(co[1]) ||
942             !finite(co[2]))
943         {
944                 return NULL;
945         }
946         
947         vec[0] = -co[0];
948         vec[1] = co[1];
949         vec[2] = co[2];
950         
951         poinval = mesh_octree_table(ob, em, vec, 'u');
952         if (poinval != -1)
953                 return (BMVert *)(poinval);
954         return NULL;
955 }
956
957 static BMVert *editbmesh_get_x_mirror_vert_topo(Object *ob, struct BMEditMesh *em, BMVert *eve, int index)
958 {
959         intptr_t poinval;
960         if (mesh_mirrtopo_table(ob, 'u') == -1)
961                 return NULL;
962
963         if (index == -1) {
964                 BMIter iter;
965                 BMVert *v;
966                 
967                 index = 0;
968                 BM_ITER_MESH (v, &iter, em->bm, BM_VERTS_OF_MESH) {
969                         if (v == eve)
970                                 break;
971                         index++;
972                 }
973                 
974                 if (index == em->bm->totvert) {
975                         return NULL;
976                 }
977         }
978
979         poinval = mesh_topo_store.index_lookup[index];
980
981         if (poinval != -1)
982                 return (BMVert *)(poinval);
983         return NULL;
984 }       
985
986 BMVert *editbmesh_get_x_mirror_vert(Object *ob, struct BMEditMesh *em, BMVert *eve, const float co[3], int index)
987 {
988         if (((Mesh *)ob->data)->editflag & ME_EDIT_MIRROR_TOPO) {
989                 return editbmesh_get_x_mirror_vert_topo(ob, em, eve, index);
990         }
991         else {
992                 return editbmesh_get_x_mirror_vert_spatial(ob, em, co);
993         }
994 }
995
996 #if 0
997
998 static float *editmesh_get_mirror_uv(BMEditMesh *em, int axis, float *uv, float *mirrCent, float *face_cent)
999 {
1000         float vec[2];
1001         float cent_vec[2];
1002         float cent[2];
1003
1004         /* ignore nan verts */
1005         if (isnan(uv[0]) || !finite(uv[0]) ||
1006             isnan(uv[1]) || !finite(uv[1])
1007             )
1008                 return NULL;
1009
1010         if (axis) {
1011                 vec[0] = uv[0];
1012                 vec[1] = -((uv[1]) - mirrCent[1]) + mirrCent[1];
1013
1014                 cent_vec[0] = face_cent[0];
1015                 cent_vec[1] = -((face_cent[1]) - mirrCent[1]) + mirrCent[1];
1016         }
1017         else {
1018                 vec[0] = -((uv[0]) - mirrCent[0]) + mirrCent[0];
1019                 vec[1] = uv[1];
1020
1021                 cent_vec[0] = -((face_cent[0]) - mirrCent[0]) + mirrCent[0];
1022                 cent_vec[1] = face_cent[1];
1023         }
1024
1025         /* TODO - Optimize */
1026         {
1027                 BMIter iter;
1028                 BMFace *efa;
1029                 
1030                 BM_ITER_MESH (efa, &iter, em->bm, BM_FACES_OF_MESH) {
1031                         uv_poly_center(em, efa, cent);
1032                         
1033                         if ( (fabs(cent[0] - cent_vec[0]) < 0.001) && (fabs(cent[1] - cent_vec[1]) < 0.001) ) {
1034                                 BMIter liter;
1035                                 BMLoop *l;
1036                                 
1037                                 BM_ITER_ELEM (l, &liter, efa, BM_LOOPS_OF_FACE) {
1038                                         MLoopUV *luv = CustomData_bmesh_get(&em->bm->ldata, l->head.data, CD_MLOOPUV);
1039                                         if ( (fabs(luv->uv[0] - vec[0]) < 0.001) && (fabs(luv->uv[1] - vec[1]) < 0.001) ) {
1040                                                 return luv->uv;
1041                                                                 
1042                                         }
1043                                 }
1044                         }
1045                 }
1046         }
1047
1048         return NULL;
1049 }
1050
1051 #endif
1052
1053 static unsigned int mirror_facehash(const void *ptr)
1054 {
1055         const MFace *mf = ptr;
1056         int v0, v1;
1057
1058         if (mf->v4) {
1059                 v0 = MIN4(mf->v1, mf->v2, mf->v3, mf->v4);
1060                 v1 = MAX4(mf->v1, mf->v2, mf->v3, mf->v4);
1061         }
1062         else {
1063                 v0 = MIN3(mf->v1, mf->v2, mf->v3);
1064                 v1 = MAX3(mf->v1, mf->v2, mf->v3);
1065         }
1066
1067         return ((v0 * 39) ^ (v1 * 31));
1068 }
1069
1070 static int mirror_facerotation(MFace *a, MFace *b)
1071 {
1072         if (b->v4) {
1073                 if (a->v1 == b->v1 && a->v2 == b->v2 && a->v3 == b->v3 && a->v4 == b->v4)
1074                         return 0;
1075                 else if (a->v4 == b->v1 && a->v1 == b->v2 && a->v2 == b->v3 && a->v3 == b->v4)
1076                         return 1;
1077                 else if (a->v3 == b->v1 && a->v4 == b->v2 && a->v1 == b->v3 && a->v2 == b->v4)
1078                         return 2;
1079                 else if (a->v2 == b->v1 && a->v3 == b->v2 && a->v4 == b->v3 && a->v1 == b->v4)
1080                         return 3;
1081         }
1082         else {
1083                 if (a->v1 == b->v1 && a->v2 == b->v2 && a->v3 == b->v3)
1084                         return 0;
1085                 else if (a->v3 == b->v1 && a->v1 == b->v2 && a->v2 == b->v3)
1086                         return 1;
1087                 else if (a->v2 == b->v1 && a->v3 == b->v2 && a->v1 == b->v3)
1088                         return 2;
1089         }
1090         
1091         return -1;
1092 }
1093
1094 static int mirror_facecmp(const void *a, const void *b)
1095 {
1096         return (mirror_facerotation((MFace *)a, (MFace *)b) == -1);
1097 }
1098
1099 /* BMESH_TODO, convert to MPoly (functions above also) */
1100 int *mesh_get_x_mirror_faces(Object *ob, BMEditMesh *em)
1101 {
1102         Mesh *me = ob->data;
1103         MVert *mv, *mvert = me->mvert;
1104         MFace mirrormf, *mf, *hashmf, *mface = me->mface;
1105         GHash *fhash;
1106         int *mirrorverts, *mirrorfaces;
1107         int a;
1108
1109         mirrorverts = MEM_callocN(sizeof(int) * me->totvert, "MirrorVerts");
1110         mirrorfaces = MEM_callocN(sizeof(int) * 2 * me->totface, "MirrorFaces");
1111
1112         mesh_octree_table(ob, em, NULL, 's');
1113
1114         for (a = 0, mv = mvert; a < me->totvert; a++, mv++)
1115                 mirrorverts[a] = mesh_get_x_mirror_vert(ob, a);
1116
1117         mesh_octree_table(ob, em, NULL, 'e');
1118
1119         fhash = BLI_ghash_new(mirror_facehash, mirror_facecmp, "mirror_facehash gh");
1120         for (a = 0, mf = mface; a < me->totface; a++, mf++)
1121                 BLI_ghash_insert(fhash, mf, mf);
1122
1123         for (a = 0, mf = mface; a < me->totface; a++, mf++) {
1124                 mirrormf.v1 = mirrorverts[mf->v3];
1125                 mirrormf.v2 = mirrorverts[mf->v2];
1126                 mirrormf.v3 = mirrorverts[mf->v1];
1127                 mirrormf.v4 = (mf->v4) ? mirrorverts[mf->v4] : 0;
1128
1129                 /* make sure v4 is not 0 if a quad */
1130                 if (mf->v4 && mirrormf.v4 == 0) {
1131                         SWAP(unsigned int, mirrormf.v1, mirrormf.v3);
1132                         SWAP(unsigned int, mirrormf.v2, mirrormf.v4);
1133                 }
1134
1135                 hashmf = BLI_ghash_lookup(fhash, &mirrormf);
1136                 if (hashmf) {
1137                         mirrorfaces[a * 2] = hashmf - mface;
1138                         mirrorfaces[a * 2 + 1] = mirror_facerotation(&mirrormf, hashmf);
1139                 }
1140                 else
1141                         mirrorfaces[a * 2] = -1;
1142         }
1143
1144         BLI_ghash_free(fhash, NULL, NULL);
1145         MEM_freeN(mirrorverts);
1146         
1147         return mirrorfaces;
1148 }