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