diff options
author | Cyrille Henry <nusmuk@users.sourceforge.net> | 2013-04-08 22:54:39 +0000 |
---|---|---|
committer | Cyrille Henry <nusmuk@users.sourceforge.net> | 2013-04-08 22:54:39 +0000 |
commit | 308f704b85f93da39a8a3559c11aab1dac9ae2f9 (patch) | |
tree | de6cf8ced157a77d342e4e38119a375df7ef5e6d /pmpd3d_stat.c | |
parent | f4e258a9d70ea1355e3e4bcce676e35c511ff3f3 (diff) |
bugfix and code refactorisation
svn path=/trunk/externals/pmpd/; revision=17071
Diffstat (limited to 'pmpd3d_stat.c')
-rw-r--r-- | pmpd3d_stat.c | 13 |
1 files changed, 6 insertions, 7 deletions
diff --git a/pmpd3d_stat.c b/pmpd3d_stat.c index 18ab6ff..b88a93b 100644 --- a/pmpd3d_stat.c +++ b/pmpd3d_stat.c @@ -443,7 +443,7 @@ void pmpd3d_linkPosMean(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) sommeX += (x->link[i].mass1->posX + x->link[i].mass2->posX)/2; sommeY += (x->link[i].mass1->posY + x->link[i].mass2->posY)/2; sommeZ += (x->link[i].mass1->posZ + x->link[i].mass2->posZ)/2; - somme += sqrt(sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posX + x->link[i].mass1->posX)/2)); + somme += sqrt(sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posY + x->link[i].mass2->posY)/2) + sqr((x->link[i].mass1->posZ + x->link[i].mass1->posZ)/2)); j++; } } @@ -455,7 +455,7 @@ void pmpd3d_linkPosMean(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) sommeX += (x->link[i].mass1->posX + x->link[i].mass2->posX)/2; sommeY += (x->link[i].mass1->posY + x->link[i].mass2->posY)/2; sommeZ += (x->link[i].mass1->posZ + x->link[i].mass2->posZ)/2; - somme += sqrt(sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posX + x->link[i].mass1->posX)/2)); + somme += sqrt(sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posY + x->link[i].mass2->posY)/2) + sqr((x->link[i].mass1->posZ + x->link[i].mass1->posZ)/2)); j++; } } @@ -668,7 +668,7 @@ void pmpd3d_linkPosStd(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) sommeX += (x->link[i].mass1->posX + x->link[i].mass2->posX)/2; sommeY += (x->link[i].mass1->posY + x->link[i].mass2->posY)/2; sommeZ += (x->link[i].mass1->posZ + x->link[i].mass2->posZ)/2; - somme += sqrt(sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posX + x->link[i].mass1->posX)/2)); + somme += sqrt(sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posY + x->link[i].mass2->posY)/2) + sqr((x->link[i].mass1->posZ + x->link[i].mass1->posZ)/2)); j++; } } @@ -680,7 +680,7 @@ void pmpd3d_linkPosStd(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) sommeX += (x->link[i].mass1->posX + x->link[i].mass2->posX)/2; sommeY += (x->link[i].mass1->posY + x->link[i].mass2->posY)/2; sommeZ += (x->link[i].mass1->posZ + x->link[i].mass2->posZ)/2; - somme += sqrt(sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posX + x->link[i].mass1->posX)/2)); + somme += sqrt(sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posY + x->link[i].mass2->posY)/2) + sqr((x->link[i].mass1->posZ + x->link[i].mass1->posZ)/2)); j++; } } @@ -702,7 +702,7 @@ void pmpd3d_linkPosStd(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) stdX += sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2 - sommeX); stdY += sqr((x->link[i].mass1->posY + x->link[i].mass2->posY)/2 - sommeY); stdZ += sqr((x->link[i].mass1->posZ + x->link[i].mass2->posZ)/2 - sommeZ); - std += sqr(sqrt(sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posX + x->link[i].mass1->posX)/2)) - somme); + std += sqr(sqrt(sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posY + x->link[i].mass2->posY)/2) + sqr((x->link[i].mass1->posZ + x->link[i].mass1->posZ)/2)) - somme); j+=1; } } @@ -714,7 +714,7 @@ void pmpd3d_linkPosStd(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) stdX += sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2 - sommeX); stdY += sqr((x->link[i].mass1->posY + x->link[i].mass2->posY)/2 - sommeY); stdZ += sqr((x->link[i].mass1->posZ + x->link[i].mass2->posZ)/2 - sommeZ); - std += sqr(sqrt(sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posX + x->link[i].mass1->posX)/2)) - somme); + std += sqr(sqrt(sqr((x->link[i].mass1->posX + x->link[i].mass2->posX)/2) + sqr((x->link[i].mass1->posY + x->link[i].mass2->posY)/2) + sqr((x->link[i].mass1->posZ + x->link[i].mass1->posZ)/2)) - somme); j+=1; } } @@ -1016,7 +1016,6 @@ void pmpd3d_linkLengthSpeedStd(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) std = sqrt(std /j); } - SETFLOAT(&(std_out[0]),stdX); SETFLOAT(&(std_out[1]),stdY); SETFLOAT(&(std_out[2]),stdZ); |