svn merge https://svn.blender.org/svnroot/bf-blender/trunk/blender -r22704:22717
[blender.git] / source / blender / blenkernel / intern / action.c
1 /**
2  * $Id$
3  *
4  * ***** BEGIN GPL LICENSE BLOCK *****
5  *
6  * This program is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU General Public License
8  * as published by the Free Software Foundation; either version 2
9  * of the License, or (at your option) any later version.
10  *
11  * This program is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14  * GNU General Public License for more details.
15  *
16  * You should have received a copy of the GNU General Public License
17  * along with this program; if not, write to the Free Software Foundation,
18  * Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
19  *
20  * The Original Code is Copyright (C) 2001-2002 by NaN Holding BV.
21  * All rights reserved.
22  *
23  * Contributor(s): Full recode, Ton Roosendaal, Crete 2005
24  *                               Full recode, Joshua Leung, 2009
25  *
26  * ***** END GPL LICENSE BLOCK *****
27  */
28
29 #ifdef HAVE_CONFIG_H
30 #include <config.h>
31 #endif
32
33 #include <string.h>
34 #include <math.h>
35 #include <stdlib.h>
36 #include <stddef.h>     
37
38 #include "MEM_guardedalloc.h"
39
40 #include "DNA_anim_types.h"
41 #include "DNA_action_types.h"
42 #include "DNA_armature_types.h"
43 #include "DNA_constraint_types.h"
44 #include "DNA_curve_types.h"
45 #include "DNA_key_types.h"
46 #include "DNA_nla_types.h"
47 #include "DNA_object_types.h"
48 #include "DNA_scene_types.h"
49
50 #include "BKE_animsys.h"
51 #include "BKE_action.h"
52 #include "BKE_anim.h"
53 #include "BKE_armature.h"
54 #include "BKE_blender.h"
55 #include "BKE_constraint.h"
56 #include "BKE_displist.h"
57 #include "BKE_global.h"
58 #include "BKE_fcurve.h"
59 #include "BKE_key.h"
60 #include "BKE_lattice.h"
61 #include "BKE_library.h"
62 #include "BKE_main.h"
63 #include "BKE_object.h"
64 #include "BKE_utildefines.h"
65
66 #include "BLI_arithb.h"
67 #include "BLI_blenlib.h"
68 #include "BLI_ghash.h"
69
70 #include "RNA_access.h"
71 #include "RNA_types.h"
72
73 /* *********************** NOTE ON POSE AND ACTION **********************
74
75   - Pose is the local (object level) component of armature. The current
76     object pose is saved in files, and (will be) is presorted for dependency
77   - Actions have fewer (or other) channels, and write data to a Pose
78   - Currently ob->pose data is controlled in where_is_pose only. The (recalc)
79     event system takes care of calling that
80   - The NLA system (here too) uses Poses as interpolation format for Actions
81   - Therefore we assume poses to be static, and duplicates of poses have channels in
82     same order, for quick interpolation reasons
83
84   ****************************** (ton) ************************************ */
85
86 /* ***************** Library data level operations on action ************** */
87
88 bAction *add_empty_action(const char name[])
89 {
90         bAction *act;
91         
92         act= alloc_libblock(&G.main->action, ID_AC, name);
93         act->id.flag |= LIB_FAKEUSER; // XXX this is nasty for new users... maybe we don't want this anymore
94         act->id.us++;
95         
96         return act;
97 }       
98
99 // does copy_fcurve...
100 void make_local_action(bAction *act)
101 {
102         // Object *ob;
103         bAction *actn;
104         int local=0, lib=0;
105         
106         if (act->id.lib==0) return;
107         if (act->id.us==1) {
108                 act->id.lib= 0;
109                 act->id.flag= LIB_LOCAL;
110                 //make_local_action_channels(act);
111                 new_id(0, (ID *)act, 0);
112                 return;
113         }
114         
115 #if 0   // XXX old animation system
116         ob= G.main->object.first;
117         while(ob) {
118                 if(ob->action==act) {
119                         if(ob->id.lib) lib= 1;
120                         else local= 1;
121                 }
122                 ob= ob->id.next;
123         }
124 #endif
125         
126         if(local && lib==0) {
127                 act->id.lib= 0;
128                 act->id.flag= LIB_LOCAL;
129                 //make_local_action_channels(act);
130                 new_id(0, (ID *)act, 0);
131         }
132         else if(local && lib) {
133                 actn= copy_action(act);
134                 actn->id.us= 0;
135                 
136 #if 0   // XXX old animation system
137                 ob= G.main->object.first;
138                 while(ob) {
139                         if(ob->action==act) {
140                                 
141                                 if(ob->id.lib==0) {
142                                         ob->action = actn;
143                                         actn->id.us++;
144                                         act->id.us--;
145                                 }
146                         }
147                         ob= ob->id.next;
148                 }
149 #endif  // XXX old animation system
150         }
151 }
152
153
154 void free_action (bAction *act)
155 {
156         /* sanity check */
157         if (act == NULL)
158                 return;
159         
160         /* Free F-Curves */
161         free_fcurves(&act->curves);
162         
163         /* Free groups */
164         if (act->groups.first)
165                 BLI_freelistN(&act->groups);
166                 
167         /* Free pose-references (aka local markers) */
168         if (act->markers.first)
169                 BLI_freelistN(&act->markers);
170 }
171
172 bAction *copy_action (bAction *src)
173 {
174         bAction *dst = NULL;
175         bActionGroup *dgrp, *sgrp;
176         FCurve *dfcu, *sfcu;
177         
178         if (src == NULL) 
179                 return NULL;
180         dst= copy_libblock(src);
181         
182         /* duplicate the lists of groups and markers */
183         BLI_duplicatelist(&dst->groups, &src->groups);
184         BLI_duplicatelist(&dst->markers, &src->markers);
185         
186         /* copy F-Curves, fixing up the links as we go */
187         dst->curves.first= dst->curves.last= NULL;
188         
189         for (sfcu= src->curves.first; sfcu; sfcu= sfcu->next) {
190                 /* duplicate F-Curve */
191                 dfcu= copy_fcurve(sfcu);
192                 BLI_addtail(&dst->curves, dfcu);
193                 
194                 /* fix group links (kindof bad list-in-list search, but this is the most reliable way) */
195                 for (dgrp=dst->groups.first, sgrp=src->groups.first; dgrp && sgrp; dgrp=dgrp->next, sgrp=sgrp->next) {
196                         if (sfcu->grp == sgrp) {
197                                 dfcu->grp= dgrp;
198                                 
199                                 if (dgrp->channels.first == sfcu)
200                                         dgrp->channels.first= dfcu;
201                                 if (dgrp->channels.last == sfcu)
202                                         dgrp->channels.last= dfcu;
203                                         
204                                 break;
205                         }
206                 }
207         }
208         
209         dst->id.flag |= LIB_FAKEUSER; // XXX this is nasty for new users... maybe we don't want this anymore
210         dst->id.us++;
211         
212         return dst;
213 }
214
215
216 /* Get the active action-group for an Action */
217 bActionGroup *get_active_actiongroup (bAction *act)
218 {
219         bActionGroup *agrp= NULL;
220         
221         if (act && act->groups.first) { 
222                 for (agrp= act->groups.first; agrp; agrp= agrp->next) {
223                         if (agrp->flag & AGRP_ACTIVE)
224                                 break;
225                 }
226         }
227         
228         return agrp;
229 }
230
231 /* Make the given Action-Group the active one */
232 void set_active_action_group (bAction *act, bActionGroup *agrp, short select)
233 {
234         bActionGroup *grp;
235         
236         /* sanity checks */
237         if (act == NULL)
238                 return;
239         
240         /* Deactive all others */
241         for (grp= act->groups.first; grp; grp= grp->next) {
242                 if ((grp==agrp) && (select))
243                         grp->flag |= AGRP_ACTIVE;
244                 else    
245                         grp->flag &= ~AGRP_ACTIVE;
246         }
247 }
248
249 /* Add given channel into (active) group 
250  *      - assumes that channel is not linked to anything anymore
251  *      - always adds at the end of the group 
252  */
253 void action_groups_add_channel (bAction *act, bActionGroup *agrp, FCurve *fcurve)
254 {
255         FCurve *fcu;
256         short done=0;
257         
258         /* sanity checks */
259         if (ELEM3(NULL, act, agrp, fcurve))
260                 return;
261         
262         /* if no channels, just add to two lists at the same time */
263         if (act->curves.first == NULL) {
264                 fcurve->next = fcurve->prev = NULL;
265                 
266                 agrp->channels.first = agrp->channels.last = fcurve;
267                 act->curves.first = act->curves.last = fcurve;
268                 
269                 fcurve->grp= agrp;
270                 return;
271         }
272         
273         /* try to find a channel to slot this in before/after */
274         for (fcu= act->curves.first; fcu; fcu= fcu->next) {
275                 /* if channel has no group, then we have ungrouped channels, which should always occur after groups */
276                 if (fcu->grp == NULL) {
277                         BLI_insertlinkbefore(&act->curves, fcu, fcurve);
278                         
279                         if (agrp->channels.first == NULL)
280                                 agrp->channels.first= fcurve;
281                         agrp->channels.last= fcurve;
282                         
283                         done= 1;
284                         break;
285                 }
286                 
287                 /* if channel has group after current, we can now insert (otherwise we have gone too far) */
288                 else if (fcu->grp == agrp->next) {
289                         BLI_insertlinkbefore(&act->curves, fcu, fcurve);
290                         
291                         if (agrp->channels.first == NULL)
292                                 agrp->channels.first= fcurve;
293                         agrp->channels.last= fcurve;
294                         
295                         done= 1;
296                         break;
297                 }
298                 
299                 /* if channel has group we're targeting, check whether it is the last one of these */
300                 else if (fcu->grp == agrp) {
301                         if ((fcu->next) && (fcu->next->grp != agrp)) {
302                                 BLI_insertlinkafter(&act->curves, fcu, fcurve);
303                                 agrp->channels.last= fcurve;
304                                 done= 1;
305                                 break;
306                         }
307                         else if (fcu->next == NULL) {
308                                 BLI_addtail(&act->curves, fcurve);
309                                 agrp->channels.last= fcurve;
310                                 done= 1;
311                                 break;
312                         }
313                 }
314                 
315                 /* if channel has group before target, check whether the next one is something after target */
316                 else if (fcu->grp == agrp->prev) {
317                         if (fcu->next) {
318                                 if ((fcu->next->grp != fcu->grp) && (fcu->next->grp != agrp)) {
319                                         BLI_insertlinkafter(&act->curves, fcu, fcurve);
320                                         
321                                         agrp->channels.first= fcurve;
322                                         agrp->channels.last= fcurve;
323                                         
324                                         done= 1;
325                                         break;
326                                 }
327                         }
328                         else {
329                                 BLI_insertlinkafter(&act->curves, fcu, fcurve);
330                                 
331                                 agrp->channels.first= fcurve;
332                                 agrp->channels.last= fcurve;
333                                 
334                                 done= 1;
335                                 break;
336                         }
337                 }
338         }
339         
340         /* only if added, set channel as belonging to this group */
341         if (done) {
342                 //printf("FCurve added to group \n");
343                 fcurve->grp= agrp;
344         }
345         else {
346                 printf("Error: FCurve '%s' couldn't be added to Group '%s' \n", fcurve->rna_path, agrp->name);
347                 BLI_addtail(&act->curves, fcurve);
348         }
349 }       
350
351 /* Remove the given channel from all groups */
352 void action_groups_remove_channel (bAction *act, FCurve *fcu)
353 {
354         /* sanity checks */
355         if (ELEM(NULL, act, fcu))       
356                 return;
357         
358         /* check if any group used this directly */
359         if (fcu->grp) {
360                 bActionGroup *agrp= fcu->grp;
361                 
362                 if (agrp->channels.first == agrp->channels.last) {
363                         if (agrp->channels.first == fcu) {
364                                 agrp->channels.first= NULL;
365                                 agrp->channels.last= NULL;
366                         }
367                 }
368                 else if (agrp->channels.first == fcu) {
369                         if ((fcu->next) && (fcu->next->grp==agrp))
370                                 agrp->channels.first= fcu->next;
371                         else
372                                 agrp->channels.first= NULL;
373                 }
374                 else if (agrp->channels.last == fcu) {
375                         if ((fcu->prev) && (fcu->prev->grp==agrp))
376                                 agrp->channels.last= fcu->prev;
377                         else
378                                 agrp->channels.last= NULL;
379                 }
380                 
381                 fcu->grp= NULL;
382         }
383         
384         /* now just remove from list */
385         BLI_remlink(&act->curves, fcu);
386 }
387
388 /* Find a group with the given name */
389 bActionGroup *action_groups_find_named (bAction *act, const char name[])
390 {
391         bActionGroup *grp;
392         
393         /* sanity checks */
394         if (ELEM3(NULL, act, act->groups.first, name) || (name[0] == 0))
395                 return NULL;
396                 
397         /* do string comparisons */
398         for (grp= act->groups.first; grp; grp= grp->next) {
399                 if (strcmp(grp->name, name) == 0)
400                         return grp;
401         }
402         
403         /* not found */
404         return NULL;
405 }
406
407 /* ************************ Pose channels *************** */
408
409 /* usually used within a loop, so we got a N^2 slowdown */
410 bPoseChannel *get_pose_channel(const bPose *pose, const char *name)
411 {
412         bPoseChannel *chan;
413         
414         if (ELEM(NULL, pose, name) || (name[0] == 0))
415                 return NULL;
416         
417         for (chan=pose->chanbase.first; chan; chan=chan->next) {
418                 if (chan->name[0] == name[0]) {
419                         if (!strcmp (chan->name, name))
420                                 return chan;
421                 }
422         }
423
424         return NULL;
425 }
426
427 /* Use with care, not on Armature poses but for temporal ones */
428 /* (currently used for action constraints and in rebuild_pose) */
429 bPoseChannel *verify_pose_channel(bPose* pose, const char* name)
430 {
431         bPoseChannel *chan;
432         
433         if (pose == NULL)
434                 return NULL;
435         
436         /* See if this channel exists */
437         for (chan=pose->chanbase.first; chan; chan=chan->next) {
438                 if (!strcmp (name, chan->name))
439                         return chan;
440         }
441         
442         /* If not, create it and add it */
443         chan = MEM_callocN(sizeof(bPoseChannel), "verifyPoseChannel");
444         
445         strncpy(chan->name, name, 31);
446         /* init vars to prevent math errors */
447         chan->quat[0] = 1.0f;
448         chan->size[0] = chan->size[1] = chan->size[2] = 1.0f;
449         
450         chan->limitmin[0]= chan->limitmin[1]= chan->limitmin[2]= -180.0f;
451         chan->limitmax[0]= chan->limitmax[1]= chan->limitmax[2]= 180.0f;
452         chan->stiffness[0]= chan->stiffness[1]= chan->stiffness[2]= 0.0f;
453         
454         Mat4One(chan->constinv);
455         
456         BLI_addtail(&pose->chanbase, chan);
457         
458         return chan;
459 }
460
461 /* Find the active posechannel for an object (we can't just use pose, as layer info is in armature) */
462 bPoseChannel *get_active_posechannel (Object *ob)
463 {
464         bArmature *arm= (ob) ? ob->data : NULL;
465         bPoseChannel *pchan;
466         
467         if ELEM3(NULL, ob, ob->pose, arm)
468                 return NULL;
469         
470         /* find active */
471         for (pchan= ob->pose->chanbase.first; pchan; pchan= pchan->next) {
472                 if ((pchan->bone) && (pchan->bone->flag & BONE_ACTIVE) && (pchan->bone->layer & arm->layer))
473                         return pchan;
474         }
475         
476         return NULL;
477 }
478
479
480 /* dst should be freed already, makes entire duplicate */
481 void copy_pose (bPose **dst, bPose *src, int copycon)
482 {
483         bPose *outPose;
484         bPoseChannel    *pchan;
485         ListBase listb;
486         
487         if (!src){
488                 *dst=NULL;
489                 return;
490         }
491         
492         if (*dst==src) {
493                 printf("copy_pose source and target are the same\n");
494                 *dst=NULL;
495                 return;
496         }
497         
498         outPose= MEM_callocN(sizeof(bPose), "pose");
499         
500         BLI_duplicatelist(&outPose->chanbase, &src->chanbase);
501         
502         if (copycon) {
503                 for (pchan=outPose->chanbase.first; pchan; pchan=pchan->next) {
504                         copy_constraints(&listb, &pchan->constraints);  // copy_constraints NULLs listb
505                         pchan->constraints= listb;
506                         pchan->path= NULL;
507                 }
508         }
509         
510         *dst=outPose;
511 }
512
513 void free_pose_channels(bPose *pose) 
514 {
515         bPoseChannel *pchan;
516         
517         if (pose->chanbase.first) {
518                 for (pchan = pose->chanbase.first; pchan; pchan=pchan->next){
519                         if (pchan->path)
520                                 MEM_freeN(pchan->path);
521                         free_constraints(&pchan->constraints);
522                 }
523                 BLI_freelistN(&pose->chanbase);
524         }
525 }
526
527 void free_pose(bPose *pose)
528 {
529         if (pose) {
530                 /* free pose-channels */
531                 free_pose_channels(pose);
532                 
533                 /* free pose-groups */
534                 if (pose->agroups.first)
535                         BLI_freelistN(&pose->agroups);
536                 
537                 /* free pose */
538                 MEM_freeN(pose);
539         }
540 }
541
542 void game_copy_pose(bPose **dst, bPose *src)
543 {
544         bPose *out;
545         bPoseChannel *pchan, *outpchan;
546         GHash *ghash;
547         
548         /* the game engine copies the current armature pose and then swaps
549          * the object pose pointer. this makes it possible to change poses
550          * without affecting the original blender data. */
551
552         if (!src) {
553                 *dst=NULL;
554                 return;
555         }
556         else if (*dst==src) {
557                 printf("copy_pose source and target are the same\n");
558                 *dst=NULL;
559                 return;
560         }
561         
562         out= MEM_dupallocN(src);
563         out->agroups.first= out->agroups.last= NULL;
564         BLI_duplicatelist(&out->chanbase, &src->chanbase);
565
566         /* remap pointers */
567         ghash= BLI_ghash_new(BLI_ghashutil_ptrhash, BLI_ghashutil_ptrcmp);
568
569         pchan= src->chanbase.first;
570         outpchan= out->chanbase.first;
571         for (; pchan; pchan=pchan->next, outpchan=outpchan->next)
572                 BLI_ghash_insert(ghash, pchan, outpchan);
573
574         for (pchan=out->chanbase.first; pchan; pchan=pchan->next) {
575                 pchan->parent= BLI_ghash_lookup(ghash, pchan->parent);
576                 pchan->child= BLI_ghash_lookup(ghash, pchan->child);
577                 pchan->path= NULL;
578         }
579
580         BLI_ghash_free(ghash, NULL, NULL);
581         
582         *dst=out;
583 }
584
585
586 /* Only allowed for Poses with identical channels */
587 void game_blend_poses(bPose *dst, bPose *src, float srcweight/*, short mode*/)
588 {
589         short mode= ACTSTRIPMODE_BLEND;
590         
591         bPoseChannel *dchan;
592         const bPoseChannel *schan;
593         bConstraint *dcon, *scon;
594         float dstweight;
595         int i;
596
597         switch (mode){
598         case ACTSTRIPMODE_BLEND:
599                 dstweight = 1.0F - srcweight;
600                 break;
601         case ACTSTRIPMODE_ADD:
602                 dstweight = 1.0F;
603                 break;
604         default :
605                 dstweight = 1.0F;
606         }
607         
608         schan= src->chanbase.first;
609         for (dchan = dst->chanbase.first; dchan; dchan=dchan->next, schan= schan->next){
610                 if (schan->flag & (POSE_ROT|POSE_LOC|POSE_SIZE)) {
611                         /* replaced quat->matrix->quat conversion with decent quaternion interpol (ton) */
612                         
613                         /* Do the transformation blend */
614                         if (schan->flag & POSE_ROT) {
615                                 /* quat interpolation done separate */
616                                 if (schan->rotmode == PCHAN_ROT_QUAT) {
617                                         float dquat[4], squat[4];
618                                         
619                                         QUATCOPY(dquat, dchan->quat);
620                                         QUATCOPY(squat, schan->quat);
621                                         if (mode==ACTSTRIPMODE_BLEND)
622                                                 QuatInterpol(dchan->quat, dquat, squat, srcweight);
623                                         else {
624                                                 QuatMulFac(squat, srcweight);
625                                                 QuatMul(dchan->quat, dquat, squat);
626                                         }
627                                         
628                                         NormalQuat(dchan->quat);
629                                 }
630                         }
631
632                         for (i=0; i<3; i++) {
633                                 /* blending for loc and scale are pretty self-explanatory... */
634                                 if (schan->flag & POSE_LOC)
635                                         dchan->loc[i] = (dchan->loc[i]*dstweight) + (schan->loc[i]*srcweight);
636                                 if (schan->flag & POSE_SIZE)
637                                         dchan->size[i] = 1.0f + ((dchan->size[i]-1.0f)*dstweight) + ((schan->size[i]-1.0f)*srcweight);
638                                 
639                                 /* euler-rotation interpolation done here instead... */
640                                 // FIXME: are these results decent?
641                                 if ((schan->flag & POSE_ROT) && (schan->rotmode))
642                                         dchan->eul[i] = (dchan->eul[i]*dstweight) + (schan->eul[i]*srcweight);
643                         }
644                         dchan->flag |= schan->flag;
645                 }
646                 for(dcon= dchan->constraints.first, scon= schan->constraints.first; dcon && scon; dcon= dcon->next, scon= scon->next) {
647                         /* no 'add' option for constraint blending */
648                         dcon->enforce= dcon->enforce*(1.0f-srcweight) + scon->enforce*srcweight;
649                 }
650         }
651         
652         /* this pose is now in src time */
653         dst->ctime= src->ctime;
654 }
655
656 void game_free_pose(bPose *pose)
657 {
658         if (pose) {
659                 /* we don't free constraints, those are owned by the original pose */
660                 if(pose->chanbase.first)
661                         BLI_freelistN(&pose->chanbase);
662                 
663                 MEM_freeN(pose);
664         }
665 }
666
667 static void copy_pose_channel_data(bPoseChannel *pchan, const bPoseChannel *chan)
668 {
669         bConstraint *pcon, *con;
670         
671         VECCOPY(pchan->loc, chan->loc);
672         VECCOPY(pchan->size, chan->size);
673         VECCOPY(pchan->eul, chan->eul);
674         QUATCOPY(pchan->quat, chan->quat);
675         pchan->rotmode= chan->rotmode;
676         Mat4CpyMat4(pchan->chan_mat, (float(*)[4])chan->chan_mat);
677         Mat4CpyMat4(pchan->pose_mat, (float(*)[4])chan->pose_mat);
678         pchan->flag= chan->flag;
679         
680         con= chan->constraints.first;
681         for(pcon= pchan->constraints.first; pcon; pcon= pcon->next, con= con->next) {
682                 pcon->enforce= con->enforce;
683                 pcon->headtail= con->headtail;
684         }
685 }
686
687 /* checks for IK constraint, and also for Follow-Path constraint.
688  * can do more constraints flags later 
689  */
690 /* pose should be entirely OK */
691 void update_pose_constraint_flags(bPose *pose)
692 {
693         bPoseChannel *pchan, *parchan;
694         bConstraint *con;
695         
696         /* clear */
697         for (pchan= pose->chanbase.first; pchan; pchan= pchan->next) {
698                 pchan->constflag= 0;
699         }
700         pose->flag &= ~POSE_CONSTRAINTS_TIMEDEPEND;
701         
702         /* detect */
703         for (pchan= pose->chanbase.first; pchan; pchan=pchan->next) {
704                 for (con= pchan->constraints.first; con; con= con->next) {
705                         if (con->type==CONSTRAINT_TYPE_KINEMATIC) {
706                                 bKinematicConstraint *data = (bKinematicConstraint*)con->data;
707                                 
708                                 pchan->constflag |= PCHAN_HAS_IK;
709                                 
710                                 if(data->tar==NULL || (data->tar->type==OB_ARMATURE && data->subtarget[0]==0))
711                                         pchan->constflag |= PCHAN_HAS_TARGET;
712                                 
713                                 /* negative rootbone = recalc rootbone index. used in do_versions */
714                                 if(data->rootbone<0) {
715                                         data->rootbone= 0;
716                                         
717                                         if(data->flag & CONSTRAINT_IK_TIP) parchan= pchan;
718                                         else parchan= pchan->parent;
719                                         
720                                         while(parchan) {
721                                                 data->rootbone++;
722                                                 if((parchan->bone->flag & BONE_CONNECTED)==0)
723                                                         break;
724                                                 parchan= parchan->parent;
725                                         }
726                                 }
727                         }
728                         else if (con->type == CONSTRAINT_TYPE_FOLLOWPATH) {
729                                 bFollowPathConstraint *data= (bFollowPathConstraint *)con->data;
730                                 
731                                 /* for drawing constraint colors when color set allows this */
732                                 pchan->constflag |= PCHAN_HAS_CONST;
733                                 
734                                 /* if we have a valid target, make sure that this will get updated on frame-change
735                                  * (needed for when there is no anim-data for this pose)
736                                  */
737                                 if ((data->tar) && (data->tar->type==OB_CURVE))
738                                         pose->flag |= POSE_CONSTRAINTS_TIMEDEPEND;
739                         }
740                         else 
741                                 pchan->constflag |= PCHAN_HAS_CONST;
742                 }
743         }
744 }
745
746 /* Clears all BONE_UNKEYED flags for every pose channel in every pose 
747  * This should only be called on frame changing, when it is acceptable to
748  * do this. Otherwise, these flags should not get cleared as poses may get lost.
749  */
750 void framechange_poses_clear_unkeyed(void)
751 {
752         Object *ob;
753         bPose *pose;
754         bPoseChannel *pchan;
755         
756         /* This needs to be done for each object that has a pose */
757         // TODO: proxies may/may not be correctly handled here... (this needs checking) 
758         for (ob= G.main->object.first; ob; ob= ob->id.next) {
759                 /* we only need to do this on objects with a pose */
760                 if ( (pose= ob->pose) ) {
761                         for (pchan= pose->chanbase.first; pchan; pchan= pchan->next) {
762                                 if (pchan->bone) 
763                                         pchan->bone->flag &= ~BONE_UNKEYED;
764                         }
765                 }
766         }
767 }
768
769 /* ************************** Bone Groups ************************** */
770
771 /* Adds a new bone-group */
772 void pose_add_group (Object *ob)
773 {
774         bPose *pose= (ob) ? ob->pose : NULL;
775         bActionGroup *grp;
776         
777         if (ELEM(NULL, ob, ob->pose))
778                 return;
779         
780         grp= MEM_callocN(sizeof(bActionGroup), "PoseGroup");
781         strcpy(grp->name, "Group");
782         BLI_addtail(&pose->agroups, grp);
783         BLI_uniquename(&pose->agroups, grp, "Group", '.', offsetof(bActionGroup, name), 32);
784         
785         pose->active_group= BLI_countlist(&pose->agroups);
786 }
787
788 /* Remove the active bone-group */
789 void pose_remove_group (Object *ob)
790 {
791         bPose *pose= (ob) ? ob->pose : NULL;
792         bActionGroup *grp = NULL;
793         bPoseChannel *pchan;
794         
795         /* sanity checks */
796         if (ELEM(NULL, ob, pose))
797                 return;
798         if (pose->active_group <= 0)
799                 return;
800         
801         /* get group to remove */
802         grp= BLI_findlink(&pose->agroups, pose->active_group-1);
803         if (grp) {
804                 /* adjust group references (the trouble of using indices!):
805                  *      - firstly, make sure nothing references it 
806                  *      - also, make sure that those after this item get corrected
807                  */
808                 for (pchan= pose->chanbase.first; pchan; pchan= pchan->next) {
809                         if (pchan->agrp_index == pose->active_group)
810                                 pchan->agrp_index= 0;
811                         else if (pchan->agrp_index > pose->active_group)
812                                 pchan->agrp_index--;
813                 }
814                 
815                 /* now, remove it from the pose */
816                 BLI_freelinkN(&pose->agroups, grp);
817                 pose->active_group= 0;
818         }
819 }
820
821 /* ************** time ****************** */
822
823 /* Check if the given action has any keyframes */
824 short action_has_motion(const bAction *act)
825 {
826         FCurve *fcu;
827         
828         /* return on the first F-Curve that has some keyframes/samples defined */
829         if (act) {
830                 for (fcu= act->curves.first; fcu; fcu= fcu->next) {
831                         if (fcu->totvert)
832                                 return 1;
833                 }
834         }
835         
836         /* nothing found */
837         return 0;
838 }
839
840 /* Calculate the extents of given action */
841 void calc_action_range(const bAction *act, float *start, float *end, int incl_hidden)
842 {
843         FCurve *fcu;
844         float min=999999999.0f, max=-999999999.0f;
845         short foundvert=0;
846
847         if (act) {
848                 for (fcu= act->curves.first; fcu; fcu= fcu->next) {
849                         if (fcu->totvert) {
850                                 float nmin, nmax;
851                                 
852                                 /* get extents for this curve */
853                                 calc_fcurve_range(fcu, &nmin, &nmax);
854                                 
855                                 /* compare to the running tally */
856                                 min= MIN2(min, nmin);
857                                 max= MAX2(max, nmax);
858                                 
859                                 foundvert= 1;
860                         }
861                 }
862         }       
863         
864         if (foundvert) {
865                 if(min==max) max+= 1.0f;
866                 *start= min;
867                 *end= max;
868         }
869         else {
870                 *start= 0.0f;
871                 *end= 1.0f;
872         }
873 }
874
875 /* Copy the data from the action-pose (src) into the pose */
876 /* both args are assumed to be valid */
877 /* exported to game engine */
878 /* Note! this assumes both poses are aligned, this isnt always true when dealing with user poses */
879 void extract_pose_from_pose(bPose *pose, const bPose *src)
880 {
881         const bPoseChannel *schan;
882         bPoseChannel *pchan= pose->chanbase.first;
883
884         if (pose==src) {
885                 printf("extract_pose_from_pose source and target are the same\n");
886                 return;
887         }
888
889         for (schan=src->chanbase.first; (schan && pchan); schan=schan->next, pchan= pchan->next) {
890                 copy_pose_channel_data(pchan, schan);
891         }
892 }
893
894 /* for do_all_pose_actions, clears the pose. Now also exported for proxy and tools */
895 void rest_pose(bPose *pose)
896 {
897         bPoseChannel *pchan;
898         int i;
899         
900         if (!pose)
901                 return;
902         
903         memset(pose->stride_offset, 0, sizeof(pose->stride_offset));
904         memset(pose->cyclic_offset, 0, sizeof(pose->cyclic_offset));
905         
906         for (pchan=pose->chanbase.first; pchan; pchan= pchan->next){
907                 for (i=0; i<3; i++) {
908                         pchan->loc[i]= 0.0f;
909                         pchan->quat[i+1]= 0.0f;
910                         pchan->eul[i]= 0.0f;
911                         pchan->size[i]= 1.0f;
912                 }
913                 pchan->quat[0]= 1.0f;
914                 
915                 pchan->flag &= ~(POSE_LOC|POSE_ROT|POSE_SIZE);
916         }
917 }
918
919 /* both poses should be in sync */
920 void copy_pose_result(bPose *to, bPose *from)
921 {
922         bPoseChannel *pchanto, *pchanfrom;
923         
924         if(to==NULL || from==NULL) {
925                 printf("pose result copy error\n"); // debug temp
926                 return;
927         }
928
929         if (to==from) {
930                 printf("copy_pose_result source and target are the same\n");
931                 return;
932         }
933
934
935         for(pchanfrom= from->chanbase.first; pchanfrom; pchanfrom= pchanfrom->next) {
936                 pchanto= get_pose_channel(to, pchanfrom->name);
937                 if(pchanto) {
938                         Mat4CpyMat4(pchanto->pose_mat, pchanfrom->pose_mat);
939                         Mat4CpyMat4(pchanto->chan_mat, pchanfrom->chan_mat);
940                         
941                         /* used for local constraints */
942                         VECCOPY(pchanto->loc, pchanfrom->loc);
943                         QUATCOPY(pchanto->quat, pchanfrom->quat);
944                         VECCOPY(pchanto->eul, pchanfrom->eul);
945                         VECCOPY(pchanto->size, pchanfrom->size);
946                         
947                         VECCOPY(pchanto->pose_head, pchanfrom->pose_head);
948                         VECCOPY(pchanto->pose_tail, pchanfrom->pose_tail);
949                         pchanto->flag= pchanfrom->flag;
950                 }
951         }
952 }
953
954 /* For the calculation of the effects of an Action at the given frame on an object 
955  * This is currently only used for the Action Constraint 
956  */
957 void what_does_obaction (Scene *scene, Object *ob, Object *workob, bPose *pose, bAction *act, char groupname[], float cframe)
958 {
959         bActionGroup *agrp= action_groups_find_named(act, groupname);
960         
961         /* clear workob */
962         clear_workob(workob);
963         
964         /* init workob */
965         Mat4CpyMat4(workob->obmat, ob->obmat);
966         Mat4CpyMat4(workob->parentinv, ob->parentinv);
967         Mat4CpyMat4(workob->constinv, ob->constinv);
968         workob->parent= ob->parent;
969         workob->track= ob->track;
970
971         workob->trackflag= ob->trackflag;
972         workob->upflag= ob->upflag;
973         
974         workob->partype= ob->partype;
975         workob->par1= ob->par1;
976         workob->par2= ob->par2;
977         workob->par3= ob->par3;
978
979         workob->constraints.first = ob->constraints.first;
980         workob->constraints.last = ob->constraints.last;
981         
982         workob->pose= pose;     /* need to set pose too, since this is used for both types of Action Constraint */
983
984         strcpy(workob->parsubstr, ob->parsubstr);
985         strcpy(workob->id.name, "OB<ConstrWorkOb>"); /* we don't use real object name, otherwise RNA screws with the real thing */
986         
987         /* if we're given a group to use, it's likely to be more efficient (though a bit more dangerous) */
988         if (agrp) {
989                 /* specifically evaluate this group only */
990                 PointerRNA id_ptr;
991                 
992                 /* get RNA-pointer for the workob's ID */
993                 RNA_id_pointer_create(&workob->id, &id_ptr);
994                 
995                 /* execute action for this group only */
996                 animsys_evaluate_action_group(&id_ptr, act, agrp, NULL, cframe);
997         }
998         else {
999                 AnimData adt;
1000                 
1001                 /* init animdata, and attach to workob */
1002                 memset(&adt, 0, sizeof(AnimData));
1003                 workob->adt= &adt;
1004                 
1005                 adt.recalc= ADT_RECALC_ANIM;
1006                 adt.action= act;
1007                 
1008                 /* execute effects of Action on to workob (or it's PoseChannels) */
1009                 BKE_animsys_evaluate_animdata(&workob->id, &adt, cframe, ADT_RECALC_ANIM);
1010         }
1011 }
1012
1013 /* ********** NLA with non-poses works with ipo channels ********** */
1014
1015 #if 0 // XXX OLD ANIMATION SYSTEM (TO BE REMOVED)
1016
1017 /* ************************ Blending with NLA *************** */
1018
1019 static void blend_pose_strides(bPose *dst, bPose *src, float srcweight, short mode)
1020 {
1021         float dstweight;
1022         
1023         switch (mode){
1024                 case ACTSTRIPMODE_BLEND:
1025                         dstweight = 1.0F - srcweight;
1026                         break;
1027                 case ACTSTRIPMODE_ADD:
1028                         dstweight = 1.0F;
1029                         break;
1030                 default :
1031                         dstweight = 1.0F;
1032         }
1033         
1034         VecLerpf(dst->stride_offset, dst->stride_offset, src->stride_offset, srcweight);
1035 }
1036
1037
1038 /* 
1039
1040 bone matching diagram, strips A and B
1041
1042                  .------------------------.
1043                  |         A              |
1044                  '------------------------'
1045                                  .          .             b2
1046                  .          .-------------v----------.
1047                  .              |         B   .          |
1048                  .          '------------------------'
1049                  .          .             .
1050                  .          .             .
1051 offset:          .    0     .    A-B      .  A-b2+B     
1052                  .          .             .
1053
1054 */
1055
1056
1057 static void blend_pose_offset_bone(bActionStrip *strip, bPose *dst, bPose *src, float srcweight, short mode)
1058 {
1059         /* matching offset bones */
1060         /* take dst offset, and put src on on that location */
1061         
1062         if(strip->offs_bone[0]==0)
1063                 return;
1064         
1065         /* are we also blending with matching bones? */
1066         if(strip->prev && strip->start>=strip->prev->start) {
1067                 bPoseChannel *dpchan= get_pose_channel(dst, strip->offs_bone);
1068                 if(dpchan) {
1069                         bPoseChannel *spchan= get_pose_channel(src, strip->offs_bone);
1070                         if(spchan) {
1071                                 float vec[3];
1072                                 
1073                                 /* dst->ctime has the internal strip->prev action time */
1074                                 /* map this time to nla time */
1075                                 
1076                                 float ctime= get_actionstrip_frame(strip, src->ctime, 1);
1077                                 
1078                                 if( ctime > strip->prev->end) {
1079                                         bActionChannel *achan;
1080                                         
1081                                         /* add src to dest, minus the position of src on strip->prev->end */
1082                                         
1083                                         ctime= get_actionstrip_frame(strip, strip->prev->end, 0);
1084                                         
1085                                         achan= get_action_channel(strip->act, strip->offs_bone);
1086                                         if(achan && achan->ipo) {
1087                                                 bPoseChannel pchan;
1088                                                 /* Evaluates and sets the internal ipo value */
1089                                                 calc_ipo(achan->ipo, ctime);
1090                                                 /* This call also sets the pchan flags */
1091                                                 execute_action_ipo(achan, &pchan);
1092                                                 
1093                                                 /* store offset that moves src to location of pchan */
1094                                                 VecSubf(vec, dpchan->loc, pchan.loc);
1095                                                 
1096                                                 Mat4Mul3Vecfl(dpchan->bone->arm_mat, vec);
1097                                         }
1098                                 }
1099                                 else {
1100                                         /* store offset that moves src to location of dst */
1101                                         
1102                                         VecSubf(vec, dpchan->loc, spchan->loc);
1103                                         Mat4Mul3Vecfl(dpchan->bone->arm_mat, vec);
1104                                 }
1105                                 
1106                                 /* if blending, we only add with factor scrweight */
1107                                 VecMulf(vec, srcweight);
1108                                 
1109                                 VecAddf(dst->cyclic_offset, dst->cyclic_offset, vec);
1110                         }
1111                 }
1112         }
1113         
1114         VecAddf(dst->cyclic_offset, dst->cyclic_offset, src->cyclic_offset);
1115 }
1116
1117 typedef struct NlaIpoChannel {
1118         struct NlaIpoChannel *next, *prev;
1119         float val;
1120         void *poin;
1121         int type;
1122 } NlaIpoChannel;
1123
1124 void extract_ipochannels_from_action(ListBase *lb, ID *id, bAction *act, const char *name, float ctime)
1125 {
1126         bActionChannel *achan= get_action_channel(act, name);
1127         IpoCurve *icu;
1128         NlaIpoChannel *nic;
1129         
1130         if(achan==NULL) return;
1131         
1132         if(achan->ipo) {
1133                 calc_ipo(achan->ipo, ctime);
1134                 
1135                 for(icu= achan->ipo->curve.first; icu; icu= icu->next) {
1136                         /* skip IPO_BITS, is for layers and cannot be blended */
1137                         if(icu->vartype != IPO_BITS) {
1138                                 nic= MEM_callocN(sizeof(NlaIpoChannel), "NlaIpoChannel");
1139                                 BLI_addtail(lb, nic);
1140                                 nic->val= icu->curval;
1141                                 nic->poin= get_ipo_poin(id, icu, &nic->type);
1142                         }
1143                 }
1144         }
1145         
1146         /* constraint channels only for objects */
1147         if(GS(id->name)==ID_OB) {
1148                 Object *ob= (Object *)id;
1149                 bConstraint *con;
1150                 bConstraintChannel *conchan;
1151                 
1152                 for (con=ob->constraints.first; con; con=con->next) {
1153                         conchan = get_constraint_channel(&achan->constraintChannels, con->name);
1154                         
1155                         if(conchan && conchan->ipo) {
1156                                 calc_ipo(conchan->ipo, ctime);
1157                                 
1158                                 icu= conchan->ipo->curve.first; // only one ipo now
1159                                 if(icu) {
1160                                         nic= MEM_callocN(sizeof(NlaIpoChannel), "NlaIpoChannel constr");
1161                                         BLI_addtail(lb, nic);
1162                                         nic->val= icu->curval;
1163                                         nic->poin= &con->enforce;
1164                                         nic->type= IPO_FLOAT;
1165                                 }
1166                         }
1167                 }
1168         }
1169 }
1170
1171 static NlaIpoChannel *find_nla_ipochannel(ListBase *lb, void *poin)
1172 {
1173         NlaIpoChannel *nic;
1174         
1175         if(poin) {
1176                 for(nic= lb->first; nic; nic= nic->next) {
1177                         if(nic->poin==poin)
1178                                 return nic;
1179                 }
1180         }
1181         return NULL;
1182 }
1183
1184
1185 static void blend_ipochannels(ListBase *dst, ListBase *src, float srcweight, int mode)
1186 {
1187         NlaIpoChannel *snic, *dnic, *next;
1188         float dstweight;
1189         
1190         switch (mode){
1191                 case ACTSTRIPMODE_BLEND:
1192                         dstweight = 1.0F - srcweight;
1193                         break;
1194                 case ACTSTRIPMODE_ADD:
1195                         dstweight = 1.0F;
1196                         break;
1197                 default :
1198                         dstweight = 1.0F;
1199         }
1200         
1201         for(snic= src->first; snic; snic= next) {
1202                 next= snic->next;
1203                 
1204                 dnic= find_nla_ipochannel(dst, snic->poin);
1205                 if(dnic==NULL) {
1206                         /* remove from src list, and insert in dest */
1207                         BLI_remlink(src, snic);
1208                         BLI_addtail(dst, snic);
1209                 }
1210                 else {
1211                         /* we do the blend */
1212                         dnic->val= dstweight*dnic->val + srcweight*snic->val;
1213                 }
1214         }
1215 }
1216
1217 int execute_ipochannels(ListBase *lb)
1218 {
1219         NlaIpoChannel *nic;
1220         int count = 0;
1221         
1222         for(nic= lb->first; nic; nic= nic->next) {
1223                 if(nic->poin) {
1224                         write_ipo_poin(nic->poin, nic->type, nic->val);
1225                         count++;
1226                 }
1227         }
1228         return count;
1229 }
1230
1231 /* nla timing */
1232
1233 /* this now only used for repeating cycles, to enable fields and blur. */
1234 /* the whole time control in blender needs serious thinking... */
1235 static float nla_time(Scene *scene, float cfra, float unit)
1236 {
1237         extern float bluroffs;  // bad construct, borrowed from object.c for now
1238         extern float fieldoffs;
1239         
1240         /* motion blur & fields */
1241         cfra+= unit*(bluroffs+fieldoffs);
1242         
1243         /* global time */
1244         cfra*= scene->r.framelen;       
1245         
1246         return cfra;
1247 }
1248
1249 /* added "sizecorr" here, to allow armatures to be scaled and still have striding.
1250    Only works for uniform scaling. In general I'd advise against scaling armatures ever though! (ton)
1251 */
1252 static float stridechannel_frame(Object *ob, float sizecorr, bActionStrip *strip, Path *path, float pathdist, float *stride_offset)
1253 {
1254         bAction *act= strip->act;
1255         const char *name= strip->stridechannel;
1256         bActionChannel *achan= get_action_channel(act, name);
1257         int stride_axis= strip->stride_axis;
1258
1259         if(achan && achan->ipo) {
1260                 IpoCurve *icu= NULL;
1261                 float minx=0.0f, maxx=0.0f, miny=0.0f, maxy=0.0f;
1262                 int foundvert= 0;
1263                 
1264                 if(stride_axis==0) stride_axis= AC_LOC_X;
1265                 else if(stride_axis==1) stride_axis= AC_LOC_Y;
1266                 else stride_axis= AC_LOC_Z;
1267                 
1268                 /* calculate the min/max */
1269                 for (icu=achan->ipo->curve.first; icu; icu=icu->next) {
1270                         if(icu->adrcode==stride_axis) {
1271                                 if(icu->totvert>1) {
1272                                         foundvert= 1;
1273                                         minx= icu->bezt[0].vec[1][0];
1274                                         maxx= icu->bezt[icu->totvert-1].vec[1][0];
1275                                         
1276                                         miny= icu->bezt[0].vec[1][1];
1277                                         maxy= icu->bezt[icu->totvert-1].vec[1][1];
1278                                 }
1279                                 break;
1280                         }
1281                 }
1282                 
1283                 if(foundvert && miny!=maxy) {
1284                         float stridelen= sizecorr*fabs(maxy-miny), striptime;
1285                         float actiondist, pdist, pdistNewNormalized, offs;
1286                         float vec1[4], vec2[4], dir[3];
1287                         
1288                         /* internal cycling, actoffs is in frames */
1289                         offs= stridelen*strip->actoffs/(maxx-minx);
1290                         
1291                         /* amount path moves object */
1292                         pdist = (float)fmod (pathdist+offs, stridelen);
1293                         striptime= pdist/stridelen;
1294                         
1295                         /* amount stride bone moves */
1296                         actiondist= sizecorr*eval_icu(icu, minx + striptime*(maxx-minx)) - miny;
1297                         
1298                         pdist = fabs(actiondist) - pdist;
1299                         pdistNewNormalized = (pathdist+pdist)/path->totdist;
1300                         
1301                         /* now we need to go pdist further (or less) on cu path */
1302                         where_on_path(ob, (pathdist)/path->totdist, vec1, dir); /* vec needs size 4 */
1303                         if (pdistNewNormalized <= 1) {
1304                                 // search for correction in positive path-direction
1305                                 where_on_path(ob, pdistNewNormalized, vec2, dir);       /* vec needs size 4 */
1306                                 VecSubf(stride_offset, vec2, vec1);
1307                         }
1308                         else {
1309                                 // we reached the end of the path, search backwards instead
1310                                 where_on_path(ob, (pathdist-pdist)/path->totdist, vec2, dir);   /* vec needs size 4 */
1311                                 VecSubf(stride_offset, vec1, vec2);
1312                         }
1313                         Mat4Mul3Vecfl(ob->obmat, stride_offset);
1314                         return striptime;
1315                 }
1316         }
1317         return 0.0f;
1318 }
1319
1320 static void cyclic_offs_bone(Object *ob, bPose *pose, bActionStrip *strip, float time)
1321 {
1322         /* only called when strip has cyclic, so >= 1.0f works... */
1323         if(time >= 1.0f) {
1324                 bActionChannel *achan= get_action_channel(strip->act, strip->offs_bone);
1325
1326                 if(achan && achan->ipo) {
1327                         IpoCurve *icu= NULL;
1328                         Bone *bone;
1329                         float min[3]={0.0f, 0.0f, 0.0f}, max[3]={0.0f, 0.0f, 0.0f};
1330                         int index=0, foundvert= 0;
1331                         
1332                         /* calculate the min/max */
1333                         for (icu=achan->ipo->curve.first; icu; icu=icu->next) {
1334                                 if(icu->totvert>1) {
1335                                         
1336                                         if(icu->adrcode==AC_LOC_X)
1337                                                 index= 0;
1338                                         else if(icu->adrcode==AC_LOC_Y)
1339                                                 index= 1;
1340                                         else if(icu->adrcode==AC_LOC_Z)
1341                                                 index= 2;
1342                                         else
1343                                                 continue;
1344                                 
1345                                         foundvert= 1;
1346                                         min[index]= icu->bezt[0].vec[1][1];
1347                                         max[index]= icu->bezt[icu->totvert-1].vec[1][1];
1348                                 }
1349                         }
1350                         if(foundvert) {
1351                                 /* bring it into armature space */
1352                                 VecSubf(min, max, min);
1353                                 bone= get_named_bone(ob->data, strip->offs_bone);       /* weak */
1354                                 if(bone) {
1355                                         Mat4Mul3Vecfl(bone->arm_mat, min);
1356                                         
1357                                         /* dominant motion, cyclic_offset was cleared in rest_pose */
1358                                         if (strip->flag & (ACTSTRIP_CYCLIC_USEX | ACTSTRIP_CYCLIC_USEY | ACTSTRIP_CYCLIC_USEZ)) {
1359                                                 if (strip->flag & ACTSTRIP_CYCLIC_USEX) pose->cyclic_offset[0]= time*min[0];
1360                                                 if (strip->flag & ACTSTRIP_CYCLIC_USEY) pose->cyclic_offset[1]= time*min[1];
1361                                                 if (strip->flag & ACTSTRIP_CYCLIC_USEZ) pose->cyclic_offset[2]= time*min[2];
1362                                         } else {
1363                                                 if( fabs(min[0]) >= fabs(min[1]) && fabs(min[0]) >= fabs(min[2]))
1364                                                         pose->cyclic_offset[0]= time*min[0];
1365                                                 else if( fabs(min[1]) >= fabs(min[0]) && fabs(min[1]) >= fabs(min[2]))
1366                                                         pose->cyclic_offset[1]= time*min[1];
1367                                                 else
1368                                                         pose->cyclic_offset[2]= time*min[2];
1369                                         }
1370                                 }
1371                         }
1372                 }
1373         }
1374 }
1375
1376 /* simple case for now; only the curve path with constraint value > 0.5 */
1377 /* blending we might do later... */
1378 static Object *get_parent_path(Object *ob)
1379 {
1380         bConstraint *con;
1381         
1382         if(ob->parent && ob->parent->type==OB_CURVE)
1383                 return ob->parent;
1384         
1385         for (con = ob->constraints.first; con; con=con->next) {
1386                 if(con->type==CONSTRAINT_TYPE_FOLLOWPATH) {
1387                         if(con->enforce>0.5f) {
1388                                 bFollowPathConstraint *data= con->data;
1389                                 return data->tar;
1390                         }
1391                 }
1392         }
1393         return NULL;
1394 }
1395
1396 /* ************** do the action ************ */
1397
1398 /* ----- nla, etc. --------- */
1399
1400 static void do_nla(Scene *scene, Object *ob, int blocktype)
1401 {
1402         bPose *tpose= NULL;
1403         Key *key= NULL;
1404         ListBase tchanbase={NULL, NULL}, chanbase={NULL, NULL};
1405         bActionStrip *strip, *striplast=NULL, *stripfirst=NULL;
1406         float striptime, frametime, length, actlength;
1407         float blendfac, stripframe;
1408         float scene_cfra= frame_to_float(scene, scene->r.cfra); 
1409         int     doit, dostride;
1410         
1411         if(blocktype==ID_AR) {
1412                 copy_pose(&tpose, ob->pose, 1);
1413                 rest_pose(ob->pose);            // potentially destroying current not-keyed pose
1414         }
1415         else {
1416                 key= ob_get_key(ob);
1417         }
1418         
1419         /* check on extend to left or right, when no strip is hit by 'cfra' */
1420         for (strip=ob->nlastrips.first; strip; strip=strip->next) {
1421                 /* escape loop on a hit */
1422                 if( scene_cfra >= strip->start && scene_cfra <= strip->end + 0.1f)      /* note 0.1 comes back below */
1423                         break;
1424                 if(scene_cfra < strip->start) {
1425                         if(stripfirst==NULL)
1426                                 stripfirst= strip;
1427                         else if(stripfirst->start > strip->start)
1428                                 stripfirst= strip;
1429                 }
1430                 else if(scene_cfra > strip->end) {
1431                         if(striplast==NULL)
1432                                 striplast= strip;
1433                         else if(striplast->end < strip->end)
1434                                 striplast= strip;
1435                 }
1436         }
1437         if(strip==NULL) {       /* extend */
1438                 if(striplast)
1439                         scene_cfra= striplast->end;
1440                 else if(stripfirst)
1441                         scene_cfra= stripfirst->start;
1442         }
1443         
1444         /* and now go over all strips */
1445         for (strip=ob->nlastrips.first; strip; strip=strip->next){
1446                 doit=dostride= 0;
1447                 
1448                 if (strip->act && !(strip->flag & ACTSTRIP_MUTE)) {     /* so theres an action */
1449                         
1450                         /* Determine if the current frame is within the strip's range */
1451                         length = strip->end-strip->start;
1452                         actlength = strip->actend-strip->actstart;
1453                         striptime = (scene_cfra-(strip->start)) / length;
1454                         stripframe = (scene_cfra-(strip->start)) ;
1455
1456                         if (striptime>=0.0){
1457                                 
1458                                 if(blocktype==ID_AR) 
1459                                         rest_pose(tpose);
1460                                 
1461                                 /* To handle repeat, we add 0.1 frame extra to make sure the last frame is included */
1462                                 if (striptime < 1.0f + 0.1f/length) {
1463                                         
1464                                         /* Handle path */
1465                                         if ((strip->flag & ACTSTRIP_USESTRIDE) && (blocktype==ID_AR) && (ob->ipoflag & OB_DISABLE_PATH)==0){
1466                                                 Object *parent= get_parent_path(ob);
1467                                                 
1468                                                 if (parent) {
1469                                                         Curve *cu = parent->data;
1470                                                         float ctime, pdist;
1471                                                         
1472                                                         if (cu->flag & CU_PATH){
1473                                                                 /* Ensure we have a valid path */
1474                                                                 if(cu->path==NULL || cu->path->data==NULL) makeDispListCurveTypes(scene, parent, 0);
1475                                                                 if(cu->path) {
1476                                                                         
1477                                                                         /* Find the position on the path */
1478                                                                         ctime= bsystem_time(scene, ob, scene_cfra, 0.0);
1479                                                                         
1480                                                                         if(calc_ipo_spec(cu->ipo, CU_SPEED, &ctime)==0) {
1481                                                                                 /* correct for actions not starting on zero */
1482                                                                                 ctime= (ctime - strip->actstart)/cu->pathlen;
1483                                                                                 CLAMP(ctime, 0.0, 1.0);
1484                                                                         }
1485                                                                         pdist = ctime*cu->path->totdist;
1486                                                                         
1487                                                                         if(tpose && strip->stridechannel[0]) {
1488                                                                                 striptime= stridechannel_frame(parent, ob->size[0], strip, cu->path, pdist, tpose->stride_offset);
1489                                                                         }                                                                       
1490                                                                         else {
1491                                                                                 if (strip->stridelen) {
1492                                                                                         striptime = pdist / strip->stridelen;
1493                                                                                         striptime = (float)fmod (striptime+strip->actoffs, 1.0);
1494                                                                                 }
1495                                                                                 else
1496                                                                                         striptime = 0;
1497                                                                         }
1498                                                                         
1499                                                                         frametime = (striptime * actlength) + strip->actstart;
1500                                                                         frametime= bsystem_time(scene, ob, frametime, 0.0);
1501                                                                         
1502                                                                         if(blocktype==ID_AR) {
1503                                                                                 extract_pose_from_action (tpose, strip->act, frametime);
1504                                                                         }
1505                                                                         else if(blocktype==ID_OB) {
1506                                                                                 extract_ipochannels_from_action(&tchanbase, &ob->id, strip->act, "Object", frametime);
1507                                                                                 if(key)
1508                                                                                         extract_ipochannels_from_action(&tchanbase, &key->id, strip->act, "Shape", frametime);
1509                                                                         }
1510                                                                         doit=dostride= 1;
1511                                                                 }
1512                                                         }
1513                                                 }
1514                                         }
1515                                         /* To handle repeat, we add 0.1 frame extra to make sure the last frame is included */
1516                                         else  {
1517                                                 
1518                                                 /* Mod to repeat */
1519                                                 if(strip->repeat!=1.0f) {
1520                                                         float cycle= striptime*strip->repeat;
1521                                                         
1522                                                         striptime = (float)fmod (cycle, 1.0f + 0.1f/length);
1523                                                         cycle-= striptime;
1524                                                         
1525                                                         if(blocktype==ID_AR)
1526                                                                 cyclic_offs_bone(ob, tpose, strip, cycle);
1527                                                 }
1528
1529                                                 frametime = (striptime * actlength) + strip->actstart;
1530                                                 frametime= nla_time(scene, frametime, (float)strip->repeat);
1531                                                         
1532                                                 if(blocktype==ID_AR) {
1533                                                         extract_pose_from_action (tpose, strip->act, frametime);
1534                                                 }
1535                                                 else if(blocktype==ID_OB) {
1536                                                         extract_ipochannels_from_action(&tchanbase, &ob->id, strip->act, "Object", frametime);
1537                                                         if(key)
1538                                                                 extract_ipochannels_from_action(&tchanbase, &key->id, strip->act, "Shape", frametime);
1539                                                 }
1540                                                 
1541                                                 doit=1;
1542                                         }
1543                                 }
1544                                 /* Handle extend */
1545                                 else {
1546                                         if (strip->flag & ACTSTRIP_HOLDLASTFRAME){
1547                                                 /* we want the strip to hold on the exact fraction of the repeat value */
1548                                                 
1549                                                 frametime = actlength * (strip->repeat-(int)strip->repeat);
1550                                                 if(frametime<=0.000001f) frametime= actlength;  /* rounding errors... */
1551                                                 frametime= bsystem_time(scene, ob, frametime+strip->actstart, 0.0);
1552                                                 
1553                                                 if(blocktype==ID_AR)
1554                                                         extract_pose_from_action (tpose, strip->act, frametime);
1555                                                 else if(blocktype==ID_OB) {
1556                                                         extract_ipochannels_from_action(&tchanbase, &ob->id, strip->act, "Object", frametime);
1557                                                         if(key)
1558                                                                 extract_ipochannels_from_action(&tchanbase, &key->id, strip->act, "Shape", frametime);
1559                                                 }
1560                                                 
1561                                                 /* handle cycle hold */
1562                                                 if(strip->repeat!=1.0f) {
1563                                                         if(blocktype==ID_AR)
1564                                                                 cyclic_offs_bone(ob, tpose, strip, strip->repeat-1.0f);
1565                                                 }
1566                                                 
1567                                                 doit=1;
1568                                         }
1569                                 }
1570                                 
1571                                 /* Handle blendin & blendout */
1572                                 if (doit){
1573                                         /* Handle blendin */
1574                                         
1575                                         if (strip->blendin>0.0 && stripframe<=strip->blendin && scene_cfra>=strip->start){
1576                                                 blendfac = stripframe/strip->blendin;
1577                                         }
1578                                         else if (strip->blendout>0.0 && stripframe>=(length-strip->blendout) && scene_cfra<=strip->end){
1579                                                 blendfac = (length-stripframe)/(strip->blendout);
1580                                         }
1581                                         else
1582                                                 blendfac = 1;
1583                                         
1584                                         if(blocktype==ID_AR) {/* Blend this pose with the accumulated pose */
1585                                                 /* offset bone, for matching cycles */
1586                                                 blend_pose_offset_bone (strip, ob->pose, tpose, blendfac, strip->mode);
1587                                                 
1588                                                 blend_poses (ob->pose, tpose, blendfac, strip->mode);
1589                                                 if(dostride)
1590                                                         blend_pose_strides (ob->pose, tpose, blendfac, strip->mode);
1591                                         }
1592                                         else {
1593                                                 blend_ipochannels(&chanbase, &tchanbase, blendfac, strip->mode);
1594                                                 BLI_freelistN(&tchanbase);
1595                                         }
1596                                 }
1597                         }                                       
1598                 }
1599         }
1600         
1601         if(blocktype==ID_OB) {
1602                 execute_ipochannels(&chanbase);
1603         }
1604         else if(blocktype==ID_AR) {
1605                 /* apply stride offset to object */
1606                 VecAddf(ob->obmat[3], ob->obmat[3], ob->pose->stride_offset);
1607         }
1608         
1609         /* free */
1610         if (tpose)
1611                 free_pose(tpose);
1612         if(chanbase.first)
1613                 BLI_freelistN(&chanbase);
1614 }
1615
1616 #endif // XXX OLD ANIMATION SYSTEM (TO BE REMOVED)