f7817171220a97ac454173c6ef744ccd690b8a3e
[blender.git] / source / blender / blenkernel / intern / object.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) 2001-2002 by NaN Holding BV.
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/blenkernel/intern/object.c
29  *  \ingroup bke
30  */
31
32
33 #include <string.h>
34 #include <math.h>
35 #include <stdio.h>
36
37 #include "MEM_guardedalloc.h"
38
39 #include "DNA_anim_types.h"
40 #include "DNA_armature_types.h"
41 #include "DNA_camera_types.h"
42 #include "DNA_constraint_types.h"
43 #include "DNA_gpencil_types.h"
44 #include "DNA_group_types.h"
45 #include "DNA_key_types.h"
46 #include "DNA_lamp_types.h"
47 #include "DNA_lattice_types.h"
48 #include "DNA_material_types.h"
49 #include "DNA_meta_types.h"
50 #include "DNA_mesh_types.h"
51 #include "DNA_meshdata_types.h"
52 #include "DNA_movieclip_types.h"
53 #include "DNA_scene_types.h"
54 #include "DNA_screen_types.h"
55 #include "DNA_sequence_types.h"
56 #include "DNA_smoke_types.h"
57 #include "DNA_space_types.h"
58 #include "DNA_view3d_types.h"
59 #include "DNA_world_types.h"
60 #include "DNA_object_types.h"
61 #include "DNA_property_types.h"
62 #include "DNA_rigidbody_types.h"
63
64 #include "BLI_blenlib.h"
65 #include "BLI_math.h"
66 #include "BLI_threads.h"
67 #include "BLI_utildefines.h"
68 #include "BLI_linklist.h"
69 #include "BLI_kdtree.h"
70
71 #include "BLT_translation.h"
72
73 #include "BKE_pbvh.h"
74 #include "BKE_main.h"
75 #include "BKE_global.h"
76 #include "BKE_idprop.h"
77 #include "BKE_armature.h"
78 #include "BKE_action.h"
79 #include "BKE_bullet.h"
80 #include "BKE_deform.h"
81 #include "BKE_depsgraph.h"
82 #include "BKE_DerivedMesh.h"
83 #include "BKE_animsys.h"
84 #include "BKE_anim.h"
85 #include "BKE_constraint.h"
86 #include "BKE_curve.h"
87 #include "BKE_displist.h"
88 #include "BKE_effect.h"
89 #include "BKE_fcurve.h"
90 #include "BKE_group.h"
91 #include "BKE_icons.h"
92 #include "BKE_key.h"
93 #include "BKE_lamp.h"
94 #include "BKE_layer.h"
95 #include "BKE_lattice.h"
96 #include "BKE_library.h"
97 #include "BKE_library_query.h"
98 #include "BKE_library_remap.h"
99 #include "BKE_linestyle.h"
100 #include "BKE_mesh.h"
101 #include "BKE_editmesh.h"
102 #include "BKE_mball.h"
103 #include "BKE_modifier.h"
104 #include "BKE_multires.h"
105 #include "BKE_node.h"
106 #include "BKE_object.h"
107 #include "BKE_paint.h"
108 #include "BKE_particle.h"
109 #include "BKE_pointcache.h"
110 #include "BKE_property.h"
111 #include "BKE_rigidbody.h"
112 #include "BKE_sca.h"
113 #include "BKE_scene.h"
114 #include "BKE_sequencer.h"
115 #include "BKE_speaker.h"
116 #include "BKE_softbody.h"
117 #include "BKE_subsurf.h"
118 #include "BKE_material.h"
119 #include "BKE_camera.h"
120 #include "BKE_image.h"
121
122 #ifdef WITH_MOD_FLUID
123 #include "LBM_fluidsim.h"
124 #endif
125
126 #ifdef WITH_PYTHON
127 #include "BPY_extern.h"
128 #endif
129
130 #include "CCGSubSurf.h"
131
132 #include "GPU_material.h"
133
134 /* Vertex parent modifies original BMesh which is not safe for threading.
135  * Ideally such a modification should be handled as a separate DAG update
136  * callback for mesh datablock, but for until it is actually supported use
137  * simpler solution with a mutex lock.
138  *                                               - sergey -
139  */
140 #define VPARENT_THREADING_HACK
141
142 #ifdef VPARENT_THREADING_HACK
143 static ThreadMutex vparent_lock = BLI_MUTEX_INITIALIZER;
144 #endif
145
146 void BKE_object_workob_clear(Object *workob)
147 {
148         memset(workob, 0, sizeof(Object));
149         
150         workob->size[0] = workob->size[1] = workob->size[2] = 1.0f;
151         workob->dscale[0] = workob->dscale[1] = workob->dscale[2] = 1.0f;
152         workob->rotmode = ROT_MODE_EUL;
153 }
154
155 void BKE_object_update_base_layer(struct Scene *scene, Object *ob)
156 {
157         Base *base = scene->base.first;
158
159         while (base) {
160                 if (base->object == ob) base->lay = ob->lay;
161                 base = base->next;
162         }
163 }
164
165 void BKE_object_free_particlesystems(Object *ob)
166 {
167         ParticleSystem *psys;
168
169         while ((psys = BLI_pophead(&ob->particlesystem))) {
170                 psys_free(ob, psys);
171         }
172 }
173
174 void BKE_object_free_softbody(Object *ob)
175 {
176         if (ob->soft) {
177                 sbFree(ob->soft);
178                 ob->soft = NULL;
179         }
180 }
181
182 void BKE_object_free_bulletsoftbody(Object *ob)
183 {
184         if (ob->bsoft) {
185                 bsbFree(ob->bsoft);
186                 ob->bsoft = NULL;
187         }
188 }
189
190 void BKE_object_free_curve_cache(Object *ob)
191 {
192         if (ob->curve_cache) {
193                 BKE_displist_free(&ob->curve_cache->disp);
194                 BKE_curve_bevelList_free(&ob->curve_cache->bev);
195                 if (ob->curve_cache->path) {
196                         free_path(ob->curve_cache->path);
197                 }
198                 BKE_nurbList_free(&ob->curve_cache->deformed_nurbs);
199                 MEM_freeN(ob->curve_cache);
200                 ob->curve_cache = NULL;
201         }
202 }
203
204 void BKE_object_free_modifiers(Object *ob)
205 {
206         ModifierData *md;
207
208         while ((md = BLI_pophead(&ob->modifiers))) {
209                 modifier_free(md);
210         }
211
212         /* particle modifiers were freed, so free the particlesystems as well */
213         BKE_object_free_particlesystems(ob);
214
215         /* same for softbody */
216         BKE_object_free_softbody(ob);
217
218         /* modifiers may have stored data in the DM cache */
219         BKE_object_free_derived_caches(ob);
220 }
221
222 void BKE_object_modifier_hook_reset(Object *ob, HookModifierData *hmd)
223 {
224         /* reset functionality */
225         if (hmd->object) {
226                 bPoseChannel *pchan = BKE_pose_channel_find_name(hmd->object->pose, hmd->subtarget);
227
228                 if (hmd->subtarget[0] && pchan) {
229                         float imat[4][4], mat[4][4];
230
231                         /* calculate the world-space matrix for the pose-channel target first, then carry on as usual */
232                         mul_m4_m4m4(mat, hmd->object->obmat, pchan->pose_mat);
233
234                         invert_m4_m4(imat, mat);
235                         mul_m4_m4m4(hmd->parentinv, imat, ob->obmat);
236                 }
237                 else {
238                         invert_m4_m4(hmd->object->imat, hmd->object->obmat);
239                         mul_m4_m4m4(hmd->parentinv, hmd->object->imat, ob->obmat);
240                 }
241         }
242 }
243
244 bool BKE_object_support_modifier_type_check(Object *ob, int modifier_type)
245 {
246         const ModifierTypeInfo *mti;
247
248         mti = modifierType_getInfo(modifier_type);
249
250         /* only geometry objects should be able to get modifiers [#25291] */
251         if (!ELEM(ob->type, OB_MESH, OB_CURVE, OB_SURF, OB_FONT, OB_LATTICE)) {
252                 return false;
253         }
254
255         if (ob->type == OB_LATTICE && (mti->flags & eModifierTypeFlag_AcceptsLattice) == 0) {
256                 return false;
257         }
258
259         if (!((mti->flags & eModifierTypeFlag_AcceptsCVs) ||
260               (ob->type == OB_MESH && (mti->flags & eModifierTypeFlag_AcceptsMesh))))
261         {
262                 return false;
263         }
264
265         return true;
266 }
267
268 void BKE_object_link_modifiers(struct Object *ob_dst, const struct Object *ob_src)
269 {
270         ModifierData *md;
271         BKE_object_free_modifiers(ob_dst);
272
273         if (!ELEM(ob_dst->type, OB_MESH, OB_CURVE, OB_SURF, OB_FONT, OB_LATTICE)) {
274                 /* only objects listed above can have modifiers and linking them to objects
275                  * which doesn't have modifiers stack is quite silly */
276                 return;
277         }
278
279         for (md = ob_src->modifiers.first; md; md = md->next) {
280                 ModifierData *nmd = NULL;
281
282                 if (ELEM(md->type,
283                          eModifierType_Hook,
284                          eModifierType_Collision))
285                 {
286                         continue;
287                 }
288
289                 if (!BKE_object_support_modifier_type_check(ob_dst, md->type))
290                         continue;
291
292                 switch (md->type) {
293                         case eModifierType_Softbody:
294                                 BKE_object_copy_softbody(ob_dst, ob_src);
295                                 break;
296                         case eModifierType_Skin:
297                                 /* ensure skin-node customdata exists */
298                                 BKE_mesh_ensure_skin_customdata(ob_dst->data);
299                                 break;
300                 }
301
302                 nmd = modifier_new(md->type);
303                 BLI_strncpy(nmd->name, md->name, sizeof(nmd->name));
304
305                 if (md->type == eModifierType_Multires) {
306                         /* Has to be done after mod creation, but *before* we actually copy its settings! */
307                         multiresModifier_sync_levels_ex(ob_dst, (MultiresModifierData *)md, (MultiresModifierData *)nmd);
308                 }
309
310                 modifier_copyData(md, nmd);
311                 BLI_addtail(&ob_dst->modifiers, nmd);
312                 modifier_unique_name(&ob_dst->modifiers, nmd);
313         }
314
315         BKE_object_copy_particlesystems(ob_dst, ob_src);
316
317         /* TODO: smoke?, cloth? */
318 }
319
320 /* free data derived from mesh, called when mesh changes or is freed */
321 void BKE_object_free_derived_caches(Object *ob)
322 {
323         /* also serves as signal to remake texspace */
324         if (ob->type == OB_MESH) {
325                 Mesh *me = ob->data;
326
327                 if (me && me->bb) {
328                         me->bb->flag |= BOUNDBOX_DIRTY;
329                 }
330         }
331         else if (ELEM(ob->type, OB_SURF, OB_CURVE, OB_FONT)) {
332                 Curve *cu = ob->data;
333
334                 if (cu && cu->bb) {
335                         cu->bb->flag |= BOUNDBOX_DIRTY;
336                 }
337         }
338
339         if (ob->bb) {
340                 MEM_freeN(ob->bb);
341                 ob->bb = NULL;
342         }
343
344         if (ob->derivedFinal) {
345                 ob->derivedFinal->needsFree = 1;
346                 ob->derivedFinal->release(ob->derivedFinal);
347                 ob->derivedFinal = NULL;
348         }
349         if (ob->derivedDeform) {
350                 ob->derivedDeform->needsFree = 1;
351                 ob->derivedDeform->release(ob->derivedDeform);
352                 ob->derivedDeform = NULL;
353         }
354         
355         BKE_object_free_curve_cache(ob);
356 }
357
358 void BKE_object_free_caches(Object *object)
359 {
360         ModifierData *md;
361         short update_flag = 0;
362
363         /* Free particle system caches holding paths. */
364         if (object->particlesystem.first) {
365                 ParticleSystem *psys;
366                 for (psys = object->particlesystem.first;
367                      psys != NULL;
368                      psys = psys->next)
369                 {
370                         psys_free_path_cache(psys, psys->edit);
371                         update_flag |= PSYS_RECALC_REDO;
372                 }
373         }
374
375         /* Free memory used by cached derived meshes in the particle system modifiers. */
376         for (md = object->modifiers.first; md != NULL; md = md->next) {
377                 if (md->type == eModifierType_ParticleSystem) {
378                         ParticleSystemModifierData *psmd = (ParticleSystemModifierData *) md;
379                         if (psmd->dm_final != NULL) {
380                                 psmd->dm_final->needsFree = 1;
381                                 psmd->dm_final->release(psmd->dm_final);
382                                 psmd->dm_final = NULL;
383                                 if (psmd->dm_deformed != NULL) {
384                                         psmd->dm_deformed->needsFree = 1;
385                                         psmd->dm_deformed->release(psmd->dm_deformed);
386                                         psmd->dm_deformed = NULL;
387                                 }
388                                 psmd->flag |= eParticleSystemFlag_file_loaded;
389                                 update_flag |= OB_RECALC_DATA;
390                         }
391                 }
392         }
393
394         /* Tag object for update, so once memory critical operation is over and
395          * scene update routines are back to it's business the object will be
396          * guaranteed to be in a known state.
397          */
398         if (update_flag != 0) {
399                 DAG_id_tag_update(&object->id, update_flag);
400         }
401 }
402
403 /** Free (or release) any data used by this object (does not free the object itself). */
404 void BKE_object_free(Object *ob)
405 {
406         BKE_animdata_free((ID *)ob, false);
407
408         BKE_object_free_modifiers(ob);
409
410         MEM_SAFE_FREE(ob->mat);
411         MEM_SAFE_FREE(ob->matbits);
412         MEM_SAFE_FREE(ob->iuser);
413         MEM_SAFE_FREE(ob->bb);
414
415         BLI_freelistN(&ob->defbase);
416         if (ob->pose) {
417                 BKE_pose_free_ex(ob->pose, false);
418                 ob->pose = NULL;
419         }
420         if (ob->mpath) {
421                 animviz_free_motionpath(ob->mpath);
422                 ob->mpath = NULL;
423         }
424         BKE_bproperty_free_list(&ob->prop);
425
426         free_sensors(&ob->sensors);
427         free_controllers(&ob->controllers);
428         free_actuators(&ob->actuators);
429         
430         BKE_constraints_free_ex(&ob->constraints, false);
431         
432         free_partdeflect(ob->pd);
433         BKE_rigidbody_free_object(ob);
434         BKE_rigidbody_free_constraint(ob);
435
436         if (ob->soft) {
437                 sbFree(ob->soft);
438                 ob->soft = NULL;
439         }
440         if (ob->bsoft) {
441                 bsbFree(ob->bsoft);
442                 ob->bsoft = NULL;
443         }
444         GPU_lamp_free(ob);
445
446         BKE_sculptsession_free(ob);
447
448         BLI_freelistN(&ob->pc_ids);
449
450         BLI_freelistN(&ob->lodlevels);
451
452         /* Free runtime curves data. */
453         if (ob->curve_cache) {
454                 BKE_curve_bevelList_free(&ob->curve_cache->bev);
455                 if (ob->curve_cache->path)
456                         free_path(ob->curve_cache->path);
457                 MEM_freeN(ob->curve_cache);
458                 ob->curve_cache = NULL;
459         }
460
461         BKE_previewimg_free(&ob->preview);
462 }
463
464 /* actual check for internal data, not context or flags */
465 bool BKE_object_is_in_editmode(Object *ob)
466 {
467         if (ob->data == NULL)
468                 return false;
469         
470         if (ob->type == OB_MESH) {
471                 Mesh *me = ob->data;
472                 if (me->edit_btmesh)
473                         return true;
474         }
475         else if (ob->type == OB_ARMATURE) {
476                 bArmature *arm = ob->data;
477                 
478                 if (arm->edbo)
479                         return true;
480         }
481         else if (ob->type == OB_FONT) {
482                 Curve *cu = ob->data;
483                 
484                 if (cu->editfont)
485                         return true;
486         }
487         else if (ob->type == OB_MBALL) {
488                 MetaBall *mb = ob->data;
489                 
490                 if (mb->editelems)
491                         return true;
492         }
493         else if (ob->type == OB_LATTICE) {
494                 Lattice *lt = ob->data;
495                 
496                 if (lt->editlatt)
497                         return true;
498         }
499         else if (ob->type == OB_SURF || ob->type == OB_CURVE) {
500                 Curve *cu = ob->data;
501
502                 if (cu->editnurb)
503                         return true;
504         }
505         return false;
506 }
507
508 bool BKE_object_is_in_editmode_vgroup(Object *ob)
509 {
510         return (OB_TYPE_SUPPORT_VGROUP(ob->type) &&
511                 BKE_object_is_in_editmode(ob));
512 }
513
514 bool BKE_object_is_in_wpaint_select_vert(Object *ob)
515 {
516         if (ob->type == OB_MESH) {
517                 Mesh *me = ob->data;
518                 return ((ob->mode & OB_MODE_WEIGHT_PAINT) &&
519                         (me->edit_btmesh == NULL) &&
520                         (ME_EDIT_PAINT_SEL_MODE(me) == SCE_SELECT_VERTEX));
521         }
522
523         return false;
524 }
525
526 bool BKE_object_exists_check(Object *obtest)
527 {
528         Object *ob;
529         
530         if (obtest == NULL) return false;
531         
532         ob = G.main->object.first;
533         while (ob) {
534                 if (ob == obtest) return true;
535                 ob = ob->id.next;
536         }
537         return false;
538 }
539
540 /* *************************************************** */
541
542 static const char *get_obdata_defname(int type)
543 {
544         switch (type) {
545                 case OB_MESH: return DATA_("Mesh");
546                 case OB_CURVE: return DATA_("Curve");
547                 case OB_SURF: return DATA_("Surf");
548                 case OB_FONT: return DATA_("Text");
549                 case OB_MBALL: return DATA_("Mball");
550                 case OB_CAMERA: return DATA_("Camera");
551                 case OB_LAMP: return DATA_("Lamp");
552                 case OB_LATTICE: return DATA_("Lattice");
553                 case OB_ARMATURE: return DATA_("Armature");
554                 case OB_SPEAKER: return DATA_("Speaker");
555                 case OB_EMPTY: return DATA_("Empty");
556                 default:
557                         printf("get_obdata_defname: Internal error, bad type: %d\n", type);
558                         return DATA_("Empty");
559         }
560 }
561
562 void *BKE_object_obdata_add_from_type(Main *bmain, int type, const char *name)
563 {
564         if (name == NULL) {
565                 name = get_obdata_defname(type);
566         }
567
568         switch (type) {
569                 case OB_MESH:      return BKE_mesh_add(bmain, name);
570                 case OB_CURVE:     return BKE_curve_add(bmain, name, OB_CURVE);
571                 case OB_SURF:      return BKE_curve_add(bmain, name, OB_SURF);
572                 case OB_FONT:      return BKE_curve_add(bmain, name, OB_FONT);
573                 case OB_MBALL:     return BKE_mball_add(bmain, name);
574                 case OB_CAMERA:    return BKE_camera_add(bmain, name);
575                 case OB_LAMP:      return BKE_lamp_add(bmain, name);
576                 case OB_LATTICE:   return BKE_lattice_add(bmain, name);
577                 case OB_ARMATURE:  return BKE_armature_add(bmain, name);
578                 case OB_SPEAKER:   return BKE_speaker_add(bmain, name);
579                 case OB_EMPTY:     return NULL;
580                 default:
581                         printf("%s: Internal error, bad type: %d\n", __func__, type);
582                         return NULL;
583         }
584 }
585
586 void BKE_object_init(Object *ob)
587 {
588         /* BLI_assert(MEMCMP_STRUCT_OFS_IS_ZERO(ob, id)); */  /* ob->type is already initialized... */
589
590         ob->col[0] = ob->col[1] = ob->col[2] = 1.0;
591         ob->col[3] = 1.0;
592         
593         ob->size[0] = ob->size[1] = ob->size[2] = 1.0;
594         ob->dscale[0] = ob->dscale[1] = ob->dscale[2] = 1.0;
595         
596         /* objects should default to having Euler XYZ rotations, 
597          * but rotations default to quaternions 
598          */
599         ob->rotmode = ROT_MODE_EUL;
600
601         unit_axis_angle(ob->rotAxis, &ob->rotAngle);
602         unit_axis_angle(ob->drotAxis, &ob->drotAngle);
603
604         unit_qt(ob->quat);
605         unit_qt(ob->dquat);
606
607         /* rotation locks should be 4D for 4 component rotations by default... */
608         ob->protectflag = OB_LOCK_ROT4D;
609         
610         unit_m4(ob->constinv);
611         unit_m4(ob->parentinv);
612         unit_m4(ob->obmat);
613         ob->dt = OB_TEXTURE;
614         ob->empty_drawtype = OB_PLAINAXES;
615         ob->empty_drawsize = 1.0;
616
617         if (ELEM(ob->type, OB_LAMP, OB_CAMERA, OB_SPEAKER)) {
618                 ob->trackflag = OB_NEGZ;
619                 ob->upflag = OB_POSY;
620         }
621         else {
622                 ob->trackflag = OB_POSY;
623                 ob->upflag = OB_POSZ;
624         }
625         
626         ob->dupon = 1; ob->dupoff = 0;
627         ob->dupsta = 1; ob->dupend = 100;
628         ob->dupfacesca = 1.0;
629
630         /* Game engine defaults*/
631         ob->mass = ob->inertia = 1.0f;
632         ob->formfactor = 0.4f;
633         ob->damping = 0.04f;
634         ob->rdamping = 0.1f;
635         ob->anisotropicFriction[0] = 1.0f;
636         ob->anisotropicFriction[1] = 1.0f;
637         ob->anisotropicFriction[2] = 1.0f;
638         ob->gameflag = OB_PROP | OB_COLLISION;
639         ob->margin = 0.04f;
640         ob->init_state = 1;
641         ob->state = 1;
642         ob->obstacleRad = 1.0f;
643         ob->step_height = 0.15f;
644         ob->jump_speed = 10.0f;
645         ob->fall_speed = 55.0f;
646         ob->max_jumps = 1;
647         ob->col_group = 0x01;
648         ob->col_mask = 0xffff;
649         ob->preview = NULL;
650
651         /* NT fluid sim defaults */
652         ob->fluidsimSettings = NULL;
653
654         BLI_listbase_clear(&ob->pc_ids);
655         
656         /* Animation Visualization defaults */
657         animviz_settings_init(&ob->avs);
658 }
659
660 /* more general add: creates minimum required data, but without vertices etc. */
661 Object *BKE_object_add_only_object(Main *bmain, int type, const char *name)
662 {
663         Object *ob;
664
665         if (!name)
666                 name = get_obdata_defname(type);
667
668         ob = BKE_libblock_alloc(bmain, ID_OB, name);
669
670         /* default object vars */
671         ob->type = type;
672
673         BKE_object_init(ob);
674
675         return ob;
676 }
677
678 /* general add: to scene, with layer from area and default name */
679 /* creates minimum required data, but without vertices etc. */
680 Object *BKE_object_add(
681         Main *bmain, Scene *scene, SceneLayer *sl,
682         int type, const char *name)
683 {
684         Object *ob;
685         ObjectBase *base;
686         LayerCollection *lc;
687
688         ob = BKE_object_add_only_object(bmain, type, name);
689
690         ob->data = BKE_object_obdata_add_from_type(bmain, type, name);
691
692         lc = BKE_layer_collection_active(sl);
693         BKE_collection_object_add(scene, lc->scene_collection, ob);
694
695         base = BKE_scene_layer_base_find(sl, ob);
696         BKE_scene_layer_base_deselect_all(sl);
697         BKE_scene_layer_base_select(sl, base);
698
699         DAG_id_tag_update_ex(bmain, &ob->id, OB_RECALC_OB | OB_RECALC_DATA | OB_RECALC_TIME);
700         return ob;
701 }
702
703
704 #ifdef WITH_GAMEENGINE
705
706 void BKE_object_lod_add(Object *ob)
707 {
708         LodLevel *lod = MEM_callocN(sizeof(LodLevel), "LoD Level");
709         LodLevel *last = ob->lodlevels.last;
710
711         /* If the lod list is empty, initialize it with the base lod level */
712         if (!last) {
713                 LodLevel *base = MEM_callocN(sizeof(LodLevel), "Base LoD Level");
714                 BLI_addtail(&ob->lodlevels, base);
715                 base->flags = OB_LOD_USE_MESH | OB_LOD_USE_MAT;
716                 base->source = ob;
717                 base->obhysteresis = 10;
718                 last = ob->currentlod = base;
719         }
720         
721         lod->distance = last->distance + 25.0f;
722         lod->obhysteresis = 10;
723         lod->flags = OB_LOD_USE_MESH | OB_LOD_USE_MAT;
724
725         BLI_addtail(&ob->lodlevels, lod);
726 }
727
728 static int lod_cmp(const void *a, const void *b)
729 {
730         const LodLevel *loda = a;
731         const LodLevel *lodb = b;
732
733         if (loda->distance < lodb->distance) return -1;
734         return loda->distance > lodb->distance;
735 }
736
737 void BKE_object_lod_sort(Object *ob)
738 {
739         BLI_listbase_sort(&ob->lodlevels, lod_cmp);
740 }
741
742 bool BKE_object_lod_remove(Object *ob, int level)
743 {
744         LodLevel *rem;
745
746         if (level < 1 || level > BLI_listbase_count(&ob->lodlevels) - 1)
747                 return false;
748
749         rem = BLI_findlink(&ob->lodlevels, level);
750
751         if (rem == ob->currentlod) {
752                 ob->currentlod = rem->prev;
753         }
754
755         BLI_remlink(&ob->lodlevels, rem);
756         MEM_freeN(rem);
757
758         /* If there are no user defined lods, remove the base lod as well */
759         if (BLI_listbase_is_single(&ob->lodlevels)) {
760                 LodLevel *base = ob->lodlevels.first;
761                 BLI_remlink(&ob->lodlevels, base);
762                 MEM_freeN(base);
763                 ob->currentlod = NULL;
764         }
765
766         return true;
767 }
768
769 static LodLevel *lod_level_select(Object *ob, const float camera_position[3])
770 {
771         LodLevel *current = ob->currentlod;
772         float dist_sq;
773
774         if (!current) return NULL;
775
776         dist_sq = len_squared_v3v3(ob->obmat[3], camera_position);
777
778         if (dist_sq < SQUARE(current->distance)) {
779                 /* check for higher LoD */
780                 while (current->prev && dist_sq < SQUARE(current->distance)) {
781                         current = current->prev;
782                 }
783         }
784         else {
785                 /* check for lower LoD */
786                 while (current->next && dist_sq > SQUARE(current->next->distance)) {
787                         current = current->next;
788                 }
789         }
790
791         return current;
792 }
793
794 bool BKE_object_lod_is_usable(Object *ob, Scene *scene)
795 {
796         bool active = (scene) ? ob == OBACT : false;
797         return (ob->mode == OB_MODE_OBJECT || !active);
798 }
799
800 void BKE_object_lod_update(Object *ob, const float camera_position[3])
801 {
802         LodLevel *cur_level = ob->currentlod;
803         LodLevel *new_level = lod_level_select(ob, camera_position);
804
805         if (new_level != cur_level) {
806                 ob->currentlod = new_level;
807         }
808 }
809
810 static Object *lod_ob_get(Object *ob, Scene *scene, int flag)
811 {
812         LodLevel *current = ob->currentlod;
813
814         if (!current || !BKE_object_lod_is_usable(ob, scene))
815                 return ob;
816
817         while (current->prev && (!(current->flags & flag) || !current->source || current->source->type != OB_MESH)) {
818                 current = current->prev;
819         }
820
821         return current->source;
822 }
823
824 struct Object *BKE_object_lod_meshob_get(Object *ob, Scene *scene)
825 {
826         return lod_ob_get(ob, scene, OB_LOD_USE_MESH);
827 }
828
829 struct Object *BKE_object_lod_matob_get(Object *ob, Scene *scene)
830 {
831         return lod_ob_get(ob, scene, OB_LOD_USE_MAT);
832 }
833
834 #endif  /* WITH_GAMEENGINE */
835
836
837 SoftBody *copy_softbody(const SoftBody *sb, bool copy_caches)
838 {
839         SoftBody *sbn;
840         
841         if (sb == NULL) return(NULL);
842         
843         sbn = MEM_dupallocN(sb);
844
845         if (copy_caches == false) {
846                 sbn->totspring = sbn->totpoint = 0;
847                 sbn->bpoint = NULL;
848                 sbn->bspring = NULL;
849         }
850         else {
851                 sbn->totspring = sb->totspring;
852                 sbn->totpoint = sb->totpoint;
853
854                 if (sbn->bpoint) {
855                         int i;
856
857                         sbn->bpoint = MEM_dupallocN(sbn->bpoint);
858
859                         for (i = 0; i < sbn->totpoint; i++) {
860                                 if (sbn->bpoint[i].springs)
861                                         sbn->bpoint[i].springs = MEM_dupallocN(sbn->bpoint[i].springs);
862                         }
863                 }
864
865                 if (sb->bspring)
866                         sbn->bspring = MEM_dupallocN(sb->bspring);
867         }
868         
869         sbn->keys = NULL;
870         sbn->totkey = sbn->totpointkey = 0;
871         
872         sbn->scratch = NULL;
873
874         sbn->pointcache = BKE_ptcache_copy_list(&sbn->ptcaches, &sb->ptcaches, copy_caches);
875
876         if (sb->effector_weights)
877                 sbn->effector_weights = MEM_dupallocN(sb->effector_weights);
878
879         return sbn;
880 }
881
882 BulletSoftBody *copy_bulletsoftbody(BulletSoftBody *bsb)
883 {
884         BulletSoftBody *bsbn;
885
886         if (bsb == NULL)
887                 return NULL;
888         bsbn = MEM_dupallocN(bsb);
889         /* no pointer in this structure yet */
890         return bsbn;
891 }
892
893 ParticleSystem *BKE_object_copy_particlesystem(ParticleSystem *psys)
894 {
895         ParticleSystem *psysn;
896         ParticleData *pa;
897         int p;
898
899         psysn = MEM_dupallocN(psys);
900         psysn->particles = MEM_dupallocN(psys->particles);
901         psysn->child = MEM_dupallocN(psys->child);
902
903         if (psys->part->type == PART_HAIR) {
904                 for (p = 0, pa = psysn->particles; p < psysn->totpart; p++, pa++)
905                         pa->hair = MEM_dupallocN(pa->hair);
906         }
907
908         if (psysn->particles && (psysn->particles->keys || psysn->particles->boid)) {
909                 ParticleKey *key = psysn->particles->keys;
910                 BoidParticle *boid = psysn->particles->boid;
911
912                 if (key)
913                         key = MEM_dupallocN(key);
914                 
915                 if (boid)
916                         boid = MEM_dupallocN(boid);
917                 
918                 for (p = 0, pa = psysn->particles; p < psysn->totpart; p++, pa++) {
919                         if (boid)
920                                 pa->boid = boid++;
921                         if (key) {
922                                 pa->keys = key;
923                                 key += pa->totkey;
924                         }
925                 }
926         }
927
928         if (psys->clmd) {
929                 psysn->clmd = (ClothModifierData *)modifier_new(eModifierType_Cloth);
930                 modifier_copyData((ModifierData *)psys->clmd, (ModifierData *)psysn->clmd);
931                 psys->hair_in_dm = psys->hair_out_dm = NULL;
932         }
933
934         BLI_duplicatelist(&psysn->targets, &psys->targets);
935
936         psysn->pathcache = NULL;
937         psysn->childcache = NULL;
938         psysn->edit = NULL;
939         psysn->pdd = NULL;
940         psysn->effectors = NULL;
941         psysn->tree = NULL;
942         psysn->bvhtree = NULL;
943         
944         BLI_listbase_clear(&psysn->pathcachebufs);
945         BLI_listbase_clear(&psysn->childcachebufs);
946         psysn->renderdata = NULL;
947         
948         psysn->pointcache = BKE_ptcache_copy_list(&psysn->ptcaches, &psys->ptcaches, false);
949
950         /* XXX - from reading existing code this seems correct but intended usage of
951          * pointcache should /w cloth should be added in 'ParticleSystem' - campbell */
952         if (psysn->clmd) {
953                 psysn->clmd->point_cache = psysn->pointcache;
954         }
955
956         id_us_plus((ID *)psysn->part);
957
958         return psysn;
959 }
960
961 void BKE_object_copy_particlesystems(Object *ob_dst, const Object *ob_src)
962 {
963         ParticleSystem *psys, *npsys;
964         ModifierData *md;
965
966         if (ob_dst->type != OB_MESH) {
967                 /* currently only mesh objects can have soft body */
968                 return;
969         }
970
971         BLI_listbase_clear(&ob_dst->particlesystem);
972         for (psys = ob_src->particlesystem.first; psys; psys = psys->next) {
973                 npsys = BKE_object_copy_particlesystem(psys);
974
975                 BLI_addtail(&ob_dst->particlesystem, npsys);
976
977                 /* need to update particle modifiers too */
978                 for (md = ob_dst->modifiers.first; md; md = md->next) {
979                         if (md->type == eModifierType_ParticleSystem) {
980                                 ParticleSystemModifierData *psmd = (ParticleSystemModifierData *)md;
981                                 if (psmd->psys == psys)
982                                         psmd->psys = npsys;
983                         }
984                         else if (md->type == eModifierType_DynamicPaint) {
985                                 DynamicPaintModifierData *pmd = (DynamicPaintModifierData *)md;
986                                 if (pmd->brush) {
987                                         if (pmd->brush->psys == psys) {
988                                                 pmd->brush->psys = npsys;
989                                         }
990                                 }
991                         }
992                         else if (md->type == eModifierType_Smoke) {
993                                 SmokeModifierData *smd = (SmokeModifierData *) md;
994                                 
995                                 if (smd->type == MOD_SMOKE_TYPE_FLOW) {
996                                         if (smd->flow) {
997                                                 if (smd->flow->psys == psys)
998                                                         smd->flow->psys = npsys;
999                                         }
1000                                 }
1001                         }
1002                 }
1003         }
1004 }
1005
1006 void BKE_object_copy_softbody(Object *ob_dst, const Object *ob_src)
1007 {
1008         if (ob_src->soft) {
1009                 ob_dst->softflag = ob_src->softflag;
1010                 ob_dst->soft = copy_softbody(ob_src->soft, false);
1011         }
1012 }
1013
1014 static void copy_object_pose(Object *obn, Object *ob)
1015 {
1016         bPoseChannel *chan;
1017         
1018         /* note: need to clear obn->pose pointer first, so that BKE_pose_copy_data works (otherwise there's a crash) */
1019         obn->pose = NULL;
1020         BKE_pose_copy_data(&obn->pose, ob->pose, 1);    /* 1 = copy constraints */
1021
1022         for (chan = obn->pose->chanbase.first; chan; chan = chan->next) {
1023                 bConstraint *con;
1024                 
1025                 chan->flag &= ~(POSE_LOC | POSE_ROT | POSE_SIZE);
1026                 
1027                 for (con = chan->constraints.first; con; con = con->next) {
1028                         const bConstraintTypeInfo *cti = BKE_constraint_typeinfo_get(con);
1029                         ListBase targets = {NULL, NULL};
1030                         bConstraintTarget *ct;
1031                         
1032                         if (cti && cti->get_constraint_targets) {
1033                                 cti->get_constraint_targets(con, &targets);
1034                                 
1035                                 for (ct = targets.first; ct; ct = ct->next) {
1036                                         if (ct->tar == ob)
1037                                                 ct->tar = obn;
1038                                 }
1039                                 
1040                                 if (cti->flush_constraint_targets)
1041                                         cti->flush_constraint_targets(con, &targets, 0);
1042                         }
1043                 }
1044         }
1045 }
1046
1047 static void copy_object_lod(Object *obn, Object *ob)
1048 {
1049         BLI_duplicatelist(&obn->lodlevels, &ob->lodlevels);
1050
1051         if (obn->lodlevels.first)
1052                 ((LodLevel *)obn->lodlevels.first)->source = obn;
1053
1054         obn->currentlod = (LodLevel *)obn->lodlevels.first;
1055 }
1056
1057 bool BKE_object_pose_context_check(Object *ob)
1058 {
1059         if ((ob) &&
1060             (ob->type == OB_ARMATURE) &&
1061             (ob->pose) &&
1062             (ob->mode & OB_MODE_POSE))
1063         {
1064                 return true;
1065         }
1066         else {
1067                 return false;
1068         }
1069 }
1070
1071 Object *BKE_object_pose_armature_get(Object *ob)
1072 {
1073         if (ob == NULL)
1074                 return NULL;
1075
1076         if (BKE_object_pose_context_check(ob))
1077                 return ob;
1078
1079         ob = modifiers_isDeformedByArmature(ob);
1080
1081         if (BKE_object_pose_context_check(ob))
1082                 return ob;
1083
1084         return NULL;
1085 }
1086
1087 void BKE_object_transform_copy(Object *ob_tar, const Object *ob_src)
1088 {
1089         copy_v3_v3(ob_tar->loc, ob_src->loc);
1090         copy_v3_v3(ob_tar->rot, ob_src->rot);
1091         copy_v3_v3(ob_tar->quat, ob_src->quat);
1092         copy_v3_v3(ob_tar->rotAxis, ob_src->rotAxis);
1093         ob_tar->rotAngle = ob_src->rotAngle;
1094         ob_tar->rotmode = ob_src->rotmode;
1095         copy_v3_v3(ob_tar->size, ob_src->size);
1096 }
1097
1098 Object *BKE_object_copy_ex(Main *bmain, Object *ob, bool copy_caches)
1099 {
1100         Object *obn;
1101         ModifierData *md;
1102         int a;
1103
1104         obn = BKE_libblock_copy(bmain, &ob->id);
1105         
1106         if (ob->totcol) {
1107                 obn->mat = MEM_dupallocN(ob->mat);
1108                 obn->matbits = MEM_dupallocN(ob->matbits);
1109                 obn->totcol = ob->totcol;
1110         }
1111
1112         if (ob->iuser) obn->iuser = MEM_dupallocN(ob->iuser);
1113         
1114         if (ob->bb) obn->bb = MEM_dupallocN(ob->bb);
1115         obn->flag &= ~OB_FROMGROUP;
1116         
1117         BLI_listbase_clear(&obn->modifiers);
1118         
1119         for (md = ob->modifiers.first; md; md = md->next) {
1120                 ModifierData *nmd = modifier_new(md->type);
1121                 BLI_strncpy(nmd->name, md->name, sizeof(nmd->name));
1122                 modifier_copyData(md, nmd);
1123                 BLI_addtail(&obn->modifiers, nmd);
1124         }
1125
1126         BLI_listbase_clear(&obn->prop);
1127         BKE_bproperty_copy_list(&obn->prop, &ob->prop);
1128
1129         BKE_sca_logic_copy(obn, ob);
1130
1131         if (ob->pose) {
1132                 copy_object_pose(obn, ob);
1133                 /* backwards compat... non-armatures can get poses in older files? */
1134                 if (ob->type == OB_ARMATURE)
1135                         BKE_pose_rebuild(obn, obn->data);
1136         }
1137         defgroup_copy_list(&obn->defbase, &ob->defbase);
1138         BKE_constraints_copy(&obn->constraints, &ob->constraints, true);
1139
1140         obn->mode = OB_MODE_OBJECT;
1141         obn->sculpt = NULL;
1142
1143         /* increase user numbers */
1144         id_us_plus((ID *)obn->data);
1145         id_us_plus((ID *)obn->gpd);
1146         id_us_plus((ID *)obn->dup_group);
1147
1148         for (a = 0; a < obn->totcol; a++) id_us_plus((ID *)obn->mat[a]);
1149         
1150         if (ob->pd) {
1151                 obn->pd = MEM_dupallocN(ob->pd);
1152                 if (obn->pd->tex)
1153                         id_us_plus(&(obn->pd->tex->id));
1154                 if (obn->pd->rng)
1155                         obn->pd->rng = MEM_dupallocN(ob->pd->rng);
1156         }
1157         obn->soft = copy_softbody(ob->soft, copy_caches);
1158         obn->bsoft = copy_bulletsoftbody(ob->bsoft);
1159         obn->rigidbody_object = BKE_rigidbody_copy_object(ob);
1160         obn->rigidbody_constraint = BKE_rigidbody_copy_constraint(ob);
1161
1162         BKE_object_copy_particlesystems(obn, ob);
1163         
1164         obn->derivedDeform = NULL;
1165         obn->derivedFinal = NULL;
1166
1167         BLI_listbase_clear(&obn->gpulamp);
1168         BLI_listbase_clear(&obn->pc_ids);
1169
1170         obn->mpath = NULL;
1171
1172         copy_object_lod(obn, ob);
1173         
1174         /* Copy runtime surve data. */
1175         obn->curve_cache = NULL;
1176
1177         BKE_id_copy_ensure_local(bmain, &ob->id, &obn->id);
1178
1179         /* Do not copy object's preview (mostly due to the fact renderers create temp copy of objects). */
1180         obn->preview = NULL;
1181
1182         return obn;
1183 }
1184
1185 /* copy objects, will re-initialize cached simulation data */
1186 Object *BKE_object_copy(Main *bmain, Object *ob)
1187 {
1188         return BKE_object_copy_ex(bmain, ob, false);
1189 }
1190
1191 void BKE_object_make_local_ex(Main *bmain, Object *ob, const bool lib_local, const bool clear_proxy)
1192 {
1193         bool is_local = false, is_lib = false;
1194
1195         /* - only lib users: do nothing (unless force_local is set)
1196          * - only local users: set flag
1197          * - mixed: make copy
1198          * In case we make a whole lib's content local, we always want to localize, and we skip remapping (done later).
1199          */
1200
1201         if (!ID_IS_LINKED_DATABLOCK(ob)) {
1202                 return;
1203         }
1204
1205         BKE_library_ID_test_usages(bmain, ob, &is_local, &is_lib);
1206
1207         if (lib_local || is_local) {
1208                 if (!is_lib) {
1209                         id_clear_lib_data(bmain, &ob->id);
1210                         BKE_id_expand_local(bmain, &ob->id);
1211                         if (clear_proxy) {
1212                                 if (ob->proxy_from != NULL) {
1213                                         ob->proxy_from->proxy = NULL;
1214                                         ob->proxy_from->proxy_group = NULL;
1215                                 }
1216                                 ob->proxy = ob->proxy_from = ob->proxy_group = NULL;
1217                         }
1218                 }
1219                 else {
1220                         Object *ob_new = BKE_object_copy(bmain, ob);
1221
1222                         ob_new->id.us = 0;
1223                         ob_new->proxy = ob_new->proxy_from = ob_new->proxy_group = NULL;
1224
1225                         /* setting newid is mandatory for complex make_lib_local logic... */
1226                         ID_NEW_SET(ob, ob_new);
1227
1228                         if (!lib_local) {
1229                                 BKE_libblock_remap(bmain, ob, ob_new, ID_REMAP_SKIP_INDIRECT_USAGE);
1230                         }
1231                 }
1232         }
1233 }
1234
1235 void BKE_object_make_local(Main *bmain, Object *ob, const bool lib_local)
1236 {
1237         BKE_object_make_local_ex(bmain, ob, lib_local, true);
1238 }
1239
1240 /* Returns true if the Object is from an external blend file (libdata) */
1241 bool BKE_object_is_libdata(Object *ob)
1242 {
1243         return (ob && ID_IS_LINKED_DATABLOCK(ob));
1244 }
1245
1246 /* Returns true if the Object data is from an external blend file (libdata) */
1247 bool BKE_object_obdata_is_libdata(Object *ob)
1248 {
1249         /* Linked objects with local obdata are forbidden! */
1250         BLI_assert(!ob || !ob->data || (ID_IS_LINKED_DATABLOCK(ob) ? ID_IS_LINKED_DATABLOCK(ob->data) : true));
1251         return (ob && ob->data && ID_IS_LINKED_DATABLOCK(ob->data));
1252 }
1253
1254 /* *************** PROXY **************** */
1255
1256 /* when you make proxy, ensure the exposed layers are extern */
1257 static void armature_set_id_extern(Object *ob)
1258 {
1259         bArmature *arm = ob->data;
1260         bPoseChannel *pchan;
1261         unsigned int lay = arm->layer_protected;
1262         
1263         for (pchan = ob->pose->chanbase.first; pchan; pchan = pchan->next) {
1264                 if (!(pchan->bone->layer & lay))
1265                         id_lib_extern((ID *)pchan->custom);
1266         }
1267                         
1268 }
1269
1270 void BKE_object_copy_proxy_drivers(Object *ob, Object *target)
1271 {
1272         if ((target->adt) && (target->adt->drivers.first)) {
1273                 FCurve *fcu;
1274                 
1275                 /* add new animdata block */
1276                 if (!ob->adt)
1277                         ob->adt = BKE_animdata_add_id(&ob->id);
1278                 
1279                 /* make a copy of all the drivers (for now), then correct any links that need fixing */
1280                 free_fcurves(&ob->adt->drivers);
1281                 copy_fcurves(&ob->adt->drivers, &target->adt->drivers);
1282                 
1283                 for (fcu = ob->adt->drivers.first; fcu; fcu = fcu->next) {
1284                         ChannelDriver *driver = fcu->driver;
1285                         DriverVar *dvar;
1286                         
1287                         for (dvar = driver->variables.first; dvar; dvar = dvar->next) {
1288                                 /* all drivers */
1289                                 DRIVER_TARGETS_LOOPER(dvar) 
1290                                 {
1291                                         if (dtar->id) {
1292                                                 if ((Object *)dtar->id == target)
1293                                                         dtar->id = (ID *)ob;
1294                                                 else {
1295                                                         /* only on local objects because this causes indirect links
1296                                                          * 'a -> b -> c', blend to point directly to a.blend
1297                                                          * when a.blend has a proxy thats linked into c.blend  */
1298                                                         if (!ID_IS_LINKED_DATABLOCK(ob))
1299                                                                 id_lib_extern((ID *)dtar->id);
1300                                                 }
1301                                         }
1302                                 }
1303                                 DRIVER_TARGETS_LOOPER_END
1304                         }
1305                 }
1306         }
1307 }
1308
1309 /* proxy rule: lib_object->proxy_from == the one we borrow from, set temporally while object_update */
1310 /*             local_object->proxy == pointer to library object, saved in files and read */
1311 /*             local_object->proxy_group == pointer to group dupli-object, saved in files and read */
1312
1313 void BKE_object_make_proxy(Object *ob, Object *target, Object *gob)
1314 {
1315         /* paranoia checks */
1316         if (ID_IS_LINKED_DATABLOCK(ob) || !ID_IS_LINKED_DATABLOCK(target)) {
1317                 printf("cannot make proxy\n");
1318                 return;
1319         }
1320         
1321         ob->proxy = target;
1322         ob->proxy_group = gob;
1323         id_lib_extern(&target->id);
1324         
1325         DAG_id_tag_update(&ob->id, OB_RECALC_OB | OB_RECALC_DATA | OB_RECALC_TIME);
1326         DAG_id_tag_update(&target->id, OB_RECALC_OB | OB_RECALC_DATA | OB_RECALC_TIME);
1327         
1328         /* copy transform
1329          * - gob means this proxy comes from a group, just apply the matrix
1330          *   so the object wont move from its dupli-transform.
1331          *
1332          * - no gob means this is being made from a linked object,
1333          *   this is closer to making a copy of the object - in-place. */
1334         if (gob) {
1335                 ob->rotmode = target->rotmode;
1336                 mul_m4_m4m4(ob->obmat, gob->obmat, target->obmat);
1337                 if (gob->dup_group) { /* should always be true */
1338                         float tvec[3];
1339                         mul_v3_mat3_m4v3(tvec, ob->obmat, gob->dup_group->dupli_ofs);
1340                         sub_v3_v3(ob->obmat[3], tvec);
1341                 }
1342                 BKE_object_apply_mat4(ob, ob->obmat, false, true);
1343         }
1344         else {
1345                 BKE_object_transform_copy(ob, target);
1346                 ob->parent = target->parent; /* libdata */
1347                 copy_m4_m4(ob->parentinv, target->parentinv);
1348         }
1349         
1350         /* copy animdata stuff - drivers only for now... */
1351         BKE_object_copy_proxy_drivers(ob, target);
1352
1353         /* skip constraints? */
1354         /* FIXME: this is considered by many as a bug */
1355         
1356         /* set object type and link to data */
1357         ob->type = target->type;
1358         ob->data = target->data;
1359         id_us_plus((ID *)ob->data);     /* ensures lib data becomes LIB_TAG_EXTERN */
1360
1361         /* copy vertex groups */
1362         defgroup_copy_list(&ob->defbase, &target->defbase);
1363
1364         /* copy material and index information */
1365         ob->actcol = ob->totcol = 0;
1366         if (ob->mat) MEM_freeN(ob->mat);
1367         if (ob->matbits) MEM_freeN(ob->matbits);
1368         ob->mat = NULL;
1369         ob->matbits = NULL;
1370         if ((target->totcol) && (target->mat) && OB_TYPE_SUPPORT_MATERIAL(ob->type)) {
1371                 int i;
1372                 
1373                 ob->actcol = target->actcol;
1374                 ob->totcol = target->totcol;
1375                 
1376                 ob->mat = MEM_dupallocN(target->mat);
1377                 ob->matbits = MEM_dupallocN(target->matbits);
1378                 for (i = 0; i < target->totcol; i++) {
1379                         /* don't need to run test_object_materials since we know this object is new and not used elsewhere */
1380                         id_us_plus((ID *)ob->mat[i]); 
1381                 }
1382         }
1383         
1384         /* type conversions */
1385         if (target->type == OB_ARMATURE) {
1386                 copy_object_pose(ob, target);   /* data copy, object pointers in constraints */
1387                 BKE_pose_rest(ob->pose);            /* clear all transforms in channels */
1388                 BKE_pose_rebuild(ob, ob->data); /* set all internal links */
1389                 
1390                 armature_set_id_extern(ob);
1391         }
1392         else if (target->type == OB_EMPTY) {
1393                 ob->empty_drawtype = target->empty_drawtype;
1394                 ob->empty_drawsize = target->empty_drawsize;
1395         }
1396
1397         /* copy IDProperties */
1398         if (ob->id.properties) {
1399                 IDP_FreeProperty(ob->id.properties);
1400                 MEM_freeN(ob->id.properties);
1401                 ob->id.properties = NULL;
1402         }
1403         if (target->id.properties) {
1404                 ob->id.properties = IDP_CopyProperty(target->id.properties);
1405         }
1406
1407         /* copy drawtype info */
1408         ob->dt = target->dt;
1409 }
1410
1411 /**
1412  * Use with newly created objects to set their size
1413  * (used to apply scene-scale).
1414  */
1415 void BKE_object_obdata_size_init(struct Object *ob, const float size)
1416 {
1417         /* apply radius as a scale to types that support it */
1418         switch (ob->type) {
1419                 case OB_EMPTY:
1420                 {
1421                         ob->empty_drawsize *= size;
1422                         break;
1423                 }
1424                 case OB_FONT:
1425                 {
1426                         Curve *cu = ob->data;
1427                         cu->fsize *= size;
1428                         break;
1429                 }
1430                 case OB_CAMERA:
1431                 {
1432                         Camera *cam = ob->data;
1433                         cam->drawsize *= size;
1434                         break;
1435                 }
1436                 case OB_LAMP:
1437                 {
1438                         Lamp *lamp = ob->data;
1439                         lamp->dist *= size;
1440                         lamp->area_size  *= size;
1441                         lamp->area_sizey *= size;
1442                         lamp->area_sizez *= size;
1443                         break;
1444                 }
1445                 /* Only lattice (not mesh, curve, mball...),
1446                  * because its got data when newly added */
1447                 case OB_LATTICE:
1448                 {
1449                         struct Lattice *lt = ob->data;
1450                         float mat[4][4];
1451
1452                         unit_m4(mat);
1453                         scale_m4_fl(mat, size);
1454
1455                         BKE_lattice_transform(lt, (float (*)[4])mat, false);
1456                         break;
1457                 }
1458         }
1459 }
1460
1461 /* *************** CALC ****************** */
1462
1463 void BKE_object_scale_to_mat3(Object *ob, float mat[3][3])
1464 {
1465         float vec[3];
1466         mul_v3_v3v3(vec, ob->size, ob->dscale);
1467         size_to_mat3(mat, vec);
1468 }
1469
1470 void BKE_object_rot_to_mat3(Object *ob, float mat[3][3], bool use_drot)
1471 {
1472         float rmat[3][3], dmat[3][3];
1473         
1474         /* 'dmat' is the delta-rotation matrix, which will get (pre)multiplied
1475          * with the rotation matrix to yield the appropriate rotation
1476          */
1477
1478         /* rotations may either be quats, eulers (with various rotation orders), or axis-angle */
1479         if (ob->rotmode > 0) {
1480                 /* euler rotations (will cause gimble lock, but this can be alleviated a bit with rotation orders) */
1481                 eulO_to_mat3(rmat, ob->rot, ob->rotmode);
1482                 eulO_to_mat3(dmat, ob->drot, ob->rotmode);
1483         }
1484         else if (ob->rotmode == ROT_MODE_AXISANGLE) {
1485                 /* axis-angle - not really that great for 3D-changing orientations */
1486                 axis_angle_to_mat3(rmat, ob->rotAxis, ob->rotAngle);
1487                 axis_angle_to_mat3(dmat, ob->drotAxis, ob->drotAngle);
1488         }
1489         else {
1490                 /* quats are normalized before use to eliminate scaling issues */
1491                 float tquat[4];
1492                 
1493                 normalize_qt_qt(tquat, ob->quat);
1494                 quat_to_mat3(rmat, tquat);
1495                 
1496                 normalize_qt_qt(tquat, ob->dquat);
1497                 quat_to_mat3(dmat, tquat);
1498         }
1499         
1500         /* combine these rotations */
1501         if (use_drot)
1502                 mul_m3_m3m3(mat, dmat, rmat);
1503         else
1504                 copy_m3_m3(mat, rmat);
1505 }
1506
1507 void BKE_object_mat3_to_rot(Object *ob, float mat[3][3], bool use_compat)
1508 {
1509         BLI_ASSERT_UNIT_M3(mat);
1510
1511         switch (ob->rotmode) {
1512                 case ROT_MODE_QUAT:
1513                 {
1514                         float dquat[4];
1515                         mat3_normalized_to_quat(ob->quat, mat);
1516                         normalize_qt_qt(dquat, ob->dquat);
1517                         invert_qt_normalized(dquat);
1518                         mul_qt_qtqt(ob->quat, dquat, ob->quat);
1519                         break;
1520                 }
1521                 case ROT_MODE_AXISANGLE:
1522                 {
1523                         float quat[4];
1524                         float dquat[4];
1525
1526                         /* without drot we could apply 'mat' directly */
1527                         mat3_normalized_to_quat(quat, mat);
1528                         axis_angle_to_quat(dquat, ob->drotAxis, ob->drotAngle);
1529                         invert_qt_normalized(dquat);
1530                         mul_qt_qtqt(quat, dquat, quat);
1531                         quat_to_axis_angle(ob->rotAxis, &ob->rotAngle, quat);
1532                         break;
1533                 }
1534                 default: /* euler */
1535                 {
1536                         float quat[4];
1537                         float dquat[4];
1538
1539                         /* without drot we could apply 'mat' directly */
1540                         mat3_normalized_to_quat(quat, mat);
1541                         eulO_to_quat(dquat, ob->drot, ob->rotmode);
1542                         invert_qt_normalized(dquat);
1543                         mul_qt_qtqt(quat, dquat, quat);
1544                         /* end drot correction */
1545
1546                         if (use_compat) quat_to_compatible_eulO(ob->rot, ob->rot, ob->rotmode, quat);
1547                         else            quat_to_eulO(ob->rot, ob->rotmode, quat);
1548                         break;
1549                 }
1550         }
1551 }
1552
1553 void BKE_object_tfm_protected_backup(const Object *ob,
1554                                      ObjectTfmProtectedChannels *obtfm)
1555 {
1556
1557 #define TFMCPY(_v) (obtfm->_v = ob->_v)
1558 #define TFMCPY3D(_v) copy_v3_v3(obtfm->_v, ob->_v)
1559 #define TFMCPY4D(_v) copy_v4_v4(obtfm->_v, ob->_v)
1560
1561         TFMCPY3D(loc);
1562         TFMCPY3D(dloc);
1563         TFMCPY3D(size);
1564         TFMCPY3D(dscale);
1565         TFMCPY3D(rot);
1566         TFMCPY3D(drot);
1567         TFMCPY4D(quat);
1568         TFMCPY4D(dquat);
1569         TFMCPY3D(rotAxis);
1570         TFMCPY3D(drotAxis);
1571         TFMCPY(rotAngle);
1572         TFMCPY(drotAngle);
1573
1574 #undef TFMCPY
1575 #undef TFMCPY3D
1576 #undef TFMCPY4D
1577
1578 }
1579
1580 void BKE_object_tfm_protected_restore(Object *ob,
1581                                       const ObjectTfmProtectedChannels *obtfm,
1582                                       const short protectflag)
1583 {
1584         unsigned int i;
1585
1586         for (i = 0; i < 3; i++) {
1587                 if (protectflag & (OB_LOCK_LOCX << i)) {
1588                         ob->loc[i] =  obtfm->loc[i];
1589                         ob->dloc[i] = obtfm->dloc[i];
1590                 }
1591
1592                 if (protectflag & (OB_LOCK_SCALEX << i)) {
1593                         ob->size[i] =  obtfm->size[i];
1594                         ob->dscale[i] = obtfm->dscale[i];
1595                 }
1596
1597                 if (protectflag & (OB_LOCK_ROTX << i)) {
1598                         ob->rot[i] =  obtfm->rot[i];
1599                         ob->drot[i] = obtfm->drot[i];
1600
1601                         ob->quat[i + 1] =  obtfm->quat[i + 1];
1602                         ob->dquat[i + 1] = obtfm->dquat[i + 1];
1603
1604                         ob->rotAxis[i] =  obtfm->rotAxis[i];
1605                         ob->drotAxis[i] = obtfm->drotAxis[i];
1606                 }
1607         }
1608
1609         if ((protectflag & OB_LOCK_ROT4D) && (protectflag & OB_LOCK_ROTW)) {
1610                 ob->quat[0] =  obtfm->quat[0];
1611                 ob->dquat[0] = obtfm->dquat[0];
1612
1613                 ob->rotAngle =  obtfm->rotAngle;
1614                 ob->drotAngle = obtfm->drotAngle;
1615         }
1616 }
1617
1618 void BKE_object_to_mat3(Object *ob, float mat[3][3]) /* no parent */
1619 {
1620         float smat[3][3];
1621         float rmat[3][3];
1622         /*float q1[4];*/
1623         
1624         /* size */
1625         BKE_object_scale_to_mat3(ob, smat);
1626
1627         /* rot */
1628         BKE_object_rot_to_mat3(ob, rmat, true);
1629         mul_m3_m3m3(mat, rmat, smat);
1630 }
1631
1632 void BKE_object_to_mat4(Object *ob, float mat[4][4])
1633 {
1634         float tmat[3][3];
1635         
1636         BKE_object_to_mat3(ob, tmat);
1637         
1638         copy_m4_m3(mat, tmat);
1639
1640         add_v3_v3v3(mat[3], ob->loc, ob->dloc);
1641 }
1642
1643 void BKE_object_matrix_local_get(struct Object *ob, float mat[4][4])
1644 {
1645         if (ob->parent) {
1646                 float par_imat[4][4];
1647
1648                 BKE_object_get_parent_matrix(NULL, ob, ob->parent, par_imat);
1649                 invert_m4(par_imat);
1650                 mul_m4_m4m4(mat, par_imat, ob->obmat);
1651         }
1652         else {
1653                 copy_m4_m4(mat, ob->obmat);
1654         }
1655 }
1656
1657 /* extern */
1658 int enable_cu_speed = 1;
1659
1660 /**
1661  * \param scene: Used when curve cache needs to be calculated, or for dupli-frame time.
1662  * \return success if \a mat is set.
1663  */
1664 static bool ob_parcurve(Scene *scene, Object *ob, Object *par, float mat[4][4])
1665 {
1666         Curve *cu = par->data;
1667         float vec[4], dir[3], quat[4], radius, ctime;
1668
1669         /* only happens on reload file, but violates depsgraph still... fix! */
1670         if (par->curve_cache == NULL) {
1671                 if (scene == NULL) {
1672                         return false;
1673                 }
1674                 BKE_displist_make_curveTypes(scene, par, 0);
1675         }
1676
1677         if (par->curve_cache->path == NULL) {
1678                 return false;
1679         }
1680
1681         /* catch exceptions: curve paths used as a duplicator */
1682         if (enable_cu_speed) {
1683                 /* ctime is now a proper var setting of Curve which gets set by Animato like any other var that's animated,
1684                  * but this will only work if it actually is animated... 
1685                  *
1686                  * we divide the curvetime calculated in the previous step by the length of the path, to get a time
1687                  * factor, which then gets clamped to lie within 0.0 - 1.0 range
1688                  */
1689                 if (cu->pathlen) {
1690                         ctime = cu->ctime / cu->pathlen;
1691                 }
1692                 else {
1693                         ctime = cu->ctime;
1694                 }
1695
1696                 CLAMP(ctime, 0.0f, 1.0f);
1697         }
1698         else {
1699                 /* For dupli-frames only */
1700                 if (scene == NULL) {
1701                         return false;
1702                 }
1703
1704                 ctime = BKE_scene_frame_get(scene);
1705                 if (cu->pathlen) {
1706                         ctime /= cu->pathlen;
1707                 }
1708                 
1709                 CLAMP(ctime, 0.0f, 1.0f);
1710         }
1711         
1712         unit_m4(mat);
1713
1714         /* vec: 4 items! */
1715         if (where_on_path(par, ctime, vec, dir, (cu->flag & CU_FOLLOW) ? quat : NULL, &radius, NULL)) {
1716
1717                 if (cu->flag & CU_FOLLOW) {
1718 #if 0
1719                         float si, q[4];
1720                         vec_to_quat(quat, dir, ob->trackflag, ob->upflag);
1721                         
1722                         /* the tilt */
1723                         normalize_v3(dir);
1724                         q[0] = cosf(0.5 * vec[3]);
1725                         si = sinf(0.5 * vec[3]);
1726                         q[1] = -si * dir[0];
1727                         q[2] = -si * dir[1];
1728                         q[3] = -si * dir[2];
1729                         mul_qt_qtqt(quat, q, quat);
1730 #else
1731                         quat_apply_track(quat, ob->trackflag, ob->upflag);
1732 #endif
1733                         normalize_qt(quat);
1734                         quat_to_mat4(mat, quat);
1735                 }
1736                 
1737                 if (cu->flag & CU_PATH_RADIUS) {
1738                         float tmat[4][4], rmat[4][4];
1739                         scale_m4_fl(tmat, radius);
1740                         mul_m4_m4m4(rmat, tmat, mat);
1741                         copy_m4_m4(mat, rmat);
1742                 }
1743
1744                 copy_v3_v3(mat[3], vec);
1745                 
1746         }
1747
1748         return true;
1749 }
1750
1751 static void ob_parbone(Object *ob, Object *par, float mat[4][4])
1752 {       
1753         bPoseChannel *pchan;
1754         float vec[3];
1755         
1756         if (par->type != OB_ARMATURE) {
1757                 unit_m4(mat);
1758                 return;
1759         }
1760         
1761         /* Make sure the bone is still valid */
1762         pchan = BKE_pose_channel_find_name(par->pose, ob->parsubstr);
1763         if (!pchan || !pchan->bone) {
1764                 printf("Object %s with Bone parent: bone %s doesn't exist\n", ob->id.name + 2, ob->parsubstr);
1765                 unit_m4(mat);
1766                 return;
1767         }
1768
1769         /* get bone transform */
1770         if (pchan->bone->flag & BONE_RELATIVE_PARENTING) {
1771                 /* the new option uses the root - expected bahaviour, but differs from old... */
1772                 /* XXX check on version patching? */
1773                 copy_m4_m4(mat, pchan->chan_mat);
1774         }
1775         else {
1776                 copy_m4_m4(mat, pchan->pose_mat);
1777
1778                 /* but for backwards compatibility, the child has to move to the tail */
1779                 copy_v3_v3(vec, mat[1]);
1780                 mul_v3_fl(vec, pchan->bone->length);
1781                 add_v3_v3(mat[3], vec);
1782         }
1783 }
1784
1785 static void give_parvert(Object *par, int nr, float vec[3])
1786 {
1787         zero_v3(vec);
1788         
1789         if (par->type == OB_MESH) {
1790                 Mesh *me = par->data;
1791                 BMEditMesh *em = me->edit_btmesh;
1792                 DerivedMesh *dm;
1793
1794                 dm = (em) ? em->derivedFinal : par->derivedFinal;
1795                         
1796                 if (dm) {
1797                         int count = 0;
1798                         int numVerts = dm->getNumVerts(dm);
1799
1800                         if (nr < numVerts) {
1801                                 bool use_special_ss_case = false;
1802
1803                                 if (dm->type == DM_TYPE_CCGDM) {
1804                                         ModifierData *md;
1805                                         VirtualModifierData virtualModifierData;
1806                                         use_special_ss_case = true;
1807                                         for (md = modifiers_getVirtualModifierList(par, &virtualModifierData);
1808                                              md != NULL;
1809                                              md = md->next)
1810                                         {
1811                                                 const ModifierTypeInfo *mti = modifierType_getInfo(md->type);
1812                                                 /* TODO(sergey): Check for disabled modifiers. */
1813                                                 if (mti->type != eModifierTypeType_OnlyDeform && md->next != NULL) {
1814                                                         use_special_ss_case = false;
1815                                                         break;
1816                                                 }
1817                                         }
1818                                 }
1819
1820                                 if (!use_special_ss_case) {
1821                                         /* avoid dm->getVertDataArray() since it allocates arrays in the dm (not thread safe) */
1822                                         if (em && dm->type == DM_TYPE_EDITBMESH) {
1823                                                 if (em->bm->elem_table_dirty & BM_VERT) {
1824 #ifdef VPARENT_THREADING_HACK
1825                                                         BLI_mutex_lock(&vparent_lock);
1826                                                         if (em->bm->elem_table_dirty & BM_VERT) {
1827                                                                 BM_mesh_elem_table_ensure(em->bm, BM_VERT);
1828                                                         }
1829                                                         BLI_mutex_unlock(&vparent_lock);
1830 #else
1831                                                         BLI_assert(!"Not safe for threading");
1832                                                         BM_mesh_elem_table_ensure(em->bm, BM_VERT);
1833 #endif
1834                                                 }
1835                                         }
1836                                 }
1837
1838                                 if (use_special_ss_case) {
1839                                         /* Special case if the last modifier is SS and no constructive modifier are in front of it. */
1840                                         CCGDerivedMesh *ccgdm = (CCGDerivedMesh *)dm;
1841                                         CCGVert *ccg_vert = ccgSubSurf_getVert(ccgdm->ss, SET_INT_IN_POINTER(nr));
1842                                         /* In case we deleted some verts, nr may refer to inexistent one now, see T42557. */
1843                                         if (ccg_vert) {
1844                                                 float *co = ccgSubSurf_getVertData(ccgdm->ss, ccg_vert);
1845                                                 add_v3_v3(vec, co);
1846                                                 count++;
1847                                         }
1848                                 }
1849                                 else if (CustomData_has_layer(&dm->vertData, CD_ORIGINDEX) &&
1850                                          !(em && dm->type == DM_TYPE_EDITBMESH))
1851                                 {
1852                                         int i;
1853
1854                                         /* Get the average of all verts with (original index == nr). */
1855                                         for (i = 0; i < numVerts; i++) {
1856                                                 const int *index = dm->getVertData(dm, i, CD_ORIGINDEX);
1857                                                 if (*index == nr) {
1858                                                         float co[3];
1859                                                         dm->getVertCo(dm, i, co);
1860                                                         add_v3_v3(vec, co);
1861                                                         count++;
1862                                                 }
1863                                         }
1864                                 }
1865                                 else {
1866                                         if (nr < numVerts) {
1867                                                 float co[3];
1868                                                 dm->getVertCo(dm, nr, co);
1869                                                 add_v3_v3(vec, co);
1870                                                 count++;
1871                                         }
1872                                 }
1873                         }
1874
1875                         if (count == 0) {
1876                                 /* keep as 0, 0, 0 */
1877                         }
1878                         else if (count > 0) {
1879                                 mul_v3_fl(vec, 1.0f / count);
1880                         }
1881                         else {
1882                                 /* use first index if its out of range */
1883                                 dm->getVertCo(dm, 0, vec);
1884                         }
1885                 }
1886                 else {
1887                         fprintf(stderr,
1888                                 "%s: DerivedMesh is needed to solve parenting, "
1889                                 "object position can be wrong now\n", __func__);
1890                 }
1891         }
1892         else if (ELEM(par->type, OB_CURVE, OB_SURF)) {
1893                 ListBase *nurb;
1894
1895                 /* Unless there's some weird depsgraph failure the cache should exist. */
1896                 BLI_assert(par->curve_cache != NULL);
1897
1898                 if (par->curve_cache->deformed_nurbs.first != NULL) {
1899                         nurb = &par->curve_cache->deformed_nurbs;
1900                 }
1901                 else {
1902                         Curve *cu = par->data;
1903                         nurb = BKE_curve_nurbs_get(cu);
1904                 }
1905
1906                 BKE_nurbList_index_get_co(nurb, nr, vec);
1907         }
1908         else if (par->type == OB_LATTICE) {
1909                 Lattice *latt  = par->data;
1910                 DispList *dl   = par->curve_cache ? BKE_displist_find(&par->curve_cache->disp, DL_VERTS) : NULL;
1911                 float (*co)[3] = dl ? (float (*)[3])dl->verts : NULL;
1912                 int tot;
1913
1914                 if (latt->editlatt) latt = latt->editlatt->latt;
1915
1916                 tot = latt->pntsu * latt->pntsv * latt->pntsw;
1917
1918                 /* ensure dl is correct size */
1919                 BLI_assert(dl == NULL || dl->nr == tot);
1920
1921                 if (nr < tot) {
1922                         if (co) {
1923                                 copy_v3_v3(vec, co[nr]);
1924                         }
1925                         else {
1926                                 copy_v3_v3(vec, latt->def[nr].vec);
1927                         }
1928                 }
1929         }
1930 }
1931
1932 static void ob_parvert3(Object *ob, Object *par, float mat[4][4])
1933 {
1934
1935         /* in local ob space */
1936         if (OB_TYPE_SUPPORT_PARVERT(par->type)) {
1937                 float cmat[3][3], v1[3], v2[3], v3[3], q[4];
1938
1939                 give_parvert(par, ob->par1, v1);
1940                 give_parvert(par, ob->par2, v2);
1941                 give_parvert(par, ob->par3, v3);
1942
1943                 tri_to_quat(q, v1, v2, v3);
1944                 quat_to_mat3(cmat, q);
1945                 copy_m4_m3(mat, cmat);
1946
1947                 mid_v3_v3v3v3(mat[3], v1, v2, v3);
1948         }
1949         else {
1950                 unit_m4(mat);
1951         }
1952 }
1953
1954
1955 void BKE_object_get_parent_matrix(Scene *scene, Object *ob, Object *par, float parentmat[4][4])
1956 {
1957         float tmat[4][4];
1958         float vec[3];
1959         bool ok;
1960
1961         switch (ob->partype & PARTYPE) {
1962                 case PAROBJECT:
1963                         ok = 0;
1964                         if (par->type == OB_CURVE) {
1965                                 if ((((Curve *)par->data)->flag & CU_PATH) &&
1966                                     (ob_parcurve(scene, ob, par, tmat)))
1967                                 {
1968                                         ok = 1;
1969                                 }
1970                         }
1971                         
1972                         if (ok) mul_m4_m4m4(parentmat, par->obmat, tmat);
1973                         else copy_m4_m4(parentmat, par->obmat);
1974                         
1975                         break;
1976                 case PARBONE:
1977                         ob_parbone(ob, par, tmat);
1978                         mul_m4_m4m4(parentmat, par->obmat, tmat);
1979                         break;
1980                 
1981                 case PARVERT1:
1982                         unit_m4(parentmat);
1983                         give_parvert(par, ob->par1, vec);
1984                         mul_v3_m4v3(parentmat[3], par->obmat, vec);
1985                         break;
1986                 case PARVERT3:
1987                         ob_parvert3(ob, par, tmat);
1988                         
1989                         mul_m4_m4m4(parentmat, par->obmat, tmat);
1990                         break;
1991                 
1992                 case PARSKEL:
1993                         copy_m4_m4(parentmat, par->obmat);
1994                         break;
1995         }
1996
1997 }
1998
1999 /**
2000  * \param r_originmat  Optional matrix that stores the space the object is in (without its own matrix applied)
2001  */
2002 static void solve_parenting(Scene *scene, Object *ob, Object *par, float obmat[4][4], float slowmat[4][4],
2003                             float r_originmat[3][3], const bool set_origin)
2004 {
2005         float totmat[4][4];
2006         float tmat[4][4];
2007         float locmat[4][4];
2008         
2009         BKE_object_to_mat4(ob, locmat);
2010         
2011         if (ob->partype & PARSLOW) copy_m4_m4(slowmat, obmat);
2012
2013         BKE_object_get_parent_matrix(scene, ob, par, totmat);
2014         
2015         /* total */
2016         mul_m4_m4m4(tmat, totmat, ob->parentinv);
2017         mul_m4_m4m4(obmat, tmat, locmat);
2018         
2019         if (r_originmat) {
2020                 /* usable originmat */
2021                 copy_m3_m4(r_originmat, tmat);
2022         }
2023         
2024         /* origin, for help line */
2025         if (set_origin) {
2026                 if ((ob->partype & PARTYPE) == PARSKEL) {
2027                         copy_v3_v3(ob->orig, par->obmat[3]);
2028                 }
2029                 else {
2030                         copy_v3_v3(ob->orig, totmat[3]);
2031                 }
2032         }
2033 }
2034
2035 static bool where_is_object_parslow(Object *ob, float obmat[4][4], float slowmat[4][4])
2036 {
2037         float *fp1, *fp2;
2038         float fac1, fac2;
2039         int a;
2040
2041         /* include framerate */
2042         fac1 = (1.0f / (1.0f + fabsf(ob->sf)));
2043         if (fac1 >= 1.0f) return false;
2044         fac2 = 1.0f - fac1;
2045
2046         fp1 = obmat[0];
2047         fp2 = slowmat[0];
2048         for (a = 0; a < 16; a++, fp1++, fp2++) {
2049                 fp1[0] = fac1 * fp1[0] + fac2 * fp2[0];
2050         }
2051
2052         return true;
2053 }
2054
2055 /* note, scene is the active scene while actual_scene is the scene the object resides in */
2056 void BKE_object_where_is_calc_time_ex(Scene *scene, Object *ob, float ctime,
2057                                       RigidBodyWorld *rbw, float r_originmat[3][3])
2058 {
2059         if (ob == NULL) return;
2060         
2061         /* execute drivers only, as animation has already been done */
2062         BKE_animsys_evaluate_animdata(scene, &ob->id, ob->adt, ctime, ADT_RECALC_DRIVERS);
2063         
2064         if (ob->parent) {
2065                 Object *par = ob->parent;
2066                 float slowmat[4][4];
2067                 
2068                 /* calculate parent matrix */
2069                 solve_parenting(scene, ob, par, ob->obmat, slowmat, r_originmat, true);
2070                 
2071                 /* "slow parent" is definitely not threadsafe, and may also give bad results jumping around 
2072                  * An old-fashioned hack which probably doesn't really cut it anymore
2073                  */
2074                 if (ob->partype & PARSLOW) {
2075                         if (!where_is_object_parslow(ob, ob->obmat, slowmat))
2076                                 return;
2077                 }
2078         }
2079         else {
2080                 BKE_object_to_mat4(ob, ob->obmat);
2081         }
2082
2083         /* try to fall back to the scene rigid body world if none given */
2084         rbw = rbw ? rbw : scene->rigidbody_world;
2085         /* read values pushed into RBO from sim/cache... */
2086         BKE_rigidbody_sync_transforms(rbw, ob, ctime);
2087         
2088         /* solve constraints */
2089         if (ob->constraints.first && !(ob->transflag & OB_NO_CONSTRAINTS)) {
2090                 bConstraintOb *cob;
2091                 cob = BKE_constraints_make_evalob(scene, ob, NULL, CONSTRAINT_OBTYPE_OBJECT);
2092                 BKE_constraints_solve(&ob->constraints, cob, ctime);
2093                 BKE_constraints_clear_evalob(cob);
2094         }
2095         
2096         /* set negative scale flag in object */
2097         if (is_negative_m4(ob->obmat)) ob->transflag |= OB_NEG_SCALE;
2098         else ob->transflag &= ~OB_NEG_SCALE;
2099 }
2100
2101 void BKE_object_where_is_calc_time(Scene *scene, Object *ob, float ctime)
2102 {
2103         BKE_object_where_is_calc_time_ex(scene, ob, ctime, NULL, NULL);
2104 }
2105
2106 /* get object transformation matrix without recalculating dependencies and
2107  * constraints -- assume dependencies are already solved by depsgraph.
2108  * no changes to object and it's parent would be done.
2109  * used for bundles orientation in 3d space relative to parented blender camera */
2110 void BKE_object_where_is_calc_mat4(Scene *scene, Object *ob, float obmat[4][4])
2111 {
2112
2113         if (ob->parent) {
2114                 float slowmat[4][4];
2115
2116                 Object *par = ob->parent;
2117                 
2118                 solve_parenting(scene, ob, par, obmat, slowmat, NULL, false);
2119                 
2120                 if (ob->partype & PARSLOW)
2121                         where_is_object_parslow(ob, obmat, slowmat);
2122         }
2123         else {
2124                 BKE_object_to_mat4(ob, obmat);
2125         }
2126 }
2127
2128 void BKE_object_where_is_calc_ex(Scene *scene, RigidBodyWorld *rbw, Object *ob, float r_originmat[3][3])
2129 {
2130         BKE_object_where_is_calc_time_ex(scene, ob, BKE_scene_frame_get(scene), rbw, r_originmat);
2131 }
2132 void BKE_object_where_is_calc(Scene *scene, Object *ob)
2133 {
2134         BKE_object_where_is_calc_time_ex(scene, ob, BKE_scene_frame_get(scene), NULL, NULL);
2135 }
2136
2137 /* for calculation of the inverse parent transform, only used for editor */
2138 void BKE_object_workob_calc_parent(Scene *scene, Object *ob, Object *workob)
2139 {
2140         BKE_object_workob_clear(workob);
2141         
2142         unit_m4(workob->obmat);
2143         unit_m4(workob->parentinv);
2144         unit_m4(workob->constinv);
2145         workob->parent = ob->parent;
2146
2147         workob->trackflag = ob->trackflag;
2148         workob->upflag = ob->upflag;
2149         
2150         workob->partype = ob->partype;
2151         workob->par1 = ob->par1;
2152         workob->par2 = ob->par2;
2153         workob->par3 = ob->par3;
2154
2155         workob->constraints.first = ob->constraints.first;
2156         workob->constraints.last = ob->constraints.last;
2157
2158         BLI_strncpy(workob->parsubstr, ob->parsubstr, sizeof(workob->parsubstr));
2159
2160         BKE_object_where_is_calc(scene, workob);
2161 }
2162
2163 /* see BKE_pchan_apply_mat4() for the equivalent 'pchan' function */
2164 void BKE_object_apply_mat4(Object *ob, float mat[4][4], const bool use_compat, const bool use_parent)
2165 {
2166         float rot[3][3];
2167
2168         if (use_parent && ob->parent) {
2169                 float rmat[4][4], diff_mat[4][4], imat[4][4], parent_mat[4][4];
2170
2171                 BKE_object_get_parent_matrix(NULL, ob, ob->parent, parent_mat);
2172
2173                 mul_m4_m4m4(diff_mat, parent_mat, ob->parentinv);
2174                 invert_m4_m4(imat, diff_mat);
2175                 mul_m4_m4m4(rmat, imat, mat); /* get the parent relative matrix */
2176
2177                 /* same as below, use rmat rather than mat */
2178                 mat4_to_loc_rot_size(ob->loc, rot, ob->size, rmat);
2179         }
2180         else {
2181                 mat4_to_loc_rot_size(ob->loc, rot, ob->size, mat);
2182         }
2183
2184         BKE_object_mat3_to_rot(ob, rot, use_compat);
2185
2186         sub_v3_v3(ob->loc, ob->dloc);
2187
2188         if (ob->dscale[0] != 0.0f) ob->size[0] /= ob->dscale[0];
2189         if (ob->dscale[1] != 0.0f) ob->size[1] /= ob->dscale[1];
2190         if (ob->dscale[2] != 0.0f) ob->size[2] /= ob->dscale[2];
2191
2192         /* BKE_object_mat3_to_rot handles delta rotations */
2193 }
2194
2195 BoundBox *BKE_boundbox_alloc_unit(void)
2196 {
2197         BoundBox *bb;
2198         const float min[3] = {-1.0f, -1.0f, -1.0f}, max[3] = {-1.0f, -1.0f, -1.0f};
2199
2200         bb = MEM_callocN(sizeof(BoundBox), "OB-BoundBox");
2201         BKE_boundbox_init_from_minmax(bb, min, max);
2202         
2203         return bb;
2204 }
2205
2206 void BKE_boundbox_init_from_minmax(BoundBox *bb, const float min[3], const float max[3])
2207 {
2208         bb->vec[0][0] = bb->vec[1][0] = bb->vec[2][0] = bb->vec[3][0] = min[0];
2209         bb->vec[4][0] = bb->vec[5][0] = bb->vec[6][0] = bb->vec[7][0] = max[0];
2210         
2211         bb->vec[0][1] = bb->vec[1][1] = bb->vec[4][1] = bb->vec[5][1] = min[1];
2212         bb->vec[2][1] = bb->vec[3][1] = bb->vec[6][1] = bb->vec[7][1] = max[1];
2213
2214         bb->vec[0][2] = bb->vec[3][2] = bb->vec[4][2] = bb->vec[7][2] = min[2];
2215         bb->vec[1][2] = bb->vec[2][2] = bb->vec[5][2] = bb->vec[6][2] = max[2];
2216 }
2217
2218 void BKE_boundbox_calc_center_aabb(const BoundBox *bb, float r_cent[3])
2219 {
2220         r_cent[0] = 0.5f * (bb->vec[0][0] + bb->vec[4][0]);
2221         r_cent[1] = 0.5f * (bb->vec[0][1] + bb->vec[2][1]);
2222         r_cent[2] = 0.5f * (bb->vec[0][2] + bb->vec[1][2]);
2223 }
2224
2225 void BKE_boundbox_calc_size_aabb(const BoundBox *bb, float r_size[3])
2226 {
2227         r_size[0] = 0.5f * fabsf(bb->vec[0][0] - bb->vec[4][0]);
2228         r_size[1] = 0.5f * fabsf(bb->vec[0][1] - bb->vec[2][1]);
2229         r_size[2] = 0.5f * fabsf(bb->vec[0][2] - bb->vec[1][2]);
2230 }
2231
2232 void BKE_boundbox_minmax(const BoundBox *bb, float obmat[4][4], float r_min[3], float r_max[3])
2233 {
2234         int i;
2235         for (i = 0; i < 8; i++) {
2236                 float vec[3];
2237                 mul_v3_m4v3(vec, obmat, bb->vec[i]);
2238                 minmax_v3v3_v3(r_min, r_max, vec);
2239         }
2240 }
2241
2242 /**
2243  * Returns a BBox which each dimensions are at least epsilon.
2244  * \note In case a given dimension needs to be enlarged, its final value will be in [epsilon, 3 * epsilon] range.
2245  *
2246  * \param bb the input bbox to check.
2247  * \param bb_temp the temp bbox to modify (\a bb content is never changed).
2248  * \param epsilon the minimum dimension to ensure.
2249  * \return either bb (if nothing needed to be changed) or bb_temp.
2250  */
2251 BoundBox *BKE_boundbox_ensure_minimum_dimensions(BoundBox *bb, BoundBox *bb_temp, const float epsilon)
2252 {
2253         if (fabsf(bb->vec[0][0] - bb->vec[4][0]) < epsilon) {
2254                 /* Flat along X axis... */
2255                 *bb_temp = *bb;
2256                 bb = bb_temp;
2257                 bb->vec[0][0] -= epsilon;
2258                 bb->vec[1][0] -= epsilon;
2259                 bb->vec[2][0] -= epsilon;
2260                 bb->vec[3][0] -= epsilon;
2261                 bb->vec[4][0] += epsilon;
2262                 bb->vec[5][0] += epsilon;
2263                 bb->vec[6][0] += epsilon;
2264                 bb->vec[7][0] += epsilon;
2265         }
2266
2267         if (fabsf(bb->vec[0][1] - bb->vec[3][1]) < epsilon) {
2268                 /* Flat along Y axis... */
2269                 if (bb != bb_temp) {
2270                         *bb_temp = *bb;
2271                         bb = bb_temp;
2272                 }
2273                 bb->vec[0][1] -= epsilon;
2274                 bb->vec[1][1] -= epsilon;
2275                 bb->vec[4][1] -= epsilon;
2276                 bb->vec[5][1] -= epsilon;
2277                 bb->vec[2][1] += epsilon;
2278                 bb->vec[3][1] += epsilon;
2279                 bb->vec[6][1] += epsilon;
2280                 bb->vec[7][1] += epsilon;
2281         }
2282
2283         if (fabsf(bb->vec[0][2] - bb->vec[1][2]) < epsilon) {
2284                 /* Flat along Z axis... */
2285                 if (bb != bb_temp) {
2286                         *bb_temp = *bb;
2287                         bb = bb_temp;
2288                 }
2289                 bb->vec[0][2] -= epsilon;
2290                 bb->vec[3][2] -= epsilon;
2291                 bb->vec[4][2] -= epsilon;
2292                 bb->vec[7][2] -= epsilon;
2293                 bb->vec[1][2] += epsilon;
2294                 bb->vec[2][2] += epsilon;
2295                 bb->vec[5][2] += epsilon;
2296                 bb->vec[6][2] += epsilon;
2297         }
2298
2299         return bb;
2300 }
2301
2302 BoundBox *BKE_object_boundbox_get(Object *ob)
2303 {
2304         BoundBox *bb = NULL;
2305         
2306         if (ob->type == OB_MESH) {
2307                 bb = BKE_mesh_boundbox_get(ob);
2308         }
2309         else if (ELEM(ob->type, OB_CURVE, OB_SURF, OB_FONT)) {
2310                 bb = BKE_curve_boundbox_get(ob);
2311         }
2312         else if (ob->type == OB_MBALL) {
2313                 bb = ob->bb;
2314         }
2315         else if (ob->type == OB_LATTICE) {
2316                 bb = BKE_lattice_boundbox_get(ob);
2317         }
2318         else if (ob->type == OB_ARMATURE) {
2319                 bb = BKE_armature_boundbox_get(ob);
2320         }
2321         return bb;
2322 }
2323
2324 /* used to temporally disable/enable boundbox */
2325 void BKE_object_boundbox_flag(Object *ob, int flag, const bool set)
2326 {
2327         BoundBox *bb = BKE_object_boundbox_get(ob);
2328         if (bb) {
2329                 if (set) bb->flag |= flag;
2330                 else bb->flag &= ~flag;
2331         }
2332 }
2333
2334 void BKE_object_dimensions_get(Object *ob, float vec[3])
2335 {
2336         BoundBox *bb = NULL;
2337         
2338         bb = BKE_object_boundbox_get(ob);
2339         if (bb) {
2340                 float scale[3];
2341                 
2342                 mat4_to_size(scale, ob->obmat);
2343                 
2344                 vec[0] = fabsf(scale[0]) * (bb->vec[4][0] - bb->vec[0][0]);
2345                 vec[1] = fabsf(scale[1]) * (bb->vec[2][1] - bb->vec[0][1]);
2346                 vec[2] = fabsf(scale[2]) * (bb->vec[1][2] - bb->vec[0][2]);
2347         }
2348         else {
2349                 zero_v3(vec);
2350         }
2351 }
2352
2353 void BKE_object_dimensions_set(Object *ob, const float value[3])
2354 {
2355         BoundBox *bb = NULL;
2356         
2357         bb = BKE_object_boundbox_get(ob);
2358         if (bb) {
2359                 float scale[3], len[3];
2360                 
2361                 mat4_to_size(scale, ob->obmat);
2362                 
2363                 len[0] = bb->vec[4][0] - bb->vec[0][0];
2364                 len[1] = bb->vec[2][1] - bb->vec[0][1];
2365                 len[2] = bb->vec[1][2] - bb->vec[0][2];
2366                 
2367                 if (len[0] > 0.f) ob->size[0] = value[0] / len[0];
2368                 if (len[1] > 0.f) ob->size[1] = value[1] / len[1];
2369                 if (len[2] > 0.f) ob->size[2] = value[2] / len[2];
2370         }
2371 }
2372
2373 void BKE_object_minmax(Object *ob, float min_r[3], float max_r[3], const bool use_hidden)
2374 {
2375         BoundBox bb;
2376         float vec[3];
2377         bool changed = false;
2378         
2379         switch (ob->type) {
2380                 case OB_CURVE:
2381                 case OB_FONT:
2382                 case OB_SURF:
2383                 {
2384                         bb = *BKE_curve_boundbox_get(ob);
2385                         BKE_boundbox_minmax(&bb, ob->obmat, min_r, max_r);
2386                         changed = true;
2387                         break;
2388                 }
2389                 case OB_LATTICE:
2390                 {
2391                         Lattice *lt = ob->data;
2392                         BPoint *bp = lt->def;
2393                         int u, v, w;
2394
2395                         for (w = 0; w < lt->pntsw; w++) {
2396                                 for (v = 0; v < lt->pntsv; v++) {
2397                                         for (u = 0; u < lt->pntsu; u++, bp++) {
2398                                                 mul_v3_m4v3(vec, ob->obmat, bp->vec);
2399                                                 minmax_v3v3_v3(min_r, max_r, vec);
2400                                         }
2401                                 }
2402                         }
2403                         changed = true;
2404                         break;
2405                 }
2406                 case OB_ARMATURE:
2407                 {
2408                         changed = BKE_pose_minmax(ob, min_r, max_r, use_hidden, false);
2409                         break;
2410                 }
2411                 case OB_MESH:
2412                 {
2413                         Mesh *me = BKE_mesh_from_object(ob);
2414
2415                         if (me) {
2416                                 bb = *BKE_mesh_boundbox_get(ob);
2417                                 BKE_boundbox_minmax(&bb, ob->obmat, min_r, max_r);
2418                                 changed = true;
2419                         }
2420                         break;
2421                 }
2422                 case OB_MBALL:
2423                 {
2424                         float ob_min[3], ob_max[3];
2425
2426                         changed = BKE_mball_minmax_ex(ob->data, ob_min, ob_max, ob->obmat, 0);
2427                         if (changed) {
2428                                 minmax_v3v3_v3(min_r, max_r, ob_min);
2429                                 minmax_v3v3_v3(min_r, max_r, ob_max);
2430                         }
2431                         break;
2432                 }
2433         }
2434
2435         if (changed == false) {
2436                 float size[3];
2437
2438                 copy_v3_v3(size, ob->size);
2439                 if (ob->type == OB_EMPTY) {
2440                         mul_v3_fl(size, ob->empty_drawsize);
2441                 }
2442
2443                 minmax_v3v3_v3(min_r, max_r, ob->obmat[3]);
2444
2445                 copy_v3_v3(vec, ob->obmat[3]);
2446                 add_v3_v3(vec, size);
2447                 minmax_v3v3_v3(min_r, max_r, vec);
2448
2449                 copy_v3_v3(vec, ob->obmat[3]);
2450                 sub_v3_v3(vec, size);
2451                 minmax_v3v3_v3(min_r, max_r, vec);
2452         }
2453 }
2454
2455 void BKE_object_empty_draw_type_set(Object *ob, const int value)
2456 {
2457         ob->empty_drawtype = value;
2458
2459         if (ob->type == OB_EMPTY && ob->empty_drawtype == OB_EMPTY_IMAGE) {
2460                 if (!ob->iuser) {
2461                         ob->iuser = MEM_callocN(sizeof(ImageUser), "image user");
2462                         ob->iuser->ok = 1;
2463                         ob->iuser->frames = 100;
2464                         ob->iuser->sfra = 1;
2465                         ob->iuser->fie_ima = 2;
2466                 }
2467         }
2468         else {
2469                 if (ob->iuser) {
2470                         MEM_freeN(ob->iuser);
2471                         ob->iuser = NULL;
2472                 }
2473         }
2474 }
2475
2476 bool BKE_object_minmax_dupli(Scene *scene, Object *ob, float r_min[3], float r_max[3], const bool use_hidden)
2477 {
2478         bool ok = false;
2479         if ((ob->transflag & OB_DUPLI) == 0) {
2480                 return ok;
2481         }
2482         else {
2483                 ListBase *lb;
2484                 DupliObject *dob;
2485                 lb = object_duplilist(G.main->eval_ctx, scene, ob);
2486                 for (dob = lb->first; dob; dob = dob->next) {
2487                         if ((use_hidden == false) && (dob->no_draw != 0)) {
2488                                 /* pass */
2489                         }
2490                         else {
2491                                 BoundBox *bb = BKE_object_boundbox_get(dob->ob);
2492
2493                                 if (bb) {
2494                                         int i;
2495                                         for (i = 0; i < 8; i++) {
2496                                                 float vec[3];
2497                                                 mul_v3_m4v3(vec, dob->mat, bb->vec[i]);
2498                                                 minmax_v3v3_v3(r_min, r_max, vec);
2499                                         }
2500
2501                                         ok = true;
2502                                 }
2503                         }
2504                 }
2505                 free_object_duplilist(lb);  /* does restore */
2506         }
2507
2508         return ok;
2509 }
2510
2511 void BKE_object_foreach_display_point(
2512         Object *ob, float obmat[4][4],
2513         void (*func_cb)(const float[3], void *), void *user_data)
2514 {
2515         float co[3];
2516
2517         if (ob->derivedFinal) {
2518                 DerivedMesh *dm = ob->derivedFinal;
2519                 MVert *mv = dm->getVertArray(dm);
2520                 int totvert = dm->getNumVerts(dm);
2521                 int i;
2522
2523                 for (i = 0; i < totvert; i++, mv++) {
2524                         mul_v3_m4v3(co, obmat, mv->co);
2525                         func_cb(co, user_data);
2526                 }
2527         }
2528         else if (ob->curve_cache && ob->curve_cache->disp.first) {
2529                 DispList *dl;
2530
2531                 for (dl = ob->curve_cache->disp.first; dl; dl = dl->next) {
2532                         const float *v3 = dl->verts;
2533                         int totvert = dl->nr;
2534                         int i;
2535
2536                         for (i = 0; i < totvert; i++, v3 += 3) {
2537                                 mul_v3_m4v3(co, obmat, v3);
2538                                 func_cb(co, user_data);
2539                         }
2540                 }
2541         }
2542 }
2543
2544 void BKE_scene_foreach_display_point(
2545         Scene *scene, View3D *v3d, const short flag,
2546         void (*func_cb)(const float[3], void *), void *user_data)
2547 {
2548         Base *base;
2549         Object *ob;
2550
2551         for (base = FIRSTBASE; base; base = base->next) {
2552                 if (BASE_VISIBLE_BGMODE(v3d, scene, base) && (base->flag & flag) == flag) {
2553                         ob = base->object;
2554
2555                         if ((ob->transflag & OB_DUPLI) == 0) {
2556                                 BKE_object_foreach_display_point(ob, ob->obmat, func_cb, user_data);
2557                         }
2558                         else {
2559                                 ListBase *lb;
2560                                 DupliObject *dob;
2561
2562                                 lb = object_duplilist(G.main->eval_ctx, scene, ob);
2563                                 for (dob = lb->first; dob; dob = dob->next) {
2564                                         if (dob->no_draw == 0) {
2565                                                 BKE_object_foreach_display_point(dob->ob, dob->mat, func_cb, user_data);
2566                                         }
2567                                 }
2568                                 free_object_duplilist(lb);  /* does restore */
2569                         }
2570                 }
2571         }
2572 }
2573
2574 /* copied from DNA_object_types.h */
2575 typedef struct ObTfmBack {
2576         float loc[3], dloc[3], orig[3];
2577         float size[3], dscale[3];   /* scale and delta scale */
2578         float rot[3], drot[3];      /* euler rotation */
2579         float quat[4], dquat[4];    /* quaternion rotation */
2580         float rotAxis[3], drotAxis[3];  /* axis angle rotation - axis part */
2581         float rotAngle, drotAngle;  /* axis angle rotation - angle part */
2582         float obmat[4][4];      /* final worldspace matrix with constraints & animsys applied */
2583         float parentinv[4][4]; /* inverse result of parent, so that object doesn't 'stick' to parent */
2584         float constinv[4][4]; /* inverse result of constraints. doesn't include effect of parent or object local transform */
2585         float imat[4][4];   /* inverse matrix of 'obmat' for during render, old game engine, temporally: ipokeys of transform  */
2586 } ObTfmBack;
2587
2588 void *BKE_object_tfm_backup(Object *ob)
2589 {
2590         ObTfmBack *obtfm = MEM_mallocN(sizeof(ObTfmBack), "ObTfmBack");
2591         copy_v3_v3(obtfm->loc, ob->loc);
2592         copy_v3_v3(obtfm->dloc, ob->dloc);
2593         copy_v3_v3(obtfm->orig, ob->orig);
2594         copy_v3_v3(obtfm->size, ob->size);
2595         copy_v3_v3(obtfm->dscale, ob->dscale);
2596         copy_v3_v3(obtfm->rot, ob->rot);
2597         copy_v3_v3(obtfm->drot, ob->drot);
2598         copy_qt_qt(obtfm->quat, ob->quat);
2599         copy_qt_qt(obtfm->dquat, ob->dquat);
2600         copy_v3_v3(obtfm->rotAxis, ob->rotAxis);
2601         copy_v3_v3(obtfm->drotAxis, ob->drotAxis);
2602         obtfm->rotAngle = ob->rotAngle;
2603         obtfm->drotAngle = ob->drotAngle;
2604         copy_m4_m4(obtfm->obmat, ob->obmat);
2605         copy_m4_m4(obtfm->parentinv, ob->parentinv);
2606         copy_m4_m4(obtfm->constinv, ob->constinv);
2607         copy_m4_m4(obtfm->imat, ob->imat);
2608
2609         return (void *)obtfm;
2610 }
2611
2612 void BKE_object_tfm_restore(Object *ob, void *obtfm_pt)
2613 {
2614         ObTfmBack *obtfm = (ObTfmBack *)obtfm_pt;
2615         copy_v3_v3(ob->loc, obtfm->loc);
2616         copy_v3_v3(ob->dloc, obtfm->dloc);
2617         copy_v3_v3(ob->orig, obtfm->orig);
2618         copy_v3_v3(ob->size, obtfm->size);
2619         copy_v3_v3(ob->dscale, obtfm->dscale);
2620         copy_v3_v3(ob->rot, obtfm->rot);
2621         copy_v3_v3(ob->drot, obtfm->drot);
2622         copy_qt_qt(ob->quat, obtfm->quat);
2623         copy_qt_qt(ob->dquat, obtfm->dquat);
2624         copy_v3_v3(ob->rotAxis, obtfm->rotAxis);
2625         copy_v3_v3(ob->drotAxis, obtfm->drotAxis);
2626         ob->rotAngle = obtfm->rotAngle;
2627         ob->drotAngle = obtfm->drotAngle;
2628         copy_m4_m4(ob->obmat, obtfm->obmat);
2629         copy_m4_m4(ob->parentinv, obtfm->parentinv);
2630         copy_m4_m4(ob->constinv, obtfm->constinv);
2631         copy_m4_m4(ob->imat, obtfm->imat);
2632 }
2633
2634 bool BKE_object_parent_loop_check(const Object *par, const Object *ob)
2635 {
2636         /* test if 'ob' is a parent somewhere in par's parents */
2637         if (par == NULL) return false;
2638         if (ob == par) return true;
2639         return BKE_object_parent_loop_check(par->parent, ob);
2640 }
2641
2642 /* proxy rule: lib_object->proxy_from == the one we borrow from, only set temporal and cleared here */
2643 /*           local_object->proxy      == pointer to library object, saved in files and read */
2644
2645 /* function below is polluted with proxy exceptions, cleanup will follow! */
2646
2647 /* the main object update call, for object matrix, constraints, keys and displist (modifiers) */
2648 /* requires flags to be set! */
2649 /* Ideally we shouldn't have to pass the rigid body world, but need bigger restructuring to avoid id */
2650 void BKE_object_handle_update_ex(EvaluationContext *eval_ctx,
2651                                  Scene *scene, Object *ob,
2652                                  RigidBodyWorld *rbw,
2653                                  const bool do_proxy_update)
2654 {
2655         if (ob->recalc & OB_RECALC_ALL) {
2656                 /* speed optimization for animation lookups */
2657                 if (ob->pose) {
2658                         BKE_pose_channels_hash_make(ob->pose);
2659                         if (ob->pose->flag & POSE_CONSTRAINTS_NEED_UPDATE_FLAGS) {
2660                                 BKE_pose_update_constraint_flags(ob->pose);
2661                         }
2662                 }
2663
2664                 if (ob->recalc & OB_RECALC_DATA) {
2665                         if (ob->type == OB_ARMATURE) {
2666                                 /* this happens for reading old files and to match library armatures
2667                                  * with poses we do it ahead of BKE_object_where_is_calc to ensure animation
2668                                  * is evaluated on the rebuilt pose, otherwise we get incorrect poses
2669                                  * on file load */
2670                                 if (ob->pose == NULL || (ob->pose->flag & POSE_RECALC))
2671                                         BKE_pose_rebuild(ob, ob->data);
2672                         }
2673                 }
2674
2675                 /* XXX new animsys warning: depsgraph tag OB_RECALC_DATA should not skip drivers, 
2676                  * which is only in BKE_object_where_is_calc now */
2677                 /* XXX: should this case be OB_RECALC_OB instead? */
2678                 if (ob->recalc & OB_RECALC_ALL) {
2679                         
2680                         if (G.debug & G_DEBUG_DEPSGRAPH)
2681                                 printf("recalcob %s\n", ob->id.name + 2);
2682                         
2683                         /* handle proxy copy for target */
2684                         if (ID_IS_LINKED_DATABLOCK(ob) && ob->proxy_from) {
2685                                 // printf("ob proxy copy, lib ob %s proxy %s\n", ob->id.name, ob->proxy_from->id.name);
2686                                 if (ob->proxy_from->proxy_group) { /* transform proxy into group space */
2687                                         Object *obg = ob->proxy_from->proxy_group;
2688                                         float imat[4][4];
2689                                         invert_m4_m4(imat, obg->obmat);
2690                                         mul_m4_m4m4(ob->obmat, imat, ob->proxy_from->obmat);
2691                                         if (obg->dup_group) { /* should always be true */
2692                                                 add_v3_v3(ob->obmat[3], obg->dup_group->dupli_ofs);
2693                                         }
2694                                 }
2695                                 else
2696                                         copy_m4_m4(ob->obmat, ob->proxy_from->obmat);
2697                         }
2698                         else
2699                                 BKE_object_where_is_calc_ex(scene, rbw, ob, NULL);
2700                 }
2701                 
2702                 if (ob->recalc & OB_RECALC_DATA) {
2703                         BKE_object_handle_data_update(eval_ctx, scene, ob);
2704                 }
2705
2706                 ob->recalc &= ~OB_RECALC_ALL;
2707         }
2708
2709         /* the case when this is a group proxy, object_update is called in group.c */
2710         if (ob->proxy) {
2711                 /* set pointer in library proxy target, for copying, but restore it */
2712                 ob->proxy->proxy_from = ob;
2713                 // printf("set proxy pointer for later group stuff %s\n", ob->id.name);
2714
2715                 /* the no-group proxy case, we call update */
2716                 if (ob->proxy_group == NULL) {
2717                         if (do_proxy_update) {
2718                                 // printf("call update, lib ob %s proxy %s\n", ob->proxy->id.name, ob->id.name);
2719                                 BKE_object_handle_update(eval_ctx, scene, ob->proxy);
2720                         }
2721                 }
2722         }
2723 }
2724 /* WARNING: "scene" here may not be the scene object actually resides in. 
2725  * When dealing with background-sets, "scene" is actually the active scene.
2726  * e.g. "scene" <-- set 1 <-- set 2 ("ob" lives here) <-- set 3 <-- ... <-- set n
2727  * rigid bodies depend on their world so use BKE_object_handle_update_ex() to also pass along the corrent rigid body world
2728  */
2729 void BKE_object_handle_update(EvaluationContext *eval_ctx, Scene *scene, Object *ob)
2730 {
2731         BKE_object_handle_update_ex(eval_ctx, scene, ob, NULL, true);
2732 }
2733
2734 void BKE_object_sculpt_modifiers_changed(Object *ob)
2735 {
2736         SculptSession *ss = ob->sculpt;
2737
2738         if (ss) {
2739                 if (!ss->cache) {
2740                         /* we free pbvh on changes, except during sculpt since it can't deal with
2741                          * changing PVBH node organization, we hope topology does not change in
2742                          * the meantime .. weak */
2743                         if (ss->pbvh) {
2744                                 BKE_pbvh_free(ss->pbvh);
2745                                 ss->pbvh = NULL;
2746                         }
2747
2748                         BKE_sculptsession_free_deformMats(ob->sculpt);
2749                 }
2750                 else {
2751                         PBVHNode **nodes;
2752                         int n, totnode;
2753
2754                         BKE_pbvh_search_gather(ss->pbvh, NULL, NULL, &nodes, &totnode);
2755
2756                         for (n = 0; n < totnode; n++)
2757                                 BKE_pbvh_node_mark_update(nodes[n]);
2758
2759                         MEM_freeN(nodes);
2760                 }
2761         }
2762 }
2763
2764 int BKE_object_obdata_texspace_get(Object *ob, short **r_texflag, float **r_loc, float **r_size, float **r_rot)
2765 {
2766         
2767         if (ob->data == NULL)
2768                 return 0;
2769         
2770         switch (GS(((ID *)ob->data)->name)) {
2771                 case ID_ME:
2772                 {
2773                         Mesh *me = ob->data;
2774                         if (me->bb == NULL || (me->bb->flag & BOUNDBOX_DIRTY)) {
2775                                 BKE_mesh_texspace_calc(me);
2776                         }
2777                         if (r_texflag) *r_texflag = &me->texflag;
2778                         if (r_loc) *r_loc = me->loc;
2779                         if (r_size) *r_size = me->size;
2780                         if (r_rot) *r_rot = me->rot;
2781                         break;
2782                 }
2783                 case ID_CU:
2784                 {
2785                         Curve *cu = ob->data;
2786                         if (cu->bb == NULL || (cu->bb->flag & BOUNDBOX_DIRTY)) {
2787                                 BKE_curve_texspace_calc(cu);
2788                         }
2789                         if (r_texflag) *r_texflag = &cu->texflag;
2790                         if (r_loc) *r_loc = cu->loc;
2791                         if (r_size) *r_size = cu->size;
2792                         if (r_rot) *r_rot = cu->rot;
2793                         break;
2794                 }
2795                 case ID_MB:
2796                 {
2797                         MetaBall *mb = ob->data;
2798                         if (r_texflag) *r_texflag = &mb->texflag;
2799                         if (r_loc) *r_loc = mb->loc;
2800                         if (r_size) *r_size = mb->size;
2801                         if (r_rot) *r_rot = mb->rot;
2802                         break;
2803                 }
2804                 default:
2805                         return 0;
2806         }
2807         return 1;
2808 }
2809
2810 static int pc_cmp(const void *a, const void *b)
2811 {
2812         const LinkData *ad = a, *bd = b;
2813         if (GET_INT_FROM_POINTER(ad->data) > GET_INT_FROM_POINTER(bd->data))
2814                 return 1;
2815         else return 0;
2816 }
2817
2818 int BKE_object_insert_ptcache(Object *ob) 
2819 {
2820         LinkData *link = NULL;
2821         int i = 0;
2822
2823         BLI_listbase_sort(&ob->pc_ids, pc_cmp);
2824
2825         for (link = ob->pc_ids.first, i = 0; link; link = link->next, i++) {
2826                 int index = GET_INT_FROM_POINTER(link->data);
2827
2828                 if (i < index)
2829                         break;
2830         }
2831
2832         link = MEM_callocN(sizeof(LinkData), "PCLink");
2833         link->data = SET_INT_IN_POINTER(i);
2834         BLI_addtail(&ob->pc_ids, link);
2835
2836         return i;
2837 }
2838
2839 static int pc_findindex(ListBase *listbase, int index)
2840 {
2841         LinkData *link = NULL;
2842         int number = 0;
2843         
2844         if (listbase == NULL) return -1;
2845         
2846         link = listbase->first;
2847         while (link) {
2848                 if (GET_INT_FROM_POINTER(link->data) == index)
2849                         return number;
2850                 
2851                 number++;
2852                 link = link->next;
2853         }
2854         
2855         return -1;
2856 }
2857
2858 void BKE_object_delete_ptcache(Object *ob, int index)
2859 {
2860         int list_index = pc_findindex(&ob->pc_ids, index);
2861         LinkData *link = BLI_findlink(&ob->pc_ids, list_index);
2862         BLI_freelinkN(&ob->pc_ids, link);
2863 }
2864
2865 /* shape key utility function */
2866
2867 /************************* Mesh ************************/
2868 static KeyBlock *insert_meshkey(Object *ob, const char *name, const bool from_mix)
2869 {
2870         Mesh *me = ob->data;
2871         Key *key = me->key;
2872         KeyBlock *kb;
2873         int newkey = 0;
2874
2875         if (key == NULL) {
2876                 key = me->key = BKE_key_add((ID *)me);
2877                 key->type = KEY_RELATIVE;
2878                 newkey = 1;
2879         }
2880
2881         if (newkey || from_mix == false) {
2882                 /* create from mesh */
2883                 kb = BKE_keyblock_add_ctime(key, name, false);
2884                 BKE_keyblock_convert_from_mesh(me, kb);
2885         }
2886         else {
2887                 /* copy from current values */
2888                 int totelem;
2889                 float *data = BKE_key_evaluate_object(ob, &totelem);
2890
2891                 /* create new block with prepared data */
2892                 kb = BKE_keyblock_add_ctime(key, name, false);
2893                 kb->data = data;
2894                 kb->totelem = totelem;
2895         }
2896
2897         return kb;
2898 }
2899 /************************* Lattice ************************/
2900 static KeyBlock *insert_lattkey(Object *ob, const char *name, const bool from_mix)
2901 {
2902         Lattice *lt = ob->data;
2903         Key *key = lt->key;
2904         KeyBlock *kb;
2905         int newkey = 0;
2906
2907         if (key == NULL) {
2908                 key = lt->key = BKE_key_add((ID *)lt);
2909                 key->type = KEY_RELATIVE;
2910                 newkey = 1;
2911         }
2912
2913         if (newkey || from_mix == false) {
2914                 kb = BKE_keyblock_add_ctime(key, name, false);
2915                 if (!newkey) {
2916                         KeyBlock *basekb = (KeyBlock *)key->block.first;
2917                         kb->data = MEM_dupallocN(basekb->data);
2918                         kb->totelem = basekb->totelem;
2919                 }
2920                 else {
2921                         BKE_keyblock_convert_from_lattice(lt, kb);
2922                 }
2923         }
2924         else {
2925                 /* copy from current values */
2926                 int totelem;
2927                 float *data = BKE_key_evaluate_object(ob, &totelem);
2928
2929                 /* create new block with prepared data */
2930                 kb = BKE_keyblock_add_ctime(key, name, false);
2931                 kb->totelem = totelem;
2932                 kb->data = data;
2933         }
2934
2935         return kb;
2936 }
2937 /************************* Curve ************************/
2938 static KeyBlock *insert_curvekey(Object *ob, const char *name, const bool from_mix)
2939 {
2940         Curve *cu = ob->data;
2941         Key *key = cu->key;
2942         KeyBlock *kb;
2943         ListBase *lb = BKE_curve_nurbs_get(cu);
2944         int newkey = 0;
2945
2946         if (key == NULL) {
2947                 key = cu->key = BKE_key_add((ID *)cu);
2948                 key->type = KEY_RELATIVE;
2949                 newkey = 1;
2950         }
2951
2952         if (newkey || from_mix == false) {
2953                 /* create from curve */
2954                 kb = BKE_keyblock_add_ctime(key, name, false);
2955                 if (!newkey) {
2956                         KeyBlock *basekb = (KeyBlock *)key->block.first;
2957                         kb->data = MEM_dupallocN(basekb->data);
2958                         kb->totelem = basekb->totelem;
2959                 }
2960                 else {
2961                         BKE_keyblock_convert_from_curve(cu, kb, lb);
2962                 }
2963         }
2964         else {
2965                 /* copy from current values */
2966                 int totelem;
2967                 float *data = BKE_key_evaluate_object(ob, &totelem);
2968
2969                 /* create new block with prepared data */
2970                 kb = BKE_keyblock_add_ctime(key, name, false);
2971                 kb->totelem = totelem;
2972                 kb->data = data;
2973         }
2974
2975         return kb;
2976 }
2977
2978 KeyBlock *BKE_object_shapekey_insert(Object *ob, const char *name, const bool from_mix)
2979 {       
2980         switch (ob->type) {
2981                 case OB_MESH:
2982                         return insert_meshkey(ob, name, from_mix);
2983                 case OB_CURVE:
2984                 case OB_SURF:
2985                         return insert_curvekey(ob, name, from_mix);
2986                 case OB_LATTICE:
2987                         return insert_lattkey(ob, name, from_mix);
2988                 default:
2989                         return NULL;
2990         }
2991
2992 }
2993
2994 bool BKE_object_shapekey_free(Main *bmain, Object *ob)
2995 {
2996         Key **key_p, *key;
2997
2998         key_p = BKE_key_from_object_p(ob);
2999         if (ELEM(NULL, key_p, *key_p)) {
3000                 return false;
3001         }
3002
3003         key = *key_p;
3004         *key_p = NULL;
3005
3006         BKE_libblock_free_us(bmain, key);
3007
3008         return false;
3009 }
3010
3011 bool BKE_object_shapekey_remove(Main *bmain, Object *ob, KeyBlock *kb)
3012 {
3013         KeyBlock *rkb;
3014         Key *key = BKE_key_from_object(ob);
3015         short kb_index;
3016
3017         if (key == NULL) {
3018                 return false;
3019         }
3020
3021         kb_index = BLI_findindex(&key->block, kb);
3022         BLI_assert(kb_index != -1);
3023
3024         for (rkb = key->block.first; rkb; rkb = rkb->next) {
3025                 if (rkb->relative == kb_index) {
3026                         /* remap to the 'Basis' */
3027                         rkb->relative = 0;
3028                 }
3029                 else if (rkb->relative >= kb_index) {
3030                         /* Fix positional shift of the keys when kb is deleted from the list */
3031                         rkb->relative -= 1;
3032                 }
3033         }
3034
3035         BLI_remlink(&key->block, kb);
3036         key->totkey--;
3037         if (key->refkey == kb) {
3038                 key->refkey = key->block.first;
3039
3040                 if (key->refkey) {
3041                         /* apply new basis key on original data */
3042                         switch (ob->type) {
3043                                 case OB_MESH:
3044                                         BKE_keyblock_convert_to_mesh(key->refkey, ob->data);
3045                                         break;
3046                                 case OB_CURVE:
3047                                 case OB_SURF:
3048                                         BKE_keyblock_convert_to_curve(key->refkey, ob->data, BKE_curve_nurbs_get(ob->data));
3049                                         break;
3050                                 case OB_LATTICE:
3051                                         BKE_keyblock_convert_to_lattice(key->refkey, ob->data);
3052                                         break;
3053                         }
3054                 }
3055         }
3056
3057         if (kb->data) {
3058                 MEM_freeN(kb->data);
3059         }
3060         MEM_freeN(kb);
3061
3062         if (ob->shapenr > 1) {
3063                 ob->shapenr--;
3064         }
3065
3066         if (key->totkey == 0) {
3067                 BKE_object_shapekey_free(bmain, ob);
3068         }
3069
3070         return true;
3071 }
3072
3073 bool BKE_object_flag_test_recursive(const Object *ob, short flag)
3074 {
3075         if (ob->flag & flag) {
3076                 return true;
3077         }
3078         else if (ob->parent) {
3079                 return BKE_object_flag_test_recursive(ob->parent, flag);
3080         }
3081         else {
3082                 return false;
3083         }
3084 }
3085
3086 bool BKE_object_is_child_recursive(Object *ob_parent, Object *ob_child)
3087 {
3088         for (ob_child = ob_child->parent; ob_child; ob_child = ob_child->parent) {
3089                 if (ob_child == ob_parent) {
3090                         return true;
3091                 }
3092         }
3093         return false;
3094 }
3095
3096 /* most important if this is modified it should _always_ return True, in certain
3097  * cases false positives are hard to avoid (shape keys for example) */
3098 int BKE_object_is_modified(Scene *scene, Object *ob)
3099 {
3100         int flag = 0;
3101
3102         if (BKE_key_from_object(ob)) {
3103                 flag |= eModifierMode_Render | eModifierMode_Realtime;
3104         }
3105         else {
3106                 ModifierData *md;
3107                 VirtualModifierData virtualModifierData;
3108                 /* cloth */
3109                 for (md = modifiers_getVirtualModifierList(ob, &virtualModifierData);
3110                      md && (flag != (eModifierMode_Render | eModifierMode_Realtime));
3111                      md = md->next)
3112                 {
3113                         if ((flag & eModifierMode_Render) == 0 && modifier_isEnabled(scene, md, eModifierMode_Render))
3114                                 flag |= eModifierMode_Render;
3115
3116                         if ((flag & eModifierMode_Realtime) == 0 && modifier_isEnabled(scene, md, eModifierMode_Realtime))
3117                                 flag |= eModifierMode_Realtime;
3118                 }
3119         }
3120
3121         return flag;
3122 }
3123
3124 /* Check of objects moves in time. */
3125 /* NOTE: This function is currently optimized for usage in combination
3126  * with mti->canDeform, so modifiers can quickly check if their target
3127  * objects moves (causing deformation motion blur) or not.
3128  *
3129  * This makes it possible to give some degree of false-positives here,
3130  * but it's currently an acceptable tradeoff between complexity and check
3131  * speed. In combination with checks of modifier stack and real life usage
3132  * percentage of false-positives shouldn't be that hight.
3133  */
3134 static bool object_moves_in_time(Object *object)
3135 {
3136         AnimData *adt = object->adt;
3137         if (adt != NULL) {
3138                 /* If object has any sort of animation data assume it is moving. */
3139                 if (adt->action != NULL ||
3140                     !BLI_listbase_is_empty(&adt->nla_tracks) ||
3141                     !BLI_listbase_is_empty(&adt->drivers) ||
3142                     !BLI_listbase_is_empty(&adt->overrides))
3143                 {
3144                         return true;
3145                 }
3146         }
3147         if (!BLI_listbase_is_empty(&object->constraints)) {
3148                 return true;
3149         }
3150         if (object->parent != NULL) {
3151                 /* TODO(sergey): Do recursive check here? */
3152                 return true;
3153         }
3154         return false;
3155 }
3156
3157 static bool object_deforms_in_time(Object *object)
3158 {
3159         if (BKE_key_from_object(object) != NULL) {
3160                 return true;
3161         }
3162         if (!BLI_listbase_is_empty(&object->modifiers)) {
3163                 return true;
3164         }
3165         return object_moves_in_time(object);
3166 }
3167
3168 static bool constructive_modifier_is_deform_modified(ModifierData *md)
3169 {
3170         /* TODO(sergey): Consider generalizing this a bit so all modifier logic
3171          * is concentrated in MOD_{modifier}.c file,
3172          */
3173         if (md->type == eModifierType_Array) {
3174                 ArrayModifierData *amd = (ArrayModifierData *)md;
3175                 /* TODO(sergey): Check if curve is deformed. */
3176                 return (amd->start_cap != NULL && object_moves_in_time(amd->start_cap)) ||
3177                        (amd->end_cap != NULL && object_moves_in_time(amd->end_cap)) ||
3178                        (amd->curve_ob != NULL && object_moves_in_time(amd->curve_ob)) ||
3179                        (amd->offset_ob != NULL && object_moves_in_time(amd->offset_ob));
3180         }
3181         else if (md->type == eModifierType_Mirror) {
3182                 MirrorModifierData *mmd = (MirrorModifierData *)md;
3183                 return mmd->mirror_ob != NULL && object_moves_in_time(mmd->mirror_ob);
3184         }
3185         else if (md->type == eModifierType_Screw) {
3186                 ScrewModifierData *smd = (ScrewModifierData *)md;
3187                 return smd->ob_axis != NULL && object_moves_in_time(smd->ob_axis);
3188         }
3189         else if (md->type == eModifierType_MeshSequenceCache) {
3190                 /* NOTE: Not ideal because it's unknown whether topology changes or not.
3191                  * This will be detected later, so by assuming it's only deformation
3192                  * going on here we allow to bake deform-only mesh to Alembic and have
3193                  * proper motion blur after that.
3194                  */
3195                 return true;
3196         }
3197         return false;
3198 }
3199
3200 static bool modifiers_has_animation_check(Object *ob)
3201 {
3202         /* TODO(sergey): This is a bit code duplication with depsgraph, but
3203          * would be nicer to solve this as a part of new dependency graph
3204          * work, so we avoid conflicts and so.
3205          */
3206         if (ob->adt != NULL) {
3207                 AnimData *adt = ob->adt;
3208                 FCurve *fcu;
3209                 if (adt->action != NULL) {
3210                         for (fcu = adt->action->curves.first; fcu; fcu = fcu->next) {
3211                                 if (fcu->rna_path && strstr(fcu->rna_path, "modifiers[")) {
3212                                         return true;
3213                                 }
3214                         }
3215                 }
3216                 for (fcu = adt->drivers.first; fcu; fcu = fcu->next) {
3217                         if (fcu->rna_path && strstr(fcu->rna_path, "modifiers[")) {
3218                                 return true;
3219                         }
3220                 }
3221         }
3222         return false;
3223 }
3224
3225 /* test if object is affected by deforming modifiers (for motion blur). again
3226  * most important is to avoid false positives, this is to skip computations
3227  * and we can still if there was actual deformation afterwards */
3228 int BKE_object_is_deform_modified(Scene *scene, Object *ob)
3229 {
3230         ModifierData *md;
3231         VirtualModifierData virtualModifierData;
3232         int flag = 0;
3233         const bool is_modifier_animated = modifiers_has_animation_check(ob);
3234
3235         if (BKE_key_from_object(ob)) {
3236                 flag |= eModifierMode_Realtime | eModifierMode_Render;
3237         }
3238
3239         if (ob->type == OB_CURVE) {
3240                 Curve *cu = (Curve *)ob->data;
3241                 if (cu->taperobj != NULL && object_deforms_in_time(cu->taperobj)) {
3242                         flag |= eModifierMode_Realtime | eModifierMode_Render;
3243                 }
3244         }
3245
3246         /* cloth */
3247         for (md = modifiers_getVirtualModifierList(ob, &virtualModifierDat