Cleanup: comment line length (editors)
[blender.git] / source / blender / editors / mesh / meshtools.c
1 /*
2  * ***** BEGIN GPL LICENSE BLOCK *****
3  *
4  * This program is free software; you can redistribute it and/or
5  * modify it under the terms of the GNU General Public License
6  * as published by the Free Software Foundation; either version 2
7  * of the License, or (at your option) any later version.
8  *
9  * This program is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12  * GNU General Public License for more details.
13  *
14  * You should have received a copy of the GNU General Public License
15  * along with this program; if not, write to the Free Software Foundation,
16  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
17  *
18  * The Original Code is Copyright (C) 2004 by Blender Foundation
19  * All rights reserved.
20  *
21  * The Original Code is: all of this file.
22  *
23  * Contributor(s): none yet.
24  *
25  * ***** END GPL LICENSE BLOCK *****
26  */
27
28 /** \file blender/editors/mesh/meshtools.c
29  *  \ingroup edmesh
30  *
31  * meshtools.c: no editmode (violated already :), mirror & join),
32  * tools operating on meshes
33  */
34
35 #include "MEM_guardedalloc.h"
36
37 #include "DNA_mesh_types.h"
38 #include "DNA_key_types.h"
39 #include "DNA_material_types.h"
40 #include "DNA_modifier_types.h"
41 #include "DNA_object_types.h"
42 #include "DNA_scene_types.h"
43 #include "DNA_screen_types.h"
44 #include "DNA_view3d_types.h"
45
46 #include "BLI_math.h"
47 #include "BLI_blenlib.h"
48
49 #include "BKE_context.h"
50 #include "BKE_depsgraph.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
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 static void join_mesh_single(
79         Main *bmain, Scene *scene,
80         Object *ob_dst, Object *ob_src, float imat[4][4],
81         MVert **mvert_pp, MEdge **medge_pp, MLoop **mloop_pp, MPoly **mpoly_pp,
82         CustomData *vdata, CustomData *edata, CustomData *ldata, CustomData *pdata,
83         int totvert, int totedge, int totloop, int totpoly,
84         Key *key, Key *nkey,
85         Material **matar, int *matmap, int totcol,
86         int *vertofs, int *edgeofs, int *loopofs, int *polyofs)
87 {
88         int a, b;
89
90         Mesh *me = ob_src->data;
91         MVert *mvert = *mvert_pp;
92         MEdge *medge = *medge_pp;
93         MLoop *mloop = *mloop_pp;
94         MPoly *mpoly = *mpoly_pp;
95
96         if (me->totvert) {
97                 /* merge customdata flag */
98                 ((Mesh *)ob_dst->data)->cd_flag |= me->cd_flag;
99
100                 /* standard data */
101                 CustomData_merge(&me->vdata, vdata, CD_MASK_MESH, CD_DEFAULT, totvert);
102                 CustomData_copy_data_named(&me->vdata, vdata, 0, *vertofs, me->totvert);
103
104                 /* vertex groups */
105                 MDeformVert *dvert = CustomData_get(vdata, *vertofs, CD_MDEFORMVERT);
106                 MDeformVert *dvert_src = CustomData_get(&me->vdata, 0, CD_MDEFORMVERT);
107
108                 /* Remap to correct new vgroup indices, if needed. */
109                 if (dvert_src) {
110                         BLI_assert(dvert != NULL);
111
112                         /* Build src to merged mapping of vgroup indices. */
113                         int *vgroup_index_map;
114                         int vgroup_index_map_len;
115                         vgroup_index_map = BKE_object_defgroup_index_map_create(ob_src, ob_dst, &vgroup_index_map_len);
116                         BKE_object_defgroup_index_map_apply(dvert, me->totvert, vgroup_index_map, vgroup_index_map_len);
117                         if (vgroup_index_map != NULL) {
118                                 MEM_freeN(vgroup_index_map);
119                         }
120                 }
121
122                 /* if this is the object we're merging into, no need to do anything */
123                 if (ob_src != ob_dst) {
124                         float cmat[4][4];
125
126                         /* watch this: switch matmul order really goes wrong */
127                         mul_m4_m4m4(cmat, imat, ob_src->obmat);
128
129                         /* transform vertex coordinates into new space */
130                         for (a = 0, mvert = *mvert_pp; a < me->totvert; a++, mvert++) {
131                                 mul_m4_v3(cmat, mvert->co);
132                         }
133
134                         /* for each shapekey in destination mesh:
135                          * - if there's a matching one, copy it across (will need to transform vertices into new space...)
136                          * - otherwise, just copy own coordinates of mesh (no need to transform vertex coordinates into new space)
137                          */
138                         if (key) {
139                                 /* if this mesh has any shapekeys, check first, otherwise just copy coordinates */
140                                 for (KeyBlock *kb = key->block.first; kb; kb = kb->next) {
141                                         /* get pointer to where to write data for this mesh in shapekey's data array */
142                                         float (*cos)[3] = ((float (*)[3])kb->data) + *vertofs;
143
144                                         /* check if this mesh has such a shapekey */
145                                         KeyBlock *okb = me->key ? BKE_keyblock_find_name(me->key, kb->name) : NULL;
146                                         if (okb) {
147                                                 /* copy this mesh's shapekey to the destination shapekey
148                                                  * (need to transform first) */
149                                                 float (*ocos)[3] = okb->data;
150                                                 for (a = 0; a < me->totvert; a++, cos++, ocos++) {
151                                                         copy_v3_v3(*cos, *ocos);
152                                                         mul_m4_v3(cmat, *cos);
153                                                 }
154                                         }
155                                         else {
156                                                 /* copy this mesh's vertex coordinates to the destination shapekey */
157                                                 for (a = 0, mvert = *mvert_pp; a < me->totvert; a++, cos++, mvert++) {
158                                                         copy_v3_v3(*cos, mvert->co);
159                                                 }
160                                         }
161                                 }
162                         }
163                 }
164                 else {
165                         /* for each shapekey in destination mesh:
166                          * - if it was an 'original', copy the appropriate data from nkey
167                          * - otherwise, copy across plain coordinates (no need to transform coordinates)
168                          */
169                         if (key) {
170                                 for (KeyBlock *kb = key->block.first; kb; kb = kb->next) {
171                                         /* get pointer to where to write data for this mesh in shapekey's data array */
172                                         float (*cos)[3] = ((float (*)[3])kb->data) + *vertofs;
173
174                                         /* check if this was one of the original shapekeys */
175                                         KeyBlock *okb = nkey ? BKE_keyblock_find_name(nkey, kb->name) : NULL;
176                                         if (okb) {
177                                                 /* copy this mesh's shapekey to the destination shapekey */
178                                                 float (*ocos)[3] = okb->data;
179                                                 for (a = 0; a < me->totvert; a++, cos++, ocos++) {
180                                                         copy_v3_v3(*cos, *ocos);
181                                                 }
182                                         }
183                                         else {
184                                                 /* copy base-coordinates to the destination shapekey */
185                                                 for (a = 0, mvert = *mvert_pp; a < me->totvert; a++, cos++, mvert++) {
186                                                         copy_v3_v3(*cos, mvert->co);
187                                                 }
188                                         }
189                                 }
190                         }
191                 }
192         }
193
194         if (me->totedge) {
195                 CustomData_merge(&me->edata, edata, CD_MASK_MESH, CD_DEFAULT, totedge);
196                 CustomData_copy_data_named(&me->edata, edata, 0, *edgeofs, me->totedge);
197
198                 for (a = 0; a < me->totedge; a++, medge++) {
199                         medge->v1 += *vertofs;
200                         medge->v2 += *vertofs;
201                 }
202         }
203
204         if (me->totloop) {
205                 if (ob_src != ob_dst) {
206                         MultiresModifierData *mmd;
207
208                         multiresModifier_prepare_join(scene, ob_src, ob_dst);
209
210                         if ((mmd = get_multires_modifier(scene, ob_src, true))) {
211                                 ED_object_iter_other(bmain, ob_src, true,
212                                                      ED_object_multires_update_totlevels_cb,
213                                                      &mmd->totlvl);
214                         }
215                 }
216
217                 CustomData_merge(&me->ldata, ldata, CD_MASK_MESH, CD_DEFAULT, totloop);
218                 CustomData_copy_data_named(&me->ldata, ldata, 0, *loopofs, me->totloop);
219
220                 for (a = 0; a < me->totloop; a++, mloop++) {
221                         mloop->v += *vertofs;
222                         mloop->e += *edgeofs;
223                 }
224         }
225
226         if (me->totpoly) {
227                 if (matmap) {
228                         /* make mapping for materials */
229                         for (a = 1; a <= ob_src->totcol; a++) {
230                                 Material *ma = give_current_material(ob_src, a);
231
232                                 for (b = 0; b < totcol; b++) {
233                                         if (ma == matar[b]) {
234                                                 matmap[a - 1] = b;
235                                                 break;
236                                         }
237                                 }
238                         }
239                 }
240
241                 CustomData_merge(&me->pdata, pdata, CD_MASK_MESH, CD_DEFAULT, totpoly);
242                 CustomData_copy_data_named(&me->pdata, pdata, 0, *polyofs, me->totpoly);
243
244                 for (a = 0; a < me->totpoly; a++, mpoly++) {
245                         mpoly->loopstart += *loopofs;
246                         mpoly->mat_nr = matmap ? matmap[mpoly->mat_nr] : 0;
247                 }
248         }
249
250         /* these are used for relinking (cannot be set earlier, or else reattaching goes wrong) */
251         *vertofs += me->totvert;
252         *mvert_pp += me->totvert;
253         *edgeofs += me->totedge;
254         *medge_pp += me->totedge;
255         *loopofs += me->totloop;
256         *mloop_pp += me->totloop;
257         *polyofs += me->totpoly;
258         *mpoly_pp += me->totpoly;
259 }
260
261 int join_mesh_exec(bContext *C, wmOperator *op)
262 {
263         Main *bmain = CTX_data_main(C);
264         Scene *scene = CTX_data_scene(C);
265         Object *ob = CTX_data_active_object(C);
266         Material **matar = NULL, *ma;
267         Mesh *me;
268         MVert *mvert = NULL;
269         MEdge *medge = NULL;
270         MPoly *mpoly = NULL;
271         MLoop *mloop = NULL;
272         Key *key, *nkey = NULL;
273         KeyBlock *kb, *kbn;
274         float imat[4][4];
275         int a, b, totcol, totmat = 0, totedge = 0, totvert = 0;
276         int totloop = 0, totpoly = 0, vertofs, *matmap = NULL;
277         int i, haskey = 0, edgeofs, loopofs, polyofs;
278         bool ok = false;
279         bDeformGroup *dg, *odg;
280         CustomData vdata, edata, fdata, ldata, pdata;
281
282         if (scene->obedit) {
283                 BKE_report(op->reports, RPT_WARNING, "Cannot join while in edit mode");
284                 return OPERATOR_CANCELLED;
285         }
286
287         /* ob is the object we are adding geometry to */
288         if (!ob || ob->type != OB_MESH) {
289                 BKE_report(op->reports, RPT_WARNING, "Active object is not a mesh");
290                 return OPERATOR_CANCELLED;
291         }
292
293         /* count & check */
294         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
295         {
296                 if (base->object->type == OB_MESH) {
297                         me = base->object->data;
298
299                         totvert += me->totvert;
300                         totedge += me->totedge;
301                         totloop += me->totloop;
302                         totpoly += me->totpoly;
303                         totmat += base->object->totcol;
304
305                         if (base->object == ob)
306                                 ok = true;
307
308                         /* check for shapekeys */
309                         if (me->key)
310                                 haskey++;
311                 }
312         }
313         CTX_DATA_END;
314
315         /* that way the active object is always selected */
316         if (ok == false) {
317                 BKE_report(op->reports, RPT_WARNING, "Active object is not a selected mesh");
318                 return OPERATOR_CANCELLED;
319         }
320
321         /* only join meshes if there are verts to join, there aren't too many, and we only had one mesh selected */
322         me = (Mesh *)ob->data;
323         key = me->key;
324
325         if (totvert == 0 || totvert == me->totvert) {
326                 BKE_report(op->reports, RPT_WARNING, "No mesh data to join");
327                 return OPERATOR_CANCELLED;
328         }
329
330         if (totvert > MESH_MAX_VERTS) {
331                 BKE_reportf(op->reports, RPT_WARNING, "Joining results in %d vertices, limit is %ld", totvert, MESH_MAX_VERTS);
332                 return OPERATOR_CANCELLED;
333         }
334
335         /* remove tessface to ensure we don't hold references to invalid faces */
336         BKE_mesh_tessface_clear(me);
337
338         /* new material indices and material array */
339         if (totmat) {
340                 matar = MEM_callocN(sizeof(*matar) * totmat, "join_mesh matar");
341                 matmap = MEM_callocN(sizeof(*matmap) * totmat, "join_mesh matmap");
342         }
343         totcol = ob->totcol;
344
345         /* obact materials in new main array, is nicer start! */
346         for (a = 0; a < ob->totcol; a++) {
347                 matar[a] = give_current_material(ob, a + 1);
348                 id_us_plus((ID *)matar[a]);
349                 /* increase id->us : will be lowered later */
350         }
351
352         /* - if destination mesh had shapekeys, move them somewhere safe, and set up placeholders
353          *   with arrays that are large enough to hold shapekey data for all meshes
354          * - if destination mesh didn't have shapekeys, but we encountered some in the meshes we're
355          *   joining, set up a new keyblock and assign to the mesh
356          */
357         if (key) {
358                 /* make a duplicate copy that will only be used here... (must remember to free it!) */
359                 nkey = BKE_key_copy(bmain, key);
360
361                 /* for all keys in old block, clear data-arrays */
362                 for (kb = key->block.first; kb; kb = kb->next) {
363                         if (kb->data) MEM_freeN(kb->data);
364                         kb->data = MEM_callocN(sizeof(float) * 3 * totvert, "join_shapekey");
365                         kb->totelem = totvert;
366                 }
367         }
368         else if (haskey) {
369                 /* add a new key-block and add to the mesh */
370                 key = me->key = BKE_key_add(bmain, (ID *)me);
371                 key->type = KEY_RELATIVE;
372         }
373
374         /* first pass over objects - copying materials and vertexgroups across */
375         CTX_DATA_BEGIN (C, Base *, base, selected_editable_bases)
376         {
377                 /* only act if a mesh, and not the one we're joining to */
378                 if ((ob != base->object) && (base->object->type == OB_MESH)) {
379                         me = base->object->data;
380
381                         /* Join this object's vertex groups to the base one's */
382                         for (dg = base->object->defbase.first; dg; dg = dg->next) {
383                                 /* See if this group exists in the object (if it doesn't, add it to the end) */
384                                 if (!defgroup_find_name(ob, dg->name)) {
385                                         odg = MEM_callocN(sizeof(bDeformGroup), "join deformGroup");
386                                         memcpy(odg, dg, sizeof(bDeformGroup));
387                                         BLI_addtail(&ob->defbase, odg);
388                                 }
389                         }
390                         if (ob->defbase.first && ob->actdef == 0)
391                                 ob->actdef = 1;
392
393
394                         if (me->totvert) {
395                                 /* Add this object's materials to the base one's if they don't exist already
396                                  * (but only if limits not exceeded yet) */
397                                 if (totcol < MAXMAT) {
398                                         for (a = 1; a <= base->object->totcol; a++) {
399                                                 ma = give_current_material(base->object, a);
400
401                                                 for (b = 0; b < totcol; b++) {
402                                                         if (ma == matar[b]) {
403                                                                 break;
404                                                         }
405                                                 }
406                                                 if (b == totcol) {
407                                                         matar[b] = ma;
408                                                         if (ma) {
409                                                                 id_us_plus(&ma->id);
410                                                         }
411                                                         totcol++;
412                                                 }
413                                                 if (totcol >= MAXMAT) {
414                                                         break;
415                                                 }
416                                         }
417                                 }
418
419                                 /* if this mesh has shapekeys,
420                                  * check if destination mesh already has matching entries too */
421                                 if (me->key && key) {
422                                         /* for remapping KeyBlock.relative */
423                                         int      *index_map = MEM_mallocN(sizeof(int)        * me->key->totkey, __func__);
424                                         KeyBlock **kb_map   = MEM_mallocN(sizeof(KeyBlock *) * me->key->totkey, __func__);
425
426                                         for (kb = me->key->block.first, i = 0; kb; kb = kb->next, i++) {
427                                                 BLI_assert(i < me->key->totkey);
428
429                                                 kbn = BKE_keyblock_find_name(key, kb->name);
430                                                 /* if key doesn't exist in destination mesh, add it */
431                                                 if (kbn) {
432                                                         index_map[i] = BLI_findindex(&key->block, kbn);
433                                                 }
434                                                 else {
435                                                         index_map[i] = key->totkey;
436
437                                                         kbn = BKE_keyblock_add(key, kb->name);
438
439                                                         BKE_keyblock_copy_settings(kbn, kb);
440
441                                                         /* adjust settings to fit (allocate a new data-array) */
442                                                         kbn->data = MEM_callocN(sizeof(float) * 3 * totvert, "joined_shapekey");
443                                                         kbn->totelem = totvert;
444                                                 }
445
446                                                 kb_map[i] = kbn;
447                                         }
448
449                                         /* remap relative index values */
450                                         for (kb = me->key->block.first, i = 0; kb; kb = kb->next, i++) {
451                                                 /* sanity check, should always be true */
452                                                 if (LIKELY(kb->relative < me->key->totkey)) {
453                                                         kb_map[i]->relative = index_map[kb->relative];
454                                                 }
455                                         }
456
457                                         MEM_freeN(index_map);
458                                         MEM_freeN(kb_map);
459                                 }
460                         }
461                 }
462         }
463         CTX_DATA_END;
464
465
466         /* setup new data for destination mesh */
467         CustomData_reset(&vdata);
468         CustomData_reset(&edata);
469         CustomData_reset(&fdata);
470         CustomData_reset(&ldata);
471         CustomData_reset(&pdata);
472
473         mvert = CustomData_add_layer(&vdata, CD_MVERT, CD_CALLOC, NULL, totvert);
474         medge = CustomData_add_layer(&edata, CD_MEDGE, CD_CALLOC, NULL, totedge);
475         mloop = CustomData_add_layer(&ldata, CD_MLOOP, CD_CALLOC, NULL, totloop);
476         mpoly = CustomData_add_layer(&pdata, CD_MPOLY, CD_CALLOC, NULL, totpoly);
477
478         vertofs = 0;
479         edgeofs = 0;
480         loopofs = 0;
481         polyofs = 0;
482
483         /* inverse transform for all selected meshes in this object */
484         invert_m4_m4(imat, ob->obmat);
485
486         /* Add back active mesh first.
487          * This allows to keep things similar as they were, as much as possible
488          * (i.e. data from active mesh will remain first ones in new result of the merge,
489          * in same order for CD layers, etc). See also T50084.
490          */
491         join_mesh_single(
492                     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                                     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_base_object_free_and_unlink(bmain, scene, base);
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         DAG_relations_tag_update(bmain);   /* removed objects, need to rebuild dag */
592
593         DAG_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         Main *bmain = CTX_data_main(C);
608         Scene *scene = CTX_data_scene(C);
609         Object *ob = CTX_data_active_object(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(bmain, (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(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(Object *ob, DerivedMesh *dm, char mode)
690 {
691         if (mode == 'u') {        /* use table */
692                 if (ED_mesh_mirrtopo_recalc_check(ob->data, dm, ob->mode, &mesh_topo_store)) {
693                         ED_mesh_mirror_topo_table(ob, dm, 's');
694                 }
695         }
696         else if (mode == 's') { /* start table */
697                 ED_mesh_mirrtopo_init(ob->data, dm, ob->mode, &mesh_topo_store, false);
698         }
699         else if (mode == 'e') { /* end table */
700                 ED_mesh_mirrtopo_free(&mesh_topo_store);
701         }
702         else {
703                 BLI_assert(0);
704         }
705
706         return 0;
707 }
708
709 /** \} */
710
711
712 static int mesh_get_x_mirror_vert_spatial(Object *ob, DerivedMesh *dm, int index)
713 {
714         Mesh *me = ob->data;
715         MVert *mvert = dm ? dm->getVertArray(dm) : me->mvert;
716         float vec[3];
717
718         mvert = &mvert[index];
719         vec[0] = -mvert->co[0];
720         vec[1] = mvert->co[1];
721         vec[2] = mvert->co[2];
722
723         return ED_mesh_mirror_spatial_table(ob, NULL, dm, vec, 'u');
724 }
725
726 static int mesh_get_x_mirror_vert_topo(Object *ob, DerivedMesh *dm, int index)
727 {
728         if (ED_mesh_mirror_topo_table(ob, dm, 'u') == -1)
729                 return -1;
730
731         return mesh_topo_store.index_lookup[index];
732 }
733
734 int mesh_get_x_mirror_vert(Object *ob, DerivedMesh *dm, int index, const bool use_topology)
735 {
736         if (use_topology) {
737                 return mesh_get_x_mirror_vert_topo(ob, dm, index);
738         }
739         else {
740                 return mesh_get_x_mirror_vert_spatial(ob, dm, index);
741         }
742 }
743
744 static BMVert *editbmesh_get_x_mirror_vert_spatial(Object *ob, BMEditMesh *em, const float co[3])
745 {
746         float vec[3];
747         int i;
748
749         /* ignore nan verts */
750         if ((isfinite(co[0]) == false) ||
751             (isfinite(co[1]) == false) ||
752             (isfinite(co[2]) == false))
753         {
754                 return NULL;
755         }
756
757         vec[0] = -co[0];
758         vec[1] = co[1];
759         vec[2] = co[2];
760
761         i = ED_mesh_mirror_spatial_table(ob, em, NULL, vec, 'u');
762         if (i != -1) {
763                 return BM_vert_at_index(em->bm, i);
764         }
765         return NULL;
766 }
767
768 static BMVert *editbmesh_get_x_mirror_vert_topo(Object *ob, struct BMEditMesh *em, BMVert *eve, int index)
769 {
770         intptr_t poinval;
771         if (ED_mesh_mirror_topo_table(ob, NULL, 'u') == -1)
772                 return NULL;
773
774         if (index == -1) {
775                 BMIter iter;
776                 BMVert *v;
777
778                 index = 0;
779                 BM_ITER_MESH (v, &iter, em->bm, BM_VERTS_OF_MESH) {
780                         if (v == eve)
781                                 break;
782                         index++;
783                 }
784
785                 if (index == em->bm->totvert) {
786                         return NULL;
787                 }
788         }
789
790         poinval = mesh_topo_store.index_lookup[index];
791
792         if (poinval != -1)
793                 return (BMVert *)(poinval);
794         return NULL;
795 }
796
797 BMVert *editbmesh_get_x_mirror_vert(Object *ob, struct BMEditMesh *em, BMVert *eve, const float co[3], int index, const bool use_topology)
798 {
799         if (use_topology) {
800                 return editbmesh_get_x_mirror_vert_topo(ob, em, eve, index);
801         }
802         else {
803                 return editbmesh_get_x_mirror_vert_spatial(ob, em, co);
804         }
805 }
806
807 /**
808  * Wrapper for objectmode/editmode.
809  *
810  * call #BM_mesh_elem_table_ensure first for editmesh.
811  */
812 int ED_mesh_mirror_get_vert(Object *ob, int index)
813 {
814         Mesh *me = ob->data;
815         BMEditMesh *em = me->edit_btmesh;
816         bool use_topology = (me->editflag & ME_EDIT_MIRROR_TOPO) != 0;
817         int index_mirr;
818
819         if (em) {
820                 BMVert *eve, *eve_mirr;
821                 eve = BM_vert_at_index(em->bm, index);
822                 eve_mirr = editbmesh_get_x_mirror_vert(ob, em, eve, eve->co, index, use_topology);
823                 index_mirr = eve_mirr ? BM_elem_index_get(eve_mirr) : -1;
824         }
825         else {
826                 index_mirr = mesh_get_x_mirror_vert(ob, NULL, index, use_topology);
827         }
828
829         return index_mirr;
830 }
831
832 #if 0
833
834 static float *editmesh_get_mirror_uv(BMEditMesh *em, int axis, float *uv, float *mirrCent, float *face_cent)
835 {
836         float vec[2];
837         float cent_vec[2];
838         float cent[2];
839
840         /* ignore nan verts */
841         if (isnan(uv[0]) || !isfinite(uv[0]) ||
842             isnan(uv[1]) || !isfinite(uv[1])
843             )
844         {
845                 return NULL;
846         }
847
848         if (axis) {
849                 vec[0] = uv[0];
850                 vec[1] = -((uv[1]) - mirrCent[1]) + mirrCent[1];
851
852                 cent_vec[0] = face_cent[0];
853                 cent_vec[1] = -((face_cent[1]) - mirrCent[1]) + mirrCent[1];
854         }
855         else {
856                 vec[0] = -((uv[0]) - mirrCent[0]) + mirrCent[0];
857                 vec[1] = uv[1];
858
859                 cent_vec[0] = -((face_cent[0]) - mirrCent[0]) + mirrCent[0];
860                 cent_vec[1] = face_cent[1];
861         }
862
863         /* TODO - Optimize */
864         {
865                 BMIter iter;
866                 BMFace *efa;
867
868                 BM_ITER_MESH (efa, &iter, em->bm, BM_FACES_OF_MESH) {
869                         uv_poly_center(efa, cent, cd_loop_uv_offset);
870
871                         if ( (fabsf(cent[0] - cent_vec[0]) < 0.001f) && (fabsf(cent[1] - cent_vec[1]) < 0.001f) ) {
872                                 BMIter liter;
873                                 BMLoop *l;
874
875                                 BM_ITER_ELEM (l, &liter, efa, BM_LOOPS_OF_FACE) {
876                                         MLoopUV *luv = CustomData_bmesh_get(&em->bm->ldata, l->head.data, CD_MLOOPUV);
877                                         if ( (fabsf(luv->uv[0] - vec[0]) < 0.001f) && (fabsf(luv->uv[1] - vec[1]) < 0.001f) ) {
878                                                 return luv->uv;
879
880                                         }
881                                 }
882                         }
883                 }
884         }
885
886         return NULL;
887 }
888
889 #endif
890
891 static unsigned int mirror_facehash(const void *ptr)
892 {
893         const MFace *mf = ptr;
894         unsigned int v0, v1;
895
896         if (mf->v4) {
897                 v0 = MIN4(mf->v1, mf->v2, mf->v3, mf->v4);
898                 v1 = MAX4(mf->v1, mf->v2, mf->v3, mf->v4);
899         }
900         else {
901                 v0 = MIN3(mf->v1, mf->v2, mf->v3);
902                 v1 = MAX3(mf->v1, mf->v2, mf->v3);
903         }
904
905         return ((v0 * 39) ^ (v1 * 31));
906 }
907
908 static int mirror_facerotation(MFace *a, MFace *b)
909 {
910         if (b->v4) {
911                 if (a->v1 == b->v1 && a->v2 == b->v2 && a->v3 == b->v3 && a->v4 == b->v4)
912                         return 0;
913                 else if (a->v4 == b->v1 && a->v1 == b->v2 && a->v2 == b->v3 && a->v3 == b->v4)
914                         return 1;
915                 else if (a->v3 == b->v1 && a->v4 == b->v2 && a->v1 == b->v3 && a->v2 == b->v4)
916                         return 2;
917                 else if (a->v2 == b->v1 && a->v3 == b->v2 && a->v4 == b->v3 && a->v1 == b->v4)
918                         return 3;
919         }
920         else {
921                 if (a->v1 == b->v1 && a->v2 == b->v2 && a->v3 == b->v3)
922                         return 0;
923                 else if (a->v3 == b->v1 && a->v1 == b->v2 && a->v2 == b->v3)
924                         return 1;
925                 else if (a->v2 == b->v1 && a->v3 == b->v2 && a->v1 == b->v3)
926                         return 2;
927         }
928
929         return -1;
930 }
931
932 static bool mirror_facecmp(const void *a, const void *b)
933 {
934         return (mirror_facerotation((MFace *)a, (MFace *)b) == -1);
935 }
936
937 /* BMESH_TODO, convert to MPoly (functions above also) */
938 int *mesh_get_x_mirror_faces(Object *ob, BMEditMesh *em, DerivedMesh *dm)
939 {
940         Mesh *me = ob->data;
941         MVert *mv, *mvert;
942         MFace mirrormf, *mf, *hashmf, *mface;
943         GHash *fhash;
944         int *mirrorverts, *mirrorfaces;
945
946         BLI_assert(em == NULL);  /* Does not work otherwise, currently... */
947
948         const bool use_topology = (me->editflag & ME_EDIT_MIRROR_TOPO) != 0;
949         const int totvert = dm ? dm->getNumVerts(dm) : me->totvert;
950         const int totface = dm ? dm->getNumTessFaces(dm) : me->totface;
951         int a;
952
953         mirrorverts = MEM_callocN(sizeof(int) * totvert, "MirrorVerts");
954         mirrorfaces = MEM_callocN(sizeof(int) * 2 * totface, "MirrorFaces");
955
956         mvert = dm ? dm->getVertArray(dm) : me->mvert;
957         mface = dm ? dm->getTessFaceArray(dm) : me->mface;
958
959         ED_mesh_mirror_spatial_table(ob, em, dm, NULL, 's');
960
961         for (a = 0, mv = mvert; a < totvert; a++, mv++)
962                 mirrorverts[a] = mesh_get_x_mirror_vert(ob, dm, a, use_topology);
963
964         ED_mesh_mirror_spatial_table(ob, em, dm, NULL, 'e');
965
966         fhash = BLI_ghash_new_ex(mirror_facehash, mirror_facecmp, "mirror_facehash gh", me->totface);
967         for (a = 0, mf = mface; a < totface; a++, mf++)
968                 BLI_ghash_insert(fhash, mf, mf);
969
970         for (a = 0, mf = mface; a < totface; a++, mf++) {
971                 mirrormf.v1 = mirrorverts[mf->v3];
972                 mirrormf.v2 = mirrorverts[mf->v2];
973                 mirrormf.v3 = mirrorverts[mf->v1];
974                 mirrormf.v4 = (mf->v4) ? mirrorverts[mf->v4] : 0;
975
976                 /* make sure v4 is not 0 if a quad */
977                 if (mf->v4 && mirrormf.v4 == 0) {
978                         SWAP(unsigned int, mirrormf.v1, mirrormf.v3);
979                         SWAP(unsigned int, mirrormf.v2, mirrormf.v4);
980                 }
981
982                 hashmf = BLI_ghash_lookup(fhash, &mirrormf);
983                 if (hashmf) {
984                         mirrorfaces[a * 2] = hashmf - mface;
985                         mirrorfaces[a * 2 + 1] = mirror_facerotation(&mirrormf, hashmf);
986                 }
987                 else
988                         mirrorfaces[a * 2] = -1;
989         }
990
991         BLI_ghash_free(fhash, NULL, NULL);
992         MEM_freeN(mirrorverts);
993
994         return mirrorfaces;
995 }
996
997 /* selection, vertex and face */
998 /* returns 0 if not found, otherwise 1 */
999
1000 /**
1001  * Face selection in object mode,
1002  * currently only weight-paint and vertex-paint use this.
1003  *
1004  * \return boolean true == Found
1005  */
1006 bool ED_mesh_pick_face(bContext *C, Object *ob, const int mval[2], unsigned int *index, int size)
1007 {
1008         ViewContext vc;
1009         Mesh *me = ob->data;
1010
1011         BLI_assert(me && GS(me->id.name) == ID_ME);
1012
1013         if (!me || me->totpoly == 0)
1014                 return false;
1015
1016         ED_view3d_viewcontext_init(C, &vc);
1017
1018         if (size) {
1019                 /* sample rect to increase chances of selecting, so that when clicking
1020                  * on an edge in the backbuf, we can still select a face */
1021
1022                 float dummy_dist;
1023                 *index = ED_view3d_backbuf_sample_rect(&vc, mval, size, 1, me->totpoly + 1, &dummy_dist);
1024         }
1025         else {
1026                 /* sample only on the exact position */
1027                 *index = ED_view3d_backbuf_sample(&vc, mval[0], mval[1]);
1028         }
1029
1030         if ((*index) == 0 || (*index) > (unsigned int)me->totpoly)
1031                 return false;
1032
1033         (*index)--;
1034
1035         return true;
1036 }
1037 static void ed_mesh_pick_face_vert__mpoly_find(
1038         /* context */
1039         struct ARegion *ar, const float mval[2],
1040         /* mesh data */
1041         DerivedMesh *dm, MPoly *mp, MLoop *mloop,
1042         /* return values */
1043         float *r_len_best, int *r_v_idx_best)
1044 {
1045         const MLoop *ml;
1046         int j = mp->totloop;
1047         for (ml = &mloop[mp->loopstart]; j--; ml++) {
1048                 float co[3], sco[2], len;
1049                 const int v_idx = ml->v;
1050                 dm->getVertCo(dm, v_idx, co);
1051                 if (ED_view3d_project_float_object(ar, co, sco, V3D_PROJ_TEST_NOP) == V3D_PROJ_RET_OK) {
1052                         len = len_manhattan_v2v2(mval, sco);
1053                         if (len < *r_len_best) {
1054                                 *r_len_best = len;
1055                                 *r_v_idx_best = v_idx;
1056                         }
1057                 }
1058         }
1059 }
1060 /**
1061  * Use when the back buffer stores face index values. but we want a vert.
1062  * This gets the face then finds the closest vertex to mval.
1063  */
1064 bool ED_mesh_pick_face_vert(bContext *C, Object *ob, const int mval[2], unsigned int *index, int size)
1065 {
1066         unsigned int poly_index;
1067         Mesh *me = ob->data;
1068
1069         BLI_assert(me && GS(me->id.name) == ID_ME);
1070
1071         if (ED_mesh_pick_face(C, ob, mval, &poly_index, size)) {
1072                 Scene *scene = CTX_data_scene(C);
1073                 struct ARegion *ar = CTX_wm_region(C);
1074
1075                 /* derived mesh to find deformed locations */
1076                 DerivedMesh *dm = mesh_get_derived_final(scene, ob, CD_MASK_BAREMESH | CD_MASK_ORIGINDEX);
1077
1078                 int v_idx_best = ORIGINDEX_NONE;
1079
1080                 /* find the vert closest to 'mval' */
1081                 const float mval_f[2] = {UNPACK2(mval)};
1082                 float len_best = FLT_MAX;
1083
1084                 MPoly *dm_mpoly;
1085                 MLoop *dm_mloop;
1086                 unsigned int dm_mpoly_tot;
1087                 const int *index_mp_to_orig;
1088
1089                 dm_mpoly = dm->getPolyArray(dm);
1090                 dm_mloop = dm->getLoopArray(dm);
1091
1092                 dm_mpoly_tot = dm->getNumPolys(dm);
1093
1094                 index_mp_to_orig = dm->getPolyDataArray(dm, CD_ORIGINDEX);
1095
1096                 /* tag all verts using this face */
1097                 if (index_mp_to_orig) {
1098                         unsigned int i;
1099
1100                         for (i = 0; i < dm_mpoly_tot; i++) {
1101                                 if (index_mp_to_orig[i] == poly_index) {
1102                                         ed_mesh_pick_face_vert__mpoly_find(
1103                                                 ar, mval_f,
1104                                                 dm, &dm_mpoly[i], dm_mloop,
1105                                                 &len_best, &v_idx_best);
1106                                 }
1107                         }
1108                 }
1109                 else {
1110                         if (poly_index < dm_mpoly_tot) {
1111                                 ed_mesh_pick_face_vert__mpoly_find(
1112                                         ar, mval_f,
1113                                         dm, &dm_mpoly[poly_index], dm_mloop,
1114                                         &len_best, &v_idx_best);
1115                         }
1116                 }
1117
1118                 /* map 'dm -> me' index if possible */
1119                 if (v_idx_best != ORIGINDEX_NONE) {
1120                         const int *index_mv_to_orig;
1121
1122                         index_mv_to_orig = dm->getVertDataArray(dm, CD_ORIGINDEX);
1123                         if (index_mv_to_orig) {
1124                                 v_idx_best = index_mv_to_orig[v_idx_best];
1125                         }
1126                 }
1127
1128                 dm->release(dm);
1129
1130                 if ((v_idx_best != ORIGINDEX_NONE) && (v_idx_best < me->totvert)) {
1131                         *index = v_idx_best;
1132                         return true;
1133                 }
1134         }
1135
1136         return false;
1137 }
1138
1139 /**
1140  * Vertex selection in object mode,
1141  * currently only weight paint uses this.
1142  *
1143  * \return boolean true == Found
1144  */
1145 typedef struct VertPickData {
1146         const MVert *mvert;
1147         const float *mval_f;  /* [2] */
1148         ARegion *ar;
1149
1150         /* runtime */
1151         float len_best;
1152         int v_idx_best;
1153 } VertPickData;
1154
1155 static void ed_mesh_pick_vert__mapFunc(void *userData, int index, const float co[3],
1156                                        const float UNUSED(no_f[3]), const short UNUSED(no_s[3]))
1157 {
1158         VertPickData *data = userData;
1159         if ((data->mvert[index].flag & ME_HIDE) == 0) {
1160                 float sco[2];
1161
1162                 if (ED_view3d_project_float_object(data->ar, co, sco, V3D_PROJ_TEST_CLIP_DEFAULT) == V3D_PROJ_RET_OK) {
1163                         const float len = len_manhattan_v2v2(data->mval_f, sco);
1164                         if (len < data->len_best) {
1165                                 data->len_best = len;
1166                                 data->v_idx_best = index;
1167                         }
1168                 }
1169         }
1170 }
1171 bool ED_mesh_pick_vert(bContext *C, Object *ob, const int mval[2], unsigned int *index, int size, bool use_zbuf)
1172 {
1173         ViewContext vc;
1174         Mesh *me = ob->data;
1175
1176         BLI_assert(me && GS(me->id.name) == ID_ME);
1177
1178         if (!me || me->totvert == 0)
1179                 return false;
1180
1181         ED_view3d_viewcontext_init(C, &vc);
1182
1183         if (use_zbuf) {
1184                 if (size > 0) {
1185                         /* sample rect to increase chances of selecting, so that when clicking
1186                          * on an face in the backbuf, we can still select a vert */
1187
1188                         float dummy_dist;
1189                         *index = ED_view3d_backbuf_sample_rect(&vc, mval, size, 1, me->totvert + 1, &dummy_dist);
1190                 }
1191                 else {
1192                         /* sample only on the exact position */
1193                         *index = ED_view3d_backbuf_sample(&vc, mval[0], mval[1]);
1194                 }
1195
1196                 if ((*index) == 0 || (*index) > (unsigned int)me->totvert)
1197                         return false;
1198
1199                 (*index)--;
1200         }
1201         else {
1202                 /* derived mesh to find deformed locations */
1203                 DerivedMesh *dm = mesh_get_derived_final(vc.scene, ob, CD_MASK_BAREMESH);
1204                 ARegion *ar = vc.ar;
1205                 RegionView3D *rv3d = ar->regiondata;
1206
1207                 /* find the vert closest to 'mval' */
1208                 const float mval_f[2] = {(float)mval[0],
1209                                          (float)mval[1]};
1210
1211                 VertPickData data = {NULL};
1212
1213                 ED_view3d_init_mats_rv3d(ob, rv3d);
1214
1215                 if (dm == NULL) {
1216                         return false;
1217                 }
1218
1219                 /* setup data */
1220                 data.mvert = me->mvert;
1221                 data.ar = ar;
1222                 data.mval_f = mval_f;
1223                 data.len_best = FLT_MAX;
1224                 data.v_idx_best = -1;
1225
1226                 dm->foreachMappedVert(dm, ed_mesh_pick_vert__mapFunc, &data, DM_FOREACH_NOP);
1227
1228                 dm->release(dm);
1229
1230                 if (data.v_idx_best == -1) {
1231                         return false;
1232                 }
1233
1234                 *index = data.v_idx_best;
1235         }
1236
1237         return true;
1238 }
1239
1240
1241 MDeformVert *ED_mesh_active_dvert_get_em(Object *ob, BMVert **r_eve)
1242 {
1243         if (ob->mode & OB_MODE_EDIT && ob->type == OB_MESH && ob->defbase.first) {
1244                 Mesh *me = ob->data;
1245                 BMesh *bm = me->edit_btmesh->bm;
1246                 const int cd_dvert_offset = CustomData_get_offset(&bm->vdata, CD_MDEFORMVERT);
1247
1248                 if (cd_dvert_offset != -1) {
1249                         BMVert *eve = BM_mesh_active_vert_get(bm);
1250
1251                         if (eve) {
1252                                 if (r_eve) *r_eve = eve;
1253                                 return BM_ELEM_CD_GET_VOID_P(eve, cd_dvert_offset);
1254                         }
1255                 }
1256         }
1257
1258         if (r_eve) *r_eve = NULL;
1259         return NULL;
1260 }
1261
1262 MDeformVert *ED_mesh_active_dvert_get_ob(Object *ob, int *r_index)
1263 {
1264         Mesh *me = ob->data;
1265         int index = BKE_mesh_mselect_active_get(me, ME_VSEL);
1266         if (r_index) *r_index = index;
1267         if (index == -1 || me->dvert == NULL) {
1268                 return NULL;
1269         }
1270         else {
1271                 return me->dvert + index;
1272         }
1273 }
1274
1275 MDeformVert *ED_mesh_active_dvert_get_only(Object *ob)
1276 {
1277         if (ob->type == OB_MESH) {
1278                 if (ob->mode & OB_MODE_EDIT) {
1279                         return ED_mesh_active_dvert_get_em(ob, NULL);
1280                 }
1281                 else {
1282                         return ED_mesh_active_dvert_get_ob(ob, NULL);
1283                 }
1284         }
1285         else {
1286                 return NULL;
1287         }
1288 }