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((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         Scene *scene = CTX_data_scene(C);
609         Object *ob = CTX_data_active_object(C);
610         Depsgraph *depsgraph = CTX_data_depsgraph(C);
611         Mesh *me = (Mesh *)ob->data;
612         Mesh *selme = NULL;
613         DerivedMesh *dm = NULL;
614         Key *key = me->key;
615         KeyBlock *kb;
616         bool ok = false, nonequal_verts = false;
617
618         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
619         {
620                 if (base->object == ob) continue;
621
622                 if (base->object->type == OB_MESH) {
623                         selme = (Mesh *)base->object->data;
624
625                         if (selme->totvert == me->totvert)
626                                 ok = true;
627                         else
628                                 nonequal_verts = 1;
629                 }
630         }
631         CTX_DATA_END;
632
633         if (!ok) {
634                 if (nonequal_verts)
635                         BKE_report(op->reports, RPT_WARNING, "Selected meshes must have equal numbers of vertices");
636                 else
637                         BKE_report(op->reports, RPT_WARNING, "No additional selected meshes with equal vertex count to join");
638                 return OPERATOR_CANCELLED;
639         }
640
641         if (key == NULL) {
642                 key = me->key = BKE_key_add((ID *)me);
643                 key->type = KEY_RELATIVE;
644
645                 /* first key added, so it was the basis. initialize it with the existing mesh */
646                 kb = BKE_keyblock_add(key, NULL);
647                 BKE_keyblock_convert_from_mesh(me, kb);
648         }
649
650         /* now ready to add new keys from selected meshes */
651         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
652         {
653                 if (base->object == ob) continue;
654
655                 if (base->object->type == OB_MESH) {
656                         selme = (Mesh *)base->object->data;
657
658                         if (selme->totvert == me->totvert) {
659                                 dm = mesh_get_derived_deform(depsgraph, scene, base->object, CD_MASK_BAREMESH);
660
661                                 if (!dm) continue;
662
663                                 kb = BKE_keyblock_add(key, base->object->id.name + 2);
664
665                                 DM_to_meshkey(dm, me, kb);
666
667                                 dm->release(dm);
668                         }
669                 }
670         }
671         CTX_DATA_END;
672
673         WM_event_add_notifier(C, NC_SCENE | ND_OB_ACTIVE, scene);
674
675         return OPERATOR_FINISHED;
676 }
677
678
679 /* -------------------------------------------------------------------- */
680 /* Mesh Mirror (Topology) */
681
682 /** \name Mesh Topology Mirror API
683  * \{ */
684
685 static MirrTopoStore_t mesh_topo_store = {NULL, -1. - 1, -1};
686
687 /* mode is 's' start, or 'e' end, or 'u' use */
688 /* if end, ob can be NULL */
689 /* note, is supposed return -1 on error, which callers are currently checking for, but is not used so far */
690 int ED_mesh_mirror_topo_table(
691         Object *ob, DerivedMesh *dm, char mode)
692 {
693         if (mode == 'u') {        /* use table */
694                 if (ED_mesh_mirrtopo_recalc_check(ob->data, dm, &mesh_topo_store)) {
695                         ED_mesh_mirror_topo_table(ob, dm, 's');
696                 }
697         }
698         else if (mode == 's') { /* start table */
699                 ED_mesh_mirrtopo_init(ob->data, dm, &mesh_topo_store, false);
700         }
701         else if (mode == 'e') { /* end table */
702                 ED_mesh_mirrtopo_free(&mesh_topo_store);
703         }
704         else {
705                 BLI_assert(0);
706         }
707
708         return 0;
709 }
710
711 /* mode is 's' start, or 'e' end, or 'u' use */
712 /* if end, ob can be NULL */
713 /* note, is supposed return -1 on error, which callers are currently checking for, but is not used so far */
714 int ED_mesh_mirror_topo_table__real_mesh(
715         Object *ob, Mesh *mesh, char mode)
716 {
717         if (mode == 'u') {        /* use table */
718                 if (ED_mesh_mirrtopo_recalc_check__real_mesh(ob->data, mesh, &mesh_topo_store)) {
719                         ED_mesh_mirror_topo_table__real_mesh(ob, mesh, 's');
720                 }
721         }
722         else if (mode == 's') { /* start table */
723                 ED_mesh_mirrtopo_init__real_mesh(ob->data, mesh, &mesh_topo_store, false);
724         }
725         else if (mode == 'e') { /* end table */
726                 ED_mesh_mirrtopo_free(&mesh_topo_store);
727         }
728         else {
729                 BLI_assert(0);
730         }
731
732         return 0;
733 }
734
735 /** \} */
736
737
738 static int mesh_get_x_mirror_vert_spatial(Object *ob, DerivedMesh *dm, int index)
739 {
740         Mesh *me = ob->data;
741         MVert *mvert = dm ? dm->getVertArray(dm) : me->mvert;
742         float vec[3];
743
744         mvert = &mvert[index];
745         vec[0] = -mvert->co[0];
746         vec[1] = mvert->co[1];
747         vec[2] = mvert->co[2];
748
749         return ED_mesh_mirror_spatial_table(ob, NULL, dm, vec, 'u');
750 }
751
752 static int mesh_get_x_mirror_vert_topo(Object *ob, DerivedMesh *dm, int index)
753 {
754         if (ED_mesh_mirror_topo_table(ob, dm, 'u') == -1)
755                 return -1;
756
757         return mesh_topo_store.index_lookup[index];
758 }
759
760 int mesh_get_x_mirror_vert(Object *ob, DerivedMesh *dm, int index, const bool use_topology)
761 {
762         if (use_topology) {
763                 return mesh_get_x_mirror_vert_topo(ob, dm, index);
764         }
765         else {
766                 return mesh_get_x_mirror_vert_spatial(ob, dm, index);
767         }
768 }
769
770 static int mesh_get_x_mirror_vert_spatial__real_mesh(Object *ob, Mesh *mesh, int index)
771 {
772         Mesh *me = ob->data;
773         MVert *mvert = mesh ? mesh->mvert : me->mvert;
774         float vec[3];
775
776         mvert = &mvert[index];
777         vec[0] = -mvert->co[0];
778         vec[1] = mvert->co[1];
779         vec[2] = mvert->co[2];
780
781         return ED_mesh_mirror_spatial_table__real_mesh(ob, NULL, mesh, vec, 'u');
782 }
783
784 static int mesh_get_x_mirror_vert_topo__real_mesh(Object *ob, Mesh *mesh, int index)
785 {
786         if (ED_mesh_mirror_topo_table__real_mesh(ob, mesh, 'u') == -1)
787                 return -1;
788
789         return mesh_topo_store.index_lookup[index];
790 }
791
792 int mesh_get_x_mirror_vert__real_mesh(Object *ob, Mesh *mesh, int index, const bool use_topology)
793 {
794         if (use_topology) {
795                 return mesh_get_x_mirror_vert_topo__real_mesh(ob, mesh, index);
796         }
797         else {
798                 return mesh_get_x_mirror_vert_spatial__real_mesh(ob, mesh, index);
799         }
800 }
801
802 static BMVert *editbmesh_get_x_mirror_vert_spatial(Object *ob, BMEditMesh *em, const float co[3])
803 {
804         float vec[3];
805         int i;
806
807         /* ignore nan verts */
808         if ((isfinite(co[0]) == false) ||
809             (isfinite(co[1]) == false) ||
810             (isfinite(co[2]) == false))
811         {
812                 return NULL;
813         }
814
815         vec[0] = -co[0];
816         vec[1] = co[1];
817         vec[2] = co[2];
818
819         i = ED_mesh_mirror_spatial_table(ob, em, NULL, vec, 'u');
820         if (i != -1) {
821                 return BM_vert_at_index(em->bm, i);
822         }
823         return NULL;
824 }
825
826 static BMVert *editbmesh_get_x_mirror_vert_topo(Object *ob, struct BMEditMesh *em, BMVert *eve, int index)
827 {
828         intptr_t poinval;
829         if (ED_mesh_mirror_topo_table(ob, NULL, 'u') == -1)
830                 return NULL;
831
832         if (index == -1) {
833                 BMIter iter;
834                 BMVert *v;
835
836                 index = 0;
837                 BM_ITER_MESH (v, &iter, em->bm, BM_VERTS_OF_MESH) {
838                         if (v == eve)
839                                 break;
840                         index++;
841                 }
842
843                 if (index == em->bm->totvert) {
844                         return NULL;
845                 }
846         }
847
848         poinval = mesh_topo_store.index_lookup[index];
849
850         if (poinval != -1)
851                 return (BMVert *)(poinval);
852         return NULL;
853 }
854
855 BMVert *editbmesh_get_x_mirror_vert(Object *ob, struct BMEditMesh *em, BMVert *eve, const float co[3], int index, const bool use_topology)
856 {
857         if (use_topology) {
858                 return editbmesh_get_x_mirror_vert_topo(ob, em, eve, index);
859         }
860         else {
861                 return editbmesh_get_x_mirror_vert_spatial(ob, em, co);
862         }
863 }
864
865 /**
866  * Wrapper for objectmode/editmode.
867  *
868  * call #BM_mesh_elem_table_ensure first for editmesh.
869  */
870 int ED_mesh_mirror_get_vert(Object *ob, int index)
871 {
872         Mesh *me = ob->data;
873         BMEditMesh *em = me->edit_btmesh;
874         bool use_topology = (me->editflag & ME_EDIT_MIRROR_TOPO) != 0;
875         int index_mirr;
876
877         if (em) {
878                 BMVert *eve, *eve_mirr;
879                 eve = BM_vert_at_index(em->bm, index);
880                 eve_mirr = editbmesh_get_x_mirror_vert(ob, em, eve, eve->co, index, use_topology);
881                 index_mirr = eve_mirr ? BM_elem_index_get(eve_mirr) : -1;
882         }
883         else {
884                 index_mirr = mesh_get_x_mirror_vert(ob, NULL, index, use_topology);
885         }
886
887         return index_mirr;
888 }
889
890 #if 0
891
892 static float *editmesh_get_mirror_uv(BMEditMesh *em, int axis, float *uv, float *mirrCent, float *face_cent)
893 {
894         float vec[2];
895         float cent_vec[2];
896         float cent[2];
897
898         /* ignore nan verts */
899         if (isnan(uv[0]) || !isfinite(uv[0]) ||
900             isnan(uv[1]) || !isfinite(uv[1])
901             )
902         {
903                 return NULL;
904         }
905
906         if (axis) {
907                 vec[0] = uv[0];
908                 vec[1] = -((uv[1]) - mirrCent[1]) + mirrCent[1];
909
910                 cent_vec[0] = face_cent[0];
911                 cent_vec[1] = -((face_cent[1]) - mirrCent[1]) + mirrCent[1];
912         }
913         else {
914                 vec[0] = -((uv[0]) - mirrCent[0]) + mirrCent[0];
915                 vec[1] = uv[1];
916
917                 cent_vec[0] = -((face_cent[0]) - mirrCent[0]) + mirrCent[0];
918                 cent_vec[1] = face_cent[1];
919         }
920
921         /* TODO - Optimize */
922         {
923                 BMIter iter;
924                 BMFace *efa;
925
926                 BM_ITER_MESH (efa, &iter, em->bm, BM_FACES_OF_MESH) {
927                         uv_poly_center(efa, cent, cd_loop_uv_offset);
928
929                         if ( (fabsf(cent[0] - cent_vec[0]) < 0.001f) && (fabsf(cent[1] - cent_vec[1]) < 0.001f) ) {
930                                 BMIter liter;
931                                 BMLoop *l;
932
933                                 BM_ITER_ELEM (l, &liter, efa, BM_LOOPS_OF_FACE) {
934                                         MLoopUV *luv = CustomData_bmesh_get(&em->bm->ldata, l->head.data, CD_MLOOPUV);
935                                         if ( (fabsf(luv->uv[0] - vec[0]) < 0.001f) && (fabsf(luv->uv[1] - vec[1]) < 0.001f) ) {
936                                                 return luv->uv;
937
938                                         }
939                                 }
940                         }
941                 }
942         }
943
944         return NULL;
945 }
946
947 #endif
948
949 static unsigned int mirror_facehash(const void *ptr)
950 {
951         const MFace *mf = ptr;
952         unsigned int v0, v1;
953
954         if (mf->v4) {
955                 v0 = MIN4(mf->v1, mf->v2, mf->v3, mf->v4);
956                 v1 = MAX4(mf->v1, mf->v2, mf->v3, mf->v4);
957         }
958         else {
959                 v0 = MIN3(mf->v1, mf->v2, mf->v3);
960                 v1 = MAX3(mf->v1, mf->v2, mf->v3);
961         }
962
963         return ((v0 * 39) ^ (v1 * 31));
964 }
965
966 static int mirror_facerotation(MFace *a, MFace *b)
967 {
968         if (b->v4) {
969                 if (a->v1 == b->v1 && a->v2 == b->v2 && a->v3 == b->v3 && a->v4 == b->v4)
970                         return 0;
971                 else if (a->v4 == b->v1 && a->v1 == b->v2 && a->v2 == b->v3 && a->v3 == b->v4)
972                         return 1;
973                 else if (a->v3 == b->v1 && a->v4 == b->v2 && a->v1 == b->v3 && a->v2 == b->v4)
974                         return 2;
975                 else if (a->v2 == b->v1 && a->v3 == b->v2 && a->v4 == b->v3 && a->v1 == b->v4)
976                         return 3;
977         }
978         else {
979                 if (a->v1 == b->v1 && a->v2 == b->v2 && a->v3 == b->v3)
980                         return 0;
981                 else if (a->v3 == b->v1 && a->v1 == b->v2 && a->v2 == b->v3)
982                         return 1;
983                 else if (a->v2 == b->v1 && a->v3 == b->v2 && a->v1 == b->v3)
984                         return 2;
985         }
986
987         return -1;
988 }
989
990 static bool mirror_facecmp(const void *a, const void *b)
991 {
992         return (mirror_facerotation((MFace *)a, (MFace *)b) == -1);
993 }
994
995 /* BMESH_TODO, convert to MPoly (functions above also) */
996 int *mesh_get_x_mirror_faces(Object *ob, BMEditMesh *em, DerivedMesh *dm)
997 {
998         Mesh *me = ob->data;
999         MVert *mv, *mvert;
1000         MFace mirrormf, *mf, *hashmf, *mface;
1001         GHash *fhash;
1002         int *mirrorverts, *mirrorfaces;
1003
1004         BLI_assert(em == NULL);  /* Does not work otherwise, currently... */
1005
1006         const bool use_topology = (me->editflag & ME_EDIT_MIRROR_TOPO) != 0;
1007         const int totvert = dm ? dm->getNumVerts(dm) : me->totvert;
1008         const int totface = dm ? dm->getNumTessFaces(dm) : me->totface;
1009         int a;
1010
1011         mirrorverts = MEM_callocN(sizeof(int) * totvert, "MirrorVerts");
1012         mirrorfaces = MEM_callocN(sizeof(int) * 2 * totface, "MirrorFaces");
1013
1014         mvert = dm ? dm->getVertArray(dm) : me->mvert;
1015         mface = dm ? dm->getTessFaceArray(dm) : me->mface;
1016
1017         ED_mesh_mirror_spatial_table(ob, em, dm, NULL, 's');
1018
1019         for (a = 0, mv = mvert; a < totvert; a++, mv++)
1020                 mirrorverts[a] = mesh_get_x_mirror_vert(ob, dm, a, use_topology);
1021
1022         ED_mesh_mirror_spatial_table(ob, em, dm, NULL, 'e');
1023
1024         fhash = BLI_ghash_new_ex(mirror_facehash, mirror_facecmp, "mirror_facehash gh", me->totface);
1025         for (a = 0, mf = mface; a < totface; a++, mf++)
1026                 BLI_ghash_insert(fhash, mf, mf);
1027
1028         for (a = 0, mf = mface; a < totface; a++, mf++) {
1029                 mirrormf.v1 = mirrorverts[mf->v3];
1030                 mirrormf.v2 = mirrorverts[mf->v2];
1031                 mirrormf.v3 = mirrorverts[mf->v1];
1032                 mirrormf.v4 = (mf->v4) ? mirrorverts[mf->v4] : 0;
1033
1034                 /* make sure v4 is not 0 if a quad */
1035                 if (mf->v4 && mirrormf.v4 == 0) {
1036                         SWAP(unsigned int, mirrormf.v1, mirrormf.v3);
1037                         SWAP(unsigned int, mirrormf.v2, mirrormf.v4);
1038                 }
1039
1040                 hashmf = BLI_ghash_lookup(fhash, &mirrormf);
1041                 if (hashmf) {
1042                         mirrorfaces[a * 2] = hashmf - mface;
1043                         mirrorfaces[a * 2 + 1] = mirror_facerotation(&mirrormf, hashmf);
1044                 }
1045                 else
1046                         mirrorfaces[a * 2] = -1;
1047         }
1048
1049         BLI_ghash_free(fhash, NULL, NULL);
1050         MEM_freeN(mirrorverts);
1051
1052         return mirrorfaces;
1053 }
1054
1055 /* This is a Mesh-based copy of mesh_get_x_mirror_faces() */
1056 int *mesh_get_x_mirror_faces__real_mesh(Object *ob, BMEditMesh *em, Mesh *mesh)
1057 {
1058         Mesh *me = ob->data;
1059         MVert *mv, *mvert;
1060         MFace mirrormf, *mf, *hashmf, *mface;
1061         GHash *fhash;
1062         int *mirrorverts, *mirrorfaces;
1063
1064         BLI_assert(em == NULL);  /* Does not work otherwise, currently... */
1065
1066         const bool use_topology = (me->editflag & ME_EDIT_MIRROR_TOPO) != 0;
1067         const int totvert = mesh ? mesh->totvert : me->totvert;
1068         const int totface = mesh ? mesh->totface : me->totface;
1069         int a;
1070
1071         mirrorverts = MEM_callocN(sizeof(int) * totvert, "MirrorVerts");
1072         mirrorfaces = MEM_callocN(sizeof(int) * 2 * totface, "MirrorFaces");
1073
1074         mvert = mesh ? mesh->mvert : me->mvert;
1075         mface = mesh ? mesh->mface : me->mface;
1076
1077         ED_mesh_mirror_spatial_table__real_mesh(ob, em, mesh, NULL, 's');
1078
1079         for (a = 0, mv = mvert; a < totvert; a++, mv++)
1080                 mirrorverts[a] = mesh_get_x_mirror_vert__real_mesh(ob, mesh, a, use_topology);
1081
1082         ED_mesh_mirror_spatial_table__real_mesh(ob, em, mesh, NULL, 'e');
1083
1084         fhash = BLI_ghash_new_ex(mirror_facehash, mirror_facecmp, "mirror_facehash gh", me->totface);
1085         for (a = 0, mf = mface; a < totface; a++, mf++)
1086                 BLI_ghash_insert(fhash, mf, mf);
1087
1088         for (a = 0, mf = mface; a < totface; a++, mf++) {
1089                 mirrormf.v1 = mirrorverts[mf->v3];
1090                 mirrormf.v2 = mirrorverts[mf->v2];
1091                 mirrormf.v3 = mirrorverts[mf->v1];
1092                 mirrormf.v4 = (mf->v4) ? mirrorverts[mf->v4] : 0;
1093
1094                 /* make sure v4 is not 0 if a quad */
1095                 if (mf->v4 && mirrormf.v4 == 0) {
1096                         SWAP(unsigned int, mirrormf.v1, mirrormf.v3);
1097                         SWAP(unsigned int, mirrormf.v2, mirrormf.v4);
1098                 }
1099
1100                 hashmf = BLI_ghash_lookup(fhash, &mirrormf);
1101                 if (hashmf) {
1102                         mirrorfaces[a * 2] = hashmf - mface;
1103                         mirrorfaces[a * 2 + 1] = mirror_facerotation(&mirrormf, hashmf);
1104                 }
1105                 else
1106                         mirrorfaces[a * 2] = -1;
1107         }
1108
1109         BLI_ghash_free(fhash, NULL, NULL);
1110         MEM_freeN(mirrorverts);
1111
1112         return mirrorfaces;
1113 }
1114
1115 /* selection, vertex and face */
1116 /* returns 0 if not found, otherwise 1 */
1117
1118 /**
1119  * Face selection in object mode,
1120  * currently only weight-paint and vertex-paint use this.
1121  *
1122  * \return boolean true == Found
1123  */
1124 bool ED_mesh_pick_face(bContext *C, Object *ob, const int mval[2], unsigned int *index, int size)
1125 {
1126         ViewContext vc;
1127         Mesh *me = ob->data;
1128
1129         BLI_assert(me && GS(me->id.name) == ID_ME);
1130
1131         if (!me || me->totpoly == 0)
1132                 return false;
1133
1134         ED_view3d_viewcontext_init(C, &vc);
1135
1136         if (size) {
1137                 /* sample rect to increase chances of selecting, so that when clicking
1138                  * on an edge in the backbuf, we can still select a face */
1139
1140                 float dummy_dist;
1141                 *index = ED_view3d_backbuf_sample_rect(&vc, mval, size, 1, me->totpoly + 1, &dummy_dist);
1142         }
1143         else {
1144                 /* sample only on the exact position */
1145                 *index = ED_view3d_backbuf_sample(&vc, mval[0], mval[1]);
1146         }
1147
1148         if ((*index) == 0 || (*index) > (unsigned int)me->totpoly)
1149                 return false;
1150
1151         (*index)--;
1152
1153         return true;
1154 }
1155 static void ed_mesh_pick_face_vert__mpoly_find(
1156         /* context */
1157         struct ARegion *ar, const float mval[2],
1158         /* mesh data (evaluated) */
1159         const MPoly *mp,
1160         const MVert *mvert, const MLoop *mloop,
1161         /* return values */
1162         float *r_len_best, int *r_v_idx_best)
1163 {
1164         const MLoop *ml;
1165         int j = mp->totloop;
1166         for (ml = &mloop[mp->loopstart]; j--; ml++) {
1167                 float sco[2];
1168                 const int v_idx = ml->v;
1169                 const float *co = mvert[v_idx].co;
1170                 if (ED_view3d_project_float_object(ar, co, sco, V3D_PROJ_TEST_NOP) == V3D_PROJ_RET_OK) {
1171                         const float len_test = len_manhattan_v2v2(mval, sco);
1172                         if (len_test < *r_len_best) {
1173                                 *r_len_best = len_test;
1174                                 *r_v_idx_best = v_idx;
1175                         }
1176                 }
1177         }
1178 }
1179 /**
1180  * Use when the back buffer stores face index values. but we want a vert.
1181  * This gets the face then finds the closest vertex to mval.
1182  */
1183 bool ED_mesh_pick_face_vert(bContext *C, Object *ob, const int mval[2], unsigned int *index, int size)
1184 {
1185         Depsgraph *depsgraph = CTX_data_depsgraph(C);
1186         unsigned int poly_index;
1187         Mesh *me = ob->data;
1188
1189         BLI_assert(me && GS(me->id.name) == ID_ME);
1190
1191         if (ED_mesh_pick_face(C, ob, mval, &poly_index, size)) {
1192                 Scene *scene = CTX_data_scene(C);
1193                 struct ARegion *ar = CTX_wm_region(C);
1194
1195                 /* derived mesh to find deformed locations */
1196                 Mesh *me_eval = mesh_get_eval_final(depsgraph, scene, ob, CD_MASK_BAREMESH | CD_MASK_ORIGINDEX);
1197
1198                 int v_idx_best = ORIGINDEX_NONE;
1199
1200                 /* find the vert closest to 'mval' */
1201                 const float mval_f[2] = {UNPACK2(mval)};
1202                 float len_best = FLT_MAX;
1203
1204                 MPoly *me_eval_mpoly;
1205                 MLoop *me_eval_mloop;
1206                 MVert *me_eval_mvert;
1207                 unsigned int me_eval_mpoly_len;
1208                 const int *index_mp_to_orig;
1209
1210                 me_eval_mpoly = me_eval->mpoly;
1211                 me_eval_mloop = me_eval->mloop;
1212                 me_eval_mvert = me_eval->mvert;
1213
1214                 me_eval_mpoly_len = me_eval->totpoly;
1215
1216                 index_mp_to_orig = CustomData_get_layer(&me_eval->pdata, CD_ORIGINDEX);
1217
1218                 /* tag all verts using this face */
1219                 if (index_mp_to_orig) {
1220                         unsigned int i;
1221
1222                         for (i = 0; i < me_eval_mpoly_len; i++) {
1223                                 if (index_mp_to_orig[i] == poly_index) {
1224                                         ed_mesh_pick_face_vert__mpoly_find(
1225                                                 ar, mval_f,
1226                                                 &me_eval_mpoly[i], me_eval_mvert, me_eval_mloop,
1227                                                 &len_best, &v_idx_best);
1228                                 }
1229                         }
1230                 }
1231                 else {
1232                         if (poly_index < me_eval_mpoly_len) {
1233                                 ed_mesh_pick_face_vert__mpoly_find(
1234                                         ar, mval_f,
1235                                         &me_eval_mpoly[poly_index], me_eval_mvert, me_eval_mloop,
1236                                         &len_best, &v_idx_best);
1237                         }
1238                 }
1239
1240                 /* map 'dm -> me' index if possible */
1241                 if (v_idx_best != ORIGINDEX_NONE) {
1242                         const int *index_mv_to_orig;
1243                         index_mv_to_orig = CustomData_get_layer(&me_eval->vdata, CD_ORIGINDEX);
1244                         if (index_mv_to_orig) {
1245                                 v_idx_best = index_mv_to_orig[v_idx_best];
1246                         }
1247                 }
1248
1249                 if ((v_idx_best != ORIGINDEX_NONE) && (v_idx_best < me->totvert)) {
1250                         *index = v_idx_best;
1251                         return true;
1252                 }
1253         }
1254
1255         return false;
1256 }
1257
1258 /**
1259  * Vertex selection in object mode,
1260  * currently only weight paint uses this.
1261  *
1262  * \return boolean true == Found
1263  */
1264 typedef struct VertPickData {
1265         const MVert *mvert;
1266         const float *mval_f;  /* [2] */
1267         ARegion *ar;
1268
1269         /* runtime */
1270         float len_best;
1271         int v_idx_best;
1272 } VertPickData;
1273
1274 static void ed_mesh_pick_vert__mapFunc(void *userData, int index, const float co[3],
1275                                        const float UNUSED(no_f[3]), const short UNUSED(no_s[3]))
1276 {
1277         VertPickData *data = userData;
1278         if ((data->mvert[index].flag & ME_HIDE) == 0) {
1279                 float sco[2];
1280
1281                 if (ED_view3d_project_float_object(data->ar, co, sco, V3D_PROJ_TEST_CLIP_DEFAULT) == V3D_PROJ_RET_OK) {
1282                         const float len = len_manhattan_v2v2(data->mval_f, sco);
1283                         if (len < data->len_best) {
1284                                 data->len_best = len;
1285                                 data->v_idx_best = index;
1286                         }
1287                 }
1288         }
1289 }
1290 bool ED_mesh_pick_vert(bContext *C, Object *ob, const int mval[2], unsigned int *index, int size, bool use_zbuf)
1291 {
1292         ViewContext vc;
1293         Mesh *me = ob->data;
1294
1295         BLI_assert(me && GS(me->id.name) == ID_ME);
1296
1297         if (!me || me->totvert == 0)
1298                 return false;
1299
1300         ED_view3d_viewcontext_init(C, &vc);
1301
1302         if (use_zbuf) {
1303                 if (size > 0) {
1304                         /* sample rect to increase chances of selecting, so that when clicking
1305                          * on an face in the backbuf, we can still select a vert */
1306
1307                         float dummy_dist;
1308                         *index = ED_view3d_backbuf_sample_rect(&vc, mval, size, 1, me->totvert + 1, &dummy_dist);
1309                 }
1310                 else {
1311                         /* sample only on the exact position */
1312                         *index = ED_view3d_backbuf_sample(&vc, mval[0], mval[1]);
1313                 }
1314
1315                 if ((*index) == 0 || (*index) > (unsigned int)me->totvert)
1316                         return false;
1317
1318                 (*index)--;
1319         }
1320         else {
1321                 /* derived mesh to find deformed locations */
1322                 DerivedMesh *dm = mesh_get_derived_final(vc.depsgraph, vc.scene, ob, CD_MASK_BAREMESH);
1323                 ARegion *ar = vc.ar;
1324                 RegionView3D *rv3d = ar->regiondata;
1325
1326                 /* find the vert closest to 'mval' */
1327                 const float mval_f[2] = {(float)mval[0],
1328                                          (float)mval[1]};
1329
1330                 VertPickData data = {NULL};
1331
1332                 ED_view3d_init_mats_rv3d(ob, rv3d);
1333
1334                 if (dm == NULL) {
1335                         return false;
1336                 }
1337
1338                 /* setup data */
1339                 data.mvert = me->mvert;
1340                 data.ar = ar;
1341                 data.mval_f = mval_f;
1342                 data.len_best = FLT_MAX;
1343                 data.v_idx_best = -1;
1344
1345                 dm->foreachMappedVert(dm, ed_mesh_pick_vert__mapFunc, &data, DM_FOREACH_NOP);
1346
1347                 dm->release(dm);
1348
1349                 if (data.v_idx_best == -1) {
1350                         return false;
1351                 }
1352
1353                 *index = data.v_idx_best;
1354         }
1355
1356         return true;
1357 }
1358
1359
1360 MDeformVert *ED_mesh_active_dvert_get_em(Object *ob, BMVert **r_eve)
1361 {
1362         if (ob->mode & OB_MODE_EDIT && ob->type == OB_MESH && ob->defbase.first) {
1363                 Mesh *me = ob->data;
1364                 BMesh *bm = me->edit_btmesh->bm;
1365                 const int cd_dvert_offset = CustomData_get_offset(&bm->vdata, CD_MDEFORMVERT);
1366
1367                 if (cd_dvert_offset != -1) {
1368                         BMVert *eve = BM_mesh_active_vert_get(bm);
1369
1370                         if (eve) {
1371                                 if (r_eve) *r_eve = eve;
1372                                 return BM_ELEM_CD_GET_VOID_P(eve, cd_dvert_offset);
1373                         }
1374                 }
1375         }
1376
1377         if (r_eve) *r_eve = NULL;
1378         return NULL;
1379 }
1380
1381 MDeformVert *ED_mesh_active_dvert_get_ob(Object *ob, int *r_index)
1382 {
1383         Mesh *me = ob->data;
1384         int index = BKE_mesh_mselect_active_get(me, ME_VSEL);
1385         if (r_index) *r_index = index;
1386         if (index == -1 || me->dvert == NULL) {
1387                 return NULL;
1388         }
1389         else {
1390                 return me->dvert + index;
1391         }
1392 }
1393
1394 MDeformVert *ED_mesh_active_dvert_get_only(Object *ob)
1395 {
1396         if (ob->type == OB_MESH) {
1397                 if (ob->mode & OB_MODE_EDIT) {
1398                         return ED_mesh_active_dvert_get_em(ob, NULL);
1399                 }
1400                 else {
1401                         return ED_mesh_active_dvert_get_ob(ob, NULL);
1402                 }
1403         }
1404         else {
1405                 return NULL;
1406         }
1407 }
1408
1409 void EDBM_mesh_stats_multi(
1410         struct Object **objects, const uint objects_len,
1411         int totelem[3], int totelem_sel[3])
1412 {
1413         if (totelem) {
1414                 totelem[0] = 0;
1415                 totelem[1] = 0;
1416                 totelem[2] = 0;
1417         }
1418         if (totelem_sel) {
1419                 totelem_sel[0] = 0;
1420                 totelem_sel[1] = 0;
1421                 totelem_sel[2] = 0;
1422         }
1423
1424         for (uint ob_index = 0; ob_index < objects_len; ob_index++) {
1425                 Object *obedit = objects[ob_index];
1426                 BMEditMesh *em = BKE_editmesh_from_object(obedit);
1427                 BMesh *bm = em->bm;
1428                 if (totelem) {
1429                         totelem[0] += bm->totvert;
1430                         totelem[1] += bm->totedge;
1431                         totelem[2] += bm->totface;
1432                 }
1433                 if (totelem_sel) {
1434                         totelem_sel[0] += bm->totvertsel;
1435                         totelem_sel[1] += bm->totedgesel;
1436                         totelem_sel[2] += bm->totfacesel;
1437                 }
1438         }
1439 }
1440
1441
1442 void EDBM_mesh_elem_index_ensure_multi(Object **objects, const uint objects_len, const char htype)
1443 {
1444         int elem_offset[4] = {0, 0, 0, 0};
1445         for (uint ob_index = 0; ob_index < objects_len; ob_index++) {
1446                 Object *obedit = objects[ob_index];
1447                 BMEditMesh *em = BKE_editmesh_from_object(obedit);
1448                 BMesh *bm = em->bm;
1449                 BM_mesh_elem_index_ensure_ex(bm, htype, elem_offset);
1450         }
1451 }