Merge branch 'master' into blender2.8
[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 #include "DNA_workspace_types.h"
46
47 #include "BLI_math.h"
48 #include "BLI_blenlib.h"
49
50 #include "BKE_context.h"
51 #include "BKE_deform.h"
52 #include "BKE_DerivedMesh.h"
53 #include "BKE_key.h"
54 #include "BKE_library.h"
55 #include "BKE_main.h"
56 #include "BKE_mesh.h"
57 #include "BKE_mesh_runtime.h"
58 #include "BKE_material.h"
59 #include "BKE_object.h"
60 #include "BKE_object_deform.h"
61 #include "BKE_report.h"
62 #include "BKE_editmesh.h"
63 #include "BKE_multires.h"
64 #include "BKE_layer.h"
65
66 #include "DEG_depsgraph.h"
67 #include "DEG_depsgraph_build.h"
68
69 #include "ED_mesh.h"
70 #include "ED_object.h"
71 #include "ED_view3d.h"
72
73 #include "WM_api.h"
74 #include "WM_types.h"
75
76 /* * ********************** no editmode!!! *********** */
77
78 /*********************** JOIN ***************************/
79
80 /* join selected meshes into the active mesh, context sensitive
81  * return 0 if no join is made (error) and 1 if the join is done */
82
83 static void join_mesh_single(
84         Depsgraph *depsgraph, Main *bmain, Scene *scene,
85         Object *ob_dst, Object *ob_src, float imat[4][4],
86         MVert **mvert_pp, MEdge **medge_pp, MLoop **mloop_pp, MPoly **mpoly_pp,
87         CustomData *vdata, CustomData *edata, CustomData *ldata, CustomData *pdata,
88         int totvert, int totedge, int totloop, int totpoly,
89         Key *key, Key *nkey,
90         Material **matar, int *matmap, int totcol,
91         int *vertofs, int *edgeofs, int *loopofs, int *polyofs)
92 {
93         int a, b;
94
95         Mesh *me = ob_src->data;
96         MVert *mvert = *mvert_pp;
97         MEdge *medge = *medge_pp;
98         MLoop *mloop = *mloop_pp;
99         MPoly *mpoly = *mpoly_pp;
100
101         if (me->totvert) {
102                 /* merge customdata flag */
103                 ((Mesh *)ob_dst->data)->cd_flag |= me->cd_flag;
104
105                 /* standard data */
106                 CustomData_merge(&me->vdata, vdata, CD_MASK_MESH, CD_DEFAULT, totvert);
107                 CustomData_copy_data_named(&me->vdata, vdata, 0, *vertofs, me->totvert);
108
109                 /* vertex groups */
110                 MDeformVert *dvert = CustomData_get(vdata, *vertofs, CD_MDEFORMVERT);
111                 MDeformVert *dvert_src = CustomData_get(&me->vdata, 0, CD_MDEFORMVERT);
112
113                 /* Remap to correct new vgroup indices, if needed. */
114                 if (dvert_src) {
115                         BLI_assert(dvert != NULL);
116
117                         /* Build src to merged mapping of vgroup indices. */
118                         int *vgroup_index_map;
119                         int vgroup_index_map_len;
120                         vgroup_index_map = BKE_object_defgroup_index_map_create(ob_src, ob_dst, &vgroup_index_map_len);
121                         BKE_object_defgroup_index_map_apply(dvert, me->totvert, vgroup_index_map, vgroup_index_map_len);
122                         if (vgroup_index_map != NULL) {
123                                 MEM_freeN(vgroup_index_map);
124                         }
125                 }
126
127                 /* if this is the object we're merging into, no need to do anything */
128                 if (ob_src != ob_dst) {
129                         float cmat[4][4];
130
131                         /* watch this: switch matmul order really goes wrong */
132                         mul_m4_m4m4(cmat, imat, ob_src->obmat);
133
134                         /* transform vertex coordinates into new space */
135                         for (a = 0, mvert = *mvert_pp; a < me->totvert; a++, mvert++) {
136                                 mul_m4_v3(cmat, mvert->co);
137                         }
138
139                         /* for each shapekey in destination mesh:
140                          *      - if there's a matching one, copy it across (will need to transform vertices into new space...)
141                          *      - otherwise, just copy own coordinates of mesh (no need to transform vertex coordinates into new space)
142                          */
143                         if (key) {
144                                 /* if this mesh has any shapekeys, check first, otherwise just copy coordinates */
145                                 for (KeyBlock *kb = key->block.first; kb; kb = kb->next) {
146                                         /* get pointer to where to write data for this mesh in shapekey's data array */
147                                         float (*cos)[3] = ((float (*)[3])kb->data) + *vertofs;
148
149                                         /* check if this mesh has such a shapekey */
150                                         KeyBlock *okb = me->key ? BKE_keyblock_find_name(me->key, kb->name) : NULL;
151                                         if (okb) {
152                                                 /* copy this mesh's shapekey to the destination shapekey (need to transform first) */
153                                                 float (*ocos)[3] = okb->data;
154                                                 for (a = 0; a < me->totvert; a++, cos++, ocos++) {
155                                                         copy_v3_v3(*cos, *ocos);
156                                                         mul_m4_v3(cmat, *cos);
157                                                 }
158                                         }
159                                         else {
160                                                 /* copy this mesh's vertex coordinates to the destination shapekey */
161                                                 for (a = 0, mvert = *mvert_pp; a < me->totvert; a++, cos++, mvert++) {
162                                                         copy_v3_v3(*cos, mvert->co);
163                                                 }
164                                         }
165                                 }
166                         }
167                 }
168                 else {
169                         /* for each shapekey in destination mesh:
170                          *      - if it was an 'original', copy the appropriate data from nkey
171                          *      - otherwise, copy across plain coordinates (no need to transform coordinates)
172                          */
173                         if (key) {
174                                 for (KeyBlock *kb = key->block.first; kb; kb = kb->next) {
175                                         /* get pointer to where to write data for this mesh in shapekey's data array */
176                                         float (*cos)[3] = ((float (*)[3])kb->data) + *vertofs;
177
178                                         /* check if this was one of the original shapekeys */
179                                         KeyBlock *okb = nkey ? BKE_keyblock_find_name(nkey, kb->name) : NULL;
180                                         if (okb) {
181                                                 /* copy this mesh's shapekey to the destination shapekey */
182                                                 float (*ocos)[3] = okb->data;
183                                                 for (a = 0; a < me->totvert; a++, cos++, ocos++) {
184                                                         copy_v3_v3(*cos, *ocos);
185                                                 }
186                                         }
187                                         else {
188                                                 /* copy base-coordinates to the destination shapekey */
189                                                 for (a = 0, mvert = *mvert_pp; a < me->totvert; a++, cos++, mvert++) {
190                                                         copy_v3_v3(*cos, mvert->co);
191                                                 }
192                                         }
193                                 }
194                         }
195                 }
196         }
197
198         if (me->totedge) {
199                 CustomData_merge(&me->edata, edata, CD_MASK_MESH, CD_DEFAULT, totedge);
200                 CustomData_copy_data_named(&me->edata, edata, 0, *edgeofs, me->totedge);
201
202                 for (a = 0; a < me->totedge; a++, medge++) {
203                         medge->v1 += *vertofs;
204                         medge->v2 += *vertofs;
205                 }
206         }
207
208         if (me->totloop) {
209                 if (ob_src != ob_dst) {
210                         MultiresModifierData *mmd;
211
212                         multiresModifier_prepare_join(depsgraph, scene, ob_src, ob_dst);
213
214                         if ((mmd = get_multires_modifier(scene, ob_src, true))) {
215                                 ED_object_iter_other(bmain, ob_src, true,
216                                                      ED_object_multires_update_totlevels_cb,
217                                                      &mmd->totlvl);
218                         }
219                 }
220
221                 CustomData_merge(&me->ldata, ldata, CD_MASK_MESH, CD_DEFAULT, totloop);
222                 CustomData_copy_data_named(&me->ldata, ldata, 0, *loopofs, me->totloop);
223
224                 for (a = 0; a < me->totloop; a++, mloop++) {
225                         mloop->v += *vertofs;
226                         mloop->e += *edgeofs;
227                 }
228         }
229
230         if (me->totpoly) {
231                 if (matmap) {
232                         /* make mapping for materials */
233                         for (a = 1; a <= ob_src->totcol; a++) {
234                                 Material *ma = give_current_material(ob_src, a);
235
236                                 for (b = 0; b < totcol; b++) {
237                                         if (ma == matar[b]) {
238                                                 matmap[a - 1] = b;
239                                                 break;
240                                         }
241                                 }
242                         }
243                 }
244
245                 CustomData_merge(&me->pdata, pdata, CD_MASK_MESH, CD_DEFAULT, totpoly);
246                 CustomData_copy_data_named(&me->pdata, pdata, 0, *polyofs, me->totpoly);
247
248                 for (a = 0; a < me->totpoly; a++, mpoly++) {
249                         mpoly->loopstart += *loopofs;
250                         mpoly->mat_nr = matmap ? matmap[mpoly->mat_nr] : 0;
251                 }
252         }
253
254         /* these are used for relinking (cannot be set earlier, or else reattaching goes wrong) */
255         *vertofs += me->totvert;
256         *mvert_pp += me->totvert;
257         *edgeofs += me->totedge;
258         *medge_pp += me->totedge;
259         *loopofs += me->totloop;
260         *mloop_pp += me->totloop;
261         *polyofs += me->totpoly;
262         *mpoly_pp += me->totpoly;
263 }
264
265 int join_mesh_exec(bContext *C, wmOperator *op)
266 {
267         Main *bmain = CTX_data_main(C);
268         Scene *scene = CTX_data_scene(C);
269         Object *ob = CTX_data_active_object(C);
270         Material **matar = NULL, *ma;
271         Mesh *me;
272         MVert *mvert = NULL;
273         MEdge *medge = NULL;
274         MPoly *mpoly = NULL;
275         MLoop *mloop = NULL;
276         Key *key, *nkey = NULL;
277         KeyBlock *kb, *kbn;
278         float imat[4][4];
279         int a, b, totcol, totmat = 0, totedge = 0, totvert = 0;
280         int totloop = 0, totpoly = 0, vertofs, *matmap = NULL;
281         int i, haskey = 0, edgeofs, loopofs, polyofs;
282         bool ok = false;
283         bDeformGroup *dg, *odg;
284         CustomData vdata, edata, fdata, ldata, pdata;
285
286         if (ob->mode & OB_MODE_EDIT) {
287                 BKE_report(op->reports, RPT_WARNING, "Cannot join while in edit mode");
288                 return OPERATOR_CANCELLED;
289         }
290
291         /* ob is the object we are adding geometry to */
292         if (!ob || ob->type != OB_MESH) {
293                 BKE_report(op->reports, RPT_WARNING, "Active object is not a mesh");
294                 return OPERATOR_CANCELLED;
295         }
296
297         Depsgraph *depsgraph = CTX_data_depsgraph(C);
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(bmain, (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
450                                                 kb_map[i] = kbn;
451                                         }
452
453                                         /* remap relative index values */
454                                         for (kb = me->key->block.first, i = 0; kb; kb = kb->next, i++) {
455                                                 if (LIKELY(kb->relative < me->key->totkey)) {  /* sanity check, should always be true */
456                                                         kb_map[i]->relative = index_map[kb->relative];
457                                                 }
458                                         }
459
460                                         MEM_freeN(index_map);
461                                         MEM_freeN(kb_map);
462                                 }
463                         }
464                 }
465         }
466         CTX_DATA_END;
467
468
469         /* setup new data for destination mesh */
470         CustomData_reset(&vdata);
471         CustomData_reset(&edata);
472         CustomData_reset(&fdata);
473         CustomData_reset(&ldata);
474         CustomData_reset(&pdata);
475
476         mvert = CustomData_add_layer(&vdata, CD_MVERT, CD_CALLOC, NULL, totvert);
477         medge = CustomData_add_layer(&edata, CD_MEDGE, CD_CALLOC, NULL, totedge);
478         mloop = CustomData_add_layer(&ldata, CD_MLOOP, CD_CALLOC, NULL, totloop);
479         mpoly = CustomData_add_layer(&pdata, CD_MPOLY, CD_CALLOC, NULL, totpoly);
480
481         vertofs = 0;
482         edgeofs = 0;
483         loopofs = 0;
484         polyofs = 0;
485
486         /* inverse transform for all selected meshes in this object */
487         invert_m4_m4(imat, ob->obmat);
488
489         /* Add back active mesh first. This allows to keep things similar as they were, as much as possible (i.e. data from
490          * active mesh will remain first ones in new result of the merge, in same order for CD layers, etc. See also T50084.
491          */
492         join_mesh_single(
493                     depsgraph, bmain, scene,
494                     ob, ob, imat,
495                     &mvert, &medge, &mloop, &mpoly,
496                     &vdata, &edata, &ldata, &pdata,
497                     totvert, totedge, totloop, totpoly,
498                     key, nkey,
499                     matar, matmap, totcol,
500                     &vertofs, &edgeofs, &loopofs, &polyofs);
501
502         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
503         {
504                 if (base->object == ob) {
505                         continue;
506                 }
507                 /* only join if this is a mesh */
508                 if (base->object->type == OB_MESH) {
509                         join_mesh_single(
510                                     depsgraph, bmain, scene,
511                                     ob, base->object, imat,
512                                     &mvert, &medge, &mloop, &mpoly,
513                                     &vdata, &edata, &ldata, &pdata,
514                                     totvert, totedge, totloop, totpoly,
515                                     key, nkey,
516                                     matar, matmap, totcol,
517                                     &vertofs, &edgeofs, &loopofs, &polyofs);
518
519                         /* free base, now that data is merged */
520                         if (base->object != ob) {
521                                 ED_object_base_free_and_unlink(bmain, scene, base->object);
522                         }
523                 }
524         }
525         CTX_DATA_END;
526
527         /* return to mesh we're merging to */
528         me = ob->data;
529
530         CustomData_free(&me->vdata, me->totvert);
531         CustomData_free(&me->edata, me->totedge);
532         CustomData_free(&me->ldata, me->totloop);
533         CustomData_free(&me->pdata, me->totpoly);
534
535         me->totvert = totvert;
536         me->totedge = totedge;
537         me->totloop = totloop;
538         me->totpoly = totpoly;
539
540         me->vdata = vdata;
541         me->edata = edata;
542         me->ldata = ldata;
543         me->pdata = pdata;
544
545         /* tessface data removed above, no need to update */
546         BKE_mesh_update_customdata_pointers(me, false);
547
548         /* update normals in case objects with non-uniform scale are joined */
549         BKE_mesh_calc_normals(me);
550
551         /* old material array */
552         for (a = 1; a <= ob->totcol; a++) {
553                 ma = ob->mat[a - 1];
554                 if (ma)
555                         id_us_min(&ma->id);
556         }
557         for (a = 1; a <= me->totcol; a++) {
558                 ma = me->mat[a - 1];
559                 if (ma)
560                         id_us_min(&ma->id);
561         }
562         MEM_SAFE_FREE(ob->mat);
563         MEM_SAFE_FREE(ob->matbits);
564         MEM_SAFE_FREE(me->mat);
565
566         if (totcol) {
567                 me->mat = matar;
568                 ob->mat = MEM_callocN(sizeof(*ob->mat) * totcol, "join obmatar");
569                 ob->matbits = MEM_callocN(sizeof(*ob->matbits) * totcol, "join obmatbits");
570                 MEM_freeN(matmap);
571         }
572
573         ob->totcol = me->totcol = totcol;
574
575         /* other mesh users */
576         test_all_objects_materials(bmain, (ID *)me);
577
578         /* free temp copy of destination shapekeys (if applicable) */
579         if (nkey) {
580                 /* We can assume nobody is using that ID currently. */
581                 BKE_libblock_free_ex(bmain, nkey, false, false);
582         }
583
584         /* ensure newly inserted keys are time sorted */
585         if (key && (key->type != KEY_RELATIVE)) {
586                 BKE_key_sort(key);
587         }
588
589         /* Due to dependnecy cycle some other object might access old derived data. */
590         BKE_object_free_derived_caches(ob);
591
592         DEG_relations_tag_update(bmain);   /* removed objects, need to rebuild dag */
593
594         DEG_id_tag_update(&ob->id, OB_RECALC_OB | OB_RECALC_DATA);
595
596         WM_event_add_notifier(C, NC_SCENE | ND_OB_ACTIVE, scene);
597
598         return OPERATOR_FINISHED;
599 }
600
601 /*********************** JOIN AS SHAPES ***************************/
602
603 /* Append selected meshes vertex locations as shapes of the active mesh,
604  * return 0 if no join is made (error) and 1 of the join is done */
605
606 int join_mesh_shapes_exec(bContext *C, wmOperator *op)
607 {
608         Main *bmain = CTX_data_main(C);
609         Scene *scene = CTX_data_scene(C);
610         Object *ob = CTX_data_active_object(C);
611         Depsgraph *depsgraph = CTX_data_depsgraph(C);
612         Mesh *me = (Mesh *)ob->data;
613         Mesh *selme = NULL;
614         DerivedMesh *dm = NULL;
615         Key *key = me->key;
616         KeyBlock *kb;
617         bool ok = false, nonequal_verts = false;
618
619         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
620         {
621                 if (base->object == ob) continue;
622
623                 if (base->object->type == OB_MESH) {
624                         selme = (Mesh *)base->object->data;
625
626                         if (selme->totvert == me->totvert)
627                                 ok = true;
628                         else
629                                 nonequal_verts = 1;
630                 }
631         }
632         CTX_DATA_END;
633
634         if (!ok) {
635                 if (nonequal_verts)
636                         BKE_report(op->reports, RPT_WARNING, "Selected meshes must have equal numbers of vertices");
637                 else
638                         BKE_report(op->reports, RPT_WARNING, "No additional selected meshes with equal vertex count to join");
639                 return OPERATOR_CANCELLED;
640         }
641
642         if (key == NULL) {
643                 key = me->key = BKE_key_add(bmain, (ID *)me);
644                 key->type = KEY_RELATIVE;
645
646                 /* first key added, so it was the basis. initialize it with the existing mesh */
647                 kb = BKE_keyblock_add(key, NULL);
648                 BKE_keyblock_convert_from_mesh(me, kb);
649         }
650
651         /* now ready to add new keys from selected meshes */
652         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
653         {
654                 if (base->object == ob) continue;
655
656                 if (base->object->type == OB_MESH) {
657                         selme = (Mesh *)base->object->data;
658
659                         if (selme->totvert == me->totvert) {
660                                 dm = mesh_get_derived_deform(depsgraph, scene, base->object, CD_MASK_BAREMESH);
661
662                                 if (!dm) continue;
663
664                                 kb = BKE_keyblock_add(key, base->object->id.name + 2);
665
666                                 DM_to_meshkey(dm, me, kb);
667
668                                 dm->release(dm);
669                         }
670                 }
671         }
672         CTX_DATA_END;
673
674         WM_event_add_notifier(C, NC_SCENE | ND_OB_ACTIVE, scene);
675
676         return OPERATOR_FINISHED;
677 }
678
679
680 /* -------------------------------------------------------------------- */
681 /* Mesh Mirror (Topology) */
682
683 /** \name Mesh Topology Mirror API
684  * \{ */
685
686 static MirrTopoStore_t mesh_topo_store = {NULL, -1. - 1, -1};
687
688 /* mode is 's' start, or 'e' end, or 'u' use */
689 /* if end, ob can be NULL */
690 /* note, is supposed return -1 on error, which callers are currently checking for, but is not used so far */
691 int ED_mesh_mirror_topo_table(
692         Object *ob, DerivedMesh *dm, char mode)
693 {
694         if (mode == 'u') {        /* use table */
695                 if (ED_mesh_mirrtopo_recalc_check(ob->data, dm, &mesh_topo_store)) {
696                         ED_mesh_mirror_topo_table(ob, dm, 's');
697                 }
698         }
699         else if (mode == 's') { /* start table */
700                 ED_mesh_mirrtopo_init(ob->data, dm, &mesh_topo_store, false);
701         }
702         else if (mode == 'e') { /* end table */
703                 ED_mesh_mirrtopo_free(&mesh_topo_store);
704         }
705         else {
706                 BLI_assert(0);
707         }
708
709         return 0;
710 }
711
712 /* mode is 's' start, or 'e' end, or 'u' use */
713 /* if end, ob can be NULL */
714 /* note, is supposed return -1 on error, which callers are currently checking for, but is not used so far */
715 int ED_mesh_mirror_topo_table__real_mesh(
716         Object *ob, Mesh *mesh, char mode)
717 {
718         if (mode == 'u') {        /* use table */
719                 if (ED_mesh_mirrtopo_recalc_check__real_mesh(ob->data, mesh, &mesh_topo_store)) {
720                         ED_mesh_mirror_topo_table__real_mesh(ob, mesh, 's');
721                 }
722         }
723         else if (mode == 's') { /* start table */
724                 ED_mesh_mirrtopo_init__real_mesh(ob->data, mesh, &mesh_topo_store, false);
725         }
726         else if (mode == 'e') { /* end table */
727                 ED_mesh_mirrtopo_free(&mesh_topo_store);
728         }
729         else {
730                 BLI_assert(0);
731         }
732
733         return 0;
734 }
735
736 /** \} */
737
738
739 static int mesh_get_x_mirror_vert_spatial(Object *ob, DerivedMesh *dm, int index)
740 {
741         Mesh *me = ob->data;
742         MVert *mvert = dm ? dm->getVertArray(dm) : me->mvert;
743         float vec[3];
744
745         mvert = &mvert[index];
746         vec[0] = -mvert->co[0];
747         vec[1] = mvert->co[1];
748         vec[2] = mvert->co[2];
749
750         return ED_mesh_mirror_spatial_table(ob, NULL, dm, vec, 'u');
751 }
752
753 static int mesh_get_x_mirror_vert_topo(Object *ob, DerivedMesh *dm, int index)
754 {
755         if (ED_mesh_mirror_topo_table(ob, dm, 'u') == -1)
756                 return -1;
757
758         return mesh_topo_store.index_lookup[index];
759 }
760
761 int mesh_get_x_mirror_vert(Object *ob, DerivedMesh *dm, int index, const bool use_topology)
762 {
763         if (use_topology) {
764                 return mesh_get_x_mirror_vert_topo(ob, dm, index);
765         }
766         else {
767                 return mesh_get_x_mirror_vert_spatial(ob, dm, index);
768         }
769 }
770
771 static int mesh_get_x_mirror_vert_spatial__real_mesh(Object *ob, Mesh *mesh, int index)
772 {
773         Mesh *me = ob->data;
774         MVert *mvert = mesh ? mesh->mvert : me->mvert;
775         float vec[3];
776
777         mvert = &mvert[index];
778         vec[0] = -mvert->co[0];
779         vec[1] = mvert->co[1];
780         vec[2] = mvert->co[2];
781
782         return ED_mesh_mirror_spatial_table__real_mesh(ob, NULL, mesh, vec, 'u');
783 }
784
785 static int mesh_get_x_mirror_vert_topo__real_mesh(Object *ob, Mesh *mesh, int index)
786 {
787         if (ED_mesh_mirror_topo_table__real_mesh(ob, mesh, 'u') == -1)
788                 return -1;
789
790         return mesh_topo_store.index_lookup[index];
791 }
792
793 int mesh_get_x_mirror_vert__real_mesh(Object *ob, Mesh *mesh, int index, const bool use_topology)
794 {
795         if (use_topology) {
796                 return mesh_get_x_mirror_vert_topo__real_mesh(ob, mesh, index);
797         }
798         else {
799                 return mesh_get_x_mirror_vert_spatial__real_mesh(ob, mesh, index);
800         }
801 }
802
803 static BMVert *editbmesh_get_x_mirror_vert_spatial(Object *ob, BMEditMesh *em, const float co[3])
804 {
805         float vec[3];
806         int i;
807
808         /* ignore nan verts */
809         if ((isfinite(co[0]) == false) ||
810             (isfinite(co[1]) == false) ||
811             (isfinite(co[2]) == false))
812         {
813                 return NULL;
814         }
815
816         vec[0] = -co[0];
817         vec[1] = co[1];
818         vec[2] = co[2];
819
820         i = ED_mesh_mirror_spatial_table(ob, em, NULL, vec, 'u');
821         if (i != -1) {
822                 return BM_vert_at_index(em->bm, i);
823         }
824         return NULL;
825 }
826
827 static BMVert *editbmesh_get_x_mirror_vert_topo(Object *ob, struct BMEditMesh *em, BMVert *eve, int index)
828 {
829         intptr_t poinval;
830         if (ED_mesh_mirror_topo_table(ob, NULL, 'u') == -1)
831                 return NULL;
832
833         if (index == -1) {
834                 BMIter iter;
835                 BMVert *v;
836
837                 index = 0;
838                 BM_ITER_MESH (v, &iter, em->bm, BM_VERTS_OF_MESH) {
839                         if (v == eve)
840                                 break;
841                         index++;
842                 }
843
844                 if (index == em->bm->totvert) {
845                         return NULL;
846                 }
847         }
848
849         poinval = mesh_topo_store.index_lookup[index];
850
851         if (poinval != -1)
852                 return (BMVert *)(poinval);
853         return NULL;
854 }
855
856 BMVert *editbmesh_get_x_mirror_vert(Object *ob, struct BMEditMesh *em, BMVert *eve, const float co[3], int index, const bool use_topology)
857 {
858         if (use_topology) {
859                 return editbmesh_get_x_mirror_vert_topo(ob, em, eve, index);
860         }
861         else {
862                 return editbmesh_get_x_mirror_vert_spatial(ob, em, co);
863         }
864 }
865
866 /**
867  * Wrapper for objectmode/editmode.
868  *
869  * call #BM_mesh_elem_table_ensure first for editmesh.
870  */
871 int ED_mesh_mirror_get_vert(Object *ob, int index)
872 {
873         Mesh *me = ob->data;
874         BMEditMesh *em = me->edit_btmesh;
875         bool use_topology = (me->editflag & ME_EDIT_MIRROR_TOPO) != 0;
876         int index_mirr;
877
878         if (em) {
879                 BMVert *eve, *eve_mirr;
880                 eve = BM_vert_at_index(em->bm, index);
881                 eve_mirr = editbmesh_get_x_mirror_vert(ob, em, eve, eve->co, index, use_topology);
882                 index_mirr = eve_mirr ? BM_elem_index_get(eve_mirr) : -1;
883         }
884         else {
885                 index_mirr = mesh_get_x_mirror_vert(ob, NULL, index, use_topology);
886         }
887
888         return index_mirr;
889 }
890
891 #if 0
892
893 static float *editmesh_get_mirror_uv(BMEditMesh *em, int axis, float *uv, float *mirrCent, float *face_cent)
894 {
895         float vec[2];
896         float cent_vec[2];
897         float cent[2];
898
899         /* ignore nan verts */
900         if (isnan(uv[0]) || !isfinite(uv[0]) ||
901             isnan(uv[1]) || !isfinite(uv[1])
902             )
903         {
904                 return NULL;
905         }
906
907         if (axis) {
908                 vec[0] = uv[0];
909                 vec[1] = -((uv[1]) - mirrCent[1]) + mirrCent[1];
910
911                 cent_vec[0] = face_cent[0];
912                 cent_vec[1] = -((face_cent[1]) - mirrCent[1]) + mirrCent[1];
913         }
914         else {
915                 vec[0] = -((uv[0]) - mirrCent[0]) + mirrCent[0];
916                 vec[1] = uv[1];
917
918                 cent_vec[0] = -((face_cent[0]) - mirrCent[0]) + mirrCent[0];
919                 cent_vec[1] = face_cent[1];
920         }
921
922         /* TODO - Optimize */
923         {
924                 BMIter iter;
925                 BMFace *efa;
926
927                 BM_ITER_MESH (efa, &iter, em->bm, BM_FACES_OF_MESH) {
928                         uv_poly_center(efa, cent, cd_loop_uv_offset);
929
930                         if ( (fabsf(cent[0] - cent_vec[0]) < 0.001f) && (fabsf(cent[1] - cent_vec[1]) < 0.001f) ) {
931                                 BMIter liter;
932                                 BMLoop *l;
933
934                                 BM_ITER_ELEM (l, &liter, efa, BM_LOOPS_OF_FACE) {
935                                         MLoopUV *luv = CustomData_bmesh_get(&em->bm->ldata, l->head.data, CD_MLOOPUV);
936                                         if ( (fabsf(luv->uv[0] - vec[0]) < 0.001f) && (fabsf(luv->uv[1] - vec[1]) < 0.001f) ) {
937                                                 return luv->uv;
938
939                                         }
940                                 }
941                         }
942                 }
943         }
944
945         return NULL;
946 }
947
948 #endif
949
950 static unsigned int mirror_facehash(const void *ptr)
951 {
952         const MFace *mf = ptr;
953         unsigned int v0, v1;
954
955         if (mf->v4) {
956                 v0 = MIN4(mf->v1, mf->v2, mf->v3, mf->v4);
957                 v1 = MAX4(mf->v1, mf->v2, mf->v3, mf->v4);
958         }
959         else {
960                 v0 = MIN3(mf->v1, mf->v2, mf->v3);
961                 v1 = MAX3(mf->v1, mf->v2, mf->v3);
962         }
963
964         return ((v0 * 39) ^ (v1 * 31));
965 }
966
967 static int mirror_facerotation(MFace *a, MFace *b)
968 {
969         if (b->v4) {
970                 if (a->v1 == b->v1 && a->v2 == b->v2 && a->v3 == b->v3 && a->v4 == b->v4)
971                         return 0;
972                 else if (a->v4 == b->v1 && a->v1 == b->v2 && a->v2 == b->v3 && a->v3 == b->v4)
973                         return 1;
974                 else if (a->v3 == b->v1 && a->v4 == b->v2 && a->v1 == b->v3 && a->v2 == b->v4)
975                         return 2;
976                 else if (a->v2 == b->v1 && a->v3 == b->v2 && a->v4 == b->v3 && a->v1 == b->v4)
977                         return 3;
978         }
979         else {
980                 if (a->v1 == b->v1 && a->v2 == b->v2 && a->v3 == b->v3)
981                         return 0;
982                 else if (a->v3 == b->v1 && a->v1 == b->v2 && a->v2 == b->v3)
983                         return 1;
984                 else if (a->v2 == b->v1 && a->v3 == b->v2 && a->v1 == b->v3)
985                         return 2;
986         }
987
988         return -1;
989 }
990
991 static bool mirror_facecmp(const void *a, const void *b)
992 {
993         return (mirror_facerotation((MFace *)a, (MFace *)b) == -1);
994 }
995
996 /* BMESH_TODO, convert to MPoly (functions above also) */
997 int *mesh_get_x_mirror_faces(Object *ob, BMEditMesh *em, DerivedMesh *dm)
998 {
999         Mesh *me = ob->data;
1000         MVert *mv, *mvert;
1001         MFace mirrormf, *mf, *hashmf, *mface;
1002         GHash *fhash;
1003         int *mirrorverts, *mirrorfaces;
1004
1005         BLI_assert(em == NULL);  /* Does not work otherwise, currently... */
1006
1007         const bool use_topology = (me->editflag & ME_EDIT_MIRROR_TOPO) != 0;
1008         const int totvert = dm ? dm->getNumVerts(dm) : me->totvert;
1009         const int totface = dm ? dm->getNumTessFaces(dm) : me->totface;
1010         int a;
1011
1012         mirrorverts = MEM_callocN(sizeof(int) * totvert, "MirrorVerts");
1013         mirrorfaces = MEM_callocN(sizeof(int) * 2 * totface, "MirrorFaces");
1014
1015         mvert = dm ? dm->getVertArray(dm) : me->mvert;
1016         mface = dm ? dm->getTessFaceArray(dm) : me->mface;
1017
1018         ED_mesh_mirror_spatial_table(ob, em, dm, NULL, 's');
1019
1020         for (a = 0, mv = mvert; a < totvert; a++, mv++)
1021                 mirrorverts[a] = mesh_get_x_mirror_vert(ob, dm, a, use_topology);
1022
1023         ED_mesh_mirror_spatial_table(ob, em, dm, NULL, 'e');
1024
1025         fhash = BLI_ghash_new_ex(mirror_facehash, mirror_facecmp, "mirror_facehash gh", me->totface);
1026         for (a = 0, mf = mface; a < totface; a++, mf++)
1027                 BLI_ghash_insert(fhash, mf, mf);
1028
1029         for (a = 0, mf = mface; a < totface; a++, mf++) {
1030                 mirrormf.v1 = mirrorverts[mf->v3];
1031                 mirrormf.v2 = mirrorverts[mf->v2];
1032                 mirrormf.v3 = mirrorverts[mf->v1];
1033                 mirrormf.v4 = (mf->v4) ? mirrorverts[mf->v4] : 0;
1034
1035                 /* make sure v4 is not 0 if a quad */
1036                 if (mf->v4 && mirrormf.v4 == 0) {
1037                         SWAP(unsigned int, mirrormf.v1, mirrormf.v3);
1038                         SWAP(unsigned int, mirrormf.v2, mirrormf.v4);
1039                 }
1040
1041                 hashmf = BLI_ghash_lookup(fhash, &mirrormf);
1042                 if (hashmf) {
1043                         mirrorfaces[a * 2] = hashmf - mface;
1044                         mirrorfaces[a * 2 + 1] = mirror_facerotation(&mirrormf, hashmf);
1045                 }
1046                 else
1047                         mirrorfaces[a * 2] = -1;
1048         }
1049
1050         BLI_ghash_free(fhash, NULL, NULL);
1051         MEM_freeN(mirrorverts);
1052
1053         return mirrorfaces;
1054 }
1055
1056 /* This is a Mesh-based copy of mesh_get_x_mirror_faces() */
1057 int *mesh_get_x_mirror_faces__real_mesh(Object *ob, BMEditMesh *em, Mesh *mesh)
1058 {
1059         Mesh *me = ob->data;
1060         MVert *mv, *mvert;
1061         MFace mirrormf, *mf, *hashmf, *mface;
1062         GHash *fhash;
1063         int *mirrorverts, *mirrorfaces;
1064
1065         BLI_assert(em == NULL);  /* Does not work otherwise, currently... */
1066
1067         const bool use_topology = (me->editflag & ME_EDIT_MIRROR_TOPO) != 0;
1068         const int totvert = mesh ? mesh->totvert : me->totvert;
1069         const int totface = mesh ? mesh->totface : me->totface;
1070         int a;
1071
1072         mirrorverts = MEM_callocN(sizeof(int) * totvert, "MirrorVerts");
1073         mirrorfaces = MEM_callocN(sizeof(int) * 2 * totface, "MirrorFaces");
1074
1075         mvert = mesh ? mesh->mvert : me->mvert;
1076         mface = mesh ? mesh->mface : me->mface;
1077
1078         ED_mesh_mirror_spatial_table__real_mesh(ob, em, mesh, NULL, 's');
1079
1080         for (a = 0, mv = mvert; a < totvert; a++, mv++)
1081                 mirrorverts[a] = mesh_get_x_mirror_vert__real_mesh(ob, mesh, a, use_topology);
1082
1083         ED_mesh_mirror_spatial_table__real_mesh(ob, em, mesh, NULL, 'e');
1084
1085         fhash = BLI_ghash_new_ex(mirror_facehash, mirror_facecmp, "mirror_facehash gh", me->totface);
1086         for (a = 0, mf = mface; a < totface; a++, mf++)
1087                 BLI_ghash_insert(fhash, mf, mf);
1088
1089         for (a = 0, mf = mface; a < totface; a++, mf++) {
1090                 mirrormf.v1 = mirrorverts[mf->v3];
1091                 mirrormf.v2 = mirrorverts[mf->v2];
1092                 mirrormf.v3 = mirrorverts[mf->v1];
1093                 mirrormf.v4 = (mf->v4) ? mirrorverts[mf->v4] : 0;
1094
1095                 /* make sure v4 is not 0 if a quad */
1096                 if (mf->v4 && mirrormf.v4 == 0) {
1097                         SWAP(unsigned int, mirrormf.v1, mirrormf.v3);
1098                         SWAP(unsigned int, mirrormf.v2, mirrormf.v4);
1099                 }
1100
1101                 hashmf = BLI_ghash_lookup(fhash, &mirrormf);
1102                 if (hashmf) {
1103                         mirrorfaces[a * 2] = hashmf - mface;
1104                         mirrorfaces[a * 2 + 1] = mirror_facerotation(&mirrormf, hashmf);
1105                 }
1106                 else
1107                         mirrorfaces[a * 2] = -1;
1108         }
1109
1110         BLI_ghash_free(fhash, NULL, NULL);
1111         MEM_freeN(mirrorverts);
1112
1113         return mirrorfaces;
1114 }
1115
1116 /* selection, vertex and face */
1117 /* returns 0 if not found, otherwise 1 */
1118
1119 /**
1120  * Face selection in object mode,
1121  * currently only weight-paint and vertex-paint use this.
1122  *
1123  * \return boolean true == Found
1124  */
1125 bool ED_mesh_pick_face(bContext *C, Object *ob, const int mval[2], unsigned int *index, int size)
1126 {
1127         ViewContext vc;
1128         Mesh *me = ob->data;
1129
1130         BLI_assert(me && GS(me->id.name) == ID_ME);
1131
1132         if (!me || me->totpoly == 0)
1133                 return false;
1134
1135         ED_view3d_viewcontext_init(C, &vc);
1136
1137         if (size) {
1138                 /* sample rect to increase chances of selecting, so that when clicking
1139                  * on an edge in the backbuf, we can still select a face */
1140
1141                 float dummy_dist;
1142                 *index = ED_view3d_backbuf_sample_rect(&vc, mval, size, 1, me->totpoly + 1, &dummy_dist);
1143         }
1144         else {
1145                 /* sample only on the exact position */
1146                 *index = ED_view3d_backbuf_sample(&vc, mval[0], mval[1]);
1147         }
1148
1149         if ((*index) == 0 || (*index) > (unsigned int)me->totpoly)
1150                 return false;
1151
1152         (*index)--;
1153
1154         return true;
1155 }
1156 static void ed_mesh_pick_face_vert__mpoly_find(
1157         /* context */
1158         struct ARegion *ar, const float mval[2],
1159         /* mesh data (evaluated) */
1160         const MPoly *mp,
1161         const MVert *mvert, const MLoop *mloop,
1162         /* return values */
1163         float *r_len_best, int *r_v_idx_best)
1164 {
1165         const MLoop *ml;
1166         int j = mp->totloop;
1167         for (ml = &mloop[mp->loopstart]; j--; ml++) {
1168                 float sco[2];
1169                 const int v_idx = ml->v;
1170                 const float *co = mvert[v_idx].co;
1171                 if (ED_view3d_project_float_object(ar, co, sco, V3D_PROJ_TEST_NOP) == V3D_PROJ_RET_OK) {
1172                         const float len_test = len_manhattan_v2v2(mval, sco);
1173                         if (len_test < *r_len_best) {
1174                                 *r_len_best = len_test;
1175                                 *r_v_idx_best = v_idx;
1176                         }
1177                 }
1178         }
1179 }
1180 /**
1181  * Use when the back buffer stores face index values. but we want a vert.
1182  * This gets the face then finds the closest vertex to mval.
1183  */
1184 bool ED_mesh_pick_face_vert(bContext *C, Object *ob, const int mval[2], unsigned int *index, int size)
1185 {
1186         Depsgraph *depsgraph = CTX_data_depsgraph(C);
1187         unsigned int poly_index;
1188         Mesh *me = ob->data;
1189
1190         BLI_assert(me && GS(me->id.name) == ID_ME);
1191
1192         if (ED_mesh_pick_face(C, ob, mval, &poly_index, size)) {
1193                 Scene *scene = CTX_data_scene(C);
1194                 struct ARegion *ar = CTX_wm_region(C);
1195
1196                 /* derived mesh to find deformed locations */
1197                 Mesh *me_eval = mesh_get_eval_final(depsgraph, scene, ob, CD_MASK_BAREMESH | CD_MASK_ORIGINDEX);
1198
1199                 int v_idx_best = ORIGINDEX_NONE;
1200
1201                 /* find the vert closest to 'mval' */
1202                 const float mval_f[2] = {UNPACK2(mval)};
1203                 float len_best = FLT_MAX;
1204
1205                 MPoly *me_eval_mpoly;
1206                 MLoop *me_eval_mloop;
1207                 MVert *me_eval_mvert;
1208                 unsigned int me_eval_mpoly_len;
1209                 const int *index_mp_to_orig;
1210
1211                 me_eval_mpoly = me_eval->mpoly;
1212                 me_eval_mloop = me_eval->mloop;
1213                 me_eval_mvert = me_eval->mvert;
1214
1215                 me_eval_mpoly_len = me_eval->totpoly;
1216
1217                 index_mp_to_orig = CustomData_get_layer(&me_eval->pdata, CD_ORIGINDEX);
1218
1219                 /* tag all verts using this face */
1220                 if (index_mp_to_orig) {
1221                         unsigned int i;
1222
1223                         for (i = 0; i < me_eval_mpoly_len; i++) {
1224                                 if (index_mp_to_orig[i] == poly_index) {
1225                                         ed_mesh_pick_face_vert__mpoly_find(
1226                                                 ar, mval_f,
1227                                                 &me_eval_mpoly[i], me_eval_mvert, me_eval_mloop,
1228                                                 &len_best, &v_idx_best);
1229                                 }
1230                         }
1231                 }
1232                 else {
1233                         if (poly_index < me_eval_mpoly_len) {
1234                                 ed_mesh_pick_face_vert__mpoly_find(
1235                                         ar, mval_f,
1236                                         &me_eval_mpoly[poly_index], me_eval_mvert, me_eval_mloop,
1237                                         &len_best, &v_idx_best);
1238                         }
1239                 }
1240
1241                 /* map 'dm -> me' index if possible */
1242                 if (v_idx_best != ORIGINDEX_NONE) {
1243                         const int *index_mv_to_orig;
1244                         index_mv_to_orig = CustomData_get_layer(&me_eval->vdata, CD_ORIGINDEX);
1245                         if (index_mv_to_orig) {
1246                                 v_idx_best = index_mv_to_orig[v_idx_best];
1247                         }
1248                 }
1249
1250                 if ((v_idx_best != ORIGINDEX_NONE) && (v_idx_best < me->totvert)) {
1251                         *index = v_idx_best;
1252                         return true;
1253                 }
1254         }
1255
1256         return false;
1257 }
1258
1259 /**
1260  * Vertex selection in object mode,
1261  * currently only weight paint uses this.
1262  *
1263  * \return boolean true == Found
1264  */
1265 typedef struct VertPickData {
1266         const MVert *mvert;
1267         const float *mval_f;  /* [2] */
1268         ARegion *ar;
1269
1270         /* runtime */
1271         float len_best;
1272         int v_idx_best;
1273 } VertPickData;
1274
1275 static void ed_mesh_pick_vert__mapFunc(void *userData, int index, const float co[3],
1276                                        const float UNUSED(no_f[3]), const short UNUSED(no_s[3]))
1277 {
1278         VertPickData *data = userData;
1279         if ((data->mvert[index].flag & ME_HIDE) == 0) {
1280                 float sco[2];
1281
1282                 if (ED_view3d_project_float_object(data->ar, co, sco, V3D_PROJ_TEST_CLIP_DEFAULT) == V3D_PROJ_RET_OK) {
1283                         const float len = len_manhattan_v2v2(data->mval_f, sco);
1284                         if (len < data->len_best) {
1285                                 data->len_best = len;
1286                                 data->v_idx_best = index;
1287                         }
1288                 }
1289         }
1290 }
1291 bool ED_mesh_pick_vert(bContext *C, Object *ob, const int mval[2], unsigned int *index, int size, bool use_zbuf)
1292 {
1293         ViewContext vc;
1294         Mesh *me = ob->data;
1295
1296         BLI_assert(me && GS(me->id.name) == ID_ME);
1297
1298         if (!me || me->totvert == 0)
1299                 return false;
1300
1301         ED_view3d_viewcontext_init(C, &vc);
1302
1303         if (use_zbuf) {
1304                 if (size > 0) {
1305                         /* sample rect to increase chances of selecting, so that when clicking
1306                          * on an face in the backbuf, we can still select a vert */
1307
1308                         float dummy_dist;
1309                         *index = ED_view3d_backbuf_sample_rect(&vc, mval, size, 1, me->totvert + 1, &dummy_dist);
1310                 }
1311                 else {
1312                         /* sample only on the exact position */
1313                         *index = ED_view3d_backbuf_sample(&vc, mval[0], mval[1]);
1314                 }
1315
1316                 if ((*index) == 0 || (*index) > (unsigned int)me->totvert)
1317                         return false;
1318
1319                 (*index)--;
1320         }
1321         else {
1322                 /* derived mesh to find deformed locations */
1323                 DerivedMesh *dm = mesh_get_derived_final(vc.depsgraph, vc.scene, ob, CD_MASK_BAREMESH);
1324                 ARegion *ar = vc.ar;
1325                 RegionView3D *rv3d = ar->regiondata;
1326
1327                 /* find the vert closest to 'mval' */
1328                 const float mval_f[2] = {(float)mval[0],
1329                                          (float)mval[1]};
1330
1331                 VertPickData data = {NULL};
1332
1333                 ED_view3d_init_mats_rv3d(ob, rv3d);
1334
1335                 if (dm == NULL) {
1336                         return false;
1337                 }
1338
1339                 /* setup data */
1340                 data.mvert = me->mvert;
1341                 data.ar = ar;
1342                 data.mval_f = mval_f;
1343                 data.len_best = FLT_MAX;
1344                 data.v_idx_best = -1;
1345
1346                 dm->foreachMappedVert(dm, ed_mesh_pick_vert__mapFunc, &data, DM_FOREACH_NOP);
1347
1348                 dm->release(dm);
1349
1350                 if (data.v_idx_best == -1) {
1351                         return false;
1352                 }
1353
1354                 *index = data.v_idx_best;
1355         }
1356
1357         return true;
1358 }
1359
1360
1361 MDeformVert *ED_mesh_active_dvert_get_em(Object *ob, BMVert **r_eve)
1362 {
1363         if (ob->mode & OB_MODE_EDIT && ob->type == OB_MESH && ob->defbase.first) {
1364                 Mesh *me = ob->data;
1365                 BMesh *bm = me->edit_btmesh->bm;
1366                 const int cd_dvert_offset = CustomData_get_offset(&bm->vdata, CD_MDEFORMVERT);
1367
1368                 if (cd_dvert_offset != -1) {
1369                         BMVert *eve = BM_mesh_active_vert_get(bm);
1370
1371                         if (eve) {
1372                                 if (r_eve) *r_eve = eve;
1373                                 return BM_ELEM_CD_GET_VOID_P(eve, cd_dvert_offset);
1374                         }
1375                 }
1376         }
1377
1378         if (r_eve) *r_eve = NULL;
1379         return NULL;
1380 }
1381
1382 MDeformVert *ED_mesh_active_dvert_get_ob(Object *ob, int *r_index)
1383 {
1384         Mesh *me = ob->data;
1385         int index = BKE_mesh_mselect_active_get(me, ME_VSEL);
1386         if (r_index) *r_index = index;
1387         if (index == -1 || me->dvert == NULL) {
1388                 return NULL;
1389         }
1390         else {
1391                 return me->dvert + index;
1392         }
1393 }
1394
1395 MDeformVert *ED_mesh_active_dvert_get_only(Object *ob)
1396 {
1397         if (ob->type == OB_MESH) {
1398                 if (ob->mode & OB_MODE_EDIT) {
1399                         return ED_mesh_active_dvert_get_em(ob, NULL);
1400                 }
1401                 else {
1402                         return ED_mesh_active_dvert_get_ob(ob, NULL);
1403                 }
1404         }
1405         else {
1406                 return NULL;
1407         }
1408 }
1409
1410 void EDBM_mesh_stats_multi(
1411         struct Object **objects, const uint objects_len,
1412         int totelem[3], int totelem_sel[3])
1413 {
1414         if (totelem) {
1415                 totelem[0] = 0;
1416                 totelem[1] = 0;
1417                 totelem[2] = 0;
1418         }
1419         if (totelem_sel) {
1420                 totelem_sel[0] = 0;
1421                 totelem_sel[1] = 0;
1422                 totelem_sel[2] = 0;
1423         }
1424
1425         for (uint ob_index = 0; ob_index < objects_len; ob_index++) {
1426                 Object *obedit = objects[ob_index];
1427                 BMEditMesh *em = BKE_editmesh_from_object(obedit);
1428                 BMesh *bm = em->bm;
1429                 if (totelem) {
1430                         totelem[0] += bm->totvert;
1431                         totelem[1] += bm->totedge;
1432                         totelem[2] += bm->totface;
1433                 }
1434                 if (totelem_sel) {
1435                         totelem_sel[0] += bm->totvertsel;
1436                         totelem_sel[1] += bm->totedgesel;
1437                         totelem_sel[2] += bm->totfacesel;
1438                 }
1439         }
1440 }
1441
1442
1443 void EDBM_mesh_elem_index_ensure_multi(Object **objects, const uint objects_len, const char htype)
1444 {
1445         int elem_offset[4] = {0, 0, 0, 0};
1446         for (uint ob_index = 0; ob_index < objects_len; ob_index++) {
1447                 Object *obedit = objects[ob_index];
1448                 BMEditMesh *em = BKE_editmesh_from_object(obedit);
1449                 BMesh *bm = em->bm;
1450                 BM_mesh_elem_index_ensure_ex(bm, htype, elem_offset);
1451         }
1452 }