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