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