Fix T47038: Particles in Particle Edit Mode get added in completely wrong location.
[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  * meshtools.c: no editmode (violated already :), mirror & join),
32  * tools operating on meshes
33  */
34
35 #include "MEM_guardedalloc.h"
36
37 #include "DNA_mesh_types.h"
38 #include "DNA_key_types.h"
39 #include "DNA_material_types.h"
40 #include "DNA_modifier_types.h"
41 #include "DNA_object_types.h"
42 #include "DNA_scene_types.h"
43 #include "DNA_screen_types.h"
44 #include "DNA_view3d_types.h"
45
46 #include "BLI_math.h"
47 #include "BLI_blenlib.h"
48
49
50 #include "BLI_kdtree.h"
51 #include "BKE_context.h"
52 #include "BKE_depsgraph.h"
53 #include "BKE_deform.h"
54 #include "BKE_DerivedMesh.h"
55 #include "BKE_key.h"
56 #include "BKE_library.h"
57 #include "BKE_main.h"
58 #include "BKE_mesh.h"
59 #include "BKE_material.h"
60 #include "BKE_report.h"
61 #include "BKE_editmesh.h"
62 #include "BKE_multires.h"
63
64 #include "ED_mesh.h"
65 #include "ED_object.h"
66 #include "ED_view3d.h"
67
68 #include "WM_api.h"
69 #include "WM_types.h"
70
71 /* * ********************** no editmode!!! *********** */
72
73 /*********************** JOIN ***************************/
74
75 /* join selected meshes into the active mesh, context sensitive
76  * return 0 if no join is made (error) and 1 if the join is done */
77
78 int join_mesh_exec(bContext *C, wmOperator *op)
79 {
80         Main *bmain = CTX_data_main(C);
81         Scene *scene = CTX_data_scene(C);
82         Object *ob = CTX_data_active_object(C);
83         Material **matar, *ma;
84         Mesh *me;
85         MVert *mvert, *mv;
86         MEdge *medge = NULL;
87         MPoly *mpoly = NULL;
88         MLoop *mloop = NULL;
89         Key *key, *nkey = NULL;
90         KeyBlock *kb, *okb, *kbn;
91         float imat[4][4], cmat[4][4], *fp1, *fp2;
92         int a, b, totcol, totmat = 0, totedge = 0, totvert = 0;
93         int totloop = 0, totpoly = 0, vertofs, *matmap = NULL;
94         int i, j, index, haskey = 0, edgeofs, loopofs, polyofs;
95         bool ok = false;
96         bDeformGroup *dg, *odg;
97         MDeformVert *dvert;
98         CustomData vdata, edata, fdata, ldata, pdata;
99
100         if (scene->obedit) {
101                 BKE_report(op->reports, RPT_WARNING, "Cannot join while in edit mode");
102                 return OPERATOR_CANCELLED;
103         }
104         
105         /* ob is the object we are adding geometry to */
106         if (!ob || ob->type != OB_MESH) {
107                 BKE_report(op->reports, RPT_WARNING, "Active object is not a mesh");
108                 return OPERATOR_CANCELLED;
109         }
110         
111         /* count & check */
112         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
113         {
114                 if (base->object->type == OB_MESH) {
115                         me = base->object->data;
116
117                         totvert += me->totvert;
118                         totedge += me->totedge;
119                         totloop += me->totloop;
120                         totpoly += me->totpoly;
121                         totmat += base->object->totcol;
122                         
123                         if (base->object == ob)
124                                 ok = true;
125                         
126                         /* check for shapekeys */
127                         if (me->key)
128                                 haskey++;
129                 }
130         }
131         CTX_DATA_END;
132         
133         /* that way the active object is always selected */ 
134         if (ok == false) {
135                 BKE_report(op->reports, RPT_WARNING, "Active object is not a selected mesh");
136                 return OPERATOR_CANCELLED;
137         }
138         
139         /* only join meshes if there are verts to join, there aren't too many, and we only had one mesh selected */
140         me = (Mesh *)ob->data;
141         key = me->key;
142
143         if (totvert == 0 || totvert == me->totvert) {
144                 BKE_report(op->reports, RPT_WARNING, "No mesh data to join");
145                 return OPERATOR_CANCELLED;
146         }
147         
148         if (totvert > MESH_MAX_VERTS) {
149                 BKE_reportf(op->reports, RPT_WARNING, "Joining results in %d vertices, limit is %ld", totvert, MESH_MAX_VERTS);
150                 return OPERATOR_CANCELLED;
151         }
152
153         /* remove tessface to ensure we don't hold references to invalid faces */
154         BKE_mesh_tessface_clear(me);
155
156         /* new material indices and material array */
157         matar = MEM_callocN(sizeof(void *) * totmat, "join_mesh matar");
158         if (totmat) matmap = MEM_callocN(sizeof(int) * totmat, "join_mesh matmap");
159         totcol = ob->totcol;
160         
161         /* obact materials in new main array, is nicer start! */
162         for (a = 0; a < ob->totcol; a++) {
163                 matar[a] = give_current_material(ob, a + 1);
164                 id_us_plus((ID *)matar[a]);
165                 /* increase id->us : will be lowered later */
166         }
167         
168         /* - if destination mesh had shapekeys, move them somewhere safe, and set up placeholders
169          *  with arrays that are large enough to hold shapekey data for all meshes
170          * -    if destination mesh didn't have shapekeys, but we encountered some in the meshes we're 
171          *      joining, set up a new keyblock and assign to the mesh
172          */
173         if (key) {
174                 /* make a duplicate copy that will only be used here... (must remember to free it!) */
175                 nkey = BKE_key_copy(key);
176                 
177                 /* for all keys in old block, clear data-arrays */
178                 for (kb = key->block.first; kb; kb = kb->next) {
179                         if (kb->data) MEM_freeN(kb->data);
180                         kb->data = MEM_callocN(sizeof(float) * 3 * totvert, "join_shapekey");
181                         kb->totelem = totvert;
182                 }
183         }
184         else if (haskey) {
185                 /* add a new key-block and add to the mesh */
186                 key = me->key = BKE_key_add((ID *)me);
187                 key->type = KEY_RELATIVE;
188         }
189         
190         /* first pass over objects - copying materials and vertexgroups across */
191         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
192         {
193                 /* only act if a mesh, and not the one we're joining to */
194                 if ((ob != base->object) && (base->object->type == OB_MESH)) {
195                         me = base->object->data;
196                         
197                         /* Join this object's vertex groups to the base one's */
198                         for (dg = base->object->defbase.first; dg; dg = dg->next) {
199                                 /* See if this group exists in the object (if it doesn't, add it to the end) */
200                                 if (!defgroup_find_name(ob, dg->name)) {
201                                         odg = MEM_callocN(sizeof(bDeformGroup), "join deformGroup");
202                                         memcpy(odg, dg, sizeof(bDeformGroup));
203                                         BLI_addtail(&ob->defbase, odg);
204                                 }
205                         }
206                         if (ob->defbase.first && ob->actdef == 0)
207                                 ob->actdef = 1;
208                         
209                         
210                         if (me->totvert) {
211                                 /* Add this object's materials to the base one's if they don't exist already (but only if limits not exceeded yet) */
212                                 if (totcol < MAXMAT) {
213                                         for (a = 1; a <= base->object->totcol; a++) {
214                                                 ma = give_current_material(base->object, a);
215
216                                                 for (b = 0; b < totcol; b++) {
217                                                         if (ma == matar[b]) break;
218                                                 }
219                                                 if (b == totcol) {
220                                                         matar[b] = ma;
221                                                         if (ma) {
222                                                                 id_us_plus(&ma->id);
223                                                         }
224                                                         totcol++;
225                                                 }
226                                                 if (totcol >= MAXMAT)
227                                                         break;
228                                         }
229                                 }
230                                 
231                                 /* if this mesh has shapekeys, check if destination mesh already has matching entries too */
232                                 if (me->key && key) {
233                                         /* for remapping KeyBlock.relative */
234                                         int      *index_map = MEM_mallocN(sizeof(int)        * me->key->totkey, __func__);
235                                         KeyBlock **kb_map   = MEM_mallocN(sizeof(KeyBlock *) * me->key->totkey, __func__);
236
237                                         for (kb = me->key->block.first, i = 0; kb; kb = kb->next, i++) {
238                                                 BLI_assert(i < me->key->totkey);
239
240                                                 kbn = BKE_keyblock_find_name(key, kb->name);
241                                                 /* if key doesn't exist in destination mesh, add it */
242                                                 if (kbn) {
243                                                         index_map[i] = BLI_findindex(&key->block, kbn);
244                                                 }
245                                                 else {
246                                                         index_map[i] = key->totkey;
247
248                                                         kbn = BKE_keyblock_add(key, kb->name);
249
250                                                         BKE_keyblock_copy_settings(kbn, kb);
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                 
256                                                         /* XXX 2.5 Animato */
257 #if 0
258                                                         /* also, copy corresponding ipo-curve to ipo-block if applicable */
259                                                         if (me->key->ipo && key->ipo) {
260                                                                 /* FIXME... this is a luxury item! */
261                                                                 puts("FIXME: ignoring IPO's when joining shapekeys on Meshes for now...");
262                                                         }
263 #endif
264                                                 }
265
266                                                 kb_map[i] = kbn;
267                                         }
268
269                                         /* remap relative index values */
270                                         for (kb = me->key->block.first, i = 0; kb; kb = kb->next, i++) {
271                                                 if (LIKELY(kb->relative < me->key->totkey)) {  /* sanity check, should always be true */
272                                                         kb_map[i]->relative = index_map[kb->relative];
273                                                 }
274                                         }
275
276                                         MEM_freeN(index_map);
277                                         MEM_freeN(kb_map);
278                                 }
279                         }
280                 }
281         }
282         CTX_DATA_END;
283
284
285         /* setup new data for destination mesh */
286         CustomData_reset(&vdata);
287         CustomData_reset(&edata);
288         CustomData_reset(&fdata);
289         CustomData_reset(&ldata);
290         CustomData_reset(&pdata);
291
292         mvert = CustomData_add_layer(&vdata, CD_MVERT, CD_CALLOC, NULL, totvert);
293         medge = CustomData_add_layer(&edata, CD_MEDGE, CD_CALLOC, NULL, totedge);
294         mloop = CustomData_add_layer(&ldata, CD_MLOOP, CD_CALLOC, NULL, totloop);
295         mpoly = CustomData_add_layer(&pdata, CD_MPOLY, CD_CALLOC, NULL, totpoly);
296
297         vertofs = 0;
298         edgeofs = 0;
299         loopofs = 0;
300         polyofs = 0;
301         
302         /* inverse transform for all selected meshes in this object */
303         invert_m4_m4(imat, ob->obmat);
304         
305         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
306         {
307                 /* only join if this is a mesh */
308                 if (base->object->type == OB_MESH) {
309                         me = base->object->data;
310                         
311                         if (me->totvert) {
312
313                                 /* merge customdata flag */
314                                 ((Mesh *)ob->data)->cd_flag |= me->cd_flag;
315
316                                 /* standard data */
317                                 CustomData_merge(&me->vdata, &vdata, CD_MASK_MESH, CD_DEFAULT, totvert);
318                                 CustomData_copy_data_named(&me->vdata, &vdata, 0, vertofs, me->totvert);
319                                 
320                                 /* vertex groups */
321                                 dvert = CustomData_get(&vdata, vertofs, CD_MDEFORMVERT);
322                                 
323                                 /* NB: vertex groups here are new version */
324                                 if (dvert) {
325                                         for (i = 0; i < me->totvert; i++) {
326                                                 for (j = 0; j < dvert[i].totweight; j++) {
327                                                         /*      Find the old vertex group */
328                                                         odg = BLI_findlink(&base->object->defbase, dvert[i].dw[j].def_nr);
329                                                         if (odg) {
330                                                                 /*      Search for a match in the new object, and set new index */
331                                                                 for (dg = ob->defbase.first, index = 0; dg; dg = dg->next, index++) {
332                                                                         if (STREQ(dg->name, odg->name)) {
333                                                                                 dvert[i].dw[j].def_nr = index;
334                                                                                 break;
335                                                                         }
336                                                                 }
337                                                         }
338                                                 }
339                                         }
340                                 }
341                                 
342                                 /* if this is the object we're merging into, no need to do anything */
343                                 if (base->object != ob) {
344                                         /* watch this: switch matmul order really goes wrong */
345                                         mul_m4_m4m4(cmat, imat, base->object->obmat);
346                                         
347                                         /* transform vertex coordinates into new space */
348                                         for (a = 0, mv = mvert; a < me->totvert; a++, mv++) {
349                                                 mul_m4_v3(cmat, mv->co);
350                                         }
351                                         
352                                         /* for each shapekey in destination mesh:
353                                          *      - if there's a matching one, copy it across (will need to transform vertices into new space...)
354                                          *      - otherwise, just copy own coordinates of mesh (no need to transform vertex coordinates into new space)
355                                          */
356                                         if (key) {
357                                                 /* if this mesh has any shapekeys, check first, otherwise just copy coordinates */
358                                                 for (kb = key->block.first; kb; kb = kb->next) {
359                                                         /* get pointer to where to write data for this mesh in shapekey's data array */
360                                                         fp1 = ((float *)kb->data) + (vertofs * 3);
361                                                         
362                                                         /* check if this mesh has such a shapekey */
363                                                         okb = me->key ? BKE_keyblock_find_name(me->key, kb->name) : NULL;
364
365                                                         if (okb) {
366                                                                 /* copy this mesh's shapekey to the destination shapekey (need to transform first) */
367                                                                 fp2 = ((float *)(okb->data));
368                                                                 for (a = 0; a < me->totvert; a++, fp1 += 3, fp2 += 3) {
369                                                                         copy_v3_v3(fp1, fp2);
370                                                                         mul_m4_v3(cmat, fp1);
371                                                                 }
372                                                         }
373                                                         else {
374                                                                 /* copy this mesh's vertex coordinates to the destination shapekey */
375                                                                 mv = mvert;
376                                                                 for (a = 0; a < me->totvert; a++, fp1 += 3, mv++) {
377                                                                         copy_v3_v3(fp1, mv->co);
378                                                                 }
379                                                         }
380                                                 }
381                                         }
382                                 }
383                                 else {
384                                         /* for each shapekey in destination mesh:
385                                          *      - if it was an 'original', copy the appropriate data from nkey
386                                          *      - otherwise, copy across plain coordinates (no need to transform coordinates)
387                                          */
388                                         if (key) {
389                                                 for (kb = key->block.first; kb; kb = kb->next) {
390                                                         /* get pointer to where to write data for this mesh in shapekey's data array */
391                                                         fp1 = ((float *)kb->data) + (vertofs * 3);
392                                                         
393                                                         /* check if this was one of the original shapekeys */
394                                                         okb = nkey ? BKE_keyblock_find_name(nkey, kb->name) : NULL;
395                                                         if (okb) {
396                                                                 /* copy this mesh's shapekey to the destination shapekey */
397                                                                 fp2 = ((float *)(okb->data));
398                                                                 for (a = 0; a < me->totvert; a++, fp1 += 3, fp2 += 3) {
399                                                                         copy_v3_v3(fp1, fp2);
400                                                                 }
401                                                         }
402                                                         else {
403                                                                 /* copy base-coordinates to the destination shapekey */
404                                                                 mv = mvert;
405                                                                 for (a = 0; a < me->totvert; a++, fp1 += 3, mv++) {
406                                                                         copy_v3_v3(fp1, mv->co);
407                                                                 }
408                                                         }
409                                                 }
410                                         }
411                                 }
412                                 
413                                 /* advance mvert pointer to end of base mesh's data */
414                                 mvert += me->totvert;
415                         }
416                         
417                         if (me->totedge) {
418                                 CustomData_merge(&me->edata, &edata, CD_MASK_MESH, CD_DEFAULT, totedge);
419                                 CustomData_copy_data_named(&me->edata, &edata, 0, edgeofs, me->totedge);
420                                 
421                                 for (a = 0; a < me->totedge; a++, medge++) {
422                                         medge->v1 += vertofs;
423                                         medge->v2 += vertofs;
424                                 }
425                         }
426
427                         if (me->totloop) {
428                                 if (base->object != ob) {
429                                         MultiresModifierData *mmd;
430
431                                         multiresModifier_prepare_join(scene, base->object, ob);
432
433                                         if ((mmd = get_multires_modifier(scene, base->object, true))) {
434                                                 ED_object_iter_other(bmain, base->object, true,
435                                                                      ED_object_multires_update_totlevels_cb,
436                                                                      &mmd->totlvl);
437                                         }
438                                 }
439                                 
440                                 CustomData_merge(&me->ldata, &ldata, CD_MASK_MESH, CD_DEFAULT, totloop);
441                                 CustomData_copy_data_named(&me->ldata, &ldata, 0, loopofs, me->totloop);
442                                 
443                                 for (a = 0; a < me->totloop; a++, mloop++) {
444                                         mloop->v += vertofs;
445                                         mloop->e += edgeofs;
446                                 }
447                         }
448                         
449                         if (me->totpoly) {
450                                 if (totmat) {
451                                         /* make mapping for materials */
452                                         for (a = 1; a <= base->object->totcol; a++) {
453                                                 ma = give_current_material(base->object, a);
454
455                                                 for (b = 0; b < totcol; b++) {
456                                                         if (ma == matar[b]) {
457                                                                 matmap[a - 1] = b;
458                                                                 break;
459                                                         }
460                                                 }
461                                         }
462                                 }
463
464                                 CustomData_merge(&me->pdata, &pdata, CD_MASK_MESH, CD_DEFAULT, totpoly);
465                                 CustomData_copy_data_named(&me->pdata, &pdata, 0, polyofs, me->totpoly);
466                                 
467                                 for (a = 0; a < me->totpoly; a++, mpoly++) {
468                                         mpoly->loopstart += loopofs;
469                                         mpoly->mat_nr = matmap ? matmap[(int)mpoly->mat_nr] : 0;
470                                 }
471                                 
472                                 polyofs += me->totpoly;
473                         }
474
475                         /* these are used for relinking (cannot be set earlier, 
476                          * or else reattaching goes wrong)
477                          */
478                         vertofs += me->totvert;
479                         edgeofs += me->totedge;
480                         loopofs += me->totloop;
481                         
482                         /* free base, now that data is merged */
483                         if (base->object != ob)
484                                 ED_base_object_free_and_unlink(bmain, scene, base);
485                 }
486         }
487         CTX_DATA_END;
488         
489         /* return to mesh we're merging to */
490         me = ob->data;
491         
492         CustomData_free(&me->vdata, me->totvert);
493         CustomData_free(&me->edata, me->totedge);
494         CustomData_free(&me->ldata, me->totloop);
495         CustomData_free(&me->pdata, me->totpoly);
496
497         me->totvert = totvert;
498         me->totedge = totedge;
499         me->totloop = totloop;
500         me->totpoly = totpoly;
501
502         me->vdata = vdata;
503         me->edata = edata;
504         me->ldata = ldata;
505         me->pdata = pdata;
506
507         /* tessface data removed above, no need to update */
508         BKE_mesh_update_customdata_pointers(me, false);
509
510         /* update normals in case objects with non-uniform scale are joined */
511         BKE_mesh_calc_normals(me);
512         
513         /* old material array */
514         for (a = 1; a <= ob->totcol; a++) {
515                 ma = ob->mat[a - 1];
516                 if (ma)
517                         id_us_min(&ma->id);
518         }
519         for (a = 1; a <= me->totcol; a++) {
520                 ma = me->mat[a - 1];
521                 if (ma)
522                         id_us_min(&ma->id);
523         }
524         if (ob->mat) MEM_freeN(ob->mat);
525         if (ob->matbits) MEM_freeN(ob->matbits);
526         if (me->mat) MEM_freeN(me->mat);
527         ob->mat = me->mat = NULL;
528         ob->matbits = NULL;
529         
530         if (totcol) {
531                 me->mat = matar;
532                 ob->mat = MEM_callocN(sizeof(void *) * totcol, "join obmatar");
533                 ob->matbits = MEM_callocN(sizeof(char) * totcol, "join obmatbits");
534         }
535         else
536                 MEM_freeN(matar);
537         
538         ob->totcol = me->totcol = totcol;
539
540         if (matmap) MEM_freeN(matmap);
541         
542         /* other mesh users */
543         test_object_materials(bmain, (ID *)me);
544         
545         /* free temp copy of destination shapekeys (if applicable) */
546         if (nkey) {
547                 /* XXX 2.5 Animato */
548 #if 0
549                 /* free it's ipo too - both are not actually freed from memory yet as ID-blocks */
550                 if (nkey->ipo) {
551                         BKE_ipo_free(nkey->ipo);
552                         BLI_remlink(&bmain->ipo, nkey->ipo);
553                         MEM_freeN(nkey->ipo);
554                 }
555 #endif
556                 
557                 BKE_key_free(nkey);
558                 BLI_remlink(&bmain->key, nkey);
559                 MEM_freeN(nkey);
560         }
561         
562         /* ensure newly inserted keys are time sorted */
563         if (key && (key->type != KEY_RELATIVE)) {
564                 BKE_key_sort(key);
565         }
566
567
568         DAG_relations_tag_update(bmain);   // removed objects, need to rebuild dag
569
570 #if 0
571         ED_object_editmode_enter(C, EM_WAITCURSOR);
572         ED_object_editmode_exit(C, EM_FREEDATA | EM_WAITCURSOR | EM_DO_UNDO);
573 #else
574         /* toggle editmode using lower level functions so this can be called from python */
575         EDBM_mesh_make(scene->toolsettings, ob);
576         EDBM_mesh_load(ob);
577         EDBM_mesh_free(me->edit_btmesh);
578         MEM_freeN(me->edit_btmesh);
579         me->edit_btmesh = NULL;
580         DAG_id_tag_update(&ob->id, OB_RECALC_OB | OB_RECALC_DATA);
581 #endif
582         WM_event_add_notifier(C, NC_SCENE | ND_OB_ACTIVE, scene);
583
584         return OPERATOR_FINISHED;
585 }
586
587 /*********************** JOIN AS SHAPES ***************************/
588
589 /* Append selected meshes vertex locations as shapes of the active mesh, 
590  * return 0 if no join is made (error) and 1 of the join is done */
591
592 int join_mesh_shapes_exec(bContext *C, wmOperator *op)
593 {
594         Scene *scene = CTX_data_scene(C);
595         Object *ob = CTX_data_active_object(C);
596         Mesh *me = (Mesh *)ob->data;
597         Mesh *selme = NULL;
598         DerivedMesh *dm = NULL;
599         Key *key = me->key;
600         KeyBlock *kb;
601         bool ok = false, nonequal_verts = false;
602         
603         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
604         {
605                 if (base->object == ob) continue;
606                 
607                 if (base->object->type == OB_MESH) {
608                         selme = (Mesh *)base->object->data;
609                         
610                         if (selme->totvert == me->totvert)
611                                 ok = true;
612                         else
613                                 nonequal_verts = 1;
614                 }
615         }
616         CTX_DATA_END;
617         
618         if (!ok) {
619                 if (nonequal_verts)
620                         BKE_report(op->reports, RPT_WARNING, "Selected meshes must have equal numbers of vertices");
621                 else
622                         BKE_report(op->reports, RPT_WARNING, "No additional selected meshes with equal vertex count to join");
623                 return OPERATOR_CANCELLED;
624         }
625         
626         if (key == NULL) {
627                 key = me->key = BKE_key_add((ID *)me);
628                 key->type = KEY_RELATIVE;
629
630                 /* first key added, so it was the basis. initialize it with the existing mesh */
631                 kb = BKE_keyblock_add(key, NULL);
632                 BKE_keyblock_convert_from_mesh(me, kb);
633         }
634         
635         /* now ready to add new keys from selected meshes */
636         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
637         {
638                 if (base->object == ob) continue;
639                 
640                 if (base->object->type == OB_MESH) {
641                         selme = (Mesh *)base->object->data;
642                         
643                         if (selme->totvert == me->totvert) {
644                                 dm = mesh_get_derived_deform(scene, base->object, CD_MASK_BAREMESH);
645                                 
646                                 if (!dm) continue;
647                                         
648                                 kb = BKE_keyblock_add(key, base->object->id.name + 2);
649                                 
650                                 DM_to_meshkey(dm, me, kb);
651                                 
652                                 dm->release(dm);
653                         }
654                 }
655         }
656         CTX_DATA_END;
657         
658         WM_event_add_notifier(C, NC_SCENE | ND_OB_ACTIVE, scene);
659         
660         return OPERATOR_FINISHED;
661 }
662
663 /* -------------------------------------------------------------------- */
664 /* Mesh Mirror (Spatial) */
665
666 /** \name Mesh Spatial Mirror API
667  * \{ */
668
669 #define KD_THRESH      0.00002f
670
671 static struct { void *tree; } MirrKdStore = {NULL};
672
673 /* mode is 's' start, or 'e' end, or 'u' use */
674 /* if end, ob can be NULL */
675 int ED_mesh_mirror_spatial_table(Object *ob, BMEditMesh *em, DerivedMesh *dm, const float co[3], char mode)
676 {
677         if (mode == 'u') {        /* use table */
678                 if (MirrKdStore.tree == NULL)
679                         ED_mesh_mirror_spatial_table(ob, em, dm, NULL, 's');
680
681                 if (MirrKdStore.tree) {
682                         KDTreeNearest nearest;
683                         const int i = BLI_kdtree_find_nearest(MirrKdStore.tree, co, &nearest);
684
685                         if (i != -1) {
686                                 if (nearest.dist < KD_THRESH) {
687                                         return i;
688                                 }
689                         }
690                 }
691                 return -1;
692         }
693         else if (mode == 's') {   /* start table */
694                 Mesh *me = ob->data;
695                 const bool use_em = (!dm && em && me->edit_btmesh == em);
696                 const int totvert = use_em ? em->bm->totvert : dm ? dm->getNumVerts(dm) : me->totvert;
697
698                 if (MirrKdStore.tree) /* happens when entering this call without ending it */
699                         ED_mesh_mirror_spatial_table(ob, em, dm, co, 'e');
700
701                 MirrKdStore.tree = BLI_kdtree_new(totvert);
702
703                 if (use_em) {
704                         BMVert *eve;
705                         BMIter iter;
706                         int i;
707
708                         /* this needs to be valid for index lookups later (callers need) */
709                         BM_mesh_elem_table_ensure(em->bm, BM_VERT);
710
711                         BM_ITER_MESH_INDEX (eve, &iter, em->bm, BM_VERTS_OF_MESH, i) {
712                                 BLI_kdtree_insert(MirrKdStore.tree, i, eve->co);
713                         }
714                 }
715                 else {
716                         MVert *mvert = dm ? dm->getVertArray(dm) : me->mvert;
717                         int i;
718                         
719                         for (i = 0; i < totvert; i++, mvert++) {
720                                 BLI_kdtree_insert(MirrKdStore.tree, i, mvert->co);
721                         }
722                 }
723
724                 BLI_kdtree_balance(MirrKdStore.tree);
725         }
726         else if (mode == 'e') { /* end table */
727                 if (MirrKdStore.tree) {
728                         BLI_kdtree_free(MirrKdStore.tree);
729                         MirrKdStore.tree = NULL;
730                 }
731         }
732         else {
733                 BLI_assert(0);
734         }
735
736         return 0;
737 }
738
739 /** \} */
740
741
742 /* -------------------------------------------------------------------- */
743 /* Mesh Mirror (Topology) */
744
745 /** \name Mesh Topology Mirror API
746  * \{ */
747
748 static MirrTopoStore_t mesh_topo_store = {NULL, -1. - 1, -1};
749
750 /* mode is 's' start, or 'e' end, or 'u' use */
751 /* if end, ob can be NULL */
752 /* note, is supposed return -1 on error, which callers are currently checking for, but is not used so far */
753 int ED_mesh_mirror_topo_table(Object *ob, DerivedMesh *dm, char mode)
754 {
755         if (mode == 'u') {        /* use table */
756                 if (ED_mesh_mirrtopo_recalc_check(ob->data, dm, ob->mode, &mesh_topo_store)) {
757                         ED_mesh_mirror_topo_table(ob, dm, 's');
758                 }
759         }
760         else if (mode == 's') { /* start table */
761                 ED_mesh_mirrtopo_init(ob->data, dm, ob->mode, &mesh_topo_store, false);
762         }
763         else if (mode == 'e') { /* end table */
764                 ED_mesh_mirrtopo_free(&mesh_topo_store);
765         }
766         else {
767                 BLI_assert(0);
768         }
769
770         return 0;
771 }
772
773 /** \} */
774
775
776 static int mesh_get_x_mirror_vert_spatial(Object *ob, DerivedMesh *dm, int index)
777 {
778         Mesh *me = ob->data;
779         MVert *mvert = dm ? dm->getVertArray(dm) : me->mvert;
780         float vec[3];
781         
782         mvert = &mvert[index];
783         vec[0] = -mvert->co[0];
784         vec[1] = mvert->co[1];
785         vec[2] = mvert->co[2];
786         
787         return ED_mesh_mirror_spatial_table(ob, NULL, dm, vec, 'u');
788 }
789
790 static int mesh_get_x_mirror_vert_topo(Object *ob, DerivedMesh *dm, int index)
791 {
792         if (ED_mesh_mirror_topo_table(ob, dm, 'u') == -1)
793                 return -1;
794
795         return mesh_topo_store.index_lookup[index];
796 }
797
798 int mesh_get_x_mirror_vert(Object *ob, DerivedMesh *dm, int index, const bool use_topology)
799 {
800         if (use_topology) {
801                 return mesh_get_x_mirror_vert_topo(ob, dm, index);
802         }
803         else {
804                 return mesh_get_x_mirror_vert_spatial(ob, dm, index);
805         }
806 }
807
808 static BMVert *editbmesh_get_x_mirror_vert_spatial(Object *ob, BMEditMesh *em, const float co[3])
809 {
810         float vec[3];
811         int i;
812         
813         /* ignore nan verts */
814         if ((finite(co[0]) == false) ||
815             (finite(co[1]) == false) ||
816             (finite(co[2]) == false))
817         {
818                 return NULL;
819         }
820         
821         vec[0] = -co[0];
822         vec[1] = co[1];
823         vec[2] = co[2];
824         
825         i = ED_mesh_mirror_spatial_table(ob, em, NULL, vec, 'u');
826         if (i != -1) {
827                 return BM_vert_at_index(em->bm, i);
828         }
829         return NULL;
830 }
831
832 static BMVert *editbmesh_get_x_mirror_vert_topo(Object *ob, struct BMEditMesh *em, BMVert *eve, int index)
833 {
834         intptr_t poinval;
835         if (ED_mesh_mirror_topo_table(ob, NULL, 'u') == -1)
836                 return NULL;
837
838         if (index == -1) {
839                 BMIter iter;
840                 BMVert *v;
841                 
842                 index = 0;
843                 BM_ITER_MESH (v, &iter, em->bm, BM_VERTS_OF_MESH) {
844                         if (v == eve)
845                                 break;
846                         index++;
847                 }
848                 
849                 if (index == em->bm->totvert) {
850                         return NULL;
851                 }
852         }
853
854         poinval = mesh_topo_store.index_lookup[index];
855
856         if (poinval != -1)
857                 return (BMVert *)(poinval);
858         return NULL;
859 }       
860
861 BMVert *editbmesh_get_x_mirror_vert(Object *ob, struct BMEditMesh *em, BMVert *eve, const float co[3], int index, const bool use_topology)
862 {
863         if (use_topology) {
864                 return editbmesh_get_x_mirror_vert_topo(ob, em, eve, index);
865         }
866         else {
867                 return editbmesh_get_x_mirror_vert_spatial(ob, em, co);
868         }
869 }
870
871 /**
872  * Wrapper for objectmode/editmode.
873  *
874  * call #BM_mesh_elem_table_ensure first for editmesh.
875  */
876 int ED_mesh_mirror_get_vert(Object *ob, int index)
877 {
878         Mesh *me = ob->data;
879         BMEditMesh *em = me->edit_btmesh;
880         bool use_topology = (me->editflag & ME_EDIT_MIRROR_TOPO) != 0;
881         int index_mirr;
882
883         if (em) {
884                 BMVert *eve, *eve_mirr;
885                 eve = BM_vert_at_index(em->bm, index);
886                 eve_mirr = editbmesh_get_x_mirror_vert(ob, em, eve, eve->co, index, use_topology);
887                 index_mirr = eve_mirr ? BM_elem_index_get(eve_mirr) : -1;
888         }
889         else {
890                 index_mirr = mesh_get_x_mirror_vert(ob, NULL, index, use_topology);
891         }
892
893         return index_mirr;
894 }
895
896 #if 0
897
898 static float *editmesh_get_mirror_uv(BMEditMesh *em, int axis, float *uv, float *mirrCent, float *face_cent)
899 {
900         float vec[2];
901         float cent_vec[2];
902         float cent[2];
903
904         /* ignore nan verts */
905         if (isnan(uv[0]) || !finite(uv[0]) ||
906             isnan(uv[1]) || !finite(uv[1])
907             )
908                 return NULL;
909
910         if (axis) {
911                 vec[0] = uv[0];
912                 vec[1] = -((uv[1]) - mirrCent[1]) + mirrCent[1];
913
914                 cent_vec[0] = face_cent[0];
915                 cent_vec[1] = -((face_cent[1]) - mirrCent[1]) + mirrCent[1];
916         }
917         else {
918                 vec[0] = -((uv[0]) - mirrCent[0]) + mirrCent[0];
919                 vec[1] = uv[1];
920
921                 cent_vec[0] = -((face_cent[0]) - mirrCent[0]) + mirrCent[0];
922                 cent_vec[1] = face_cent[1];
923         }
924
925         /* TODO - Optimize */
926         {
927                 BMIter iter;
928                 BMFace *efa;
929                 
930                 BM_ITER_MESH (efa, &iter, em->bm, BM_FACES_OF_MESH) {
931                         uv_poly_center(efa, cent, cd_loop_uv_offset);
932                         
933                         if ( (fabsf(cent[0] - cent_vec[0]) < 0.001f) && (fabsf(cent[1] - cent_vec[1]) < 0.001f) ) {
934                                 BMIter liter;
935                                 BMLoop *l;
936                                 
937                                 BM_ITER_ELEM (l, &liter, efa, BM_LOOPS_OF_FACE) {
938                                         MLoopUV *luv = CustomData_bmesh_get(&em->bm->ldata, l->head.data, CD_MLOOPUV);
939                                         if ( (fabsf(luv->uv[0] - vec[0]) < 0.001f) && (fabsf(luv->uv[1] - vec[1]) < 0.001f) ) {
940                                                 return luv->uv;
941                                                                 
942                                         }
943                                 }
944                         }
945                 }
946         }
947
948         return NULL;
949 }
950
951 #endif
952
953 static unsigned int mirror_facehash(const void *ptr)
954 {
955         const MFace *mf = ptr;
956         unsigned int v0, v1;
957
958         if (mf->v4) {
959                 v0 = MIN4(mf->v1, mf->v2, mf->v3, mf->v4);
960                 v1 = MAX4(mf->v1, mf->v2, mf->v3, mf->v4);
961         }
962         else {
963                 v0 = MIN3(mf->v1, mf->v2, mf->v3);
964                 v1 = MAX3(mf->v1, mf->v2, mf->v3);
965         }
966
967         return ((v0 * 39) ^ (v1 * 31));
968 }
969
970 static int mirror_facerotation(MFace *a, MFace *b)
971 {
972         if (b->v4) {
973                 if (a->v1 == b->v1 && a->v2 == b->v2 && a->v3 == b->v3 && a->v4 == b->v4)
974                         return 0;
975                 else if (a->v4 == b->v1 && a->v1 == b->v2 && a->v2 == b->v3 && a->v3 == b->v4)
976                         return 1;
977                 else if (a->v3 == b->v1 && a->v4 == b->v2 && a->v1 == b->v3 && a->v2 == b->v4)
978                         return 2;
979                 else if (a->v2 == b->v1 && a->v3 == b->v2 && a->v4 == b->v3 && a->v1 == b->v4)
980                         return 3;
981         }
982         else {
983                 if (a->v1 == b->v1 && a->v2 == b->v2 && a->v3 == b->v3)
984                         return 0;
985                 else if (a->v3 == b->v1 && a->v1 == b->v2 && a->v2 == b->v3)
986                         return 1;
987                 else if (a->v2 == b->v1 && a->v3 == b->v2 && a->v1 == b->v3)
988                         return 2;
989         }
990         
991         return -1;
992 }
993
994 static bool mirror_facecmp(const void *a, const void *b)
995 {
996         return (mirror_facerotation((MFace *)a, (MFace *)b) == -1);
997 }
998
999 /* BMESH_TODO, convert to MPoly (functions above also) */
1000 int *mesh_get_x_mirror_faces(Object *ob, BMEditMesh *em, DerivedMesh *dm)
1001 {
1002         Mesh *me = ob->data;
1003         MVert *mv, *mvert;
1004         MFace mirrormf, *mf, *hashmf, *mface;
1005         GHash *fhash;
1006         int *mirrorverts, *mirrorfaces;
1007
1008         BLI_assert(em == NULL);  /* Does not work otherwise, currently... */
1009
1010         const bool use_topology = (me->editflag & ME_EDIT_MIRROR_TOPO) != 0;
1011         const int totvert = dm ? dm->getNumVerts(dm) : me->totvert;
1012         const int totface = dm ? dm->getNumTessFaces(dm) : me->totface;
1013         int a;
1014
1015         mirrorverts = MEM_callocN(sizeof(int) * totvert, "MirrorVerts");
1016         mirrorfaces = MEM_callocN(sizeof(int) * 2 * totface, "MirrorFaces");
1017
1018         mvert = dm ? dm->getVertArray(dm) : me->mvert;
1019         mface = dm ? dm->getTessFaceArray(dm) : me->mface;
1020
1021         ED_mesh_mirror_spatial_table(ob, em, dm, NULL, 's');
1022
1023         for (a = 0, mv = mvert; a < totvert; a++, mv++)
1024                 mirrorverts[a] = mesh_get_x_mirror_vert(ob, dm, a, use_topology);
1025
1026         ED_mesh_mirror_spatial_table(ob, em, dm, NULL, 'e');
1027
1028         fhash = BLI_ghash_new_ex(mirror_facehash, mirror_facecmp, "mirror_facehash gh", me->totface);
1029         for (a = 0, mf = mface; a < totface; a++, mf++)
1030                 BLI_ghash_insert(fhash, mf, mf);
1031
1032         for (a = 0, mf = mface; a < totface; a++, mf++) {
1033                 mirrormf.v1 = mirrorverts[mf->v3];
1034                 mirrormf.v2 = mirrorverts[mf->v2];
1035                 mirrormf.v3 = mirrorverts[mf->v1];
1036                 mirrormf.v4 = (mf->v4) ? mirrorverts[mf->v4] : 0;
1037
1038                 /* make sure v4 is not 0 if a quad */
1039                 if (mf->v4 && mirrormf.v4 == 0) {
1040                         SWAP(unsigned int, mirrormf.v1, mirrormf.v3);
1041                         SWAP(unsigned int, mirrormf.v2, mirrormf.v4);
1042                 }
1043
1044                 hashmf = BLI_ghash_lookup(fhash, &mirrormf);
1045                 if (hashmf) {
1046                         mirrorfaces[a * 2] = hashmf - mface;
1047                         mirrorfaces[a * 2 + 1] = mirror_facerotation(&mirrormf, hashmf);
1048                 }
1049                 else
1050                         mirrorfaces[a * 2] = -1;
1051         }
1052
1053         BLI_ghash_free(fhash, NULL, NULL);
1054         MEM_freeN(mirrorverts);
1055         
1056         return mirrorfaces;
1057 }
1058
1059 /* selection, vertex and face */
1060 /* returns 0 if not found, otherwise 1 */
1061
1062 /**
1063  * Face selection in object mode,
1064  * currently only weight-paint and vertex-paint use this.
1065  *
1066  * \return boolean true == Found
1067  */
1068 bool ED_mesh_pick_face(bContext *C, Object *ob, const int mval[2], unsigned int *index, int size)
1069 {
1070         ViewContext vc;
1071         Mesh *me = ob->data;
1072
1073         BLI_assert(me && GS(me->id.name) == ID_ME);
1074
1075         if (!me || me->totpoly == 0)
1076                 return false;
1077
1078         view3d_set_viewcontext(C, &vc);
1079
1080         if (size) {
1081                 /* sample rect to increase chances of selecting, so that when clicking
1082                  * on an edge in the backbuf, we can still select a face */
1083
1084                 float dummy_dist;
1085                 *index = ED_view3d_backbuf_sample_rect(&vc, mval, size, 1, me->totpoly + 1, &dummy_dist);
1086         }
1087         else {
1088                 /* sample only on the exact position */
1089                 *index = ED_view3d_backbuf_sample(&vc, mval[0], mval[1]);
1090         }
1091
1092         if ((*index) == 0 || (*index) > (unsigned int)me->totpoly)
1093                 return false;
1094
1095         (*index)--;
1096
1097         return true;
1098 }
1099 static void ed_mesh_pick_face_vert__mpoly_find(
1100         /* context */
1101         struct ARegion *ar, const float mval[2],
1102         /* mesh data */
1103         DerivedMesh *dm, MPoly *mp, MLoop *mloop,
1104         /* return values */
1105         float *r_len_best, int *r_v_idx_best)
1106 {
1107         const MLoop *ml;
1108         int j = mp->totloop;
1109         for (ml = &mloop[mp->loopstart]; j--; ml++) {
1110                 float co[3], sco[2], len;
1111                 const int v_idx = ml->v;
1112                 dm->getVertCo(dm, v_idx, co);
1113                 if (ED_view3d_project_float_object(ar, co, sco, V3D_PROJ_TEST_NOP) == V3D_PROJ_RET_OK) {
1114                         len = len_manhattan_v2v2(mval, sco);
1115                         if (len < *r_len_best) {
1116                                 *r_len_best = len;
1117                                 *r_v_idx_best = v_idx;
1118                         }
1119                 }
1120         }
1121 }
1122 /**
1123  * Use when the back buffer stores face index values. but we want a vert.
1124  * This gets the face then finds the closest vertex to mval.
1125  */
1126 bool ED_mesh_pick_face_vert(bContext *C, Object *ob, const int mval[2], unsigned int *index, int size)
1127 {
1128         unsigned int poly_index;
1129         Mesh *me = ob->data;
1130
1131         BLI_assert(me && GS(me->id.name) == ID_ME);
1132
1133         if (ED_mesh_pick_face(C, ob, mval, &poly_index, size)) {
1134                 Scene *scene = CTX_data_scene(C);
1135                 struct ARegion *ar = CTX_wm_region(C);
1136
1137                 /* derived mesh to find deformed locations */
1138                 DerivedMesh *dm = mesh_get_derived_final(scene, ob, CD_MASK_BAREMESH | CD_MASK_ORIGINDEX);
1139
1140                 int v_idx_best = ORIGINDEX_NONE;
1141
1142                 /* find the vert closest to 'mval' */
1143                 const float mval_f[2] = {UNPACK2(mval)};
1144                 float len_best = FLT_MAX;
1145
1146                 MPoly *dm_mpoly;
1147                 MLoop *dm_mloop;
1148                 unsigned int dm_mpoly_tot;
1149                 const int *index_mp_to_orig;
1150
1151                 dm_mpoly = dm->getPolyArray(dm);
1152                 dm_mloop = dm->getLoopArray(dm);
1153
1154                 dm_mpoly_tot = dm->getNumPolys(dm);
1155
1156                 index_mp_to_orig = dm->getPolyDataArray(dm, CD_ORIGINDEX);
1157
1158                 /* tag all verts using this face */
1159                 if (index_mp_to_orig) {
1160                         unsigned int i;
1161
1162                         for (i = 0; i < dm_mpoly_tot; i++) {
1163                                 if (index_mp_to_orig[i] == poly_index) {
1164                                         ed_mesh_pick_face_vert__mpoly_find(
1165                                                 ar, mval_f,
1166                                                 dm, &dm_mpoly[i], dm_mloop,
1167                                                 &len_best, &v_idx_best);
1168                                 }
1169                         }
1170                 }
1171                 else {
1172                         if (poly_index < dm_mpoly_tot) {
1173                                 ed_mesh_pick_face_vert__mpoly_find(
1174                                         ar, mval_f,
1175                                         dm, &dm_mpoly[poly_index], dm_mloop,
1176                                         &len_best, &v_idx_best);
1177                         }
1178                 }
1179
1180                 /* map 'dm -> me' index if possible */
1181                 if (v_idx_best != ORIGINDEX_NONE) {
1182                         const int *index_mv_to_orig;
1183
1184                         index_mv_to_orig = dm->getVertDataArray(dm, CD_ORIGINDEX);
1185                         if (index_mv_to_orig) {
1186                                 v_idx_best = index_mv_to_orig[v_idx_best];
1187                         }
1188                 }
1189
1190                 dm->release(dm);
1191
1192                 if ((v_idx_best != ORIGINDEX_NONE) && (v_idx_best < me->totvert)) {
1193                         *index = v_idx_best;
1194                         return true;
1195                 }
1196         }
1197
1198         return false;
1199 }
1200
1201 /**
1202  * Vertex selection in object mode,
1203  * currently only weight paint uses this.
1204  *
1205  * \return boolean true == Found
1206  */
1207 typedef struct VertPickData {
1208         const MVert *mvert;
1209         const float *mval_f;  /* [2] */
1210         ARegion *ar;
1211
1212         /* runtime */
1213         float len_best;
1214         int v_idx_best;
1215 } VertPickData;
1216
1217 static void ed_mesh_pick_vert__mapFunc(void *userData, int index, const float co[3],
1218                                        const float UNUSED(no_f[3]), const short UNUSED(no_s[3]))
1219 {
1220         VertPickData *data = userData;
1221         if ((data->mvert[index].flag & ME_HIDE) == 0) {
1222                 float sco[2];
1223
1224                 if (ED_view3d_project_float_object(data->ar, co, sco, V3D_PROJ_TEST_CLIP_DEFAULT) == V3D_PROJ_RET_OK) {
1225                         const float len = len_manhattan_v2v2(data->mval_f, sco);
1226                         if (len < data->len_best) {
1227                                 data->len_best = len;
1228                                 data->v_idx_best = index;
1229                         }
1230                 }
1231         }
1232 }
1233 bool ED_mesh_pick_vert(bContext *C, Object *ob, const int mval[2], unsigned int *index, int size, bool use_zbuf)
1234 {
1235         ViewContext vc;
1236         Mesh *me = ob->data;
1237
1238         BLI_assert(me && GS(me->id.name) == ID_ME);
1239
1240         if (!me || me->totvert == 0)
1241                 return false;
1242
1243         view3d_set_viewcontext(C, &vc);
1244
1245         if (use_zbuf) {
1246                 if (size > 0) {
1247                         /* sample rect to increase chances of selecting, so that when clicking
1248                          * on an face in the backbuf, we can still select a vert */
1249
1250                         float dummy_dist;
1251                         *index = ED_view3d_backbuf_sample_rect(&vc, mval, size, 1, me->totvert + 1, &dummy_dist);
1252                 }
1253                 else {
1254                         /* sample only on the exact position */
1255                         *index = ED_view3d_backbuf_sample(&vc, mval[0], mval[1]);
1256                 }
1257
1258                 if ((*index) == 0 || (*index) > (unsigned int)me->totvert)
1259                         return false;
1260
1261                 (*index)--;
1262         }
1263         else {
1264                 /* derived mesh to find deformed locations */
1265                 DerivedMesh *dm = mesh_get_derived_final(vc.scene, ob, CD_MASK_BAREMESH);
1266                 ARegion *ar = vc.ar;
1267                 RegionView3D *rv3d = ar->regiondata;
1268
1269                 /* find the vert closest to 'mval' */
1270                 const float mval_f[2] = {(float)mval[0],
1271                                          (float)mval[1]};
1272
1273                 VertPickData data = {NULL};
1274
1275                 ED_view3d_init_mats_rv3d(ob, rv3d);
1276
1277                 if (dm == NULL) {
1278                         return false;
1279                 }
1280
1281                 /* setup data */
1282                 data.mvert = me->mvert;
1283                 data.ar = ar;
1284                 data.mval_f = mval_f;
1285                 data.len_best = FLT_MAX;
1286                 data.v_idx_best = -1;
1287
1288                 dm->foreachMappedVert(dm, ed_mesh_pick_vert__mapFunc, &data, DM_FOREACH_NOP);
1289
1290                 dm->release(dm);
1291
1292                 if (data.v_idx_best == -1) {
1293                         return false;
1294                 }
1295
1296                 *index = data.v_idx_best;
1297         }
1298
1299         return true;
1300 }
1301
1302
1303 MDeformVert *ED_mesh_active_dvert_get_em(Object *ob, BMVert **r_eve)
1304 {
1305         if (ob->mode & OB_MODE_EDIT && ob->type == OB_MESH && ob->defbase.first) {
1306                 Mesh *me = ob->data;
1307                 BMesh *bm = me->edit_btmesh->bm;
1308                 const int cd_dvert_offset = CustomData_get_offset(&bm->vdata, CD_MDEFORMVERT);
1309
1310                 if (cd_dvert_offset != -1) {
1311                         BMVert *eve = BM_mesh_active_vert_get(bm);
1312
1313                         if (eve) {
1314                                 if (r_eve) *r_eve = eve;
1315                                 return BM_ELEM_CD_GET_VOID_P(eve, cd_dvert_offset);
1316                         }
1317                 }
1318         }
1319
1320         if (r_eve) *r_eve = NULL;
1321         return NULL;
1322 }
1323
1324 MDeformVert *ED_mesh_active_dvert_get_ob(Object *ob, int *r_index)
1325 {
1326         Mesh *me = ob->data;
1327         int index = BKE_mesh_mselect_active_get(me, ME_VSEL);
1328         if (r_index) *r_index = index;
1329         if (index == -1 || me->dvert == NULL) {
1330                 return NULL;
1331         }
1332         else {
1333                 return me->dvert + index;
1334         }
1335 }
1336
1337 MDeformVert *ED_mesh_active_dvert_get_only(Object *ob)
1338 {
1339         if (ob->type == OB_MESH) {
1340                 if (ob->mode & OB_MODE_EDIT) {
1341                         return ED_mesh_active_dvert_get_em(ob, NULL);
1342                 }
1343                 else {
1344                         return ED_mesh_active_dvert_get_ob(ob, NULL);
1345                 }
1346         }
1347         else {
1348                 return NULL;
1349         }
1350 }