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