code refactorisation, and few new fonctionalities copied from pmpd3d

svn path=/trunk/externals/pmpd/; revision=17073

1 files changed, 258 insertions, 3848 deletions

diff --git a/pmpd2d.c b/pmpd2d.c
index 4b5bef3..0c60128 100644
--- a/pmpd2d.c
+++ b/pmpd2d.c
@@ -30,3768 +30,17 @@
 #include "stdio.h"

 #include "math.h"

 

-#define nb_max_link   100000 

-#define nb_max_mass   100000

-

-#define max(a,b) ( ((a) > (b)) ? (a) : (b) ) 

-#define min(a,b) ( ((a) < (b)) ? (a) : (b) ) 

-

-t_float sign_ch(t_float v)

-{

-    return v > 0 ? 1 : -1;

-}

-

-t_float sqr(t_float x)

-{

-    return x*x ;

-}

-

-t_float pow_ch(t_float x, t_float y)

-{

-    return x > 0 ? pow(x,y) : -pow(-x,y);

-}

-    

-static t_class *pmpd2d_class;

-

-typedef struct _mass {

-    t_symbol *Id;

-    int mobile;

-    t_float invM;

-    t_float speedX;

-    t_float speedY;

-    t_float posX;

-    t_float posY;

-    t_float forceX;

-    t_float forceY;

-    t_float D2;

-    t_float D2offset;

-    int num;

-} foo;

-

-typedef struct _link {

-    t_symbol *Id;

-    int lType;

-    struct _mass *mass1;

-    struct _mass *mass2;

-    t_float K;

-    t_float D;

-    t_float L;

-    t_float Pow;

-    t_float Lmin;

-    t_float Lmax;

-    t_float distance;

-    t_float VX; // vecteur portant la liaison, si c'est le cas

-    t_float VY;

-    t_symbol *arrayK;

-    t_symbol *arrayD;

-    t_float K_L; // longeur du tabeau K

-    t_float D_L; // longeur du tabeau D

-} foo1 ;

-

-typedef struct _pmpd2d {

-     t_object  x_obj;

-    struct _link link[nb_max_link];

-    struct _mass mass[nb_max_mass];

-    t_outlet *main_outlet;

-    t_outlet *info_outlet;

-    int nb_link;

-    int nb_mass;

-    t_float minX, maxX, minY, maxY;

-    t_int grab; // si on grab une mass ou pas

-    t_int grab_nb; // la masse grabé

-} t_pmpd2d;

-

-t_float tabread2(t_pmpd2d *x, t_float pos, t_symbol *array)

-{

-    t_garray *a;

-    int npoints;

-    t_word *vec;

-    t_float posx;

-    

-    if (!(a = (t_garray *)pd_findbyclass(array, garray_class)))

-        pd_error(x, "%s: no such array", array->s_name);

-    else if (!garray_getfloatwords(a, &npoints, &vec))

-        pd_error(x, "%s: bad template for tabLink", array->s_name);

-    else

-    {

-        posx = fabs(pos)*npoints;

-        int n=posx;

-        if (n >= npoints - 1) 

-            return (sign_ch(pos)*vec[npoints-1].w_float);

-        float fract = posx-n;

-        return (sign_ch(pos) * ( fract*vec[n+1].w_float+(1-fract)*vec[n].w_float));

-    }

-    return( pos); // si il y a un pb sur le tableau, on renvoie l'identité

-}

-

-void pmpd2d_reset(t_pmpd2d *x)

-{

-    x->nb_link = 0;

-    x->nb_mass = 0;

-    x->minX = -1000000;

-    x->maxX = 1000000;

-    x->minY = -1000000;

-    x->maxY = 1000000;

-    x->grab = 0;

-}

-

-void pmpd2d_infosL(t_pmpd2d *x)

-{

-    int i;

-    post("list of mass");

-    post("number, Id, mobile, mass, damping, positionX Y, speedX Y, forcesX Y");

-    for(i=0; i < x->nb_mass; i++)

-    {

-        post("masse %i: %s, %d, %f, %f, %f, %f, %f, %f, %f, %f",i, x->mass[i].Id->s_name, \

-            x->mass[i].mobile, 1/x->mass[i].invM, x->mass[i].D2, x->mass[i].posX, x->mass[i].posY, \

-            x->mass[i].speedX, x->mass[i].speedY, x->mass[i].forceX, x->mass[i].forceY );

-    }

-

-    post("list of link");

-    post("number, Id, mass1, mass2, K, D, Pow, L, Lmin, Lmax");

-    for(i=0; i < x->nb_link; i++)

-    {

-        switch(x->link[i].lType)

-        {

-        case 0 :

-            post("link %i: %s, %i, %i, %f, %f, %f, %f, %f, %f", i, x->link[i].Id->s_name, \

-                x->link[i].mass1->num, x->link[i].mass2->num, x->link[i].K, x->link[i].D, \

-                x->link[i].Pow, x->link[i].L, x->link[i].Lmin, x->link[i].Lmax);

-            break;

-        case 1 :

-            post("tLink %i: %s, %i, %i, %f, %f, %f, %f, %f, %f, %f, %f", i, x->link[i].Id->s_name, \

-                x->link[i].mass1->num, x->link[i].mass2->num, x->link[i].K, x->link[i].D, \

-                x->link[i].Pow, x->link[i].L, x->link[i].Lmin, x->link[i].Lmax, x->link[i].VX, x->link[i].VY);

-            break;

-        case 2 :

-            post("tabLink %i: %s, %i, %i, %f, %f, %s, %f, %s, %f", i, x->link[i].Id->s_name, \

-                x->link[i].mass1->num, x->link[i].mass2->num, x->link[i].K, x->link[i].D, \

-                x->link[i].arrayK->s_name, x->link[i].K_L, x->link[i].arrayD->s_name, x->link[i].D_L);

-            break;

-        }

-    }

-}

-

-void pmpd2d_bang(t_pmpd2d *x)

-{

-// this part is doing all the PM

-    t_float F, L, Lx,Ly, Fx, Fy, tmpX, tmpY,speed;

-    t_int i;

-    // post("bang");

-

-    for (i=0; i<x->nb_mass; i++)

-    // compute new masses position

-        if (x->mass[i].mobile > 0) // only if mobile

-        {

-            x->mass[i].speedX += x->mass[i].forceX * x->mass[i].invM;

-            x->mass[i].speedY += x->mass[i].forceY * x->mass[i].invM;

-            // x->mass[i].forceX = 0;

-            // x->mass[i].forceY = 0;        

-            x->mass[i].posX += x->mass[i].speedX ;

-            x->mass[i].posY += x->mass[i].speedY ;

-            if ( (x->mass[i].posX < x->minX) || (x->mass[i].posX > x->maxX) || (x->mass[i].posY < x->minY) || (x->mass[i].posY > x->maxY) ) 

-            {

-                tmpX = min(x->maxX,max(x->minX,x->mass[i].posX));

-                tmpY = min(x->maxY,max(x->minY,x->mass[i].posY));

-                x->mass[i].speedX -= x->mass[i].posX - tmpX;

-                x->mass[i].speedY -= x->mass[i].posY - tmpY;

-                x->mass[i].posX = tmpX;

-                x->mass[i].posY = tmpY;

-            }

-            x->mass[i].forceX = -x->mass[i].D2 * x->mass[i].speedX;

-            x->mass[i].forceY = -x->mass[i].D2 * x->mass[i].speedY;

-            speed = sqrt(x->mass[i].speedX * x->mass[i].speedX + x->mass[i].speedY * x->mass[i].speedY);

-            if (speed != 0) {

-                x->mass[i].forceX += x->mass[i].D2offset * (x->mass[i].speedX/speed);

-                x->mass[i].forceY += x->mass[i].D2offset * (x->mass[i].speedY/speed);

-            }

-        }

-

-    for (i=0; i<x->nb_link; i++)

-    // compute link forces

-    {

-        Lx = x->link[i].mass1->posX - x->link[i].mass2->posX;

-        Ly = x->link[i].mass1->posY - x->link[i].mass2->posY;

-        L = sqrt( sqr(Lx) + sqr(Ly) );

-        

-        if ( (L >= x->link[i].Lmin) && (L < x->link[i].Lmax)  && (L != 0))

-        {

-            if (x->link[i].lType == 2)

-            { // K et D viennent d'une table

-                F  = x->link[i].D * tabread2(x, (L - x->link[i].distance) / x->link[i].D_L, x->link[i].arrayD);

-                F += x->link[i].K * tabread2(x, L / x->link[i].K_L, x->link[i].arrayK);

-            }

-            else

-            {            

-                F  = x->link[i].D * (L - x->link[i].distance) ;

-                F += x->link[i].K *  pow_ch( L - x->link[i].L, x->link[i].Pow);

-            }

-

-            Fx = F * Lx/L;

-            Fy = F * Ly/L;    

-                

-            if (x->link[i].lType == 1)

-            { // on projette selon 1 axe

-                // F = Fx*x->link[i].VX + Fy*x->link[i].VY; // produit scalaire de la force sur le vecteur qui la porte

-                Fx = Fx*x->link[i].VX; // V est unitaire, dc on projete sans pb

-                Fy = Fy*x->link[i].VY;                

-            }

-            

-            x->link[i].mass1->forceX -= Fx;

-            x->link[i].mass1->forceY -= Fy;

-            x->link[i].mass2->forceX += Fx;

-            x->link[i].mass2->forceY += Fy;

-        }

-        x->link[i].distance=L;

-    }

-}

-

-void pmpd2d_mass(t_pmpd2d *x, t_symbol *Id, t_float mobile, t_float M, t_float posX, t_float posY )

-{ // add a mass : Id, invM speedX posX

-

-    if (M<=0) M=1;

-    x->mass[x->nb_mass].Id = Id;

-    x->mass[x->nb_mass].mobile = (int)mobile;

-    x->mass[x->nb_mass].invM = 1/M;

-    x->mass[x->nb_mass].speedX = 0;

-    x->mass[x->nb_mass].speedY = 0;

-    x->mass[x->nb_mass].posX = posX;

-    x->mass[x->nb_mass].posY = posY;

-    x->mass[x->nb_mass].forceX = 0;

-    x->mass[x->nb_mass].forceY = 0;

-    x->mass[x->nb_mass].num = x->nb_mass;

-    x->mass[x->nb_mass].D2 = 0;

-    x->mass[x->nb_mass].D2offset = 0;

-

-    x->nb_mass++ ;

-    x->nb_mass = min ( nb_max_mass -1, x->nb_mass );

-}

-

-void pmpd2d_create_link(t_pmpd2d *x, t_symbol *Id, int mass1, int mass2, t_float K, t_float D, t_float Pow, t_float Lmin, t_float Lmax, t_int type)

-{ // create a link based on mass number

-

-    if ((x->nb_mass>1) && (mass1 != mass2) && (mass1 >= 0) && (mass2 >= 0) && (mass1 < x->nb_mass) && (mass2 < x->nb_mass) )

-    {

-        x->link[x->nb_link].lType = type;

-        x->link[x->nb_link].Id = Id;

-        x->link[x->nb_link].mass1 = &x->mass[mass1]; 

-        x->link[x->nb_link].mass2 = &x->mass[mass2];

-        x->link[x->nb_link].K = K;

-        x->link[x->nb_link].D = D;

-        x->link[x->nb_link].L = sqrt(sqr(x->mass[mass1].posX - x->mass[mass2].posX) + sqr(x->mass[mass1].posY - x->mass[mass2].posY));

-        x->link[x->nb_link].Pow = Pow;

-        x->link[x->nb_link].Lmin = Lmin;

-        x->link[x->nb_link].Lmax = Lmax;

-        x->link[x->nb_link].distance = x->link[x->nb_link].L ;

-

-        x->nb_link++ ;

-        x->nb_link = min ( nb_max_link -1, x->nb_link );

-    }

-}

-

-void pmpd2d_link(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{ // add a link : Id, *mass1, *mass2, K, D, Pow, Lmin, Lmax;

-

-    int i, j;

-

-    t_symbol *Id = atom_getsymbolarg(0,argc,argv);

-    int mass1 = atom_getfloatarg(1, argc, argv);

-    int mass2 = atom_getfloatarg(2, argc, argv);

-    t_float K = atom_getfloatarg(3, argc, argv);

-    t_float D = atom_getfloatarg(4, argc, argv);

-    t_float Pow = 1; 

-    if (argc > 5) Pow = atom_getfloatarg(5, argc, argv);

-    t_float Lmin = -1000000;

-    if (argc > 6) Lmin = atom_getfloatarg(6, argc, argv);

-    t_float Lmax =  1000000;

-    if (argc > 7) Lmax = atom_getfloatarg(7, argc, argv);

-//    post("%d,%d, %f,%f", mass1, mass2, K, D);

-

-    if ( ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        pmpd2d_create_link(x, Id, mass1, mass2, K, D, Pow, Lmin, Lmax, 0);

-    }

-    else

-    if ( ( argv[1].a_type == A_SYMBOL ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(1,argc,argv) == x->mass[i].Id)

-            {

-                pmpd2d_create_link(x, Id, i, mass2, K, D, Pow, Lmin, Lmax, 0);

-            }

-        }

-    }

-    else

-    if ( ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_SYMBOL ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(2,argc,argv) == x->mass[i].Id)

-            {

-                pmpd2d_create_link(x, Id, mass1, i, K, D, Pow, Lmin, Lmax, 0);

-            }

-        }

-    }

-    else

-    if ( ( argv[1].a_type == A_SYMBOL ) && ( argv[2].a_type == A_SYMBOL ) )

-    {

-        for (i=0; i < x->nb_mass; i++)

-        {

-            for (j=0; j < x->nb_mass; j++)

-            {

-                if ( (atom_getsymbolarg(1,argc,argv) == x->mass[i].Id)&&(atom_getsymbolarg(2,argc,argv) == x->mass[j].Id))

-                {

-                    if (!( (x->mass[i].Id == x->mass[j].Id) && (i>j) )) 

-                        // si lien entre 2 serie de masses identique entres elle, alors on ne creer qu'un lien sur 2, pour evider les redondances

-                        pmpd2d_create_link(x, Id, i, j, K, D, Pow, Lmin, Lmax, 0);

-                }

-            }   

-        }

-    }

-}

-

-void pmpd2d_tLink(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{ // add a link : Id, *mass1, *mass2, K, D, Vx, Vy, Pow, Lmin, Lmax;

-

-    int i, j;

-    

-    t_symbol *Id = atom_getsymbolarg(0,argc,argv);

-    int mass1 = atom_getfloatarg(1, argc, argv);

-    int mass2 = atom_getfloatarg(2, argc, argv);

-    t_float K = atom_getfloatarg(3, argc, argv);

-    t_float D = atom_getfloatarg(4, argc, argv);

-    t_float vecteurX = atom_getfloatarg(5, argc, argv);

-    t_float vecteurY = atom_getfloatarg(6, argc, argv);

-    t_float vecteur = sqrt( sqr(vecteurX) + sqr(vecteurY) );

-    vecteurX /= vecteur;

-    vecteurY /= vecteur;

-    t_float Pow = 1; 

-    if (argc > 7) Pow = atom_getfloatarg(7, argc, argv);

-    t_float Lmin = 0;

-    if (argc > 8) Lmin = atom_getfloatarg(8, argc, argv);

-    t_float Lmax =  1000000;

-    if (argc > 9) Lmax = atom_getfloatarg(9, argc, argv);

-

-    if ( ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        pmpd2d_create_link(x, Id, mass1, mass2, K, D, Pow, Lmin, Lmax, 1);

-        x->link[x->nb_link-1].VX = vecteurX;

-        x->link[x->nb_link-1].VY = vecteurY;

-    }

-    else

-    if ( ( argv[1].a_type == A_SYMBOL ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(1,argc,argv) == x->mass[i].Id)

-            {

-                pmpd2d_create_link(x, Id, i, mass2, K, D, Pow, Lmin, Lmax, 1);

-                x->link[x->nb_link-1].VX = vecteurX;

-                x->link[x->nb_link-1].VY = vecteurY;

-            }

-        }

-    }

-    else

-    if ( ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_SYMBOL ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(2,argc,argv) == x->mass[i].Id)

-            {

-                pmpd2d_create_link(x, Id, mass1, i, K, D, Pow, Lmin, Lmax, 1);

-                x->link[x->nb_link-1].VX = vecteurX;

-                x->link[x->nb_link-1].VY = vecteurY;

-            }

-        }

-    }

-    else

-    if ( ( argv[1].a_type == A_SYMBOL ) && ( argv[2].a_type == A_SYMBOL ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            for (j=0; j< x->nb_mass; j++)

-            {

-                if ( (atom_getsymbolarg(1,argc,argv) == x->mass[i].Id)&(atom_getsymbolarg(2,argc,argv) == x->mass[j].Id))

-                {

-                    if (!( (x->mass[i].Id == x->mass[j].Id) && (i>j) )) 

-                        // si lien entre 2 serie de masses identique entres elle, alors on ne creer qu'un lien sur 2, pour evider les redondances

-                    {

-                        pmpd2d_create_link(x, Id, i, j, K, D, Pow, Lmin, Lmax, 1);

-                        x->link[x->nb_link-1].VX = vecteurX;

-                        x->link[x->nb_link-1].VY = vecteurY;

-                    }

-

-                }

-            }   

-        }

-    }

-}

-

-void pmpd2d_tabLink(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j;

-    

-    t_symbol *Id = atom_getsymbolarg(0,argc,argv);

-    int mass1 = atom_getfloatarg(1, argc, argv);

-    int mass2 = atom_getfloatarg(2, argc, argv);

-    t_symbol *arrayK = atom_getsymbolarg(3,argc,argv);

-    t_float Kl = atom_getfloatarg(4, argc, argv);

-       if (Kl <= 0) Kl = 1;

-    t_symbol *arrayD = atom_getsymbolarg(5,argc,argv);    

-    t_float Dl = atom_getfloatarg(6, argc, argv);

-    if (Dl <= 0) Dl = 1;

-

-    if ( ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        pmpd2d_create_link(x, Id, mass1, mass2, 1, 1, 1, 0, 1000000, 2);

-        x->link[x->nb_link-1].arrayK = arrayK;

-        x->link[x->nb_link-1].arrayD = arrayD;

-        x->link[x->nb_link-1].K_L = Kl;

-        x->link[x->nb_link-1].D_L = Dl;    

-    }

-    else

-    if ( ( argv[1].a_type == A_SYMBOL ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(1,argc,argv) == x->mass[i].Id)

-            {

-                pmpd2d_create_link(x, Id, i, mass2, 1, 1, 1, 0, 1000000, 2);

-                x->link[x->nb_link-1].arrayK = arrayK;

-                x->link[x->nb_link-1].arrayD = arrayD;

-                x->link[x->nb_link-1].K_L = Kl;

-                x->link[x->nb_link-1].D_L = Dl;    

-            }

-        }

-    }

-    else

-    if ( ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_SYMBOL ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(2,argc,argv) == x->mass[i].Id)

-            {

-                pmpd2d_create_link(x, Id, mass1, i, 1, 1, 1, 0, 1000000, 2);

-                x->link[x->nb_link-1].arrayK = arrayK;

-                x->link[x->nb_link-1].arrayD = arrayD;

-                x->link[x->nb_link-1].K_L = Kl;

-                x->link[x->nb_link-1].D_L = Dl;    

-            }

-        }

-    }

-    else

-    if ( ( argv[1].a_type == A_SYMBOL ) && ( argv[2].a_type == A_SYMBOL ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            for (j=0; j< x->nb_mass; j++)

-            {

-                if ( (atom_getsymbolarg(1,argc,argv) == x->mass[i].Id)&(atom_getsymbolarg(2,argc,argv) == x->mass[j].Id))

-                {

-                    if (!( (x->mass[i].Id == x->mass[j].Id) && (i>j) )) 

-                        // si lien entre 2 serie de masses identique entres elle, alors on ne creer qu'un lien sur 2, pour evider les redondances

-                    {

-                        pmpd2d_create_link(x, Id, i, j, 1, 1, 1, 0, 1000000, 2);

-                        x->link[x->nb_link-1].arrayK = arrayK;

-                        x->link[x->nb_link-1].arrayD = arrayD;

-                        x->link[x->nb_link-1].K_L = Kl;

-                        x->link[x->nb_link-1].D_L = Dl;    

-                    }

-                }

-            }   

-        }

-    }

-

-}

-

-void pmpd2d_setK(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_link-1, tmp));

-        x->link[tmp].K = atom_getfloatarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_link; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->link[i].Id)

-            {

-                x->link[i].K = atom_getfloatarg(1, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_setD(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_link-1, tmp));

-        x->link[tmp].D = atom_getfloatarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_link; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->link[i].Id)

-            {

-                x->link[i].D = atom_getfloatarg(1, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_setL(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_link-1, tmp));

-        x->link[tmp].L = atom_getfloatarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_link; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->link[i].Id)

-            {

-                x->link[i].L = atom_getfloatarg(1, argc, argv);

-            }

-        }

-    }

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argc == 1 ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_link-1, tmp));

-        x->link[tmp].L = x->link[tmp].mass2->posX - x->link[tmp].mass1->posX;

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argc == 1 ) )

-    {

-        for (i=0; i< x->nb_link; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->link[i].Id)

-            {

-                x->link[i].L = x->link[i].mass2->posX - x->link[i].mass1->posX;

-            }

-        }

-    }

-}

-

-void pmpd2d_setLKTab(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-    t_float K_l = atom_getfloatarg(1, argc, argv);

-    if (K_l <=  0) K_l = 1;

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_link-1, tmp));

-        x->link[tmp].K_L = K_l;

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_link; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->link[i].Id)

-            {

-                x->link[i].K_L = K_l;

-            }

-        }

-    }

-}

-

-void pmpd2d_setLDTab(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-    t_float D_l = atom_getfloatarg(1, argc, argv);

-    if (D_l <=  0) D_l = 1;

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_link-1, tmp));

-        x->link[tmp].D_L = D_l;

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_link; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->link[i].Id)

-            {

-                x->link[i].D_L = D_l;

-            }

-        }

-    }

-}

-

-void pmpd2d_setLinkId(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_SYMBOL ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_link-1, tmp));

-        x->link[tmp].Id = atom_getsymbolarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_SYMBOL ) )

-    {

-        for (i=0; i< x->nb_link; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->link[i].Id)

-            {

-                x->link[i].Id = atom_getsymbolarg(1, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_setMassId(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_SYMBOL ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].Id = atom_getsymbolarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_SYMBOL ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].Id = atom_getsymbolarg(1, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_forceX(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-// add a force to a specific mass

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].forceX += atom_getfloatarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].forceX += atom_getfloatarg(1, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_forceY(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-// add a force to a specific mass

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].forceY += atom_getfloatarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].forceY += atom_getfloatarg(1, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_force(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-// add a force to a specific mass

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].forceX += atom_getfloatarg(1, argc, argv);

-        x->mass[tmp].forceY += atom_getfloatarg(2, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].forceX += atom_getfloatarg(1, argc, argv);

-                x->mass[i].forceY += atom_getfloatarg(2, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_posX(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-// displace a mass to a certain position

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].posX = atom_getfloatarg(1, argc, argv);

-           x->mass[tmp].speedX = 0; 

-        x->mass[tmp].forceX = 0; 

-        

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].posX = atom_getfloatarg(1, argc, argv);

-                x->mass[i].speedX = 0; 

-                x->mass[i].forceX = 0;

-

-            }

-        }

-    }

-}

-

-void pmpd2d_posY(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-// displace a mass to a certain position

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].posY = atom_getfloatarg(1, argc, argv);

-           x->mass[tmp].speedY = 0; 

-        x->mass[tmp].forceY = 0; 

-        

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].posY = atom_getfloatarg(1, argc, argv);

-                x->mass[i].speedY = 0; 

-                x->mass[i].forceY = 0;

-            }

-        }

-    }

-}

-

-void pmpd2d_pos(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-// displace a mass to a certain position

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].posX = atom_getfloatarg(1, argc, argv);

-           x->mass[tmp].speedX = 0; 

-        x->mass[tmp].forceX = 0; 

-           x->mass[tmp].posY = atom_getfloatarg(2, argc, argv);

-           x->mass[tmp].speedY = 0; 

-        x->mass[tmp].forceY = 0; 

-        

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].posX = atom_getfloatarg(1, argc, argv);

-                x->mass[i].speedX = 0; 

-                x->mass[i].forceX = 0;

-                x->mass[i].posY = atom_getfloatarg(2, argc, argv);

-                x->mass[i].speedY = 0; 

-                x->mass[i].forceY = 0;

-

-            }

-        }

-    }

-}

-

-void pmpd2d_minX(t_pmpd2d *x, t_float min)

-{

-    x->minX = min;

-}

-

-void pmpd2d_maxX(t_pmpd2d *x, t_float max)

-{

-    x->maxX = max;

-}

-

-void pmpd2d_minY(t_pmpd2d *x, t_float min)

-{

-    x->minY = min;

-}

-

-void pmpd2d_maxY(t_pmpd2d *x, t_float max)

-{

-    x->maxY = max;

-}

-

-void pmpd2d_min(t_pmpd2d *x, t_float minX, t_float minY)

-{

-    x->minX = minX;

-    x->minY = minY;

-}

-

-void pmpd2d_max(t_pmpd2d *x, t_float maxX, t_float maxY)

-{

-    x->maxX = maxX;

-    x->maxY = maxY;

-}

-

-void pmpd2d_setFixed(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( argv[0].a_type == A_FLOAT ) 

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].mobile = 0;

-    }

-    if ( argv[0].a_type == A_SYMBOL )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].mobile = 0;

-            }

-        }

-    }

-}

-

-void pmpd2d_setMobile(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( argv[0].a_type == A_FLOAT )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].mobile = 1;

-    }

-    if ( argv[0].a_type == A_SYMBOL ) 

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].mobile = 1;

-            }

-        }

-    }

-}

-

-void pmpd2d_setD2(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].D2 = atom_getfloatarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].D2 = atom_getfloatarg(1, argc, argv);

-            }

-        }

-    }

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argc == 1 ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            x->mass[i].D2 = atom_getfloatarg(0, argc, argv);

-        }

-    }

-}

-

-void pmpd2d_setD2offset(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].D2offset = atom_getfloatarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].D2offset = atom_getfloatarg(1, argc, argv);

-            }

-        }

-    }

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argc == 1 ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            x->mass[i].D2offset = atom_getfloatarg(0, argc, argv);

-        }

-    }

-}

-

-void pmpd2d_setSpeedX(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].speedX = atom_getfloatarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].speedX = atom_getfloatarg(1, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_setSpeedY(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].speedY = atom_getfloatarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].speedY = atom_getfloatarg(1, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_setSpeed(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].speedX = atom_getfloatarg(1, argc, argv);

-        x->mass[tmp].speedY = atom_getfloatarg(2, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].speedX = atom_getfloatarg(1, argc, argv);

-                x->mass[i].speedY = atom_getfloatarg(2, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_addPosX(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].posX += atom_getfloatarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].posX += atom_getfloatarg(1, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_addPosY(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].posY += atom_getfloatarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].posY += atom_getfloatarg(1, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_addPos(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].posX += atom_getfloatarg(1, argc, argv);

-        x->mass[tmp].posY += atom_getfloatarg(2, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].posX += atom_getfloatarg(1, argc, argv);

-                x->mass[i].posY += atom_getfloatarg(2, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_setForceX(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].forceX = atom_getfloatarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].forceX = atom_getfloatarg(1, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_setForceY(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].forceY = atom_getfloatarg(1, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].forceY = atom_getfloatarg(1, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_setForce(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int tmp, i;

-

-    if ( ( argv[0].a_type == A_FLOAT ) && ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        tmp = atom_getfloatarg(0, argc, argv);

-        tmp = max(0, min( x->nb_mass-1, tmp));

-        x->mass[tmp].forceX = atom_getfloatarg(1, argc, argv);

-        x->mass[tmp].forceY = atom_getfloatarg(2, argc, argv);

-    }

-    if ( ( argv[0].a_type == A_SYMBOL ) && ( argv[1].a_type == A_FLOAT ) && ( argv[2].a_type == A_FLOAT ) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            {

-                x->mass[i].forceX = atom_getfloatarg(1, argc, argv);

-                x->mass[i].forceY = atom_getfloatarg(2, argc, argv);

-            }

-        }

-    }

-}

-

-void pmpd2d_get(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i;

-    t_symbol *toget; 

-    t_atom  toout[5];

-    toget = atom_getsymbolarg(0, argc, argv);

-

-    if ( (toget == gensym("massesPos")) && (argv[1].a_type == A_FLOAT) )

-    {

-        i = atom_getfloatarg(1, argc, argv);

-        if ( (i>=0) && (i<x->nb_mass) )

-        {

-            SETFLOAT(&(toout[0]), i);

-            SETFLOAT(&(toout[1]), x->mass[i].posX);

-            SETFLOAT(&(toout[2]), x->mass[i].posY);

-            outlet_anything(x->main_outlet, gensym("massesPosNo"), 3, toout);

-        }  

-    }

-    else

-    if ( (toget == gensym("massesPos")) && (argv[1].a_type == A_SYMBOL) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(1,argc,argv) == x->mass[i].Id)

-            {

-                SETFLOAT(&(toout[0]), i);

-                SETFLOAT(&(toout[1]), x->mass[i].posX);

-                SETFLOAT(&(toout[2]), x->mass[i].posY);

-                outlet_anything(x->main_outlet, gensym("massesPosId"), 3, toout);

-            }

-        } 

-    }

-    else

-    if ( (toget == gensym("massesPos")) && (argc == 1) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            SETFLOAT(&(toout[0]), i);

-            SETFLOAT(&(toout[1]), x->mass[i].posX);

-            SETFLOAT(&(toout[2]), x->mass[i].posY);

-            outlet_anything(x->main_outlet, gensym("massesPos"), 3, toout);

-        } 

-    }

-    else

-    if ( (toget == gensym("massesPosName")) && (argv[1].a_type == A_FLOAT) )

-    {

-        i = atom_getfloatarg(1, argc, argv);

-        if ( (i>=0) && (i<x->nb_mass) )

-        {

-            SETSYMBOL(&(toout[0]), x->mass[i].Id);

-            SETFLOAT(&(toout[1]), x->mass[i].posX);

-            SETFLOAT(&(toout[2]), x->mass[i].posY);

-            outlet_anything(x->main_outlet, gensym("massesPosNameNo"), 3, toout);

-        }  

-    }

-    else

-    if ( (toget == gensym("massesPosName")) && (argv[1].a_type == A_SYMBOL) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(1,argc,argv) == x->mass[i].Id)

-            {

-                SETSYMBOL(&(toout[0]), x->mass[i].Id);

-                SETFLOAT(&(toout[1]), x->mass[i].posX);

-                SETFLOAT(&(toout[2]), x->mass[i].posY);

-                outlet_anything(x->main_outlet, gensym("massesPosNameId"), 3, toout);

-            }

-        } 

-    }

-    else

-    if ( (toget == gensym("massesPosName")) && (argc == 1) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            SETSYMBOL(&(toout[0]), x->mass[i].Id);

-            SETFLOAT(&(toout[1]), x->mass[i].posX);

-            SETFLOAT(&(toout[2]), x->mass[i].posY);

-            outlet_anything(x->main_outlet, gensym("massesPosName"), 3, toout);

-        } 

-    }

-    else

-    if ( (toget == gensym("massesSpeeds")) && (argv[1].a_type == A_FLOAT) )

-    {

-        i = atom_getfloatarg(1, argc, argv);

-        if ( (i>=0) && (i<x->nb_mass) )

-        {

-            SETFLOAT(&(toout[0]), i);

-            SETFLOAT(&(toout[1]), x->mass[i].speedX);

-            SETFLOAT(&(toout[2]), x->mass[i].speedY);

-            outlet_anything(x->main_outlet, gensym("massesSpeedsNo"), 3, toout);

-        }  

-    }

-    else

-    if ( (toget == gensym("massesSpeeds")) && (argv[1].a_type == A_SYMBOL) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(1,argc,argv) == x->mass[i].Id)

-            {

-                SETFLOAT(&(toout[0]), i);

-                SETFLOAT(&(toout[1]), x->mass[i].speedX);

-                SETFLOAT(&(toout[2]), x->mass[i].speedY);

-                outlet_anything(x->main_outlet, gensym("massesSpeedsId"), 3, toout);

-            }

-        } 

-    }

-    else

-    if ( (toget == gensym("massesSpeeds")) && (argc == 1) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            SETFLOAT(&(toout[0]), i);

-            SETFLOAT(&(toout[1]), x->mass[i].speedX);

-            SETFLOAT(&(toout[2]), x->mass[i].speedY);

-            outlet_anything(x->main_outlet, gensym("massesSpeeds"), 3, toout);

-        } 

-    }

-    else

-    if ( (toget == gensym("massesSpeedsName")) && (argv[1].a_type == A_FLOAT) )

-    {

-        i = atom_getfloatarg(1, argc, argv);

-        if ( (i>=0) && (i<x->nb_mass) )

-        {

-            SETSYMBOL(&(toout[0]), x->mass[i].Id);

-            SETFLOAT(&(toout[1]), x->mass[i].speedX);

-            SETFLOAT(&(toout[2]), x->mass[i].speedY);

-            outlet_anything(x->main_outlet, gensym("massesSpeedsNameNo"), 3, toout);

-        }  

-    }

-    else

-    if ( (toget == gensym("massesSpeedsName")) && (argv[1].a_type == A_SYMBOL) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(1,argc,argv) == x->mass[i].Id)

-            {

-                SETSYMBOL(&(toout[0]), x->mass[i].Id);

-                SETFLOAT(&(toout[1]), x->mass[i].speedX);

-                SETFLOAT(&(toout[2]), x->mass[i].speedY);

-                outlet_anything(x->main_outlet, gensym("massesSpeedsNameId"), 3, toout);

-            }

-        } 

-    }

-    else

-    if ( (toget == gensym("massesSpeedsName")) && (argc == 1) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            SETSYMBOL(&(toout[0]), x->mass[i].Id);

-            SETFLOAT(&(toout[1]), x->mass[i].speedX);

-            SETFLOAT(&(toout[2]), x->mass[i].speedY);

-            outlet_anything(x->main_outlet, gensym("massesSpeedsName"), 3, toout);

-        } 

-    }

-    else

-    if ( (toget == gensym("massesForces")) && (argv[1].a_type == A_FLOAT) )

-    {

-        i = atom_getfloatarg(1, argc, argv);

-        if ( (i>=0) && (i<x->nb_mass) )

-        {

-            SETFLOAT(&(toout[0]), i);

-            SETFLOAT(&(toout[1]), x->mass[i].forceX);

-            SETFLOAT(&(toout[2]), x->mass[i].forceY);

-            outlet_anything(x->main_outlet, gensym("massesForcesNo"), 3, toout);

-        }  

-    }

-    else

-    if ( (toget == gensym("massesForces")) && (argv[1].a_type == A_SYMBOL) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(1,argc,argv) == x->mass[i].Id)

-            {

-                SETFLOAT(&(toout[0]), i);

-                SETFLOAT(&(toout[1]), x->mass[i].forceX);

-                SETFLOAT(&(toout[2]), x->mass[i].forceY);

-                outlet_anything(x->main_outlet, gensym("massesForcesId"), 3, toout);

-            }

-        } 

-    }

-    else

-    if ( (toget == gensym("massesForces")) && (argc == 1) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            SETFLOAT(&(toout[0]), i);

-            SETFLOAT(&(toout[1]), x->mass[i].forceX);

-            SETFLOAT(&(toout[2]), x->mass[i].forceY);

-            outlet_anything(x->main_outlet, gensym("massesForces"), 3, toout);

-        } 

-    }

-    else

-    if ( (toget == gensym("massesForcesName")) && (argv[1].a_type == A_FLOAT) )

-    {

-        i = atom_getfloatarg(1, argc, argv);

-        if ( (i>=0) && (i<x->nb_mass) )

-        {

-            SETSYMBOL(&(toout[0]), x->mass[i].Id);

-            SETFLOAT(&(toout[1]), x->mass[i].forceX);

-            SETFLOAT(&(toout[2]), x->mass[i].forceY);

-            outlet_anything(x->main_outlet, gensym("massesForcesNameNo"), 3, toout);

-        }

-    }

-    else

-    if ( (toget == gensym("massesForcesName")) && (argv[1].a_type == A_SYMBOL) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if ( atom_getsymbolarg(1,argc,argv) == x->mass[i].Id)

-            {

-                SETSYMBOL(&(toout[0]), x->mass[i].Id);

-                SETFLOAT(&(toout[1]), x->mass[i].forceX);

-                SETFLOAT(&(toout[2]), x->mass[i].forceY);

-                outlet_anything(x->main_outlet, gensym("massesForcesNameId"), 3, toout);

-            }

-        } 

-    }

-    else

-    if ( (toget == gensym("massesForcesName")) && (argc == 1) )

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            SETSYMBOL(&(toout[0]), x->mass[i].Id);

-            SETFLOAT(&(toout[1]), x->mass[i].forceX);

-            SETFLOAT(&(toout[2]), x->mass[i].forceY);

-            outlet_anything(x->main_outlet, gensym("massesForcesName"), 3, toout);

-        } 

-    }

-    else

-    if ( (toget == gensym("linksPos")) && (argv[1].a_type == A_FLOAT) )

-    {

-        i = atom_getfloatarg(1, argc, argv);

-        if ( (i>=0) && (i<x->nb_mass) )

-        {

-            SETFLOAT(&(toout[0]), i);

-            SETFLOAT(&(toout[1]), x->link[i].mass1->posX);

-            SETFLOAT(&(toout[2]), x->link[i].mass1->posY);

-            SETFLOAT(&(toout[3]), x->link[i].mass2->posX);

-            SETFLOAT(&(toout[4]), x->link[i].mass2->posY);

-            outlet_anything(x->main_outlet, gensym("linksPosNo"), 5, toout);

-        }

-    }

-    else

-    if ( (toget == gensym("linksPos")) && (argv[1].a_type == A_SYMBOL) )

-    {

-        for (i=0; i< x->nb_link; i++)

-        {

-            if ( atom_getsymbolarg(1,argc,argv) == x->link[i].Id)

-            {

-                SETFLOAT(&(toout[0]), i);

-                SETFLOAT(&(toout[1]), x->link[i].mass1->posX);

-                SETFLOAT(&(toout[2]), x->link[i].mass1->posY);

-                SETFLOAT(&(toout[3]), x->link[i].mass2->posX);

-                SETFLOAT(&(toout[4]), x->link[i].mass2->posY);

-                outlet_anything(x->main_outlet, gensym("linksPosNo"), 5, toout);

-            }

-        } 

-    }

-    else

-    if ( (toget == gensym("linksPos")) && (argc == 1) )

-    {

-        for (i=0; i< x->nb_link; i++)

-        {

-            SETFLOAT(&(toout[0]), i);

-            SETFLOAT(&(toout[1]), x->link[i].mass1->posX);

-            SETFLOAT(&(toout[2]), x->link[i].mass1->posY);

-            SETFLOAT(&(toout[3]), x->link[i].mass2->posX);

-            SETFLOAT(&(toout[4]), x->link[i].mass2->posY);

-            outlet_anything(x->main_outlet, gensym("linksPosNo"), 5, toout);

-        } 

-    }

-    else

-    if ( (toget == gensym("linksPosName")) && (argv[1].a_type == A_FLOAT) )

-    {

-        i = atom_getfloatarg(1, argc, argv);

-        if ( (i>=0) && (i<x->nb_mass) )

-        {

-            SETFLOAT(&(toout[0]), i);

-            SETFLOAT(&(toout[1]), x->link[i].mass1->posX);

-            SETFLOAT(&(toout[2]), x->link[i].mass1->posY);

-            SETFLOAT(&(toout[3]), x->link[i].mass2->posX);

-            SETFLOAT(&(toout[4]), x->link[i].mass2->posY);

-            outlet_anything(x->main_outlet, gensym("linksPosNo"), 5, toout);

-        }

-    }

-    else

-    if ( (toget == gensym("linksPosName")) && (argv[1].a_type == A_SYMBOL) )

-    {

-        for (i=0; i< x->nb_link; i++)

-        {

-            if ( atom_getsymbolarg(1,argc,argv) == x->link[i].Id)

-            {

-            SETFLOAT(&(toout[0]), i);

-            SETFLOAT(&(toout[1]), x->link[i].mass1->posX);

-            SETFLOAT(&(toout[2]), x->link[i].mass1->posY);

-            SETFLOAT(&(toout[3]), x->link[i].mass2->posX);

-            SETFLOAT(&(toout[4]), x->link[i].mass2->posY);

-            outlet_anything(x->main_outlet, gensym("linksPosNo"), 5, toout);

-            }

-        } 

-    }

-    else

-    if ( (toget == gensym("linksPosName")) && (argc == 1) )

-    {

-        for (i=0; i< x->nb_link; i++)

-        {

-            SETFLOAT(&(toout[0]), i);

-            SETFLOAT(&(toout[1]), x->link[i].mass1->posX);

-            SETFLOAT(&(toout[2]), x->link[i].mass1->posY);

-            SETFLOAT(&(toout[3]), x->link[i].mass2->posX);

-            SETFLOAT(&(toout[4]), x->link[i].mass2->posY);

-            outlet_anything(x->main_outlet, gensym("linksPosNo"), 5, toout);

-        } 

-    }

-    else

-        error("not get attribute");

-}

-

-void pmpd2d_massesPosXL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_mass];

-

-    for (i=0; i < x->nb_mass; i++)

-    {

-        SETFLOAT(&(pos_list[i]),x->mass[i].posX);

-    }

-    outlet_anything(x->main_outlet, gensym("massesPosXL"),x->nb_mass , pos_list);

-}

-

-void pmpd2d_massesForcesXL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_mass];

-

-    for (i=0; i< x->nb_mass; i++)

-    {

-        SETFLOAT(&(pos_list[i]),x->mass[i].forceX);

-    }

-    outlet_anything(x->main_outlet, gensym("massesForcesXL"),x->nb_mass , pos_list);

-}

-

-void pmpd2d_massesSpeedsXL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_mass];

-

-    for (i=0; i< x->nb_mass; i++)

-    {

-        SETFLOAT(&(pos_list[i]),x->mass[i].speedX);

-    }

-    outlet_anything(x->main_outlet, gensym("massesSpeedsXL"),x->nb_mass , pos_list);

-}

-

-void pmpd2d_massesPosXT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_mass;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = x->mass[i].posX;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_mass))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->mass[j].Id)

-                {

-                    vec[i].w_float = x->mass[j].posX;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_massesSpeedsXT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {        

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_mass;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = x->mass[i].speedX;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_mass))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->mass[j].Id)

-                {

-                    vec[i].w_float = x->mass[j].speedX;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_massesForcesXT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_mass;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = x->mass[i].forceX;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_mass))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->mass[j].Id)

-                {

-                    vec[i].w_float = x->mass[j].forceX;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_massesPosYL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_mass];

-

-    for (i=0; i < x->nb_mass; i++)

-    {

-        SETFLOAT(&(pos_list[i]),x->mass[i].posY);

-    }

-    outlet_anything(x->main_outlet, gensym("massesPosYL"),x->nb_mass , pos_list);

-}

-

-void pmpd2d_massesForcesYL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_mass];

-

-    for (i=0; i< x->nb_mass; i++)

-    {

-        SETFLOAT(&(pos_list[i]),x->mass[i].forceY);

-    }

-    outlet_anything(x->main_outlet, gensym("massesForcesYL"),x->nb_mass , pos_list);

-}

-

-void pmpd2d_massesSpeedsYL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_mass];

-

-    for (i=0; i< x->nb_mass; i++)

-    {

-        SETFLOAT(&(pos_list[i]),x->mass[i].speedY);

-    }

-    outlet_anything(x->main_outlet, gensym("massesSpeedsYL"),x->nb_mass , pos_list);

-}

-

-void pmpd2d_massesPosYT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_mass;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = x->mass[i].posY;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_mass))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->mass[j].Id)

-                {

-                    vec[i].w_float = x->mass[j].posY;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_massesSpeedsYT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_mass;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = x->mass[i].speedY;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_mass))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->mass[j].Id)

-                {

-                    vec[i].w_float = x->mass[j].speedY;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_massesForcesYT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_mass;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = x->mass[i].forceY;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_mass))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->mass[j].Id)

-                {

-                    vec[i].w_float = x->mass[j].forceY;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_massesPosL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[2*x->nb_mass];

-

-    for (i=0; i < x->nb_mass; i++)

-    {

-        SETFLOAT(&(pos_list[2*i]  ),x->mass[i].posX);

-        SETFLOAT(&(pos_list[2*i+1]),x->mass[i].posY);

-    }

-    outlet_anything(x->main_outlet, gensym("massesPosL"),2*x->nb_mass , pos_list);

-}

-

-void pmpd2d_massesForcesL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[2*x->nb_mass];

-

-    for (i=0; i< x->nb_mass; i++)

-    {

-        SETFLOAT(&(pos_list[2*i]  ),x->mass[i].forceX);

-        SETFLOAT(&(pos_list[2*i+1]),x->mass[i].forceY);

-    }

-    outlet_anything(x->main_outlet, gensym("massesForcesL"),2*x->nb_mass , pos_list);

-}

-

-void pmpd2d_massesSpeedsL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[2*x->nb_mass];

-

-    for (i=0; i< x->nb_mass; i++)

-    {

-        SETFLOAT(&(pos_list[2*i]  ),x->mass[i].speedX);

-        SETFLOAT(&(pos_list[2*i+1]),x->mass[i].speedY);

-    }

-    outlet_anything(x->main_outlet, gensym("massesSpeedsL"),2*x->nb_mass , pos_list);

-}

-

-void pmpd2d_massesPosT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_mass;

-            taille_max = min(taille_max, (vecsize-1)/2);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[2*i  ].w_float = x->mass[i].posX;

-                vec[2*i+1].w_float = x->mass[i].posY;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize-1) && (j < x->nb_mass))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->mass[j].Id)

-                {

-                    vec[i].w_float = x->mass[j].posX;

-                    i++;

-                    vec[i].w_float = x->mass[j].posY;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_massesSpeedsT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_mass;

-            taille_max = min(taille_max, (vecsize-1)/2);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[2*i  ].w_float = x->mass[i].speedX;

-                vec[2*i+1].w_float = x->mass[i].speedY;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize-1) && (j < x->nb_mass))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->mass[j].Id)

-                {

-                    vec[i].w_float = x->mass[j].speedX;

-                    i++;

-                    vec[i].w_float = x->mass[j].speedY;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_massesForcesT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_mass;

-            taille_max = min(taille_max, (vecsize-1)/2);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[2*i  ].w_float = x->mass[i].forceX;

-                vec[2*i+1].w_float = x->mass[i].forceY;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize-1) && (j < x->nb_mass))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->mass[j].Id)

-                {

-                    vec[i].w_float = x->mass[j].forceX;

-                    i++;

-                    vec[i].w_float = x->mass[j].forceY;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_massesPosNormL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_mass];

-

-    for (i=0; i < x->nb_mass; i++)

-    {

-        SETFLOAT(&(pos_list[i]),sqrt(sqr(x->mass[i].posX)+sqr(x->mass[i].posY)));

-    }

-    outlet_anything(x->main_outlet, gensym("massesPosNormL"),x->nb_mass , pos_list);

-}

-

-void pmpd2d_massesForcesNormL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_mass];

-

-    for (i=0; i< x->nb_mass; i++)

-    {

-        SETFLOAT(&(pos_list[i]),sqrt(sqr(x->mass[i].forceX)+sqr(x->mass[i].forceY)));

-    }

-    outlet_anything(x->main_outlet, gensym("massesForcesNormL"),x->nb_mass , pos_list);

-}

-

-void pmpd2d_massesSpeedsNormL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_mass];

-

-    for (i=0; i< x->nb_mass; i++)

-    {

-        SETFLOAT(&(pos_list[i]),sqrt(sqr(x->mass[i].speedX)+sqr(x->mass[i].speedY)));

-    }

-    outlet_anything(x->main_outlet, gensym("massesSpeedsNormL"),x->nb_mass , pos_list);

-}

-

-void pmpd2d_massesPosNormT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_mass;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = sqrt(sqr(x->mass[i].posX)+sqr(x->mass[i].posY));

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_mass))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->mass[j].Id)

-                {

-                    vec[i].w_float = sqrt(sqr(x->mass[j].posX)+sqr(x->mass[j].posY));

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_massesSpeedsNormT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_mass;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = sqrt(sqr(x->mass[i].speedX)+sqr(x->mass[i].speedY));

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_mass))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->mass[j].Id)

-                {

-                    vec[i].w_float = sqrt(sqr(x->mass[j].speedX)+sqr(x->mass[j].speedY));

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_massesForcesNormT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_mass;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = sqrt(sqr(x->mass[i].forceX)+sqr(x->mass[i].forceY));

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_mass))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->mass[j].Id)

-                {

-                    vec[i].w_float = sqrt(sqr(x->mass[j].forceX)+sqr(x->mass[j].forceY));

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_massesPosMean(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    t_float sommeX, sommeY, somme;

-    t_int i,j;

-    t_atom mean[3];

-

-    sommeX = 0;

-    sommeY = 0;

-    somme = 0;

-    j = 0;

-    

-    if ( (argc >= 1) && (argv[0].a_type == A_SYMBOL) ) 

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if (atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            { 

-                sommeX += x->mass[i].posX;

-                sommeY += x->mass[i].posY;

-                somme +=  sqrt(sqr(x->mass[i].posX) + sqr(x->mass[i].posY)); // distance au centre

-                j++;

-            }

-        }

-    }

-    else

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-                sommeX += x->mass[i].posX;

-                sommeY += x->mass[i].posY;

-                somme +=  sqrt(sqr(x->mass[i].posX) + sqr(x->mass[i].posY)); // distance au centre

-                j++;

-        }

-    }    

-    

-    sommeX /= j;

-    sommeY /= j;

-    somme  /= j;    

-    

-    SETFLOAT(&(mean[0]),sommeX);

-    SETFLOAT(&(mean[1]),sommeY);

-    SETFLOAT(&(mean[2]),somme);

-    

-    outlet_anything(x->main_outlet, gensym("massesPosMean"),3 , mean);

-}

-

-void pmpd2d_massesPosStd(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    t_float sommeX, sommeY, somme;

-    t_int i,j;

-    t_float stdX, stdY,std;

-    t_atom std_out[3];

-

-    sommeX = 0;

-    sommeY = 0;

-    somme  = 0;

-    stdX = 0;

-    stdY = 0;

-    std  = 0;

-    j = 0;

-    

-    if ( (argc >= 1) && (argv[0].a_type == A_SYMBOL) ) 

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if (atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            { 

-                sommeX += x->mass[i].posX;

-                sommeY += x->mass[i].posY;

-                somme +=  sqrt(sqr(x->mass[i].posX) + sqr(x->mass[i].posY)); // distance au centre

-                j++;

-            }

-        }

-        sommeX /= j;

-        sommeY /= j;

-        somme /= j;

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if (atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            { 

-                stdX += sqr(x->mass[i].posX-sommeX);

-                stdY += sqr(x->mass[i].posY-sommeY);

-                std  +=  sqr(sqrt(sqr(x->mass[i].posX) + sqr(x->mass[i].posY))-somme);

-            }

-        }        

-    }

-    else

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            sommeX += x->mass[i].posX;

-            sommeY += x->mass[i].posY;

-            somme +=  sqrt(sqr(x->mass[i].posX) + sqr(x->mass[i].posY)); // distance au centre

-            j++;

-        }

-        sommeX /= j;

-        sommeY /= j;

-        somme /= j;

-        for (i=0; i< x->nb_mass; i++)

-        {

-            stdX += sqr(x->mass[i].posX-sommeX);

-            stdY += sqr(x->mass[i].posY-sommeY);

-            std  += sqr(sqrt(sqr(x->mass[i].posX) + sqr(x->mass[i].posY))-somme);

-        }

-    }    

-    

-    stdX = sqrt(stdX/j);

-    stdY = sqrt(stdY/j);

-    std  = sqrt(std /j);    

-

-    SETFLOAT(&(std_out[0]),stdX);

-    SETFLOAT(&(std_out[1]),stdY);

-    SETFLOAT(&(std_out[2]),std);

-    

-    outlet_anything(x->main_outlet, gensym("massesPosStd"),3 , std_out);

-}

-

-void pmpd2d_massesForcesMean(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    t_float sommeX, sommeY, somme;

-    t_int i,j;

-    t_atom mean[3];

-

-    sommeX = 0;

-    sommeY = 0;

-    somme = 0;

-    j = 0;

-    

-    if ( (argc >= 1) && (argv[0].a_type == A_SYMBOL) ) 

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if (atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            { 

-                sommeX += x->mass[i].forceX;

-                sommeY += x->mass[i].forceY;

-                somme +=  sqrt(sqr(x->mass[i].forceX) + sqr(x->mass[i].forceY)); // distance au centre

-                j++;

-            }

-        }

-    }

-    else

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-                sommeX += x->mass[i].forceX;

-                sommeY += x->mass[i].forceY;

-                somme +=  sqrt(sqr(x->mass[i].forceX) + sqr(x->mass[i].forceY)); // distance au centre

-                j++;

-        }

-    }    

-    

-    sommeX /= j;

-    sommeY /= j;

-    somme  /= j;    

-    

-    SETFLOAT(&(mean[0]),sommeX);

-    SETFLOAT(&(mean[1]),sommeY);

-    SETFLOAT(&(mean[2]),somme);

-    

-    outlet_anything(x->main_outlet, gensym("massesForcesMean"),3 , mean);

-}

-

-void pmpd2d_massesForcesStd(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    t_float sommeX, sommeY, somme;

-    t_int i,j;

-    t_float stdX, stdY,std;

-    t_atom std_out[3];

-

-    sommeX = 0;

-    sommeY = 0;

-    somme  = 0;

-    stdX = 0;

-    stdY = 0;

-    std  = 0;

-    j = 0;

-    

-    if ( (argc >= 1) && (argv[0].a_type == A_SYMBOL) ) 

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if (atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            { 

-                sommeX += x->mass[i].forceX;

-                sommeY += x->mass[i].forceY;

-                somme +=  sqrt(sqr(x->mass[i].forceX) + sqr(x->mass[i].forceY)); // distance au centre

-                j++;

-            }

-        }

-        sommeX /= j;

-        sommeY /= j;

-        somme /= j;

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if (atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            { 

-                stdX += sqr(x->mass[i].forceX-sommeX);

-                stdY += sqr(x->mass[i].forceY-sommeY);

-                std  +=  sqr(sqrt(sqr(x->mass[i].forceX) + sqr(x->mass[i].forceY))-somme);

-            }

-        }        

-    }

-    else

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            sommeX += x->mass[i].forceX;

-            sommeY += x->mass[i].forceY;

-            somme +=  sqrt(sqr(x->mass[i].forceX) + sqr(x->mass[i].forceY)); // distance au centre

-            j++;

-        }

-        sommeX /= j;

-        sommeY /= j;

-        somme /= j;

-        for (i=0; i< x->nb_mass; i++)

-        {

-            stdX += sqr(x->mass[i].forceX-sommeX);

-            stdY += sqr(x->mass[i].forceY-sommeY);

-            std  += sqr(sqrt(sqr(x->mass[i].forceX) + sqr(x->mass[i].forceY))-somme);

-        }

-    }    

-    

-    stdX = sqrt(stdX/j);

-    stdY = sqrt(stdY/j);

-    std  = sqrt(std /j);    

-

-    SETFLOAT(&(std_out[0]),stdX);

-    SETFLOAT(&(std_out[1]),stdY);

-    SETFLOAT(&(std_out[2]),std);

-    

-    outlet_anything(x->main_outlet, gensym("massesForcesStd"),3 , std_out);

-}

-

-void pmpd2d_massesSpeedsMean(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    t_float sommeX, sommeY, somme;

-    t_int i,j;

-    t_atom mean[3];

-

-    sommeX = 0;

-    sommeY = 0;

-    somme = 0;

-    j = 0;

-    

-    if ( (argc >= 1) && (argv[0].a_type == A_SYMBOL) ) 

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if (atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            { 

-                sommeX += x->mass[i].speedX;

-                sommeY += x->mass[i].speedY;

-                somme +=  sqrt(sqr(x->mass[i].speedX) + sqr(x->mass[i].speedY)); // distance au centre

-                j++;

-            }

-        }

-    }

-    else

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-                sommeX += x->mass[i].speedX;

-                sommeY += x->mass[i].speedY;

-                somme +=  sqrt(sqr(x->mass[i].speedX) + sqr(x->mass[i].speedY)); // distance au centre

-                j++;

-        }

-    }    

-    

-    sommeX /= j;

-    sommeY /= j;

-    somme  /= j;    

-    

-    SETFLOAT(&(mean[0]),sommeX);

-    SETFLOAT(&(mean[1]),sommeY);

-    SETFLOAT(&(mean[2]),somme);

-    

-    outlet_anything(x->main_outlet, gensym("massesSpeedsMean"),3 , mean);

-}

-

-void pmpd2d_massesSpeedsStd(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    t_float sommeX, sommeY, somme;

-    t_int i,j;

-    t_float stdX, stdY,std;

-    t_atom std_out[3];

-

-    sommeX = 0;

-    sommeY = 0;

-    somme  = 0;

-    stdX = 0;

-    stdY = 0;

-    std  = 0;

-    j = 0;

-    

-    if ( (argc >= 1) && (argv[0].a_type == A_SYMBOL) ) 

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if (atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            { 

-                sommeX += x->mass[i].speedX;

-                sommeY += x->mass[i].speedY;

-                somme +=  sqrt(sqr(x->mass[i].speedX) + sqr(x->mass[i].speedY)); // distance au centre

-                j++;

-            }

-        }

-        sommeX /= j;

-        sommeY /= j;

-        somme /= j;

-        for (i=0; i< x->nb_mass; i++)

-        {

-            if (atom_getsymbolarg(0,argc,argv) == x->mass[i].Id)

-            { 

-                stdX += sqr(x->mass[i].speedX-sommeX);

-                stdY += sqr(x->mass[i].speedY-sommeY);

-                std  +=  sqr(sqrt(sqr(x->mass[i].speedX) + sqr(x->mass[i].speedY))-somme);

-            }

-        }        

-    }

-    else

-    {

-        for (i=0; i< x->nb_mass; i++)

-        {

-            sommeX += x->mass[i].speedX;

-            sommeY += x->mass[i].speedY;

-            somme +=  sqrt(sqr(x->mass[i].speedX) + sqr(x->mass[i].speedY)); // distance au centre

-            j++;

-        }

-        sommeX /= j;

-        sommeY /= j;

-        somme /= j;

-        for (i=0; i< x->nb_mass; i++)

-        {

-            stdX += sqr(x->mass[i].speedX-sommeX);

-            stdY += sqr(x->mass[i].speedY-sommeY);

-            std  += sqr(sqrt(sqr(x->mass[i].speedX) + sqr(x->mass[i].speedY))-somme);

-        }

-    }    

-    

-    stdX = sqrt(stdX/j);

-    stdY = sqrt(stdY/j);

-    std  = sqrt(std /j);    

-

-    SETFLOAT(&(std_out[0]),stdX);

-    SETFLOAT(&(std_out[1]),stdY);

-    SETFLOAT(&(std_out[2]),std);

-    

-    outlet_anything(x->main_outlet, gensym("massesSpeedsStd"),3 , std_out);

-}

-

-// --------------------------------------------

-

-void pmpd2d_linksPosXL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[i]),(x->link[i].mass1->posX + x->link[i].mass2->posX)/2);

-    }

-    outlet_anything(x->main_outlet, gensym("linksPosXL"),x->nb_link , pos_list);

-}

-

-void pmpd2d_linksLengthXL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[i]),x->link[i].mass2->posX - x->link[i].mass1->posX);

-    }

-    outlet_anything(x->main_outlet, gensym("linksLengthXL"),x->nb_link , pos_list);

-}

-

-void pmpd2d_linksPosSpeedXL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[i]),(x->link[i].mass1->speedX + x->link[i].mass2->speedX)/2);

-    }

-    outlet_anything(x->main_outlet, gensym("linksPosSpeedXL"),x->nb_link , pos_list);

-}

-

-void pmpd2d_linksLengthSpeedXL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[i]),x->link[i].mass2->speedX - x->link[i].mass1->speedX);

-    }

-    outlet_anything(x->main_outlet, gensym("linksLengthSpeedXL"),x->nb_link , pos_list);

-}

-

-void pmpd2d_linksPosXT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = (x->link[i].mass1->posX + x->link[i].mass2->posX)/2;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = (x->link[j].mass1->posX + x->link[j].mass2->posX)/2;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksLengthXT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = x->link[i].mass2->posX - x->link[i].mass1->posX;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = x->link[j].mass2->posX - x->link[j].mass1->posX;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksPosSpeedXT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = (x->link[i].mass1->speedX + x->link[i].mass2->speedX)/2;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = (x->link[j].mass1->speedX + x->link[j].mass2->speedX)/2;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksLengthSpeedXT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = x->link[i].mass2->speedX - x->link[i].mass1->speedX;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = x->link[j].mass2->speedX - x->link[j].mass1->speedX;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksPosYL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[i]),(x->link[i].mass1->posY + x->link[i].mass2->posY)/2);

-    }

-    outlet_anything(x->main_outlet, gensym("linksPosYL"),x->nb_link , pos_list);

-}

-

-void pmpd2d_linksLengthYL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[i]),x->link[i].mass2->posY - x->link[i].mass1->posY);

-    }

-    outlet_anything(x->main_outlet, gensym("linksLengthYL"),x->nb_link , pos_list);

-}

-

-void pmpd2d_linksPosSpeedYL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[i]),(x->link[i].mass1->speedY + x->link[i].mass2->speedY)/2);

-    }

-    outlet_anything(x->main_outlet, gensym("linksPosSpeedYL"),x->nb_link , pos_list);

-}

-

-void pmpd2d_linksLengthSpeedYL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[i]),x->link[i].mass2->speedY - x->link[i].mass1->speedY);

-    }

-    outlet_anything(x->main_outlet, gensym("linksLengthSpeedYL"),x->nb_link , pos_list);

-}

-

-void pmpd2d_linksPosYT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = (x->link[i].mass1->posY + x->link[i].mass2->posY)/2;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = (x->link[j].mass1->posY + x->link[j].mass2->posY)/2;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksLengthYT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = x->link[i].mass2->posY - x->link[i].mass1->posY;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = x->link[j].mass2->posY - x->link[j].mass1->posY;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksPosSpeedYT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = (x->link[i].mass1->speedY + x->link[i].mass2->speedY)/2;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = (x->link[j].mass1->speedY + x->link[j].mass2->speedY)/2;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksLengthSpeedYT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = x->link[i].mass2->speedY - x->link[i].mass1->speedY;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = x->link[j].mass2->speedY - x->link[j].mass1->speedY;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksPosL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[2*x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[2*i]  ),(x->link[i].mass2->posX + x->link[i].mass1->posX)/2);

-        SETFLOAT(&(pos_list[2*i+1]),(x->link[i].mass2->posY + x->link[i].mass1->posY)/2);

-    }

-    outlet_anything(x->main_outlet, gensym("linksPosL"),2*x->nb_link , pos_list);

-}

-

-void pmpd2d_linksLengthL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[2*x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[2*i]  ),x->link[i].mass2->posX - x->link[i].mass1->posX);

-        SETFLOAT(&(pos_list[2*i+1]),x->link[i].mass2->posY - x->link[i].mass1->posY);

-    }

-    outlet_anything(x->main_outlet, gensym("linksLengthL"),2*x->nb_link , pos_list);

-}

-

-void pmpd2d_linksPosSpeedL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[2*x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[2*i]  ),(x->link[i].mass2->speedX + x->link[i].mass1->speedX)/2);

-        SETFLOAT(&(pos_list[2*i+1]),(x->link[i].mass2->speedY + x->link[i].mass1->speedY)/2);

-    }

-    outlet_anything(x->main_outlet, gensym("linksPosSpeedL"),3*x->nb_link , pos_list);

-}

-

-void pmpd2d_linksLengthSpeedL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[2*x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[2*i]  ),x->link[i].mass2->speedX - x->link[i].mass1->speedX);

-        SETFLOAT(&(pos_list[2*i+1]),x->link[i].mass2->speedY - x->link[i].mass1->speedY);

-    }

-    outlet_anything(x->main_outlet, gensym("linksLengthSpeedL"),2*x->nb_link , pos_list);

-}

-

-void pmpd2d_linksPosT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, (vecsize-2)/2);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[2*i  ].w_float = (x->link[i].mass2->posX + x->link[i].mass1->posX)/2;

-                vec[2*i+1].w_float = (x->link[i].mass2->posY + x->link[i].mass1->posY)/2;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize-2) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = (x->link[j].mass2->posX + x->link[j].mass1->posX)/2;

-                    i++;

-                    vec[i].w_float = (x->link[j].mass2->posY + x->link[j].mass1->posY)/2;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksLengthT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, (vecsize-2)/2);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[2*i  ].w_float = x->link[i].mass2->posX - x->link[i].mass1->posX;

-                vec[2*i+1].w_float = x->link[i].mass2->posY - x->link[i].mass1->posY;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize-2) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = x->link[j].mass2->posX + x->link[j].mass1->posX;

-                    i++;

-                    vec[i].w_float = x->link[j].mass2->posY + x->link[j].mass1->posY;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksPosSpeedT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, (vecsize-2)/2);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[2*i  ].w_float = (x->link[i].mass2->speedX + x->link[i].mass1->speedX)/2;

-                vec[2*i+1].w_float = (x->link[i].mass2->speedY + x->link[i].mass1->speedY)/2;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize-2) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = (x->link[j].mass2->speedX + x->link[j].mass1->speedX)/2;

-                    i++;

-                    vec[i].w_float = (x->link[j].mass2->speedY + x->link[j].mass1->speedY)/2;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksLengthSpeedT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, (vecsize-2)/2);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[2*i  ].w_float = x->link[i].mass2->speedX - x->link[i].mass1->speedX;

-                vec[2*i+1].w_float = x->link[i].mass2->speedY - x->link[i].mass1->speedY;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize-2) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = x->link[j].mass2->speedX + x->link[j].mass1->speedX;

-                    i++;

-                    vec[i].w_float = x->link[j].mass2->speedY + x->link[j].mass1->speedY;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksPosNormL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[i]),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) ));

-    }

-    outlet_anything(x->main_outlet, gensym("linksPosNormL"),x->nb_link , pos_list);

-}

-

-void pmpd2d_linksLengthNormL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[i]),sqrt( \

-                            sqr(x->link[i].mass2->posX - x->link[i].mass1->posX) + \

-                            sqr(x->link[i].mass2->posY - x->link[i].mass1->posY)  ));

-    }

-    outlet_anything(x->main_outlet, gensym("linksLengthNormL"),x->nb_link , pos_list);

-}

-

-void pmpd2d_linksPosSpeedNormL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[i]),sqrt( \

-                            sqr((x->link[i].mass1->speedX + x->link[i].mass2->speedX)/2) + \

-                            sqr((x->link[i].mass1->speedY + x->link[i].mass2->speedY)/2)  ));

-    }

-    outlet_anything(x->main_outlet, gensym("linksPosSpeedNormL"),x->nb_link , pos_list);

-}

-

-void pmpd2d_linksLengthSpeedNormL(t_pmpd2d *x)

-{

-    int i;

-    t_atom pos_list[x->nb_link];

-

-    for (i=0; i < x->nb_link; i++)

-    {

-        SETFLOAT(&(pos_list[i]),sqrt( \

-                            sqr(x->link[i].mass2->speedX - x->link[i].mass1->speedX) + \

-                            sqr(x->link[i].mass2->speedY - x->link[i].mass1->speedY) ));

-    }

-    outlet_anything(x->main_outlet, gensym("linksLengthSpeedNormL"),x->nb_link , pos_list);

-}

-

-void pmpd2d_linksPosNormT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = 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) );

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[j].w_float = sqrt( \

-                            sqr((x->link[j].mass1->posX + x->link[j].mass2->posX)/2) + \

-                            sqr((x->link[j].mass1->posY + x->link[j].mass2->posY)/2) );

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksLengthNormT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = sqrt( \

-                            sqr(x->link[i].mass2->posX - x->link[i].mass1->posX) + \

-                            sqr(x->link[i].mass2->posY - x->link[i].mass1->posY) );

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = sqrt( \

-                            sqr(x->link[j].mass2->posX - x->link[j].mass1->posX) + \

-                            sqr(x->link[j].mass2->posY - x->link[j].mass1->posY) );

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksPosSpeedNormT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = sqrt( \

-                            sqr((x->link[i].mass1->speedX + x->link[i].mass2->speedX)/2) + \

-                            sqr((x->link[i].mass1->speedY + x->link[i].mass2->speedY)/2) );

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = sqrt( \

-                            sqr((x->link[j].mass1->speedX + x->link[j].mass2->speedX)/2) + \

-                            sqr((x->link[j].mass1->speedY + x->link[j].mass2->speedY)/2) );

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-void pmpd2d_linksLengthSpeedNormT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, vecsize);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[i].w_float = sqrt( \

-                            sqr(x->link[i].mass2->speedX - x->link[i].mass1->speedX) + \

-                            sqr(x->link[i].mass2->speedY - x->link[i].mass1->speedY) );

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = sqrt( \

-                            sqr(x->link[j].mass2->speedX - x->link[j].mass1->speedX) + \

-                            sqr(x->link[j].mass2->speedY - x->link[j].mass1->speedY) );

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-

-void pmpd2d_linksEndT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)

-{

-    int i, j, vecsize;

-    t_garray *a;

-    t_word *vec;

-    

-    if ( (argc==1) && (argv[0].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {

-            int taille_max = x->nb_link;

-            taille_max = min(taille_max, (vecsize-3)/4);

-            for (i=0; i < taille_max ; i++)

-            {

-                vec[4*i  ].w_float = x->link[i].mass1->posX;

-                vec[4*i+1].w_float = x->link[i].mass1->posY;

-                vec[4*i+2].w_float = x->link[i].mass2->posX;

-                vec[4*i+3].w_float = x->link[i].mass2->posY;

-            }

-            garray_redraw(a);

-        }

-    }

-    else 

-    if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) )

-    {

-        t_symbol *tab_name = atom_getsymbolarg(0, argc, argv);

-        if (!(a = (t_garray *)pd_findbyclass(tab_name, garray_class)))

-            pd_error(x, "%s: no such array", tab_name->s_name);

-        else if (!garray_getfloatwords(a, &vecsize, &vec))

-            pd_error(x, "%s: bad template for tabwrite", tab_name->s_name);

-        else

-        {    

-            i = 0;

-            j = 0;

-            while ((i < vecsize-3) && (j < x->nb_link))

-            {

-                if (atom_getsymbolarg(1,argc,argv) == x->link[j].Id)

-                {

-                    vec[i].w_float = x->link[j].mass1->posX;

-                    i++;

-                    vec[i].w_float = x->link[j].mass1->posY;

-                    i++;

-                    vec[i].w_float = x->link[j].mass2->posX;

-                    i++;

-                    vec[i].w_float = x->link[j].mass2->posY;

-                    i++;

-                }

-                j++;

-            }

-            garray_redraw(a);

-        }

-    }

-}

-

-//----------------------------------------------

-

-void pmpd2d_grabMass(t_pmpd2d *x, t_float posX, t_float posY, t_float grab)

-{

-    t_float dist, tmp;

-    t_int i;

-    

-    if (grab == 0)

-        x->grab=0;

-    if ((x->grab == 0)&(grab == 1)&(x->nb_mass > 0))

-    {

-        x->grab=1;

-        x->grab_nb= 0;

-        dist = sqr(x->mass[0].posX - posX) + sqr(x->mass[0].posY - posY);

-        for (i=1; i<x->nb_mass; i++)

-        {

-            tmp = sqr(x->mass[i].posX - posX) + sqr(x->mass[i].posY - posY);

-            if (tmp < dist)

-            {

-                dist = tmp;

-                x->grab_nb= i;

-            }

-        }

-    }

-    if (x->grab == 1)

-    {

-        x->mass[x->grab_nb].posX = posX;

-        x->mass[x->grab_nb].posY = posY;        

-    }

-}

-

-void pmpd2d_closestMass(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv)// t_float posX, t_float posY)

-{

-    t_float dist, tmp;

-    t_int i;

-    t_atom std_out[3];

-    

-    t_float posX, posY;

-    posX = atom_getfloatarg(0, argc, argv);

-    posY = atom_getfloatarg(1, argc, argv);    

-    if ( (argc == 3)  && (argv[2].a_type == A_SYMBOL) )

-    {

-        //t_symbol *mass_name = atom_getsymbolarg(2, argc, argv);

-        if ((x->nb_mass > 0))

-        {

-            dist = 1000000000;//sqr(x->mass[0].posX - posX) + sqr(x->mass[0].posY - posY);

-            for (i=0; i<x->nb_mass; i++)

-            {

-                if (atom_getsymbolarg(2,argc,argv) == x->mass[i].Id)

-                {

-                    tmp = sqr(x->mass[i].posX - posX) + sqr(x->mass[i].posY - posY);

-                    if (tmp < dist)

-                    {

-                        dist = tmp;

-                        x->grab_nb= i;

-                    }

-                }

-            }

-        }

-    }

-    else {

-        if ((x->nb_mass > 0))

-        {

-            dist = sqr(x->mass[0].posX - posX) + sqr(x->mass[0].posY - posY);

-            for (i=1; i<x->nb_mass; i++)

-            {

-                tmp = sqr(x->mass[i].posX - posX) + sqr(x->mass[i].posY - posY);

-                if (tmp < dist)

-                {

-                    dist = tmp;

-                    x->grab_nb= i;

-                }

-            }

-        }

-    }

-    

-        

-    SETFLOAT(&(std_out[0]),x->grab_nb);

-    SETFLOAT(&(std_out[1]), x->mass[x->grab_nb].posX);

-    SETFLOAT(&(std_out[2]), x->mass[x->grab_nb].posY);

-    outlet_anything(x->main_outlet, gensym("closestMass"),3,std_out);

-

-}

-

-void *pmpd2d_new()

-{

-    t_pmpd2d *x = (t_pmpd2d *)pd_new(pmpd2d_class);

-

-    pmpd2d_reset(x);

-    

-    x->main_outlet=outlet_new(&x->x_obj, 0);

-    // x->info_outlet=outlet_new(&x->x_obj, 0); // TODO

-

-    return (void *)x;

-}

+#include "pmpd2d.h"

+#include "pmpd2d_core.c"

+#include "pmpd2d_set.c"

+#include "pmpd2d_get.c"

+#include "pmpd2d_list.c"

+#include "pmpd2d_tab.c"

+#include "pmpd2d_test.c"

+#include "pmpd2d_stat.c"

+#include "pmpd2d_interactor.c"

+#include "pmpd2d_various.c"

+#include "pmpd2d_deprecated.c"

 

 void pmpd2d_setup(void) 

 {

@@ -3799,147 +48,308 @@ void pmpd2d_setup(void)
         (t_newmethod)pmpd2d_new,

         0, sizeof(t_pmpd2d),CLASS_DEFAULT, 0);

 

+/*

+ pmpd2d_core

+ --

+ Basic functions for creation of the structure

+*/

+

     class_addbang(pmpd2d_class, pmpd2d_bang);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_reset,           gensym("reset"), 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_infosL,          gensym("infosL"), 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_infosL,          gensym("print"), 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_mass,            gensym("mass"), A_DEFSYMBOL, A_DEFFLOAT, A_DEFFLOAT, A_DEFFLOAT, A_DEFFLOAT, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_link,            gensym("link"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_tabLink,         gensym("tabLink"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_tLink,           gensym("tLink"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_tabLink,         gensym("tabLink"), A_GIMME, 0);

+

+/*    

+ pmpd2d_set

+ --

+ Functions to modify the internal state of the pmpd3d object

+*/

+

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setK,            gensym("setK"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setD,            gensym("setD"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_setDEnv,         gensym("setDEnv"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_setDEnvOffset,   gensym("setDEnvOffset"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setL,            gensym("setL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_setLCurrent,     gensym("setLCurrent"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setLKTab,        gensym("setLKTab"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setLDTab,        gensym("setLDTab"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setLinkId,       gensym("setLinkId"), A_GIMME, 0);

-

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setMassId,       gensym("setMassId"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_pos,             gensym("pos"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_posX,            gensym("posX"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_posY,            gensym("posY"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_force,           gensym("force"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_forceX,          gensym("forceX"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_forceY,          gensym("forceY"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_minX,            gensym("Xmin"), A_DEFFLOAT, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_maxX,            gensym("Xmax"), A_DEFFLOAT, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_minY,            gensym("Ymin"), A_DEFFLOAT, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_maxY,            gensym("Ymax"), A_DEFFLOAT, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_min,             gensym("min"), A_DEFFLOAT, A_DEFFLOAT, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_max,             gensym("max"), A_DEFFLOAT, A_DEFFLOAT, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setFixed,        gensym("setFixed"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setMobile,       gensym("setMobile"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_setD2,           gensym("setDEnv"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_setD2offset,     gensym("setDEnvOffset"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setSpeed,        gensym("setSpeed"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setSpeedX,       gensym("setSpeedX"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setSpeedY,       gensym("setSpeedY"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setForce,        gensym("setForce"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setForceX,       gensym("setForceX"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_setForceY,       gensym("setForceY"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_addPos,          gensym("addPos"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_addPosX,         gensym("addPosX"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_addPosY,         gensym("addPosY"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_setActive,       gensym("setActive"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_setInactive,     gensym("setInactive"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_pos,             gensym("pos"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_posX,            gensym("posX"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_posY,            gensym("posY"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_pos,             gensym("setPos"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_posX,            gensym("setPosX"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_posY,            gensym("setPosY"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_overdamp,	  	gensym("setOverdamp"), A_GIMME, 0);

+

+/*        

+ pmpd2d_get

+ --

+ Basic functions to output internal state of the object

+ Output syntax : 1 line by element (mass/link)

+*/	

 

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_get,             gensym("get"), A_GIMME, 0);

-    

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosL,      gensym("massesPosL"), 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsL,   gensym("massesSpeedsL"), 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesL,   gensym("massesForcesL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massPos,       	gensym("massPos"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massSpeed,    	gensym("massSpeed"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massForce,    	gensym("massForce"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkPos,     	gensym("linkPos"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkEnd,	       	gensym("linkEnd"), A_GIMME, 0);

+

+/*

+ pmpd2d_list 

+ --

+ Fucntions to output internal state of the object in lists

+ Output Syntax : 1 list with all elements

+*/

+

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massPosL,      gensym("massPosL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massSpeedsL,   gensym("massSpeedL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massForcesL,   gensym("massForceL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massPosXL,     gensym("massPosXL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massSpeedsXL,  gensym("massSpeedXL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massForcesXL,  gensym("massForceXL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massPosYL,     gensym("massPosYL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massSpeedsYL,  gensym("massSpeedYL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massForcesYL,  gensym("massForceYL"), 0);

+

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massPosNormL,      gensym("massPosNormL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massSpeedsNormL,   gensym("massSpeedNormL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massForcesNormL,   gensym("massForceNormL"), 0);

+

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkPosL,            gensym("linkPosL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkLengthL,         gensym("linkLengthL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkPosSpeedL,       gensym("linkPosSpeedL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkLengthSpeedL,    gensym("linkLengthSpeedL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkPosXL,           gensym("linkPosXL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkLengthXL,        gensym("linkLengthXL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkPosSpeedXL,      gensym("linkPosSpeedXL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkLengthSpeedXL,   gensym("linkLengthSpeedXL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkPosYL,           gensym("linkPosYL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkLengthYL,        gensym("linkLengthYL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkPosSpeedYL,      gensym("linkPosSpeedYL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkLengthSpeedYL,   gensym("linkLengthSpeedYL"), A_GIMME, 0);

+

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkPosNormL,               gensym("linkPosNormL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkLengthNormL,           gensym("linkLengthNormL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkPosSpeedNormL,        gensym("linkPosSpeedNormL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkLengthSpeedNormL,    gensym("linkLengthSpeedNormL"), A_GIMME, 0);

+

+/*

+ pmpd2d_tab

+ --

+ Functions to copy the internal state of the object in arrays

+ Output Syntax : none

+*/

+

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosT,      gensym("massesPosT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsT,   gensym("massesSpeedsT"),A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesT,   gensym("massesForcesT"),A_GIMME, 0);

-    

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosXL,     gensym("massesPosXL"), 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsXL,  gensym("massesSpeedsXL"), 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesXL,  gensym("massesForcesXL"), 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosXT,     gensym("massesPosXT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsXT,  gensym("massesSpeedsXT"),A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesXT,  gensym("massesForcesXT"),A_GIMME, 0);

-    

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosYL,     gensym("massesPosYL"), 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsYL,  gensym("massesSpeedsYL"), 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesYL,  gensym("massesForcesYL"), 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosYT,     gensym("massesPosYT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsYT,  gensym("massesSpeedsYT"),A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesYT,  gensym("massesForcesYT"),A_GIMME, 0);

-    

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosNormL,      gensym("massesPosNormL"), 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsNormL,   gensym("massesSpeedsNormL"), 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesNormL,   gensym("massesForcesNormL"), 0);

+

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosT,      gensym("massPosT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsT,   gensym("massSpeedT"),A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesT,   gensym("massForceT"),A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosXT,     gensym("massPosXT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsXT,  gensym("massSpeedXT"),A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesXT,  gensym("massForceXT"),A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosYT,     gensym("massPosYT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsYT,  gensym("massSpeedYT"),A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesYT,  gensym("massForceYT"),A_GIMME, 0);

+

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosNormT,      gensym("massesPosNormT"),  A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsNormT,   gensym("massesSpeedsNormT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesNormT,   gensym("massesForcesNormT"), A_GIMME, 0);

-    

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosMean,       gensym("massesPosMean"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosStd,        gensym("massesPosStd"),A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesMean,    gensym("massesForcesMean"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesStd,     gensym("massesForcesStd"),A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsMean,    gensym("massesSpeedsMean"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsStd,     gensym("massesSpeedsStd"),A_GIMME, 0);

 

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosXL,               gensym("linksPosXL"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthXL,           gensym("linksLengthXL"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedXL,        gensym("linksPosSpeedXL"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedXL,    gensym("linksLengthSpeedXL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosNormT,      gensym("massPosNormT"),  A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsNormT,   gensym("massSpeedsNormT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesNormT,   gensym("massForcesNormT"), A_GIMME, 0);

+

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosT,               gensym("linksPosT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthT,           gensym("linksLengthT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedT,        gensym("linksPosSpeedT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedT,    gensym("linksLengthSpeedT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosXT,               gensym("linksPosXT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthXT,           gensym("linksLengthXT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedXT,        gensym("linksPosSpeedXT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedXT,    gensym("linksLengthSpeedXT"), A_GIMME, 0);

-    

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosYL,               gensym("linksPosYL"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthYL,           gensym("linksLengthYL"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedYL,        gensym("linksPosSpeedYL"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedYL,    gensym("linksLengthSpeedYL"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosYT,               gensym("linksPosYT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthYT,           gensym("linksLengthYT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedYT,        gensym("linksPosSpeedYT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedYT,    gensym("linksLengthSpeedYT"), A_GIMME, 0);

 

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosL,               gensym("linksPosL"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthL,           gensym("linksLengthL"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedL,        gensym("linksPosSpeedL"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedL,    gensym("linksLengthSpeedL"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosT,               gensym("linksPosT"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthT,           gensym("linksLengthT"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedT,        gensym("linksPosSpeedT"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedT,    gensym("linksLengthSpeedT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosT,               gensym("linkPosT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthT,           gensym("linkLengthT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedT,        gensym("linkPosSpeedT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedT,    gensym("linkLengthSpeedT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosXT,               gensym("linkPosXT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthXT,           gensym("linkLengthXT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedXT,        gensym("linkPosSpeedXT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedXT,    gensym("linkLengthSpeedXT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosYT,               gensym("linkPosYT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthYT,           gensym("linkLengthYT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedYT,        gensym("linkPosSpeedYT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedYT,    gensym("linkLengthSpeedYT"), A_GIMME, 0);

 

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosNormL,               gensym("linksPosNormL"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthNormL,           gensym("linksLengthNormL"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedNormL,        gensym("linksPosSpeedNormL"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedNormL,    gensym("linksLengthSpeedNormL"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosNormT,               gensym("linksPosNormT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthNormT,           gensym("linksLengthNormT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedNormT,        gensym("linksPosSpeedNormT"), A_GIMME, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedNormT,    gensym("linksLengthSpeedNormT"), A_GIMME, 0);

-    

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksEndT,               gensym("linksEndT"),   A_GIMME, 0);

-/*    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksEnd1T,               gensym("linksEnd1T"),  A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksEnd2T,               gensym("linksEnd2T"),  A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksEndXT,               gensym("linksEndXT"),  A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksEnd1XT,               gensym("linksEnd1XT"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksEnd2XT,               gensym("linksEnd2XT"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksEndYT,               gensym("linksEndYT"),  A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksEnd1YT,               gensym("linksEnd1YT"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksEnd2YT,               gensym("linksEnd2YT"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksEndZT,               gensym("linksEndZT"),  A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksEnd1ZT,               gensym("linksEnd1ZT"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksEnd2ZT,               gensym("linksEnd2ZT"), A_GIMME, 0);

+

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosNormT,               gensym("linkPosNormT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthNormT,           gensym("linkLengthNormT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedNormT,        gensym("linkPosSpeedNormT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedNormT,    gensym("linkLengthSpeedNormT"), A_GIMME, 0);

+    

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkEndT,               gensym("linkEndT"),   A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkEnd1T,               gensym("linkEnd1T"),  A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkEnd2T,               gensym("linkEnd2T"),  A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkEndXT,               gensym("linkEndXT"),  A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkEnd1XT,               gensym("linkEnd1XT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkEnd2XT,               gensym("linkEnd2XT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkEndYT,               gensym("linkEndYT"),  A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkEnd1YT,               gensym("linkEnd1YT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkEnd2YT,               gensym("linkEnd2YT"), A_GIMME, 0);

+

+/*

+ pmpd2d_test

+ --

+ Functions to list all elements that fit specific conditions

+ Output syntax : depends of the function

+*/

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_testMass,       			gensym("testMass"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_testLink,       			gensym("testLink"), A_GIMME, 0); 

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_testMassT,       		gensym("testMassT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_testLinkT,       		gensym("testLinkT"), A_GIMME, 0); 

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_testMassL,       		gensym("testMassL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_testLinkL,       		gensym("testLinkL"), A_GIMME, 0); 

+

+/*    

+ pmpd2d_stat

+ --

+ Functions to get global stats

 */

 

-/*    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosMean,               gensym("linksPosMean"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthMean,               gensym("linksLengthMean"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedMean,        gensym("linksPosSpeedMean"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedMean,    gensym("linksLengthSpeedMean"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosStd,                   gensym("linksPosStd"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthStd,               gensym("linksLengthStd"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedStd,        gensym("linksPosSpeedStd"), A_GIMME, 0);

-    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedStd0,    gensym("linksLengthSpeedStd"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massPosMean,       gensym("massPosMean"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massPosStd,        gensym("massPosStd"),A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massForcesMean,    gensym("massForcesMean"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massForcesStd,     gensym("massForcesStd"),A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massSpeedsMean,    gensym("massSpeedsMean"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massSpeedsStd,     gensym("massSpeedsStd"),A_GIMME, 0);

+

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkPosMean,               gensym("linkPosMean"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkLengthMean,               gensym("linkLengthMean"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkPosSpeedMean,        gensym("linkPosSpeedMean"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkLengthSpeedMean,    gensym("linkLengthSpeedMean"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkPosStd,                   gensym("linkPosStd"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkLengthStd,               gensym("linkLengthStd"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkPosSpeedStd,        gensym("linkPosSpeedStd"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkLengthSpeedStd,    gensym("linkLengthSpeedStd"), A_GIMME, 0);

+

+/* 	

+ pmpd2d_interactor

+ --

+ Functions to add a global interaction with a specific shape

+*/ 

+ 

+ //   class_addmethod(pmpd2d_class, (t_method)pmpd2d_iCircle,				gensym("iCircle"), A_GIMME, 0);

+ //   class_addmethod(pmpd2d_class, (t_method)pmpd2d_iLine,					gensym("iLine"), A_GIMME, 0);

+

+/*

+ pmpd2d_various

+ --

+ Others

 */

 

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_infosL,          gensym("infosL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_infosL,          gensym("print"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_force,           gensym("force"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_forceX,          gensym("forceX"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_forceY,          gensym("forceY"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_min,             gensym("min"), A_DEFFLOAT, A_DEFFLOAT, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_max,             gensym("max"), A_DEFFLOAT, A_DEFFLOAT, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_minX,            gensym("Xmin"), A_DEFFLOAT, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_maxX,            gensym("Xmax"), A_DEFFLOAT, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_minY,            gensym("Ymin"), A_DEFFLOAT, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_maxY,            gensym("Ymax"), A_DEFFLOAT, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_minX,            gensym("minX"), A_DEFFLOAT, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_maxX,            gensym("maxX"), A_DEFFLOAT, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_minY,            gensym("minY"), A_DEFFLOAT, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_maxY,            gensym("maxY"), A_DEFFLOAT, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_addPos,          gensym("addPos"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_addPosX,         gensym("addPosX"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_addPosY,         gensym("addPosY"), A_GIMME, 0);

+

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_grabMass,        gensym("grabMass"), A_DEFFLOAT, A_DEFFLOAT, A_DEFFLOAT, 0);

     class_addmethod(pmpd2d_class, (t_method)pmpd2d_closestMass,        gensym("closestMass"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massDistances,     	gensym("massDistance"), A_GIMME, 0);

+

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_forcesXT,       			gensym("forceXT"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_forcesYT,       			gensym("forceYT"), A_GIMME, 0);

+

+/*

+ pmpd3d_deprecated

+ --

+ Functions in which the output selector has been modified

+ It is now the same as the input slector : all singular

+*/

+

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosL,      gensym("massesPosL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsL,   gensym("massesSpeedsL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesL,   gensym("massesForcesL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosXL,     gensym("massesPosXL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsXL,  gensym("massesSpeedsXL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesXL,  gensym("massesForcesXL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosYL,     gensym("massesPosYL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsYL,  gensym("massesSpeedsYL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesYL,  gensym("massesForcesYL"), 0);

+

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosNormL,      gensym("massesPosNormL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsNormL,   gensym("massesSpeedsNormL"), 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesNormL,   gensym("massesForcesNormL"), 0);

 

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosL,            gensym("linksPosL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthL,         gensym("linksLengthL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedL,       gensym("linksPosSpeedL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedL,    gensym("linksLengthSpeedL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosXL,           gensym("linksPosXL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthXL,        gensym("linksLengthXL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedXL,      gensym("linksPosSpeedXL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedXL,   gensym("linksLengthSpeedXL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosYL,           gensym("linksPosYL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthYL,        gensym("linksLengthYL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedYL,      gensym("linksPosSpeedYL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedYL,   gensym("linksLengthSpeedYL"), A_GIMME, 0);

+

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosNormL,               gensym("linksPosNormL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthNormL,           gensym("linksLengthNormL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksPosSpeedNormL,        gensym("linksPosSpeedNormL"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_linksLengthSpeedNormL,    gensym("linksLengthSpeedNormL"), A_GIMME, 0);

+

+    

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosMean,       gensym("massesPosMean"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesPosStd,        gensym("massesPosStd"),A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesMean,    gensym("massesForcesMean"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesForcesStd,     gensym("massesForcesStd"),A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsMean,    gensym("massesSpeedsMean"), A_GIMME, 0);

+    class_addmethod(pmpd2d_class, (t_method)pmpd2d_massesSpeedsStd,     gensym("massesSpeedsStd"),A_GIMME, 0);

 }