From e80dc405e548667e8c21b433dfc9030a6c7ed75b Mon Sep 17 00:00:00 2001 From: Cyrille Henry Date: Tue, 26 Jul 2011 19:29:42 +0000 Subject: adding lot's of message to get links informations svn path=/trunk/externals/pmpd/; revision=15166 --- Makefile | 2 +- pmpd-help.pd | 407 +++++++++------ pmpd.c | 290 ++++++++++- pmpd2d-help.pd | 573 +++++++++++++-------- pmpd2d.c | 1239 ++++++++++++++++++++++++++++++++++++++++++--- pmpd3d-help.pd | 629 ++++++++++++++--------- pmpd3d.c | 1548 +++++++++++++++++++++++++++++++++++++++++++++++++++++--- 7 files changed, 3926 insertions(+), 762 deletions(-) diff --git a/Makefile b/Makefile index b1e7d9c..f45387d 100644 --- a/Makefile +++ b/Makefile @@ -114,7 +114,7 @@ ifeq ($(UNAME),Linux) EXTENSION = pd_linux OS = linux PD_PATH = /usr - OPT_CFLAGS = -O6 -funroll-loops -fomit-frame-pointer -ffast-math -fno-trapping-math + OPT_CFLAGS = -O6 -funroll-loops -fomit-frame-pointer -ffast-math -fno-trapping-math -mfpmath=sse CFLAGS += -fPIC LDFLAGS += -Wl,--export-dynamic -shared -fPIC LIBS += -lc diff --git a/pmpd-help.pd b/pmpd-help.pd index 284071f..82b4235 100644 --- a/pmpd-help.pd +++ b/pmpd-help.pd @@ -1,23 +1,23 @@ -#N canvas 553 116 575 730 10; +#N canvas 37 61 575 730 10; #X obj 4 369 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 0; -#X obj 5 669 cnv 15 550 40 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 5 689 cnv 15 550 40 empty empty empty 20 12 0 14 -158509 -66577 0; #X obj 8 6 cnv 15 550 30 empty empty empty 20 12 0 14 -233017 -66577 0; #X obj 14 389 cnv 15 75 15 empty empty empty 20 12 0 14 -233017 -66577 0; -#X obj 14 584 cnv 15 75 15 empty empty empty 20 12 0 14 -233017 -66577 +#X obj 14 604 cnv 15 75 15 empty empty empty 20 12 0 14 -233017 -66577 0; -#X obj 4 527 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 4 547 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 0; -#X obj 4 562 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 4 582 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 0; #X obj 7 76 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 0; -#X obj 5 619 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 5 639 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 0; -#X obj 5 639 cnv 15 550 30 empty empty empty 20 12 0 14 -233017 -66577 +#X obj 5 659 cnv 15 550 30 empty empty empty 20 12 0 14 -233017 -66577 0; #N canvas 76 31 777 741 More_Info 0; #X text 25 65 It is designed to implement particules physical model @@ -176,13 +176,13 @@ masses at the creation.; #X text 26 46 pmpd is the 1D object of the msd objects collection. ; #X connect 2 0 3 0; -#X restore 16 675 pd More_Info; +#X restore 16 695 pd More_Info; #X text 12 76 Examples:; #X text 9 369 Inlets:; #X text 19 388 - Left:; -#X text 10 526 Arguments:; -#X text 11 562 Outlets:; -#X text 19 619 See Also:; +#X text 10 546 Arguments:; +#X text 11 582 Outlets:; +#X text 19 639 See Also:; #X text 74 48 Full Name:; #N canvas 58 22 262 70 Related_Objects 0; #X obj 3 10 cnv 15 250 15 empty empty empty 20 12 0 14 -158509 -66577 @@ -192,7 +192,7 @@ masses at the creation.; #X text 5 10 Externals and libraries; #X obj 44 37 pmpd2d; #X obj 140 37 pmpd3d; -#X restore 122 675 pd Related_Objects; +#X restore 122 695 pd Related_Objects; #N canvas 541 387 524 302 init 0; #X msg 89 187 Xmax 127 \, Xmin 0; #X obj 89 215 t a; @@ -250,149 +250,112 @@ souple 3 2 1 0.5; #X obj 261 152 unpack f f f f; #X text 101 388 Bang - A bang at the left inlet compute the new model state based on previous instant.; -#X text 158 478 To set the model parameters after creation.; -#X text 158 456 To create the model masses and links.; -#X text 158 501 To get the model parameters; -#N canvas 818 245 561 671 creation________ 0; -#X obj 11 95 cnv 15 100 35 empty empty empty 20 12 0 14 -233017 -66577 -0; -#X obj 5 75 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 -0; -#X obj 5 3 cnv 15 550 15 empty empty empty 20 12 0 14 -233017 -66577 -0; -#X text 7 2 CREATION Messages; -#X obj 10 156 cnv 15 150 140 empty empty empty 20 12 0 14 -233017 -66577 -0; -#X obj 4 137 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 -0; -#X obj 10 332 cnv 15 240 220 empty empty empty 20 12 0 14 -233017 -66577 -0; -#X obj 3 308 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 -0; -#X msg 32 104 reset; -#X text 202 167 Add a mass; -#X text 288 168 \$1 : Id (symbol); -#X text 287 188 \$2 : fixed or mobile (0/1); -#X text 287 206 \$3 : mass; -#X msg 32 167 mass \$1 \$2 \$3 \$4; -#X text 287 224 \$4 : initial position; -#X text 7 137 Masses :; -#X text 7 74 Reset :; -#X text 129 105 Delete all masses \, links and internal variables; -#X text 6 308 Links :; -#X text 271 360 \$1 : Id (symbol); -#X text 271 360 \$1 : Id (symbol); -#X text 271 413 \$4 : rigidity; -#X text 271 378 \$2 : creation No/Id of mass1; -#X text 271 396 \$3 : creation No/Id of mass2; -#X text 266 341 Add link(s); -#X text 271 431 \$5 : damping; -#X text 10 28 Creation messages are used to define the structure of -the model. Messages create links and masses or destroy them.; -#X msg 30 341 link \$1 \$2 \$3 \$4 \$5 (\$6 \$7 \$8); -#X text 271 449 \$6 : Power of the rigidity distance; -#X msg 35 519 tabLink \$1 \$2 \$3 \$4 \$5 \$6 \$7; -#X text 267 468 (\$7) : minimum lenght of link; -#X text 268 485 (\$8) : maximum lenght of link; -#X text 266 535 \$1 : Id (symbol); -#X text 266 535 \$1 : Id (symbol); -#X text 266 553 \$2 : creation No/Id of mass1; -#X text 266 571 \$3 : creation No/Id of mass2; -#X text 261 516 Add non linear link(s); -#X text 264 588 \$4 : table name for the rigidity; -#X text 264 606 \$5 : length of the table; -#X text 264 624 \$6 : table name for the damping; -#X text 263 643 \$7 : length of the table; -#X restore 12 457 pd creation________; -#X text 103 542 None; -#X text 18 583 - Left:; -#X text 101 584 Outputs the model parameters asked with the attributes +#X text 188 478 To set the model parameters after creation.; +#X text 188 456 To create the model masses and links.; +#X text 186 512 To get the model parameters; +#X text 103 562 None; +#X text 18 603 - Left:; +#X text 101 604 Outputs the model parameters asked with the attributes messages.; #X text 101 420 Messages - Different messages are used to control the msd object. They are of three types :; -#X text 9 641 CATEGORY: control; -#N canvas 841 74 1139 731 dynamic 0; -#X obj 5 3 cnv 15 550 15 empty empty empty 20 12 0 14 -233017 -66577 +#X text 9 661 CATEGORY: control; +#X text 9 671 KEYWORDS: physical model mass spring damper link; +#X text 111 639 01_msdtest.pd; +#X obj 157 48 pmpd; +#X text 12 8 HELP: pmpd; +#X text 12 20 DESCRIPTION: physical modeling for pd \, based on msd. +; +#X obj 17 319 pmpd; +#X text 266 696 - Cyrille Henry \, 2011; +#X text 283 712 based on MSD by Nicolas Montgermont; +#X obj 17 187 loadbang; +#N canvas 446 66 1236 814 links_attributes 0; +#X obj 5 3 cnv 15 590 15 empty empty empty 20 12 0 14 -233017 -66577 0; -#X obj 10 85 cnv 15 130 580 empty empty empty 20 12 0 14 -233017 -66577 +#X obj 10 108 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 0; -#X obj 4 62 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 4 87 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 0; -#X obj 575 86 cnv 15 130 380 empty empty empty 20 12 0 14 -233017 -66577 +#X text 7 2 ATTRIBUTES Messages; +#X text 7 87 Lists :; +#X obj 11 222 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 0; -#X obj 568 62 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 11 335 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 0; -#X text 7 62 Masses :; -#X text 571 62 Links :; -#X text 7 2 DYNAMIC SETTINGS Messages; -#X msg 31 96 posX \$1 \$2; -#X text 190 144 Add force on mass(es); -#X msg 30 152 forceX \$1 \$2; -#X text 192 212 \$1 : Value; -#X text 193 107 \$1 : Id (symbol) or No; -#X text 193 161 \$1 : Id (symbol) or No; -#X msg 30 195 Xmin \$1; -#X msg 30 217 Xmax \$1; -#X msg 29 244 setMobile \$1; -#X msg 29 265 setFixed \$1; -#X text 193 89 Set position of fixed mass(es); -#X text 193 125 \$2 : Value; -#X text 193 179 \$2 : Value; -#X text 189 196 Set minimimum and maximum position of all masses; -#X text 188 244 Set mass to mobile or fixed; -#X msg 594 150 setD \$1 \$2; -#X text 749 136 \$2 : New value; -#X msg 594 99 setK \$1 \$2; -#X text 749 192 \$2 : New value; -#X text 749 245 \$2 : New value; -#X text 743 100 Set rigidity of link(s); -#X text 743 156 Set damping of link(s); -#X msg 594 205 setL \$1 \$2; -#X text 743 209 Set initial lenght of link(s); -#X text 749 118 \$1 : Id (symbol) or No; -#X text 749 174 \$1 : Id (symbol) or No; -#X text 749 227 \$1 : Id (symbol) or No; -#X text 191 261 \$1 : Id (symbol) or No; -#X text 10 25 Dynamic settings messages allows the user to redefine -internal parameters of links and masses.; -#X msg 28 286 setDEnv \$1; -#X text 184 287 Change environement damping; -#X text 196 304 \$1 : damping; -#X msg 588 280 setLinkId \$1 \$2; -#X text 742 281 Change the Id of a link; -#X text 743 298 \$1 : link (number or Id); -#X text 746 315 \$2 : New Id; -#X msg 591 353 setLKTab \$1 \$2; -#X text 750 353 change the table (K) size of a non linear link; -#X text 751 373 \$1 : link (number or Id); -#X text 752 390 \$2 : New KL; -#X text 751 435 \$1 : link (number or Id); -#X msg 591 415 setLDTab \$1 \$2; -#X text 750 415 change the table (D) size of a non linear link; -#X text 752 452 \$2 : New DL; -#X text 191 373 \$2 : New Id; -#X msg 25 339 setMassId \$1 \$2; -#X text 191 339 Change the Id of a mass; -#X text 192 356 \$1 : mass (number or Id); -#X msg 23 400 setForce \$1 \$2; -#X msg 23 468 setSpeed \$1 \$2; -#X msg 24 531 addPos \$1 \$2; -#X text 186 420 \$1 : Id (symbol) or No; -#X text 186 438 \$2 : Value; -#X text 183 403 set force on mass(es); -#X text 185 487 \$1 : Id (symbol) or No; -#X text 185 505 \$2 : Value; -#X text 182 470 set speed of a mass(es); -#X text 184 548 \$1 : Id (symbol) or No; -#X text 184 566 \$2 : Value; -#X text 181 531 add pos of a mass(es); -#X msg 26 593 grabMass \$1 \$2 \$3; -#X text 183 592 grab a mass \, and move it; -#X text 185 614 \$1 : mass number; -#X text 185 630 \$2 : position; -#X text 185 646 \$3 : grab or not (0 or 1); -#X restore 12 478 pd dynamic settings; -#N canvas 888 196 1230 682 attributes______ 0; +#X obj 611 26 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 605 5 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X obj 612 140 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 +-66577 0; +#X obj 612 253 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 +-66577 0; +#X text 369 530 \$1 (optional): Masses Id; +#X text 834 500 \$2 (optional): Masses Id; +#X text 833 487 \$1 : table name; +#X obj 17 584 cnv 15 150 220 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 11 564 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X msg 22 597 massesPosMean; +#X text 177 613 \$1 : none or Id; +#X msg 23 629 massesPosStd \$1; +#X text 179 644 \$1 : none or Id; +#X text 173 597 output the average position of masses; +#X text 174 628 output the standard deviation of masses position; +#X text 177 674 \$1 : none or Id; +#X text 179 711 \$1 : none or Id; +#X text 173 658 output the average position of masses; +#X text 174 695 output the standard deviation of masses position; +#X text 177 749 \$1 : none or Id; +#X text 178 786 \$1 : none or Id; +#X text 173 733 output the average position of masses; +#X text 173 770 output the standard deviation of masses position; +#X msg 22 696 massesForcesStd \$1; +#X msg 22 658 massesForcesMean; +#X msg 23 735 massesSpeedsMean; +#X msg 22 771 massesSpeedsStd \$1; +#X text 608 5 tables :; +#X text 15 30 The attributes messages ask the object to output some +of his internal parameters. They can be output by lists for positions +and forces of links; +#X text 14 564 statistics (not implemented yet) :; +#X msg 11 119 linksPosL (\$1); +#X text 170 120 Output all links center positions in a list on outlet +No 1; +#X msg 11 138 linksPosXL (\$1); +#X msg 13 229 linksLengthL (\$1); +#X msg 13 248 linksLengthXL (\$1); +#X obj 12 451 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X msg 14 458 linksLengthSpeedL (\$1); +#X msg 14 477 linksLengthSpeedXL (\$1); +#X text 184 459 Output all link lenght deformation speed in a list +on outlet No 1; +#X msg 615 38 linksPosT \$1 (\$2); +#X msg 616 57 linksPosXT \$1 (\$2); +#X msg 616 152 linksLengthT \$1 (\$2); +#X obj 612 369 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 +-66577 0; +#X text 808 377 Output all link lenght deformation speed in a list +on outlet No 1; +#X msg 614 376 linksLengthSpeedT \$1 (\$2); +#X msg 613 395 linksLengthSpeedXT \$1 (\$2); +#X text 793 47 Output all links center positions in a table (\$1); +#X text 170 233 Output all link lenght in a list on outlet No 1; +#X text 171 349 Output all link center speed in a list on outlet No +1; +#X text 797 154 Output all link length in a list on outlet No 1; +#X msg 16 342 linksPosSpeedL (\$1); +#X msg 16 361 linksPosSpeedXL (\$1); +#X msg 616 259 linksPosSpeedT \$1 (\$2); +#X text 795 264 Output all link center speed in a list on outlet No +1; +#X msg 616 171 linksLengthXT \$1 (\$2); +#X msg 616 278 linksPosSpeedXT \$1 (\$2); +#X restore 12 519 pd links_attributes ___; +#N canvas 888 196 1230 682 attributes__________ 0; #X obj 11 95 cnv 15 100 35 empty empty empty 20 12 0 14 -233017 -66577 0; #X obj 5 75 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 @@ -489,18 +452,140 @@ name \, \$2 (optional): Masses Id); #X text 339 645 [massesPosName mas x2(; #X text 337 581 [massesPosId 1 x1(; #X text 337 596 [massesPosId 2 x2(; -#X restore 12 499 pd attributes______; -#X text 9 651 KEYWORDS: physical model mass spring damper link; -#X text 111 619 01_msdtest.pd; -#X obj 157 48 pmpd; -#X text 12 8 HELP: pmpd; -#X text 12 20 DESCRIPTION: physical modeling for pd \, based on msd. -; -#X obj 17 319 pmpd; -#X text 266 676 - Cyrille Henry \, 2011; -#X text 283 692 based on MSD by Nicolas Montgermont; -#X obj 17 187 loadbang; -#X connect 21 0 51 0; +#X restore 12 499 pd attributes__________; +#N canvas 841 74 1139 731 dynamic 0; +#X obj 5 3 cnv 15 550 15 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 10 85 cnv 15 130 580 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 4 62 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X obj 575 86 cnv 15 130 380 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 568 62 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X text 7 62 Masses :; +#X text 571 62 Links :; +#X text 7 2 DYNAMIC SETTINGS Messages; +#X msg 31 96 posX \$1 \$2; +#X text 190 144 Add force on mass(es); +#X msg 30 152 forceX \$1 \$2; +#X text 192 212 \$1 : Value; +#X text 193 107 \$1 : Id (symbol) or No; +#X text 193 161 \$1 : Id (symbol) or No; +#X msg 30 195 Xmin \$1; +#X msg 30 217 Xmax \$1; +#X msg 29 244 setMobile \$1; +#X msg 29 265 setFixed \$1; +#X text 193 89 Set position of fixed mass(es); +#X text 193 125 \$2 : Value; +#X text 193 179 \$2 : Value; +#X text 189 196 Set minimimum and maximum position of all masses; +#X text 188 244 Set mass to mobile or fixed; +#X msg 594 150 setD \$1 \$2; +#X text 749 136 \$2 : New value; +#X msg 594 99 setK \$1 \$2; +#X text 749 192 \$2 : New value; +#X text 749 245 \$2 : New value; +#X text 743 100 Set rigidity of link(s); +#X text 743 156 Set damping of link(s); +#X msg 594 205 setL \$1 \$2; +#X text 743 209 Set initial lenght of link(s); +#X text 749 118 \$1 : Id (symbol) or No; +#X text 749 174 \$1 : Id (symbol) or No; +#X text 749 227 \$1 : Id (symbol) or No; +#X text 191 261 \$1 : Id (symbol) or No; +#X text 10 25 Dynamic settings messages allows the user to redefine +internal parameters of links and masses.; +#X msg 28 286 setDEnv \$1; +#X text 184 287 Change environement damping; +#X text 196 304 \$1 : damping; +#X msg 588 280 setLinkId \$1 \$2; +#X text 742 281 Change the Id of a link; +#X text 743 298 \$1 : link (number or Id); +#X text 746 315 \$2 : New Id; +#X msg 591 353 setLKTab \$1 \$2; +#X text 750 353 change the table (K) size of a non linear link; +#X text 751 373 \$1 : link (number or Id); +#X text 752 390 \$2 : New KL; +#X text 751 435 \$1 : link (number or Id); +#X msg 591 415 setLDTab \$1 \$2; +#X text 750 415 change the table (D) size of a non linear link; +#X text 752 452 \$2 : New DL; +#X text 191 373 \$2 : New Id; +#X msg 25 339 setMassId \$1 \$2; +#X text 191 339 Change the Id of a mass; +#X text 192 356 \$1 : mass (number or Id); +#X msg 23 400 setForce \$1 \$2; +#X msg 23 468 setSpeed \$1 \$2; +#X msg 24 531 addPos \$1 \$2; +#X text 186 420 \$1 : Id (symbol) or No; +#X text 186 438 \$2 : Value; +#X text 183 403 set force on mass(es); +#X text 185 487 \$1 : Id (symbol) or No; +#X text 185 505 \$2 : Value; +#X text 182 470 set speed of a mass(es); +#X text 184 548 \$1 : Id (symbol) or No; +#X text 184 566 \$2 : Value; +#X text 181 531 add pos of a mass(es); +#X msg 26 593 grabMass \$1 \$2 \$3; +#X text 183 592 grab a mass \, and move it; +#X text 185 614 \$1 : mass number; +#X text 185 630 \$2 : position; +#X text 185 646 \$3 : grab or not (0 or 1); +#X restore 12 478 pd dynamic settings ___; +#N canvas 818 245 561 671 creation____________ 0; +#X obj 11 95 cnv 15 100 35 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 5 75 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X obj 5 3 cnv 15 550 15 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X text 7 2 CREATION Messages; +#X obj 10 156 cnv 15 150 140 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 4 137 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X obj 10 332 cnv 15 240 220 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 3 308 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X msg 32 104 reset; +#X text 202 167 Add a mass; +#X text 288 168 \$1 : Id (symbol); +#X text 287 188 \$2 : fixed or mobile (0/1); +#X text 287 206 \$3 : mass; +#X msg 32 167 mass \$1 \$2 \$3 \$4; +#X text 287 224 \$4 : initial position; +#X text 7 137 Masses :; +#X text 7 74 Reset :; +#X text 129 105 Delete all masses \, links and internal variables; +#X text 6 308 Links :; +#X text 271 360 \$1 : Id (symbol); +#X text 271 360 \$1 : Id (symbol); +#X text 271 413 \$4 : rigidity; +#X text 271 378 \$2 : creation No/Id of mass1; +#X text 271 396 \$3 : creation No/Id of mass2; +#X text 266 341 Add link(s); +#X text 271 431 \$5 : damping; +#X text 10 28 Creation messages are used to define the structure of +the model. Messages create links and masses or destroy them.; +#X msg 30 341 link \$1 \$2 \$3 \$4 \$5 (\$6 \$7 \$8); +#X text 271 449 \$6 : Power of the rigidity distance; +#X msg 35 519 tabLink \$1 \$2 \$3 \$4 \$5 \$6 \$7; +#X text 267 468 (\$7) : minimum lenght of link; +#X text 268 485 (\$8) : maximum lenght of link; +#X text 266 535 \$1 : Id (symbol); +#X text 266 535 \$1 : Id (symbol); +#X text 266 553 \$2 : creation No/Id of mass1; +#X text 266 571 \$3 : creation No/Id of mass2; +#X text 261 516 Add non linear link(s); +#X text 264 588 \$4 : table name for the rigidity; +#X text 264 606 \$5 : length of the table; +#X text 264 624 \$6 : table name for the damping; +#X text 263 643 \$7 : length of the table; +#X restore 12 457 pd creation____________; +#X connect 21 0 48 0; #X connect 22 0 25 0; #X connect 25 0 24 0; #X connect 27 0 26 0; @@ -511,5 +596,5 @@ name \, \$2 (optional): Masses Id); #X connect 33 1 23 0; #X connect 33 2 31 0; #X connect 33 3 32 0; -#X connect 51 0 20 0; -#X connect 54 0 27 0; +#X connect 48 0 20 0; +#X connect 51 0 27 0; diff --git a/pmpd.c b/pmpd.c index 74dfc7f..8e49da8 100644 --- a/pmpd.c +++ b/pmpd.c @@ -31,8 +31,8 @@ #include "stdio.h" #include "math.h" -#define nb_max_link 100000 -#define nb_max_mass 100000 +#define nb_max_link 10000000 +#define nb_max_mass 1000000 #define max(a,b) ( ((a) > (b)) ? (a) : (b) ) #define min(a,b) ( ((a) < (b)) ? (a) : (b) ) @@ -293,7 +293,8 @@ void pmpd_link(t_pmpd *x, t_symbol *s, int argc, t_atom *argv) { 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 + 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 pmpd_create_link(x, Id, i, j, K, D, Pow, Lmin, Lmax, 0); } } @@ -362,7 +363,8 @@ void pmpd_tabLink(t_pmpd *x, t_symbol *s, int argc, t_atom *argv) { 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 + 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 { pmpd_create_link(x, Id, i, j, 1, 1, 1, -1000000, 1000000, 2); x->link[x->nb_link-1].arrayK = arrayK; @@ -1484,6 +1486,258 @@ void pmpd_massesSpeedsStd(t_pmpd *x, t_symbol *s, int argc, t_atom *argv) outlet_anything(x->main_outlet, gensym("massesSpeedsStd"),1 , std_out); } +// -------------------------------------------- + +void pmpd_linksPosL(t_pmpd *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 pmpd_linksLengthL(t_pmpd *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 pmpd_linksPosSpeedL(t_pmpd *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 pmpd_linksLengthSpeedL(t_pmpd *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 pmpd_linksPosT(t_pmpd *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 pmpd_linksLengthT(t_pmpd *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 pmpd_linksPosSpeedT(t_pmpd *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 pmpd_linksLengthSpeedT(t_pmpd *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 pmpd_grabMass(t_pmpd *x, t_float posX, t_float grab) { t_float dist, tmp; @@ -1586,6 +1840,34 @@ void pmpd_setup(void) class_addmethod(pmpd_class, (t_method)pmpd_massesSpeedsMean,gensym("massesSpeedsMean"), A_GIMME, 0); class_addmethod(pmpd_class, (t_method)pmpd_massesSpeedsStd, gensym("massesSpeedsStd"),A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksPosL, gensym("linksPosXL"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksLengthL, gensym("linksLengthXL"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksPosSpeedL, gensym("linksPosSpeedXL"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksLengthSpeedL, gensym("linksLengthXL"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksPosT, gensym("linksPosXT"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksLengthT, gensym("linksLengthXT"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksPosSpeedT, gensym("linksPosSpeedXT"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksLengthSpeedT, gensym("linksLengthSpeedXT"), A_GIMME, 0); + + class_addmethod(pmpd_class, (t_method)pmpd_linksPosL, gensym("linksPosL"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksLengthL, gensym("linksLengthL"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksPosSpeedL, gensym("linksPosSpeedL"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksLengthSpeedL, gensym("linksLengthL"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksPosT, gensym("linksPosT"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksLengthT, gensym("linksLengthT"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksPosSpeedT, gensym("linksPosSpeedT"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksLengthSpeedT, gensym("linksLengthSpeedT"), A_GIMME, 0); + +/* class_addmethod(pmpd_class, (t_method)pmpd_linksPosMean, gensym("linksPosMean"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksLengthMean, gensym("linksLengthMean"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksPosSpeedMean, gensym("linksPosSpeedMean"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksLengthSpeedMean, gensym("linksLengthSpeedMean"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksPosStd, gensym("linksPosStd"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksLengthStd, gensym("linksLengthStd"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksPosSpeedStd, gensym("linksPosSpeedStd"), A_GIMME, 0); + class_addmethod(pmpd_class, (t_method)pmpd_linksLengthSpeedStd0, gensym("linksLengthSpeedStd"), A_GIMME, 0); +*/ + class_addmethod(pmpd_class, (t_method)pmpd_grabMass, gensym("grabMass"), A_DEFFLOAT, A_DEFFLOAT, 0); } diff --git a/pmpd2d-help.pd b/pmpd2d-help.pd index 3fcd199..937a5b4 100644 --- a/pmpd2d-help.pd +++ b/pmpd2d-help.pd @@ -1,23 +1,23 @@ -#N canvas 150 72 568 710 10; +#N canvas 150 72 571 767 10; #X obj 4 369 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 0; -#X obj 5 665 cnv 15 550 30 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 5 685 cnv 15 550 30 empty empty empty 20 12 0 14 -158509 -66577 0; #X obj 8 6 cnv 15 550 30 empty empty empty 20 12 0 14 -233017 -66577 0; #X obj 14 389 cnv 15 75 15 empty empty empty 20 12 0 14 -233017 -66577 0; -#X obj 14 584 cnv 15 75 15 empty empty empty 20 12 0 14 -233017 -66577 +#X obj 14 604 cnv 15 75 15 empty empty empty 20 12 0 14 -233017 -66577 0; -#X obj 4 527 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 4 547 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 0; -#X obj 4 562 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 4 582 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 0; #X obj 7 76 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 0; -#X obj 5 615 cnv 15 550 20 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 5 635 cnv 15 550 20 empty empty empty 20 12 0 14 -158509 -66577 0; -#X obj 5 635 cnv 15 550 30 empty empty empty 20 12 0 14 -233017 -66577 +#X obj 5 655 cnv 15 550 30 empty empty empty 20 12 0 14 -233017 -66577 0; #N canvas 76 31 921 714 More_Info 0; #X text 96 12 pmpd : mass - spring - damper model; @@ -226,13 +226,13 @@ their Id.; #X text 28 41 pmpd is the 2D object of the pmpd objects collection. ; #X connect 3 0 4 0; -#X restore 16 673 pd More_Info; +#X restore 16 693 pd More_Info; #X text 12 76 Examples:; #X text 9 369 Inlets:; #X text 19 388 - Left:; -#X text 10 526 Arguments:; -#X text 11 562 Outlets:; -#X text 19 617 See Also:; +#X text 10 546 Arguments:; +#X text 11 582 Outlets:; +#X text 19 637 See Also:; #X text 74 48 Full Name:; #N canvas 58 22 262 70 Related_Objects 0; #X obj 3 10 cnv 15 250 15 empty empty empty 20 12 0 14 -158509 -66577 @@ -242,7 +242,7 @@ their Id.; #X text 5 10 Externals and libraries; #X obj 44 37 pmpd; #X obj 141 39 pmpd3d; -#X restore 122 673 pd Related_Objects; +#X restore 122 693 pd Related_Objects; #N canvas 257 262 759 345 init 0; #X obj 89 215 t a; #X obj 89 33 loadbang; @@ -292,10 +292,232 @@ their Id.; -1 -1; #X text 101 388 Bang - A bang at the left inlet compute the new model state based on previous instant.; -#X text 158 478 To set the model parameters after creation.; -#X text 158 456 To create the model masses and links.; -#X text 158 501 To get the model parameters; -#N canvas 41 48 967 626 creation________ 0; +#X text 182 479 To set the model parameters after creation.; +#X text 182 457 To create the model masses and links.; +#X text 183 509 To get the model parameters; +#X text 103 562 None; +#X text 18 603 - Left:; +#X text 101 604 Outputs the model parameters asked with the attributes +messages.; +#X text 101 420 Messages - Different messages are used to control the +pmpd object. They are of three types :; +#X text 9 659 CATEGORY: control; +#X text 9 669 KEYWORDS: physical model mass spring damper link; +#X text 12 8 HELP: pmpd2d; +#X obj 157 48 pmpd2d; +#X text 112 638 01_pmpd2dtest.pd; +#X obj 18 296 pmpd2d; +#N canvas 731 296 450 300 gemmouse 0; +#X obj 189 77 gemmouse; +#X obj 189 184 pack f f; +#X obj 189 218 spigot; +#X obj 109 243 outlet; +#X obj 109 36 inlet; +#X obj 109 98 t b; +#X obj 238 131 + 4; +#X msg 57 182 posX fix \$1; +#X msg 57 203 posY fix \$2; +#X obj 189 108 / 62.5; +#X obj 237 107 / -62.5; +#X obj 189 131 - 4; +#X obj 109 125 list 0 0; +#X connect 0 0 9 0; +#X connect 0 1 10 0; +#X connect 0 2 2 1; +#X connect 1 0 2 0; +#X connect 2 0 12 1; +#X connect 4 0 5 0; +#X connect 5 0 12 0; +#X connect 6 0 1 1; +#X connect 7 0 3 0; +#X connect 8 0 3 0; +#X connect 9 0 11 0; +#X connect 10 0 6 0; +#X connect 11 0 1 0; +#X connect 12 0 7 0; +#X connect 12 0 8 0; +#X restore 393 295 pd gemmouse; +#X obj 393 274 gemhead; +#X obj 346 216 gemwin; +#X msg 357 193 0 \, destroy; +#N canvas 472 258 550 319 gemrender 0; +#X obj 48 203 translateXYZ; +#X obj 48 229 sphere 0.1; +#X obj 127 24 inlet; +#X obj 360 32 inlet; +#X obj 275 232 curve 2; +#X msg 431 81 \$4 \$5 0; +#X msg 359 77 \$2 \$3 0; +#X obj 359 105 t b a; +#X obj 127 62 unpack f f f; +#X obj 166 88 t b f; +#X msg 48 110 0; +#X obj 48 77 loadbang; +#X obj 48 137 gemhead; +#X msg 275 160 0; +#X obj 275 127 loadbang; +#X obj 275 187 gemhead; +#X connect 0 0 1 0; +#X connect 2 0 8 0; +#X connect 3 0 5 0; +#X connect 3 0 6 0; +#X connect 5 0 4 2; +#X connect 6 0 7 0; +#X connect 7 0 15 0; +#X connect 7 1 4 1; +#X connect 8 1 9 0; +#X connect 8 2 0 2; +#X connect 9 0 12 0; +#X connect 9 1 0 1; +#X connect 10 0 12 0; +#X connect 11 0 10 0; +#X connect 12 0 0 0; +#X connect 13 0 15 0; +#X connect 14 0 13 0; +#X connect 15 0 4 0; +#X restore 172 320 pd gemrender; +#X msg 346 129 reset \, create \, 1; +#X text 135 98 Sorry \, you need GEM for this example...; +#X text 169 132 1 Create window -->; +#X text 169 163 2 Drag the structure with the mouse; +#X text 169 195 3 Destroy the window -->; +#X obj 172 294 route massesPos linksPosNo; +#X text 11 19 DESCRIPTION: Physical modeling for Pd in 2d; +#X text 265 689 - Cyrille Henry 2011; +#X text 265 701 based on MSD by Nicolas Montgermont \, 2005; +#N canvas 446 66 1236 814 links_attributes 0; +#X obj 5 3 cnv 15 590 15 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 10 108 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 4 87 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X text 7 2 ATTRIBUTES Messages; +#X text 7 87 Lists :; +#X text 169 167 Output all norm of position vector (sqrt(x*x+y*y)) +in a list on outlet No 1; +#X obj 11 222 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 11 335 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 611 26 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 605 5 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X obj 612 140 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 +-66577 0; +#X obj 612 253 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 +-66577 0; +#X msg 13 286 massesSpeedsNormL (\$1); +#X text 402 546 \$1 (optional): Masses Id; +#X text 850 488 \$2 (optional): Masses Id; +#X text 849 475 \$1 : table name; +#X obj 17 584 cnv 15 150 220 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 11 564 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X msg 22 597 massesPosMean; +#X text 177 613 \$1 : none or Id; +#X msg 23 629 massesPosStd \$1; +#X text 179 644 \$1 : none or Id; +#X text 173 597 output the average position of masses; +#X text 174 628 output the standard deviation of masses position; +#X text 177 674 \$1 : none or Id; +#X text 179 711 \$1 : none or Id; +#X text 173 658 output the average position of masses; +#X text 174 695 output the standard deviation of masses position; +#X text 177 749 \$1 : none or Id; +#X text 178 786 \$1 : none or Id; +#X text 173 733 output the average position of masses; +#X text 173 770 output the standard deviation of masses position; +#X msg 22 696 massesForcesStd \$1; +#X msg 22 658 massesForcesMean; +#X msg 23 735 massesSpeedsMean; +#X msg 22 771 massesSpeedsStd \$1; +#X text 608 5 tables :; +#X text 15 30 The attributes messages ask the object to output some +of his internal parameters. They can be output by lists for positions +and forces of links; +#X text 14 564 statistics (not implemented yet) :; +#X msg 11 119 linksPosL (\$1); +#X text 170 120 Output all links center positions in a list on outlet +No 1; +#X msg 11 138 linksPosXL (\$1); +#X msg 12 157 linksPosYL (\$1); +#X msg 12 176 linksPosNormL (\$1); +#X msg 13 229 linksLengthL (\$1); +#X msg 13 248 linksLengthXL (\$1); +#X msg 13 267 linksLengthYL (\$1); +#X text 170 144 Output all links x or y in a list on outlet No 1; +#X obj 12 451 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X msg 14 458 linksLengthSpeedL (\$1); +#X msg 14 477 linksLengthSpeedXL (\$1); +#X msg 14 496 linksLengthSpeedYL (\$1); +#X msg 14 515 massesSpeedsSpeedNormL (\$1); +#X text 184 459 Output all link lenght deformation speed in a list +on outlet No 1; +#X text 184 485 Output all link lenght deformation speed in x or y +in a list on outlet No 1; +#X text 215 515 Output all norm of link lenght deformation speed vector +(sqrt(x*x+y*y)) in a list on outlet No 1; +#X msg 615 38 linksPosT \$1 (\$2); +#X msg 616 57 linksPosXT \$1 (\$2); +#X msg 615 76 linksPosYT \$1 (\$2); +#X msg 615 96 linksPosNormT \$1 (\$2); +#X msg 616 152 linksLengthT \$1 (\$2); +#X msg 616 210 linksLengthNormT \$1 (\$2); +#X obj 612 369 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 +-66577 0; +#X text 808 377 Output all link lenght deformation speed in a list +on outlet No 1; +#X text 806 401 Output all link lenght deformation speed in x or y +in a list on outlet No 1; +#X text 839 430 Output all norm of link lenght deformation speed vector +(sqrt(x*x+y*y)) in a list on outlet No 1; +#X msg 614 376 linksLengthSpeedT \$1 (\$2); +#X msg 613 395 linksLengthSpeedXT \$1 (\$2); +#X msg 614 414 linksLengthSpeedYT \$1 (\$2); +#X msg 614 434 massesSpeedsSpeedNormT \$1 (\$2); +#X text 794 35 Output all links center positions in a table (\$1); +#X text 794 67 Output all links center positions x or y in a table +(\$1); +#X text 789 90 Output all norm of links center positions vector (sqrt(x*x+y*y)) +in a table (\$1); +#X text 170 233 Output all link lenght in a list on outlet No 1; +#X text 171 255 Output all link lenght in x or y in a list on outlet +No 1; +#X text 181 280 Output all norm of link lenght vector (sqrt(x*x+y*y)) +in a list on outlet No 1; +#X text 171 349 Output all link center speed in a list on outlet No +1; +#X text 175 363 Output all link center speed in x or y in a list on +outlet No 1; +#X text 192 393 Output all norm of link center speed vector (sqrt(x*x+y*y)) +in a list on outlet No 1; +#X text 797 154 Output all link length in a list on outlet No 1; +#X text 796 181 Output all link length in x or y in a list on outlet +No 1; +#X text 799 201 Output all norm of link length vector (sqrt(x*x+y*y)) +in a list on outlet No 1; +#X msg 16 342 linksPosSpeedL (\$1); +#X msg 16 361 linksPosSpeedXL (\$1); +#X msg 15 380 linksPosSpeedYL (\$1); +#X msg 15 399 linksPosSpeedNormL (\$1); +#X msg 616 259 linksPosSpeedT \$1 (\$2); +#X msg 616 316 linksPosSpeedNormT \$1 (\$2); +#X text 795 264 Output all link center speed in a list on outlet No +1; +#X text 798 282 Output all link center speed in x or y in a list on +outlet No 1; +#X text 813 313 Output all norm of link center speed vector (sqrt(x*x+y*y)) +in a list on outlet No 1; +#X msg 616 171 linksLengthXT \$1 (\$2); +#X msg 616 190 linksLengthYT \$1 (\$2); +#X msg 616 278 linksPosSpeedXT \$1 (\$2); +#X msg 616 297 linksPosSpeedYT \$1 (\$2); +#X restore 12 519 pd links_attributes ___; +#N canvas 41 48 967 626 creation____________ 0; #X obj 5 75 cnv 15 450 15 empty empty empty 20 12 0 14 -158509 -66577 0; #X obj 5 3 cnv 15 450 15 empty empty empty 20 12 0 14 -233017 -66577 @@ -367,114 +589,8 @@ following a vector.; 0; #X text 497 351 Table links :; #X connect 29 0 31 0; -#X restore 12 457 pd creation________; -#X text 103 542 None; -#X text 18 583 - Left:; -#X text 101 584 Outputs the model parameters asked with the attributes -messages.; -#X text 101 420 Messages - Different messages are used to control the -pmpd object. They are of three types :; -#X text 9 639 CATEGORY: control; -#N canvas 100 72 1078 872 dynamic 0; -#X obj 5 3 cnv 15 550 15 empty empty empty 20 12 0 14 -233017 -66577 -0; -#X obj 10 85 cnv 15 150 750 empty empty empty 20 12 0 14 -233017 -66577 -0; -#X obj 4 62 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 -0; -#X obj 574 86 cnv 15 130 350 empty empty empty 20 12 0 14 -233017 -66577 -0; -#X obj 567 62 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 -0; -#X text 7 62 Masses :; -#X text 570 62 Links :; -#X text 7 2 DYNAMIC SETTINGS Messages; -#X msg 31 116 posX \$1 \$2; -#X text 191 178 Add force on mass(es); -#X msg 30 193 forceX \$1 \$2; -#X text 190 114 \$1 : Id (symbol) or No; -#X text 193 195 \$1 : Id (symbol) or No; -#X msg 31 285 Xmin \$1; -#X msg 90 285 Xmax \$1; -#X msg 29 334 setMobile \$1; -#X msg 29 355 setFixed \$1; -#X text 190 96 Set position of fixed mass(es); -#X text 177 268 Set minimimum and maximum position of all masses; -#X text 188 334 Set mass to mobile or fixed; -#X msg 593 150 setD \$1 \$2; -#X text 748 136 \$2 : New value; -#X msg 593 99 setK \$1 \$2; -#X text 748 192 \$2 : New value; -#X text 748 245 \$2 : New value; -#X text 742 100 Set rigidity of link(s); -#X text 742 156 Set damping of link(s); -#X msg 593 205 setL \$1 \$2; -#X text 742 209 Set initial lenght of link(s); -#X text 748 118 \$1 : Id (symbol) or No; -#X text 748 174 \$1 : Id (symbol) or No; -#X text 748 227 \$1 : Id (symbol) or No; -#X text 191 351 \$1 : Id (symbol) or No; -#X text 10 25 Dynamic settings messages allows the user to redefine -internal parameters of links and masses.; -#X msg 29 391 grabMass \$1 \$2 \$3; -#X text 186 391 Grab nearest mass; -#X text 191 407 \$1 \, \$2 : position; -#X text 190 424 \$3 : grab or not (0/1); -#X msg 31 137 posY \$1 \$2; -#X msg 30 214 forceY \$1 \$2; -#X msg 31 307 Ymin \$1; -#X msg 90 307 Ymax \$1; -#X text 191 488 \$2 : New Id; -#X msg 30 456 setMassId \$1 \$2; -#X text 191 454 Change the Id of a mass; -#X text 192 471 \$1 : mass (number or Id); -#X msg 23 583 setSpeed \$1 \$2; -#X text 186 535 \$1 : Id (symbol) or No; -#X text 183 518 set force on mass(es); -#X text 185 602 \$1 : Id (symbol) or No; -#X text 182 585 set speed of a mass(es); -#X text 184 663 \$1 : Id (symbol) or No; -#X text 181 646 add pos of a mass(es); -#X msg 23 737 grabMass \$1 \$2 \$3; -#X text 183 707 grab a mass \, and move it; -#X text 185 729 \$1 : mass number; -#X text 185 745 \$2 : position; -#X text 185 761 \$3 : grab or not (0 or 1); -#X msg 24 515 setForce \$1 \$2 \$3; -#X text 186 553 \$2 / \$3 : Value; -#X msg 24 535 setForceX \$1 \$2; -#X msg 24 554 setForceY \$1 \$2; -#X msg 31 96 pos \$1 \$2 \$3; -#X text 190 132 \$2 / \$3 : Value; -#X msg 30 173 force \$1 \$2 \$3; -#X text 191 213 \$2 / \$3 : Value; -#X text 180 284 \$1 / \$2 : Value; -#X msg 27 805 setDEnv \$1; -#X text 178 798 Change environement damping; -#X text 190 815 \$1 : damping; -#X msg 23 602 setSpeedX \$1 \$2; -#X msg 23 621 setSpeedY \$1 \$2; -#X text 185 619 \$2 / \$3 : Value; -#X text 185 681 \$2 / \$3 : Value; -#X msg 24 646 addPos \$1 \$2 \$3; -#X msg 588 280 setLinkId \$1 \$2; -#X text 742 281 Change the Id of a link; -#X text 743 298 \$1 : link (number or Id); -#X text 746 315 \$2 : New Id; -#X msg 591 353 setLKTab \$1 \$2; -#X text 750 353 change the table (K) size of a non linear link; -#X text 751 373 \$1 : link (number or Id); -#X text 752 390 \$2 : New KL; -#X text 751 435 \$1 : link (number or Id); -#X msg 591 415 setLDTab \$1 \$2; -#X text 750 415 change the table (D) size of a non linear link; -#X text 752 452 \$2 : New DL; -#X msg 31 245 min \$1 \$2; -#X msg 31 265 max \$1 \$2; -#X msg 24 665 addPosX \$1 \$2; -#X msg 24 685 addPosY \$1 \$2; -#X restore 12 478 pd dynamic settings; -#N canvas 30 26 1236 814 attributes______ 0; +#X restore 12 457 pd creation____________; +#N canvas 30 26 1236 814 attributes__________ 0; #X obj 11 95 cnv 15 100 35 empty empty empty 20 12 0 14 -233017 -66577 0; #X obj 5 75 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 @@ -618,98 +734,121 @@ a table (\$1); #X text 941 729 [massesPosId 2 x1 y2(; #X text 942 765 [massesPosName mas x1 y1(; #X text 942 780 [massesPosName mas x2 y2(; -#X restore 12 499 pd attributes______; -#X text 9 649 KEYWORDS: physical model mass spring damper link; -#X text 12 8 HELP: pmpd2d; -#X obj 157 48 pmpd2d; -#X text 112 618 01_pmpd2dtest.pd; -#X obj 18 296 pmpd2d; -#N canvas 731 296 450 300 gemmouse 0; -#X obj 189 77 gemmouse; -#X obj 189 184 pack f f; -#X obj 189 218 spigot; -#X obj 109 243 outlet; -#X obj 109 36 inlet; -#X obj 109 98 t b; -#X obj 238 131 + 4; -#X msg 57 182 posX fix \$1; -#X msg 57 203 posY fix \$2; -#X obj 189 108 / 62.5; -#X obj 237 107 / -62.5; -#X obj 189 131 - 4; -#X obj 109 125 list 0 0; -#X connect 0 0 9 0; -#X connect 0 1 10 0; -#X connect 0 2 2 1; -#X connect 1 0 2 0; -#X connect 2 0 12 1; -#X connect 4 0 5 0; -#X connect 5 0 12 0; -#X connect 6 0 1 1; -#X connect 7 0 3 0; -#X connect 8 0 3 0; -#X connect 9 0 11 0; -#X connect 10 0 6 0; -#X connect 11 0 1 0; -#X connect 12 0 7 0; -#X connect 12 0 8 0; -#X restore 393 295 pd gemmouse; -#X obj 393 274 gemhead; -#X obj 346 216 gemwin; -#X msg 357 193 0 \, destroy; -#N canvas 472 258 550 319 gemrender 0; -#X obj 48 203 translateXYZ; -#X obj 48 229 sphere 0.1; -#X obj 127 24 inlet; -#X obj 360 32 inlet; -#X obj 275 232 curve 2; -#X msg 431 81 \$4 \$5 0; -#X msg 359 77 \$2 \$3 0; -#X obj 359 105 t b a; -#X obj 127 62 unpack f f f; -#X obj 166 88 t b f; -#X msg 48 110 0; -#X obj 48 77 loadbang; -#X obj 48 137 gemhead; -#X msg 275 160 0; -#X obj 275 127 loadbang; -#X obj 275 187 gemhead; -#X connect 0 0 1 0; -#X connect 2 0 8 0; -#X connect 3 0 5 0; -#X connect 3 0 6 0; -#X connect 5 0 4 2; -#X connect 6 0 7 0; -#X connect 7 0 15 0; -#X connect 7 1 4 1; -#X connect 8 1 9 0; -#X connect 8 2 0 2; -#X connect 9 0 12 0; -#X connect 9 1 0 1; -#X connect 10 0 12 0; -#X connect 11 0 10 0; -#X connect 12 0 0 0; -#X connect 13 0 15 0; -#X connect 14 0 13 0; -#X connect 15 0 4 0; -#X restore 172 320 pd gemrender; -#X msg 346 129 reset \, create \, 1; -#X text 135 98 Sorry \, you need GEM for this example...; -#X text 169 132 1 Create window -->; -#X text 169 163 2 Drag the structure with the mouse; -#X text 169 195 3 Destroy the window -->; -#X obj 172 294 route massesPos linksPosNo; -#X text 11 19 DESCRIPTION: Physical modeling for Pd in 2d; -#X text 265 669 - Cyrille Henry 2011; -#X text 265 681 based on MSD by Nicolas Montgermont \, 2005; -#X connect 21 0 43 0; +#X restore 12 499 pd attributes__________; +#N canvas 100 72 1078 872 dynamic 0; +#X obj 5 3 cnv 15 550 15 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 10 85 cnv 15 150 750 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 4 62 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X obj 574 86 cnv 15 130 350 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 567 62 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X text 7 62 Masses :; +#X text 570 62 Links :; +#X text 7 2 DYNAMIC SETTINGS Messages; +#X msg 31 116 posX \$1 \$2; +#X text 191 168 Add force on mass(es); +#X msg 30 183 forceX \$1 \$2; +#X text 190 114 \$1 : Id (symbol) or No; +#X text 193 185 \$1 : Id (symbol) or No; +#X msg 31 265 Xmin \$1; +#X msg 90 265 Xmax \$1; +#X msg 29 311 setMobile \$1; +#X msg 29 332 setFixed \$1; +#X text 190 96 Set position of fixed mass(es); +#X text 177 248 Set minimimum and maximum position of all masses; +#X text 188 311 Set mass to mobile or fixed; +#X msg 593 150 setD \$1 \$2; +#X text 748 136 \$2 : New value; +#X msg 593 99 setK \$1 \$2; +#X text 748 192 \$2 : New value; +#X text 748 245 \$2 : New value; +#X text 742 100 Set rigidity of link(s); +#X text 742 156 Set damping of link(s); +#X msg 593 205 setL \$1 \$2; +#X text 742 209 Set initial lenght of link(s); +#X text 748 118 \$1 : Id (symbol) or No; +#X text 748 174 \$1 : Id (symbol) or No; +#X text 748 227 \$1 : Id (symbol) or No; +#X text 191 328 \$1 : Id (symbol) or No; +#X text 10 25 Dynamic settings messages allows the user to redefine +internal parameters of links and masses.; +#X msg 29 358 grabMass \$1 \$2 \$3; +#X text 186 358 Grab nearest mass; +#X text 191 374 \$1 \, \$2 : position; +#X text 190 391 \$3 : grab or not (0/1); +#X msg 31 137 posY \$1 \$2; +#X msg 30 204 forceY \$1 \$2; +#X msg 31 287 Ymin \$1; +#X msg 90 287 Ymax \$1; +#X text 191 445 \$2 : New Id; +#X msg 30 413 setMassId \$1 \$2; +#X text 191 411 Change the Id of a mass; +#X text 192 428 \$1 : mass (number or Id); +#X msg 23 530 setSpeed \$1 \$2; +#X text 186 482 \$1 : Id (symbol) or No; +#X text 183 465 set force on mass(es); +#X text 185 549 \$1 : Id (symbol) or No; +#X text 182 532 set speed of a mass(es); +#X text 184 610 \$1 : Id (symbol) or No; +#X text 181 593 add pos of a mass(es); +#X msg 23 684 grabMass \$1 \$2 \$3; +#X text 183 654 grab a mass \, and move it; +#X text 185 676 \$1 : mass number; +#X text 185 692 \$2 : position; +#X text 185 708 \$3 : grab or not (0 or 1); +#X msg 24 462 setForce \$1 \$2 \$3; +#X text 186 500 \$2 / \$3 : Value; +#X msg 24 482 setForceX \$1 \$2; +#X msg 24 501 setForceY \$1 \$2; +#X msg 31 96 pos \$1 \$2 \$3; +#X text 190 132 \$2 / \$3 : Value; +#X msg 30 163 force \$1 \$2 \$3; +#X text 191 203 \$2 / \$3 : Value; +#X text 180 264 \$1 / \$2 : Value; +#X msg 23 549 setSpeedX \$1 \$2; +#X msg 23 568 setSpeedY \$1 \$2; +#X text 185 566 \$2 / \$3 : Value; +#X text 185 628 \$2 / \$3 : Value; +#X msg 24 593 addPos \$1 \$2 \$3; +#X msg 588 280 setLinkId \$1 \$2; +#X text 742 281 Change the Id of a link; +#X text 743 298 \$1 : link (number or Id); +#X text 746 315 \$2 : New Id; +#X msg 591 353 setLKTab \$1 \$2; +#X text 750 353 change the table (K) size of a non linear link; +#X text 751 373 \$1 : link (number or Id); +#X text 752 390 \$2 : New KL; +#X text 751 435 \$1 : link (number or Id); +#X msg 591 415 setLDTab \$1 \$2; +#X text 750 415 change the table (D) size of a non linear link; +#X text 752 452 \$2 : New DL; +#X msg 31 225 min \$1 \$2; +#X msg 31 245 max \$1 \$2; +#X msg 24 612 addPosX \$1 \$2; +#X msg 24 632 addPosY \$1 \$2; +#X text 571 580 if only 1 argument is previded \, message generally +used this argument as the value applayed to all masses (or link); +#X text 186 819 \$2 : damping; +#X text 185 804 \$1 : Id (symbol) or No; +#X msg 22 788 setDEnvOffset \$1 \$2; +#X text 174 788 Change environement damping offset; +#X text 173 733 Change environement damping; +#X msg 23 733 setDEnv \$1 \$2; +#X text 185 764 \$2 : damping; +#X text 184 749 \$1 : Id (symbol) or No; +#X restore 12 478 pd dynamic settings ___; +#X connect 21 0 40 0; #X connect 24 0 23 0; -#X connect 25 0 54 0; +#X connect 25 0 51 0; #X connect 26 0 19 0; -#X connect 43 0 20 0; -#X connect 44 0 22 0; -#X connect 45 0 44 0; -#X connect 47 0 46 0; -#X connect 49 0 46 0; -#X connect 54 0 48 0; -#X connect 54 1 48 1; +#X connect 40 0 20 0; +#X connect 41 0 22 0; +#X connect 42 0 41 0; +#X connect 44 0 43 0; +#X connect 46 0 43 0; +#X connect 51 0 45 0; +#X connect 51 1 45 1; diff --git a/pmpd2d.c b/pmpd2d.c index 7920388..d5fbd50 100644 --- a/pmpd2d.c +++ b/pmpd2d.c @@ -31,8 +31,8 @@ #include "stdio.h" #include "math.h" -#define nb_max_link 100000 -#define nb_max_mass 100000 +#define nb_max_link 10000000 +#define nb_max_mass 1000000 #define max(a,b) ( ((a) > (b)) ? (a) : (b) ) #define min(a,b) ( ((a) < (b)) ? (a) : (b) ) @@ -89,7 +89,7 @@ typedef struct _link { t_float D_L; // longeur du tabeau D } foo1 ; -typedef struct _pmpd2d2d { +typedef struct _pmpd2d { t_object x_obj; struct _link link[nb_max_link]; struct _mass mass[nb_max_mass]; @@ -143,7 +143,9 @@ void pmpd2d_infosL(t_pmpd2d *x) 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("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"); @@ -153,13 +155,19 @@ void pmpd2d_infosL(t_pmpd2d *x) 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); + 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); + 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); + 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; } } @@ -329,13 +337,14 @@ void pmpd2d_link(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv) else if ( ( argv[1].a_type == A_SYMBOL ) && ( argv[2].a_type == A_SYMBOL ) ) { - for (i=0; i< x->nb_mass; i++) + for (i=0; i < x->nb_mass; i++) { - for (j=0; j< x->nb_mass; j++) + 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 ( (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 + 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); } } @@ -406,7 +415,8 @@ void pmpd2d_tLink(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv) { 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 + 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; @@ -480,7 +490,8 @@ void pmpd2d_tabLink(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv) { 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 + 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; @@ -2593,78 +2604,1114 @@ void pmpd2d_massesSpeedsStd(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv) outlet_anything(x->main_outlet, gensym("massesSpeedsStd"),3 , std_out); } -void pmpd2d_grabMass(t_pmpd2d *x, t_float posX, t_float posY, t_float grab) +// -------------------------------------------- + +void pmpd2d_linksPosXL(t_pmpd2d *x) { - t_float dist, tmp; - t_int i; - - if (grab == 0) - x->grab=0; - if ((x->grab == 0)&(grab == 1)&(x->nb_mass > 0)) + 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) ) { - x->grab=1; - x->grab_nb= 0; - dist = sqr(x->mass[0].posX - posX) + sqr(x->mass[0].posY - posY); - for (i=1; inb_mass; i++) + 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 { - tmp = sqr(x->mass[i].posX - posX) + sqr(x->mass[i].posY - posY); - if (tmp < dist) + int taille_max = x->nb_link; + taille_max = min(taille_max, vecsize); + for (i=0; i < taille_max ; i++) { - dist = tmp; - x->grab_nb= i; + vec[i].w_float = x->link[i].mass2->posX - x->link[i].mass1->posX; } + garray_redraw(a); } } - if (x->grab == 1) + else + if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) ) { - x->mass[x->grab_nb].posX = posX; - x->mass[x->grab_nb].posY = posY; + 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_new() +void pmpd2d_linksPosSpeedXT(t_pmpd2d *x, t_symbol *s, int argc, t_atom *argv) { - t_pmpd2d *x = (t_pmpd2d *)pd_new(pmpd2d_class); + 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); + } + } +} - pmpd2d_reset(x); - - x->main_outlet=outlet_new(&x->x_obj, 0); - // x->info_outlet=outlet_new(&x->x_obj, 0); // TODO +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); + } + } +} - return (void *)x; +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_setup(void) +void pmpd2d_linksLengthYL(t_pmpd2d *x) { - pmpd2d_class = class_new(gensym("pmpd2d"), - (t_newmethod)pmpd2d_new, - 0, sizeof(t_pmpd2d),CLASS_DEFAULT, 0); + int i; + t_atom pos_list[x->nb_link]; - 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_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_setL, gensym("setL"), 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); + 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); +} - 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); +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->mass[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->mass[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->mass[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->mass[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_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; inb_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_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; +} + +void pmpd2d_setup(void) +{ + pmpd2d_class = class_new(gensym("pmpd2d"), + (t_newmethod)pmpd2d_new, + 0, sizeof(t_pmpd2d),CLASS_DEFAULT, 0); + + 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_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_setL, gensym("setL"), 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); @@ -2712,15 +3759,59 @@ void pmpd2d_setup(void) 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_linkLMean, gensym("linkLMean"),A_GIMME, 0); - // class_addmethod(pmpd2d_class, (t_method)pmpd2d_linkLStd, gensym("linkLStd"),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("linksLengthXL"), 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("linksLengthYL"), 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("linksLengthL"), 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_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("linksLengthNormL"), 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_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_grabMass, gensym("grabMass"), A_DEFFLOAT, A_DEFFLOAT, A_DEFFLOAT, 0); } diff --git a/pmpd3d-help.pd b/pmpd3d-help.pd index e81aae1..4ce2282 100644 --- a/pmpd3d-help.pd +++ b/pmpd3d-help.pd @@ -1,23 +1,23 @@ -#N canvas 1045 104 613 729 10; +#N canvas 52 60 624 761 10; #X obj 4 369 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 0; -#X obj 5 665 cnv 15 550 30 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 5 695 cnv 15 550 30 empty empty empty 20 12 0 14 -158509 -66577 0; #X obj 8 6 cnv 15 550 30 empty empty empty 20 12 0 14 -233017 -66577 0; #X obj 14 389 cnv 15 75 15 empty empty empty 20 12 0 14 -233017 -66577 0; -#X obj 14 584 cnv 15 75 15 empty empty empty 20 12 0 14 -233017 -66577 +#X obj 14 614 cnv 15 75 15 empty empty empty 20 12 0 14 -233017 -66577 0; -#X obj 4 527 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 4 557 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 0; -#X obj 4 562 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 4 592 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 0; #X obj 7 76 cnv 15 550 15 empty empty empty 20 12 0 14 -158509 -66577 0; -#X obj 5 615 cnv 15 550 20 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 5 645 cnv 15 550 20 empty empty empty 20 12 0 14 -158509 -66577 0; -#X obj 5 635 cnv 15 550 30 empty empty empty 20 12 0 14 -233017 -66577 +#X obj 5 665 cnv 15 550 30 empty empty empty 20 12 0 14 -233017 -66577 0; #N canvas 76 31 921 714 More_Info 0; #X text 96 12 pmpd : mass - spring - damper model; @@ -226,13 +226,13 @@ their Id.; #X text 28 41 pmpd is the 2D object of the pmpd objects collection. ; #X connect 3 0 4 0; -#X restore 16 673 pd More_Info; +#X restore 16 703 pd More_Info; #X text 12 76 Examples:; #X text 9 369 Inlets:; #X text 19 388 - Left:; -#X text 10 526 Arguments:; -#X text 11 562 Outlets:; -#X text 19 617 See Also:; +#X text 10 556 Arguments:; +#X text 11 592 Outlets:; +#X text 19 647 See Also:; #X text 74 48 Full Name:; #N canvas 58 22 262 70 Related_Objects 0; #X obj 3 10 cnv 15 250 15 empty empty empty 20 12 0 14 -158509 -66577 @@ -242,13 +242,122 @@ their Id.; #X text 5 10 Externals and libraries; #X obj 44 37 pmpd; #X obj 141 39 pmpd3d; -#X restore 122 673 pd Related_Objects; +#X restore 122 703 pd Related_Objects; #X text 101 388 Bang - A bang at the left inlet compute the new model state based on previous instant.; -#X text 158 478 To set the model parameters after creation.; -#X text 158 456 To create the model masses and links.; -#X text 158 501 To get the model parameters; -#N canvas 556 79 967 626 creation________ 0; +#X text 188 456 To create the model masses and links.; +#X text 188 510 To get the model parameters; +#X text 103 572 None; +#X text 18 613 - Left:; +#X text 101 614 Outputs the model parameters asked with the attributes +messages.; +#X text 101 420 Messages - Different messages are used to control the +pmpd object. They are of three types :; +#X text 9 669 CATEGORY: control; +#X text 9 679 KEYWORDS: physical model mass spring damper link; +#X text 135 98 Sorry \, you need GEM for this example...; +#X text 265 699 - Cyrille Henry 2011; +#X text 265 711 based on MSD by Nicolas Montgermont \, 2005; +#X text 12 8 HELP: pmpd3d; +#X text 11 19 DESCRIPTION: Physical modeling for Pd in 3d; +#X obj 157 48 pmpd3d; +#N canvas 257 262 759 345 init 0; +#X obj 89 215 t a; +#X obj 89 33 loadbang; +#X obj 89 241 s \$0-in; +#X obj 89 59 t b b b b b; +#X msg 161 80 reset; +#X obj 44 13 inlet; +#X msg 107 161 link souple fix mob 10 10 \, link souple mob mob 10 +10; +#X msg 125 120 mass mob 1 100 0 -2 0 \, mass mob 1 100 2 0 0 \, mass +mob 1 100 0 2 0 \, mass mob 1 100 -2 0 0; +#X msg 89 187 Xmax 4 \, Xmin -4 \, Ymax 4 \, Ymin -4 \, Zmax 4 \, Zmin +-4; +#X msg 143 100 mass fix 0 10 0 0 0; +#X connect 0 0 2 0; +#X connect 1 0 3 0; +#X connect 3 0 8 0; +#X connect 3 1 6 0; +#X connect 3 2 7 0; +#X connect 3 3 9 0; +#X connect 3 4 4 0; +#X connect 4 0 0 0; +#X connect 5 0 3 0; +#X connect 6 0 0 0; +#X connect 7 0 0 0; +#X connect 8 0 0 0; +#X connect 9 0 0 0; +#X restore 17 156 pd init; +#X obj 18 321 s \$0-out; +#X obj 18 272 r \$0-in; +#X obj 393 317 s \$0-in; +#N canvas 565 515 355 193 compute 0; +#X obj 27 29 inlet; +#X obj 27 127 s \$0-in; +#X msg 27 96 bang \, get massesPos \, get linksPos; +#X obj 27 66 gemhead; +#X obj 160 36 bng 15 250 50 0 empty empty empty 0 -6 0 8 -262144 -1 +-1; +#X connect 0 0 3 0; +#X connect 2 0 1 0; +#X connect 3 0 2 0; +#X connect 4 0 2 0; +#X restore 17 209 pd compute; +#X obj 17 181 tgl 15 0 empty empty ON/OFF 25 10 1 10 -262144 -1 -1 +0 1; +#X obj 172 270 r \$0-out; +#X obj 17 126 bng 15 250 50 0 empty empty reset 25 10 1 10 -262144 +-1 -1; +#X obj 336 216 gemwin; +#X msg 336 194 0 \, destroy; +#N canvas 363 318 550 319 gemrender 0; +#X obj 48 203 translateXYZ; +#X obj 48 229 sphere 0.1; +#X obj 127 24 inlet; +#X obj 390 13 inlet; +#X obj 303 213 curve 2; +#X obj 125 136 t b f; +#X obj 127 62 unpack f f f f; +#X msg 390 58 \$2 \$3 \$4; +#X msg 462 62 \$5 \$6 \$7; +#X obj 48 167 gemhead; +#X msg 48 143 0; +#X obj 48 115 loadbang; +#X obj 303 147 gemhead; +#X msg 303 123 0; +#X obj 303 95 loadbang; +#X obj 390 86 t b l; +#X connect 0 0 1 0; +#X connect 2 0 6 0; +#X connect 3 0 8 0; +#X connect 3 0 7 0; +#X connect 5 0 9 0; +#X connect 5 1 0 1; +#X connect 6 1 5 0; +#X connect 6 2 0 2; +#X connect 6 3 0 3; +#X connect 7 0 15 0; +#X connect 8 0 4 2; +#X connect 9 0 0 0; +#X connect 10 0 9 0; +#X connect 11 0 10 0; +#X connect 12 0 4 0; +#X connect 13 0 12 0; +#X connect 14 0 13 0; +#X connect 15 0 12 0; +#X connect 15 1 4 1; +#X restore 172 321 pd gemrender; +#X msg 336 129 reset \, create \, 1; +#X obj 18 296 pmpd3d; +#X msg 393 271 forceZ mob -10; +#X msg 393 293 forceZ mob 10; +#X text 392 246 2 Send forces; +#X text 170 132 1 Create window -->; +#X text 170 195 3 And destroy it -->; +#X obj 172 294 route massesPos linksPosNo; +#X text 112 648 01_pmpd3dtest.pd; +#N canvas 556 79 967 626 creation____________ 0; #X obj 5 75 cnv 15 450 15 empty empty empty 20 12 0 14 -158509 -66577 0; #X obj 5 3 cnv 15 450 15 empty empty empty 20 12 0 14 -233017 -66577 @@ -320,15 +429,8 @@ following a vector.; #X msg 32 167 mass \$1 \$2 \$3 \$4 \$5 \$6; #X text 58 247 \$4 \, \$5 \, \$6 : initial position; #X connect 27 0 29 0; -#X restore 12 457 pd creation________; -#X text 103 542 None; -#X text 18 583 - Left:; -#X text 101 584 Outputs the model parameters asked with the attributes -messages.; -#X text 101 420 Messages - Different messages are used to control the -pmpd object. They are of three types :; -#X text 9 639 CATEGORY: control; -#N canvas 675 81 1078 872 dynamic 0; +#X restore 12 457 pd creation____________; +#N canvas 541 24 1078 872 dynamic 0; #X obj 5 3 cnv 15 550 15 empty empty empty 20 12 0 14 -233017 -66577 0; #X obj 10 85 cnv 15 150 750 empty empty empty 20 12 0 14 -233017 -66577 @@ -369,46 +471,44 @@ pmpd object. They are of three types :; #X text 191 364 \$1 : Id (symbol) or No; #X text 10 25 Dynamic settings messages allows the user to redefine internal parameters of links and masses.; -#X msg 29 404 grabMass \$1 \$2 \$3; -#X text 186 404 Grab nearest mass; -#X text 191 420 \$1 \, \$2 : position; -#X text 190 437 \$3 : grab or not (0/1); +#X msg 29 391 grabMass \$1 \$2 \$3; +#X text 186 391 Grab nearest mass; +#X text 191 407 \$1 \, \$2 : position; +#X text 190 424 \$3 : grab or not (0/1); #X msg 31 134 posY \$1 \$2; #X msg 30 210 forceY \$1 \$2; #X msg 31 307 Ymin \$1; #X msg 90 307 Ymax \$1; -#X text 191 501 \$2 : New Id; -#X msg 29 469 setMassId \$1 \$2; -#X text 191 467 Change the Id of a mass; -#X text 192 484 \$1 : mass (number or Id); -#X msg 24 594 setSpeed \$1 \$2; -#X text 186 548 \$1 : Id (symbol) or No; -#X text 183 531 set force on mass(es); -#X text 185 615 \$1 : Id (symbol) or No; -#X text 182 598 set speed of a mass(es); -#X text 178 694 \$1 : Id (symbol) or No; -#X text 175 677 add pos of a mass(es); -#X msg 27 764 grabMass \$1 \$2 \$3; -#X text 183 733 grab a mass \, and move it; -#X text 185 755 \$1 : mass number; -#X text 185 771 \$2 : position; -#X text 185 787 \$3 : grab or not (0 or 1); -#X msg 24 514 setForce \$1 \$2 \$3; -#X text 186 566 \$2 / \$3 : Value; -#X msg 24 533 setForceX \$1 \$2; -#X msg 24 552 setForceY \$1 \$2; +#X text 191 476 \$2 : New Id; +#X msg 29 444 setMassId \$1 \$2; +#X text 191 442 Change the Id of a mass; +#X text 192 459 \$1 : mass (number or Id); +#X msg 24 564 setSpeed \$1 \$2; +#X text 186 518 \$1 : Id (symbol) or No; +#X text 183 501 set force on mass(es); +#X text 185 585 \$1 : Id (symbol) or No; +#X text 182 568 set speed of a mass(es); +#X text 178 664 \$1 : Id (symbol) or No; +#X text 175 647 add pos of a mass(es); +#X msg 27 734 grabMass \$1 \$2 \$3; +#X text 183 703 grab a mass \, and move it; +#X text 185 725 \$1 : mass number; +#X text 185 741 \$2 : position; +#X text 185 757 \$3 : grab or not (0 or 1); +#X msg 24 484 setForce \$1 \$2 \$3; +#X text 186 536 \$2 / \$3 : Value; +#X msg 24 503 setForceX \$1 \$2; +#X msg 24 522 setForceY \$1 \$2; #X msg 31 96 pos \$1 \$2 \$3; #X text 190 132 \$2 / \$3 : Value; #X msg 30 172 force \$1 \$2 \$3; #X text 191 213 \$2 / \$3 : Value; -#X msg 29 811 setDEnv \$1; -#X text 178 811 Change environement damping; -#X text 190 828 \$1 : damping; -#X msg 24 613 setSpeedX \$1 \$2; -#X msg 24 632 setSpeedY \$1 \$2; -#X text 185 632 \$2 / \$3 : Value; -#X text 179 712 \$2 / \$3 : Value; -#X msg 24 675 addPos \$1 \$2 \$3; +#X text 178 781 Change environement damping; +#X msg 24 583 setSpeedX \$1 \$2; +#X msg 24 602 setSpeedY \$1 \$2; +#X text 185 602 \$2 / \$3 : Value; +#X text 179 682 \$2 / \$3 : Value; +#X msg 24 645 addPos \$1 \$2 \$3; #X msg 588 280 setLinkId \$1 \$2; #X text 742 281 Change the Id of a link; #X text 743 298 \$1 : link (number or Id); @@ -428,32 +528,42 @@ internal parameters of links and masses.; #X text 180 288 \$1 / \$2 / \$3 : Value; #X msg 31 326 Zmin \$1; #X msg 90 326 Zmax \$1; -#X msg 24 571 setForceZ \$1 \$2; -#X msg 24 651 setSpeedZ \$1 \$2; -#X msg 24 694 addPosX \$1 \$2; -#X msg 24 713 addPosY \$1 \$2; -#X msg 25 732 addPosZ \$1 \$2; -#X restore 12 478 pd dynamic settings; -#N canvas 30 26 1236 814 attributes______ 0; -#X obj 11 95 cnv 15 100 35 empty empty empty 20 12 0 14 -233017 -66577 -0; -#X obj 5 75 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 +#X msg 24 541 setForceZ \$1 \$2; +#X msg 24 621 setSpeedZ \$1 \$2; +#X msg 24 664 addPosX \$1 \$2; +#X msg 24 683 addPosY \$1 \$2; +#X msg 25 702 addPosZ \$1 \$2; +#X msg 29 781 setDEnv \$1 \$2; +#X text 190 812 \$2 : damping; +#X text 189 797 \$1 : Id (symbol) or No; +#X text 571 580 if only 1 argument is previded \, message generally +used this argument as the value applayed to all masses (or link); +#X text 191 855 \$2 : damping; +#X text 190 840 \$1 : Id (symbol) or No; +#X msg 30 824 setDEnvOffset \$1 \$2; +#X text 179 824 Change environement damping offset; +#X restore 12 478 pd dynamic settings ___; +#X text 188 478 To set the model parameters after creation.; +#N canvas 30 26 1241 841 masses_attributes 0; +#X obj 11 115 cnv 15 100 35 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 5 95 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 0; #X obj 5 3 cnv 15 590 15 empty empty empty 20 12 0 14 -233017 -66577 0; -#X obj 10 158 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +#X obj 10 178 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 0; -#X obj 4 137 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 +#X obj 4 157 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 0; #X obj 616 422 cnv 15 110 330 empty empty empty 20 12 0 14 -233017 -66577 0; #X obj 610 398 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 0; -#X text 7 74 General :; +#X text 7 94 General :; #X text 7 2 ATTRIBUTES Messages; -#X text 7 137 Lists :; -#X text 136 104 Get infos on all masses and links on right outlet; -#X text 170 170 Output all masses positions in a list on outlet No +#X text 7 157 Lists :; +#X text 136 124 Get infos on all masses and links on right outlet; +#X text 170 190 Output all masses positions in a list on outlet No 1; #X text 746 438 Get specific attribute on specific element; #X msg 626 436 get \$1 (\$2); @@ -465,29 +575,26 @@ is made of the asked parameter and the way you asked for it.; #X text 987 609 response; #X text 746 590 Examples with 3 masses numbered 0 \, 1 and 2 and named mas:; -#X text 15 30 The attributes messages ask the object to output some -of his internal parameters. They can be output by lists for positions -and forces of masses.; -#X text 172 399 Output all forces applied on masses in a list on outlet +#X text 172 419 Output all forces applied on masses in a list on outlet No 1; #X text 746 500 (\$2) : - If not defined all the attributes are send for all the elements. - Ids or/and creations No; -#X text 170 202 Output all masses x or y in a list on outlet No 1; -#X msg 33 104 print; -#X text 170 230 Output all norm of position vector (sqrt(x*x+y*y)) +#X text 170 222 Output all masses x or y in a list on outlet No 1; +#X msg 33 124 print; +#X text 170 250 Output all norm of position vector (sqrt(x*x+y*y)) in a list on outlet No 1; -#X obj 11 272 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +#X obj 11 292 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 0; -#X text 171 284 Output all masses speeds in a list on outlet No 1; -#X text 172 316 Output all masses speeds in x or y in a list on outlet +#X text 171 304 Output all masses speeds in a list on outlet No 1; +#X text 172 336 Output all masses speeds in x or y in a list on outlet No 1; -#X text 179 353 Output all norm of speed vector (sqrt(x*x+y*y)) in +#X text 179 373 Output all norm of speed vector (sqrt(x*x+y*y)) in a list on outlet No 1; -#X obj 11 385 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +#X obj 11 405 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 0; -#X text 175 429 Output all masses forces in x or y in a list on outlet +#X text 175 449 Output all masses forces in x or y in a list on outlet No 1; -#X text 175 463 Output all norm of forces vector (sqrt(x*x+y*y)) in +#X text 175 483 Output all norm of forces vector (sqrt(x*x+y*y)) in a list on outlet No 1; #X text 745 457 \$1 : Attribute type ( massesPos / massesPosName / massesSpeeds / massesSpeedsName / massesForces / massesForcesName / @@ -511,15 +618,15 @@ a list on outlet No 1; No 1; #X text 804 322 Output all norm of forces vector (sqrt(x*x+y*y)) in a list on outlet No 1; -#X msg 12 169 massesPosL (\$1); -#X msg 12 188 massesPosXL (\$1); -#X msg 12 207 massesPosYL (\$1); -#X msg 12 245 massesPosNormL (\$1); -#X msg 13 279 massesSpeedsL (\$1); -#X msg 14 355 massesSpeedsNormL (\$1); -#X msg 15 392 massesForcesL (\$1); -#X msg 12 169 massesPosL (\$1); -#X msg 14 468 massesSpeedsNormL (\$1); +#X msg 12 189 massesPosL (\$1); +#X msg 12 208 massesPosXL (\$1); +#X msg 12 227 massesPosYL (\$1); +#X msg 12 265 massesPosNormL (\$1); +#X msg 13 299 massesSpeedsL (\$1); +#X msg 14 375 massesSpeedsNormL (\$1); +#X msg 15 412 massesForcesL (\$1); +#X msg 12 189 massesPosL (\$1); +#X msg 14 488 massesSpeedsNormL (\$1); #X msg 615 38 massesPosT \$1 (\$2); #X msg 615 57 massesPosXT \$1 (\$2); #X msg 615 76 massesPosYT \$1 (\$2); @@ -531,40 +638,40 @@ a list on outlet No 1; #X text 794 67 Output all masses x or y in a table (\$1); #X text 794 95 Output all norm of position vector (sqrt(x*x+y*y)) in a table (\$1); -#X text 174 498 \$1 (optional): Masses Id; +#X text 174 518 \$1 (optional): Masses Id; #X text 815 376 \$2 (optional): Masses Id; #X text 814 363 \$1 : table name; -#X obj 17 544 cnv 15 150 270 empty empty empty 20 12 0 14 -233017 -66577 -0; -#X obj 11 524 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 -0; -#X text 14 524 statistics :; -#X msg 22 557 massesPosMean; -#X text 177 573 \$1 : none or Id; -#X msg 23 595 massesPosStd \$1; -#X text 179 610 \$1 : none or Id; -#X text 173 557 output the average position of masses; -#X text 174 594 output the standard deviation of masses position; -#X text 177 660 \$1 : none or Id; -#X text 179 697 \$1 : none or Id; -#X text 173 644 output the average position of masses; -#X text 174 681 output the standard deviation of masses position; -#X text 177 745 \$1 : none or Id; -#X text 178 782 \$1 : none or Id; -#X text 173 729 output the average position of masses; -#X text 173 766 output the standard deviation of masses position; -#X msg 22 682 massesForcesStd \$1; -#X msg 22 644 massesForcesMean; -#X msg 23 731 massesSpeedsMean; -#X msg 22 767 massesSpeedsStd \$1; +#X obj 17 564 cnv 15 150 270 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 11 544 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X text 14 544 statistics :; +#X msg 22 577 massesPosMean; +#X text 177 593 \$1 : none or Id; +#X msg 23 615 massesPosStd \$1; +#X text 179 630 \$1 : none or Id; +#X text 173 577 output the average position of masses; +#X text 174 614 output the standard deviation of masses position; +#X text 177 680 \$1 : none or Id; +#X text 179 717 \$1 : none or Id; +#X text 173 664 output the average position of masses; +#X text 174 701 output the standard deviation of masses position; +#X text 177 765 \$1 : none or Id; +#X text 178 802 \$1 : none or Id; +#X text 173 749 output the average position of masses; +#X text 173 786 output the standard deviation of masses position; +#X msg 22 702 massesForcesStd \$1; +#X msg 22 664 massesForcesMean; +#X msg 23 751 massesSpeedsMean; +#X msg 22 787 massesSpeedsStd \$1; #X text 608 5 tables :; -#X msg 12 226 massesPosZL (\$1); -#X msg 14 298 massesSpeedsXL (\$1); -#X msg 13 317 massesSpeedsYL (\$1); -#X msg 14 336 massesSpeedsZL (\$1); -#X msg 15 411 massesForcesXL (\$1); -#X msg 15 430 massesForcesYL (\$1); -#X msg 14 449 massesForcesZL (\$1); +#X msg 12 246 massesPosZL (\$1); +#X msg 14 318 massesSpeedsXL (\$1); +#X msg 13 337 massesSpeedsYL (\$1); +#X msg 14 356 massesSpeedsZL (\$1); +#X msg 15 431 massesForcesXL (\$1); +#X msg 15 450 massesForcesYL (\$1); +#X msg 14 469 massesForcesZL (\$1); #X text 780 629 [get massesPos( -----> [massesPos 0 x0 y0 z0(; #X text 940 641 [massesPos 1 x1 y1 z1(; #X text 940 656 [massesPos 2 x2 y2 z2(; @@ -585,118 +692,158 @@ z0(; #X msg 616 278 massesForcesXT \$1 (\$2); #X msg 616 297 massesForcesYT \$1 (\$2); #X msg 616 316 massesForcesZT \$1 (\$2); -#X restore 12 499 pd attributes______; -#X text 9 649 KEYWORDS: physical model mass spring damper link; -#X text 135 98 Sorry \, you need GEM for this example...; -#X text 265 669 - Cyrille Henry 2011; -#X text 265 681 based on MSD by Nicolas Montgermont \, 2005; -#X text 12 8 HELP: pmpd3d; -#X text 11 19 DESCRIPTION: Physical modeling for Pd in 3d; -#X obj 157 48 pmpd3d; -#N canvas 257 262 759 345 init 0; -#X obj 89 215 t a; -#X obj 89 33 loadbang; -#X obj 89 241 s \$0-in; -#X obj 89 59 t b b b b b; -#X msg 161 80 reset; -#X obj 44 13 inlet; -#X msg 107 161 link souple fix mob 10 10 \, link souple mob mob 10 -10; -#X msg 125 120 mass mob 1 100 0 -2 0 \, mass mob 1 100 2 0 0 \, mass -mob 1 100 0 2 0 \, mass mob 1 100 -2 0 0; -#X msg 89 187 Xmax 4 \, Xmin -4 \, Ymax 4 \, Ymin -4 \, Zmax 4 \, Zmin --4; -#X msg 143 100 mass fix 0 10 0 0 0; -#X connect 0 0 2 0; -#X connect 1 0 3 0; -#X connect 3 0 8 0; -#X connect 3 1 6 0; -#X connect 3 2 7 0; -#X connect 3 3 9 0; -#X connect 3 4 4 0; -#X connect 4 0 0 0; -#X connect 5 0 3 0; -#X connect 6 0 0 0; -#X connect 7 0 0 0; -#X connect 8 0 0 0; -#X connect 9 0 0 0; -#X restore 17 156 pd init; -#X obj 18 321 s \$0-out; -#X obj 18 272 r \$0-in; -#X obj 393 317 s \$0-in; -#N canvas 565 515 355 193 compute 0; -#X obj 27 29 inlet; -#X obj 27 127 s \$0-in; -#X msg 27 96 bang \, get massesPos \, get linksPos; -#X obj 27 66 gemhead; -#X obj 160 36 bng 15 250 50 0 empty empty empty 0 -6 0 8 -262144 -1 --1; -#X connect 0 0 3 0; -#X connect 2 0 1 0; -#X connect 3 0 2 0; -#X connect 4 0 2 0; -#X restore 17 209 pd compute; -#X obj 17 181 tgl 15 0 empty empty ON/OFF 25 10 1 10 -262144 -1 -1 -1 1; -#X obj 172 270 r \$0-out; -#X obj 17 126 bng 15 250 50 0 empty empty reset 25 10 1 10 -262144 --1 -1; -#X obj 336 216 gemwin; -#X msg 336 194 0 \, destroy; -#N canvas 363 318 550 319 gemrender 0; -#X obj 48 203 translateXYZ; -#X obj 48 229 sphere 0.1; -#X obj 127 24 inlet; -#X obj 390 13 inlet; -#X obj 303 213 curve 2; -#X obj 125 136 t b f; -#X obj 127 62 unpack f f f f; -#X msg 390 58 \$2 \$3 \$4; -#X msg 462 62 \$5 \$6 \$7; -#X obj 48 167 gemhead; -#X msg 48 143 0; -#X obj 48 115 loadbang; -#X obj 303 147 gemhead; -#X msg 303 123 0; -#X obj 303 95 loadbang; -#X obj 390 86 t b l; -#X connect 0 0 1 0; -#X connect 2 0 6 0; -#X connect 3 0 8 0; -#X connect 3 0 7 0; -#X connect 5 0 9 0; -#X connect 5 1 0 1; -#X connect 6 1 5 0; -#X connect 6 2 0 2; -#X connect 6 3 0 3; -#X connect 7 0 15 0; -#X connect 8 0 4 2; -#X connect 9 0 0 0; -#X connect 10 0 9 0; -#X connect 11 0 10 0; -#X connect 12 0 4 0; -#X connect 13 0 12 0; -#X connect 14 0 13 0; -#X connect 15 0 12 0; -#X connect 15 1 4 1; -#X restore 172 321 pd gemrender; -#X msg 336 129 reset \, create \, 1; -#X obj 18 296 pmpd3d; -#X msg 393 271 forceZ mob -10; -#X msg 393 293 forceZ mob 10; -#X text 392 246 2 Send forces; -#X text 170 132 1 Create window -->; -#X text 170 195 3 And destroy it -->; -#X obj 172 294 route massesPos linksPosNo; -#X text 112 618 01_pmpd3dtest.pd; -#X connect 40 0 50 0; +#X text 15 30 The attributes messages ask the object to output some +of his internal parameters. They can be output by lists for positions +and forces of masses \, or directly send value to a table.; +#X restore 12 499 pd masses_attributes __; +#N canvas 446 66 1236 814 links_attributes 0; +#X obj 5 3 cnv 15 590 15 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 10 108 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 4 87 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X text 7 2 ATTRIBUTES Messages; +#X text 7 87 Lists :; +#X text 170 180 Output all norm of position vector (sqrt(x*x+y*y)) +in a list on outlet No 1; +#X obj 11 222 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 11 335 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 611 26 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 605 5 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X obj 612 140 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 +-66577 0; +#X obj 612 253 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 +-66577 0; +#X msg 14 305 massesSpeedsNormL (\$1); +#X text 402 546 \$1 (optional): Masses Id; +#X text 833 502 \$2 (optional): Masses Id; +#X text 832 489 \$1 : table name; +#X obj 17 584 cnv 15 150 220 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X obj 11 564 cnv 15 590 15 empty empty empty 20 12 0 14 -158509 -66577 +0; +#X msg 22 597 massesPosMean; +#X text 177 613 \$1 : none or Id; +#X msg 23 629 massesPosStd \$1; +#X text 179 644 \$1 : none or Id; +#X text 173 597 output the average position of masses; +#X text 174 628 output the standard deviation of masses position; +#X text 177 674 \$1 : none or Id; +#X text 179 711 \$1 : none or Id; +#X text 173 658 output the average position of masses; +#X text 174 695 output the standard deviation of masses position; +#X text 177 749 \$1 : none or Id; +#X text 178 786 \$1 : none or Id; +#X text 173 733 output the average position of masses; +#X text 173 770 output the standard deviation of masses position; +#X msg 22 696 massesForcesStd \$1; +#X msg 22 658 massesForcesMean; +#X msg 23 735 massesSpeedsMean; +#X msg 22 771 massesSpeedsStd \$1; +#X text 608 5 tables :; +#X text 15 30 The attributes messages ask the object to output some +of his internal parameters. They can be output by lists for positions +and forces of links; +#X text 14 564 statistics (not implemented yet) :; +#X msg 11 119 linksPosL (\$1); +#X text 170 120 Output all links center positions in a list on outlet +No 1; +#X msg 11 138 linksPosXL (\$1); +#X msg 12 157 linksPosYL (\$1); +#X msg 12 176 linksPosZL (\$1); +#X msg 12 195 linksPosNormL (\$1); +#X msg 13 229 linksLengthL (\$1); +#X msg 13 248 linksLengthXL (\$1); +#X msg 13 267 linksLengthYL (\$1); +#X msg 14 286 linksLengthZL (\$1); +#X text 170 152 Output all links x or y in a list on outlet No 1; +#X obj 12 451 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 -66577 +0; +#X msg 14 458 linksLengthSpeedL (\$1); +#X msg 14 477 linksLengthSpeedXL (\$1); +#X msg 14 496 linksLengthSpeedYL (\$1); +#X msg 15 515 linksLengthSpeedZL (\$1); +#X msg 15 534 massesSpeedsSpeedNormL (\$1); +#X text 184 459 Output all link lenght deformation speed in a list +on outlet No 1; +#X text 184 495 Output all link lenght deformation speed in x or y +in a list on outlet No 1; +#X text 214 523 Output all norm of link lenght deformation speed vector +(sqrt(x*x+y*y)) in a list on outlet No 1; +#X msg 615 38 linksPosT \$1 (\$2); +#X msg 616 57 linksPosXT \$1 (\$2); +#X msg 615 76 linksPosYT \$1 (\$2); +#X msg 615 95 linksPosZT \$1 (\$2); +#X msg 614 114 linksPosNormT \$1 (\$2); +#X msg 616 152 linksLengthT \$1 (\$2); +#X msg 615 228 linksLengthNormT \$1 (\$2); +#X obj 612 369 cnv 15 150 110 empty empty empty 20 12 0 14 -233017 +-66577 0; +#X text 808 377 Output all link lenght deformation speed in a list +on outlet No 1; +#X text 808 413 Output all link lenght deformation speed in x or y +in a list on outlet No 1; +#X text 838 441 Output all norm of link lenght deformation speed vector +(sqrt(x*x+y*y)) in a list on outlet No 1; +#X msg 614 376 linksLengthSpeedT \$1 (\$2); +#X msg 613 395 linksLengthSpeedXT \$1 (\$2); +#X msg 614 414 linksLengthSpeedYT \$1 (\$2); +#X msg 615 433 linksLengthSpeedZT \$1 (\$2); +#X msg 615 452 massesSpeedsSpeedNormT \$1 (\$2); +#X text 794 35 Output all links center positions in a table (\$1); +#X text 794 67 Output all links center positions x or y in a table +(\$1); +#X text 793 105 Output all norm of links center positions vector (sqrt(x*x+y*y)) +in a table (\$1); +#X text 170 233 Output all link lenght in a list on outlet No 1; +#X text 172 266 Output all link lenght in x or y in a list on outlet +No 1; +#X text 181 294 Output all norm of link lenght vector (sqrt(x*x+y*y)) +in a list on outlet No 1; +#X text 171 349 Output all link center speed in a list on outlet No +1; +#X text 175 379 Output all link center speed in x or y in a list on +outlet No 1; +#X text 194 415 Output all norm of link center speed vector (sqrt(x*x+y*y)) +in a list on outlet No 1; +#X text 797 154 Output all link length in a list on outlet No 1; +#X text 796 181 Output all link length in x or y in a list on outlet +No 1; +#X text 799 218 Output all norm of link length vector (sqrt(x*x+y*y)) +in a list on outlet No 1; +#X msg 16 342 linksPosSpeedL (\$1); +#X msg 16 361 linksPosSpeedXL (\$1); +#X msg 15 380 linksPosSpeedYL (\$1); +#X msg 15 399 linksPosSpeedsZL (\$1); +#X msg 14 418 linksPosSpeedNormL (\$1); +#X msg 616 259 linksPosSpeedT \$1 (\$2); +#X msg 616 335 linksPosSpeedNormT \$1 (\$2); +#X text 795 264 Output all link center speed in a list on outlet No +1; +#X text 799 294 Output all link center speed in x or y in a list on +outlet No 1; +#X text 815 326 Output all norm of link center speed vector (sqrt(x*x+y*y)) +in a list on outlet No 1; +#X msg 616 171 linksLengthXT \$1 (\$2); +#X msg 616 190 linksLengthYT \$1 (\$2); +#X msg 616 209 linksLengthZT \$1 (\$2); +#X msg 616 278 linksPosSpeedXT \$1 (\$2); +#X msg 616 297 linksPosSpeedYT \$1 (\$2); +#X msg 616 316 linksPosSpeedZT \$1 (\$2); +#X restore 12 519 pd links_attributes ___; +#X connect 36 0 46 0; +#X connect 39 0 38 0; +#X connect 40 0 52 0; +#X connect 41 0 34 0; #X connect 43 0 42 0; -#X connect 44 0 56 0; -#X connect 45 0 38 0; -#X connect 47 0 46 0; -#X connect 49 0 46 0; -#X connect 50 0 39 0; -#X connect 51 0 41 0; -#X connect 52 0 41 0; -#X connect 56 0 48 0; -#X connect 56 1 48 1; +#X connect 45 0 42 0; +#X connect 46 0 35 0; +#X connect 47 0 37 0; +#X connect 48 0 37 0; +#X connect 52 0 44 0; +#X connect 52 1 44 1; diff --git a/pmpd3d.c b/pmpd3d.c index 5d9204c..f67b225 100644 --- a/pmpd3d.c +++ b/pmpd3d.c @@ -31,8 +31,8 @@ #include "stdio.h" #include "math.h" -#define nb_max_link 100000 -#define nb_max_mass 100000 +#define nb_max_link 10000000 +#define nb_max_mass 1000000 #define max(a,b) ( ((a) > (b)) ? (a) : (b) ) #define min(a,b) ( ((a) < (b)) ? (a) : (b) ) @@ -68,6 +68,7 @@ typedef struct _mass { t_float forceY; t_float forceZ; t_float D2; + t_float D2offset; int num; } foo; @@ -92,7 +93,7 @@ typedef struct _link { t_float D_L; // longeur du tabeau D } foo1 ; -typedef struct _pmpd3d2d { +typedef struct _pmpd3d { t_object x_obj; struct _link link[nb_max_link]; struct _mass mass[nb_max_mass]; @@ -148,7 +149,10 @@ void pmpd3d_infosL(t_pmpd3d *x) post("number, Id, mobile, mass, Damping, positionX Y Z, speedX Y Z, forcesX Y Z"); for(i=0; i < x->nb_mass; i++) { - post("masse %i: %s, %d, %f, %f, %f, %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].posZ, x->mass[i].speedX, x->mass[i].speedY, x->mass[i].speedZ, x->mass[i].forceX, x->mass[i].forceY, x->mass[i].forceZ ); + post("masse %i: %s, %d, %f, %f, %f, %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].posZ, x->mass[i].speedX, x->mass[i].speedY, x->mass[i].speedZ, \ + x->mass[i].forceX, x->mass[i].forceY, x->mass[i].forceZ ); } post("list of link"); @@ -158,13 +162,20 @@ void pmpd3d_infosL(t_pmpd3d *x) 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); + 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, %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, x->link[i].VZ); + post("tLink %i: %s, %i, %i, %f, %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, x->link[i].VZ); 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); + 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; } } @@ -173,7 +184,7 @@ void pmpd3d_infosL(t_pmpd3d *x) void pmpd3d_bang(t_pmpd3d *x) { // this part is doing all the PM - t_float F, L, Lx,Ly, Lz, Fx, Fy, Fz, tmpX, tmpY, tmpZ; + t_float F, L, Lx,Ly, Lz, Fx, Fy, Fz, tmpX, tmpY, tmpZ, speed; t_int i; // post("bang"); @@ -190,7 +201,8 @@ void pmpd3d_bang(t_pmpd3d *x) x->mass[i].posX += x->mass[i].speedX ; x->mass[i].posY += x->mass[i].speedY ; x->mass[i].posZ += x->mass[i].speedZ ; - 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) || (x->mass[i].posZ < x->minZ) || (x->mass[i].posZ > x->maxZ) ) + 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) || (x->mass[i].posZ < x->minZ) || (x->mass[i].posZ > x->maxZ) ) { tmpX = min(x->maxX,max(x->minX,x->mass[i].posX)); tmpY = min(x->maxY,max(x->minY,x->mass[i].posY)); @@ -205,6 +217,12 @@ void pmpd3d_bang(t_pmpd3d *x) x->mass[i].forceX = -x->mass[i].D2 * x->mass[i].speedX; x->mass[i].forceY = -x->mass[i].D2 * x->mass[i].speedY; x->mass[i].forceZ = -x->mass[i].D2 * x->mass[i].speedZ; + speed = sqrt(x->mass[i].speedX * x->mass[i].speedX + x->mass[i].speedY * x->mass[i].speedY + x->mass[i].speedZ * x->mass[i].speedZ); + 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); + x->mass[i].forceZ += x->mass[i].D2offset * (x->mass[i].speedZ/speed); + } } for (i=0; inb_link; i++) @@ -250,7 +268,8 @@ void pmpd3d_bang(t_pmpd3d *x) } } -void pmpd3d_mass(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv)//, t_symbol *Id, t_float mobile, t_float M, t_float posX, t_float posY, t_float posZ ) +void pmpd3d_mass(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) +//, t_symbol *Id, t_float mobile, t_float M, t_float posX, t_float posY, t_float posZ ) { // add a mass : Id, invM speedX posX if(argv[0].a_type == A_SYMBOL) { @@ -269,6 +288,9 @@ void pmpd3d_mass(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv)//, t_symbol * x->mass[x->nb_mass].forceY = 0; x->mass[x->nb_mass].forceZ = 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 ); @@ -286,7 +308,8 @@ void pmpd3d_create_link(t_pmpd3d *x, t_symbol *Id, int mass1, int mass2, t_float 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) + sqr(x->mass[mass1].posZ - x->mass[mass2].posZ)); + 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) + sqr(x->mass[mass1].posZ - x->mass[mass2].posZ)); x->link[x->nb_link].Pow = Pow; x->link[x->nb_link].Lmin = Lmin; x->link[x->nb_link].Lmax = Lmax; @@ -350,7 +373,8 @@ void pmpd3d_link(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) { 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 + 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 pmpd3d_create_link(x, Id, i, j, K, D, Pow, Lmin, Lmax, 0); } } @@ -421,7 +445,8 @@ void pmpd3d_tLink(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) { 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 + 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 { pmpd3d_create_link(x, Id, i, j, K, D, Pow, Lmin, Lmax, 1); x->link[x->nb_link-1].VX = vecteurX; @@ -494,7 +519,8 @@ void pmpd3d_tabLink(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) { 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 + 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 { pmpd3d_create_link(x, Id, i, j, 1, 1, 1, 0, 1000000, 2); x->link[x->nb_link-1].arrayK = arrayK; @@ -1018,6 +1044,35 @@ void pmpd3d_setD2(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) } } +void pmpd3d_setD2offset(t_pmpd3d *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 pmpd3d_setSpeedX(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) { int tmp, i; @@ -2208,9 +2263,9 @@ void pmpd3d_massesPosL(t_pmpd3d *x) 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); - SETFLOAT(&(pos_list[2*i+2]),x->mass[i].posZ); + SETFLOAT(&(pos_list[3*i] ),x->mass[i].posX); + SETFLOAT(&(pos_list[3*i+1]),x->mass[i].posY); + SETFLOAT(&(pos_list[3*i+2]),x->mass[i].posZ); } outlet_anything(x->main_outlet, gensym("massesPosL"),3*x->nb_mass , pos_list); } @@ -2222,9 +2277,9 @@ void pmpd3d_massesForcesL(t_pmpd3d *x) 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); - SETFLOAT(&(pos_list[2*i+2]),x->mass[i].forceZ); + SETFLOAT(&(pos_list[3*i] ),x->mass[i].forceX); + SETFLOAT(&(pos_list[3*i+1]),x->mass[i].forceY); + SETFLOAT(&(pos_list[3*i+2]),x->mass[i].forceZ); } outlet_anything(x->main_outlet, gensym("massesForcesL"),3*x->nb_mass , pos_list); } @@ -2236,9 +2291,9 @@ void pmpd3d_massesSpeedsL(t_pmpd3d *x) 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); - SETFLOAT(&(pos_list[2*i+2]),x->mass[i].speedZ); + SETFLOAT(&(pos_list[3*i] ),x->mass[i].speedX); + SETFLOAT(&(pos_list[3*i+1]),x->mass[i].speedY); + SETFLOAT(&(pos_list[3*i+2]),x->mass[i].speedZ); } outlet_anything(x->main_outlet, gensym("massesSpeedsL"),3*x->nb_mass , pos_list); } @@ -2993,59 +3048,1370 @@ void pmpd3d_massesSpeedsStd(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) outlet_anything(x->main_outlet, gensym("massesSpeedsStd"),4 , std_out); } -void pmpd3d_grabMass(t_pmpd3d *x, t_float posX, t_float posY, t_float posZ, t_float grab) +// -------------------------------------------- + +void pmpd3d_linksPosXL(t_pmpd3d *x) { - t_float dist, tmp; - t_int i; - - if (grab == 0) - x->grab=0; - if ((x->grab == 0)&(grab == 1)&(x->nb_mass > 0)) + 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 pmpd3d_linksLengthXL(t_pmpd3d *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 pmpd3d_linksPosSpeedXL(t_pmpd3d *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 pmpd3d_linksLengthSpeedXL(t_pmpd3d *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 pmpd3d_linksPosXT(t_pmpd3d *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) ) { - x->grab=1; - x->grab_nb= 0; - dist = sqr(x->mass[0].posX - posX) + sqr(x->mass[0].posY - posY) + sqr(x->mass[0].posZ - posZ); - for (i=1; inb_mass; i++) + 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 pmpd3d_linksLengthXT(t_pmpd3d *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 { - tmp = sqr(x->mass[i].posX - posX) + sqr(x->mass[i].posY - posY) + sqr(x->mass[i].posZ - posZ); - if (tmp < dist) + int taille_max = x->nb_link; + taille_max = min(taille_max, vecsize); + for (i=0; i < taille_max ; i++) { - dist = tmp; - x->grab_nb= i; + vec[i].w_float = x->link[i].mass2->posX - x->link[i].mass1->posX; } + garray_redraw(a); } } - if (x->grab == 1) + else + if ( (argc==2) && (argv[0].a_type == A_SYMBOL) && (argv[1].a_type == A_SYMBOL) ) { - x->mass[x->grab_nb].posX = posX; - x->mass[x->grab_nb].posY = posY; - x->mass[x->grab_nb].posZ = posZ; + 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 *pmpd3d_new() +void pmpd3d_linksPosSpeedXT(t_pmpd3d *x, t_symbol *s, int argc, t_atom *argv) { - t_pmpd3d *x = (t_pmpd3d *)pd_new(pmpd3d_class); + 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); + } + } +} - pmpd3d_reset(x); - - x->main_outlet=outlet_new(&x->x_obj, 0); - // x->info_outlet=outlet_new(&x->x_obj, 0); // TODO +void pmpd3d_linksLengthSpeedXT(t_pmpd3d *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); + } + } +} - return (void *)x; +void pmpd3d_linksPosYL(t_pmpd3d *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 pmpd3d_setup(void) +void pmpd3d_linksLengthYL(t_pmpd3d *x) { - pmpd3d_class = class_new(gensym("pmpd3d"), - (t_newmethod)pmpd3d_new, - 0, sizeof(t_pmpd3d),CLASS_DEFAULT, 0); + int i; + t_atom pos_list[x->nb_link]; - class_addbang(pmpd3d_class, pmpd3d_bang); - class_addmethod(pmpd3d_class, (t_method)pmpd3d_reset, gensym("reset"), 0); + 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 pmpd3d_linksPosSpeedYL(t_pmpd3d *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 pmpd3d_linksLengthSpeedYL(t_pmpd3d *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 pmpd3d_linksPosYT(t_pmpd3d *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 pmpd3d_linksLengthYT(t_pmpd3d *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 pmpd3d_linksPosSpeedYT(t_pmpd3d *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 pmpd3d_linksLengthSpeedYT(t_pmpd3d *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 pmpd3d_linksPosZL(t_pmpd3d *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->posZ + x->link[i].mass2->posZ)/2); + } + outlet_anything(x->main_outlet, gensym("linksPosZL"),x->nb_link , pos_list); +} + +void pmpd3d_linksLengthZL(t_pmpd3d *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->posZ - x->link[i].mass1->posZ); + } + outlet_anything(x->main_outlet, gensym("linksLengthZL"),x->nb_link , pos_list); +} + +void pmpd3d_linksPosSpeedZL(t_pmpd3d *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->speedZ + x->link[i].mass2->speedZ)/2); + } + outlet_anything(x->main_outlet, gensym("linksPosSpeedZL"),x->nb_link , pos_list); +} + +void pmpd3d_linksLengthSpeedZL(t_pmpd3d *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->speedZ - x->link[i].mass1->speedZ); + } + outlet_anything(x->main_outlet, gensym("linksLengthSpeedZL"),x->nb_link , pos_list); +} + +void pmpd3d_linksPosZT(t_pmpd3d *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->posZ + x->link[i].mass2->posZ)/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->posZ + x->link[j].mass2->posZ)/2; + i++; + } + j++; + } + garray_redraw(a); + } + } +} + +void pmpd3d_linksLengthZT(t_pmpd3d *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->posZ - x->link[i].mass1->posZ; + } + 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->posZ - x->link[j].mass1->posZ; + i++; + } + j++; + } + garray_redraw(a); + } + } +} + +void pmpd3d_linksPosSpeedZT(t_pmpd3d *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->speedZ + x->link[i].mass2->speedZ)/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->speedZ + x->link[j].mass2->speedZ)/2; + i++; + } + j++; + } + garray_redraw(a); + } + } +} + +void pmpd3d_linksLengthSpeedZT(t_pmpd3d *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->speedZ - x->link[i].mass1->speedZ; + } + 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->speedZ - x->link[j].mass1->speedZ; + i++; + } + j++; + } + garray_redraw(a); + } + } +} + +void pmpd3d_linksPosL(t_pmpd3d *x) +{ + int i; + t_atom pos_list[3*x->nb_link]; + + for (i=0; i < x->nb_link; i++) + { + SETFLOAT(&(pos_list[3*i] ),(x->link[i].mass2->posX + x->link[i].mass1->posX)/2); + SETFLOAT(&(pos_list[3*i+1]),(x->link[i].mass2->posY + x->link[i].mass1->posY)/2); + SETFLOAT(&(pos_list[3*i+2]),(x->link[i].mass2->posZ + x->link[i].mass1->posZ)/2); + } + outlet_anything(x->main_outlet, gensym("linksPosL"),3*x->nb_link , pos_list); +} + +void pmpd3d_linksLengthL(t_pmpd3d *x) +{ + int i; + t_atom pos_list[3*x->nb_link]; + + for (i=0; i < x->nb_link; i++) + { + SETFLOAT(&(pos_list[3*i] ),x->link[i].mass2->posX - x->link[i].mass1->posX); + SETFLOAT(&(pos_list[3*i+1]),x->link[i].mass2->posY - x->link[i].mass1->posY); + SETFLOAT(&(pos_list[3*i+2]),x->link[i].mass2->posZ - x->link[i].mass1->posZ); + } + outlet_anything(x->main_outlet, gensym("linksLengthL"),3*x->nb_link , pos_list); +} + +void pmpd3d_linksPosSpeedL(t_pmpd3d *x) +{ + int i; + t_atom pos_list[3*x->nb_link]; + + for (i=0; i < x->nb_link; i++) + { + SETFLOAT(&(pos_list[3*i] ),(x->link[i].mass2->speedX + x->link[i].mass1->speedX)/2); + SETFLOAT(&(pos_list[3*i+1]),(x->link[i].mass2->speedY + x->link[i].mass1->speedY)/2); + SETFLOAT(&(pos_list[3*i+2]),(x->link[i].mass2->speedZ + x->link[i].mass1->speedZ)/2); + } + outlet_anything(x->main_outlet, gensym("linksPosSpeedL"),3*x->nb_link , pos_list); +} + +void pmpd3d_linksLengthSpeedL(t_pmpd3d *x) +{ + int i; + t_atom pos_list[3*x->nb_link]; + + for (i=0; i < x->nb_link; i++) + { + SETFLOAT(&(pos_list[3*i] ),x->link[i].mass2->speedX - x->link[i].mass1->speedX); + SETFLOAT(&(pos_list[3*i+1]),x->link[i].mass2->speedY - x->link[i].mass1->speedY); + SETFLOAT(&(pos_list[3*i+2]),x->link[i].mass2->speedZ - x->link[i].mass1->speedZ); + } + outlet_anything(x->main_outlet, gensym("linksLengthSpeedL"),3*x->nb_link , pos_list); +} + +void pmpd3d_linksPosT(t_pmpd3d *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)/3); + for (i=0; i < taille_max ; i++) + { + vec[3*i ].w_float = (x->link[i].mass2->posX + x->link[i].mass1->posX)/2; + vec[3*i+1].w_float = (x->link[i].mass2->posY + x->link[i].mass1->posY)/2; + vec[3*i+2].w_float = (x->link[i].mass2->posZ + x->link[i].mass1->posZ)/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->mass[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++; + vec[i].w_float = (x->link[j].mass2->posZ + x->link[j].mass1->posZ)/2; + i++; + } + j++; + } + garray_redraw(a); + } + } +} + +void pmpd3d_linksLengthT(t_pmpd3d *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)/3); + for (i=0; i < taille_max ; i++) + { + vec[3*i ].w_float = x->link[i].mass2->posX - x->link[i].mass1->posX; + vec[3*i+1].w_float = x->link[i].mass2->posY - x->link[i].mass1->posY; + vec[3*i+2].w_float = x->link[i].mass2->posZ - x->link[i].mass1->posZ; + } + 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->mass[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++; + vec[i].w_float = x->link[j].mass2->posZ + x->link[j].mass1->posZ; + i++; + } + j++; + } + garray_redraw(a); + } + } +} + +void pmpd3d_linksPosSpeedT(t_pmpd3d *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)/3); + for (i=0; i < taille_max ; i++) + { + vec[3*i ].w_float = (x->link[i].mass2->speedX + x->link[i].mass1->speedX)/2; + vec[3*i+1].w_float = (x->link[i].mass2->speedY + x->link[i].mass1->speedY)/2; + vec[3*i+2].w_float = (x->link[i].mass2->speedZ + x->link[i].mass1->speedZ)/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->mass[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++; + vec[i].w_float = (x->link[j].mass2->speedZ + x->link[j].mass1->speedZ)/2; + i++; + } + j++; + } + garray_redraw(a); + } + } +} + +void pmpd3d_linksLengthSpeedT(t_pmpd3d *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)/3); + for (i=0; i < taille_max ; i++) + { + vec[3*i ].w_float = x->link[i].mass2->speedX - x->link[i].mass1->speedX; + vec[3*i+1].w_float = x->link[i].mass2->speedY - x->link[i].mass1->speedY; + vec[3*i+2].w_float = x->link[i].mass2->speedZ - x->link[i].mass1->speedZ; + } + 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->mass[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++; + vec[i].w_float = x->link[j].mass2->speedZ + x->link[j].mass1->speedZ; + i++; + } + j++; + } + garray_redraw(a); + } + } +} + +void pmpd3d_linksPosNormL(t_pmpd3d *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) + \ + sqr((x->link[i].mass1->posZ + x->link[i].mass2->posZ)/2) )); + } + outlet_anything(x->main_outlet, gensym("linksPosNormL"),x->nb_link , pos_list); +} + +void pmpd3d_linksLengthNormL(t_pmpd3d *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) + \ + sqr(x->link[i].mass2->posZ - x->link[i].mass1->posZ) )); + } + outlet_anything(x->main_outlet, gensym("linksLengthNormL"),x->nb_link , pos_list); +} + +void pmpd3d_linksPosSpeedNormL(t_pmpd3d *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) + \ + sqr((x->link[i].mass1->speedZ + x->link[i].mass2->speedZ)/2) )); + } + outlet_anything(x->main_outlet, gensym("linksPosSpeedNormL"),x->nb_link , pos_list); +} + +void pmpd3d_linksLengthSpeedNormL(t_pmpd3d *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) + \ + sqr(x->link[i].mass2->speedZ - x->link[i].mass1->speedZ) )); + } + outlet_anything(x->main_outlet, gensym("linksLengthSpeedNormL"),x->nb_link , pos_list); +} + +void pmpd3d_linksPosNormT(t_pmpd3d *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) + \ + sqr((x->link[i].mass1->posZ + x->link[i].mass2->posZ)/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) + \ + sqr((x->link[j].mass1->posZ + x->link[j].mass2->posZ)/2) ); + i++; + } + j++; + } + garray_redraw(a); + } + } +} + +void pmpd3d_linksLengthNormT(t_pmpd3d *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) + \ + sqr(x->link[i].mass2->posZ - x->link[i].mass1->posZ) ); + } + 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) + \ + sqr(x->link[j].mass2->posZ - x->link[j].mass1->posZ) ); + i++; + } + j++; + } + garray_redraw(a); + } + } +} + +void pmpd3d_linksPosSpeedNormT(t_pmpd3d *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) + \ + sqr((x->link[i].mass1->speedZ + x->link[i].mass2->speedZ)/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) + \ + sqr((x->link[j].mass1->speedZ + x->link[j].mass2->speedZ)/2) ); + i++; + } + j++; + } + garray_redraw(a); + } + } +} + +void pmpd3d_linksLengthSpeedNormT(t_pmpd3d *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) + \ + sqr(x->link[i].mass2->speedZ - x->link[i].mass1->speedZ) ); + } + 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) + \ + sqr(x->link[j].mass2->speedZ - x->link[j].mass1->speedZ) ); + i++; + } + j++; + } + garray_redraw(a); + } + } +} + +//---------------------------------------------- + +void pmpd3d_grabMass(t_pmpd3d *x, t_float posX, t_float posY, t_float posZ, 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) + sqr(x->mass[0].posZ - posZ); + for (i=1; inb_mass; i++) + { + tmp = sqr(x->mass[i].posX - posX) + sqr(x->mass[i].posY - posY) + sqr(x->mass[i].posZ - posZ); + 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; + x->mass[x->grab_nb].posZ = posZ; + } +} + +void *pmpd3d_new() +{ + t_pmpd3d *x = (t_pmpd3d *)pd_new(pmpd3d_class); + + pmpd3d_reset(x); + + x->main_outlet=outlet_new(&x->x_obj, 0); + // x->info_outlet=outlet_new(&x->x_obj, 0); // TODO + + return (void *)x; +} + +void pmpd3d_setup(void) +{ + pmpd3d_class = class_new(gensym("pmpd3d"), + (t_newmethod)pmpd3d_new, + 0, sizeof(t_pmpd3d),CLASS_DEFAULT, 0); + + class_addbang(pmpd3d_class, pmpd3d_bang); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_reset, gensym("reset"), 0); class_addmethod(pmpd3d_class, (t_method)pmpd3d_infosL, gensym("infosL"), 0); class_addmethod(pmpd3d_class, (t_method)pmpd3d_infosL, gensym("print"), 0); -// class_addmethod(pmpd3d_class, (t_method)pmpd3d_mass, gensym("mass"), A_DEFSYMBOL, A_DEFFLOAT, A_DEFFLOAT, A_DEFFLOAT, A_DEFFLOAT, A_DEFFLOAT, 0); class_addmethod(pmpd3d_class, (t_method)pmpd3d_mass, gensym("mass"), A_GIMME, 0); class_addmethod(pmpd3d_class, (t_method)pmpd3d_link, gensym("link"), A_GIMME, 0); class_addmethod(pmpd3d_class, (t_method)pmpd3d_tabLink, gensym("tabLink"), A_GIMME, 0); @@ -3077,6 +4443,7 @@ void pmpd3d_setup(void) class_addmethod(pmpd3d_class, (t_method)pmpd3d_setFixed, gensym("setFixed"), A_GIMME, 0); class_addmethod(pmpd3d_class, (t_method)pmpd3d_setMobile, gensym("setMobile"), A_GIMME, 0); class_addmethod(pmpd3d_class, (t_method)pmpd3d_setD2, gensym("setDEnv"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_setD2offset, gensym("setDEnvOffset"), A_GIMME, 0); class_addmethod(pmpd3d_class, (t_method)pmpd3d_setSpeed, gensym("setSpeed"), A_GIMME, 0); class_addmethod(pmpd3d_class, (t_method)pmpd3d_setSpeedX, gensym("setSpeedX"), A_GIMME, 0); class_addmethod(pmpd3d_class, (t_method)pmpd3d_setSpeedY, gensym("setSpeedY"), A_GIMME, 0); @@ -3127,15 +4494,68 @@ void pmpd3d_setup(void) class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesSpeedsNormT, gensym("massesSpeedsNormT"), A_GIMME, 0); class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesForcesNormT, gensym("massesForcesNormT"), A_GIMME, 0); - class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesPosMean, gensym("massesPosMean"), A_GIMME, 0); - class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesPosStd, gensym("massesPosStd"),A_GIMME, 0); - class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesForcesMean,gensym("massesForecesMean"), A_GIMME, 0); - class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesForcesStd, gensym("massesForcesStd"),A_GIMME, 0); - class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesSpeedsMean,gensym("massesSpeedsMean"), A_GIMME, 0); - class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesSpeedsStd, gensym("massesSpeedsStd"),A_GIMME, 0); - // class_addmethod(pmpd3d_class, (t_method)pmpd3d_linkLMean, gensym("linkLMean"),A_GIMME, 0); - // class_addmethod(pmpd3d_class, (t_method)pmpd3d_linkLStd, gensym("linkLStd"),A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesPosMean, gensym("massesPosMean"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesPosStd, gensym("massesPosStd"),A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesForcesMean, gensym("massesForecesMean"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesForcesStd, gensym("massesForcesStd"),A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesSpeedsMean, gensym("massesSpeedsMean"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_massesSpeedsStd, gensym("massesSpeedsStd"),A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosXL, gensym("linksPosXL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthXL, gensym("linksLengthXL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosSpeedXL, gensym("linksPosSpeedXL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthSpeedXL, gensym("linksLengthXL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosXT, gensym("linksPosXT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthXT, gensym("linksLengthXT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosSpeedXT, gensym("linksPosSpeedXT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthSpeedXT, gensym("linksLengthSpeedXT"), A_GIMME, 0); + + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosYL, gensym("linksPosYL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthYL, gensym("linksLengthYL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosSpeedYL, gensym("linksPosSpeedYL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthSpeedYL, gensym("linksLengthYL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosYT, gensym("linksPosYT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthYT, gensym("linksLengthYT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosSpeedYT, gensym("linksPosSpeedYT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthSpeedYT, gensym("linksLengthSpeedYT"), A_GIMME, 0); + + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosZL, gensym("linksPosZL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthZL, gensym("linksLengthZL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosSpeedZL, gensym("linksPosSpeedZL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthSpeedZL, gensym("linksLengthZL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosZT, gensym("linksPosZT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthZT, gensym("linksLengthZT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosSpeedZT, gensym("linksPosSpeedZT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthSpeedZT, gensym("linksLengthSpeedZT"), A_GIMME, 0); + + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosL, gensym("linksPosL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthL, gensym("linksLengthL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosSpeedL, gensym("linksPosSpeedL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthSpeedL, gensym("linksLengthL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosT, gensym("linksPosT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthT, gensym("linksLengthT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosSpeedT, gensym("linksPosSpeedT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthSpeedT, gensym("linksLengthSpeedT"), A_GIMME, 0); + + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosNormL, gensym("linksPosNormL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthNormL, gensym("linksLengthNormL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosSpeedNormL, gensym("linksPosSpeedNormL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthSpeedNormL, gensym("linksLengthNormL"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosNormT, gensym("linksPosNormT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthNormT, gensym("linksLengthNormT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosSpeedNormT, gensym("linksPosSpeedNormT"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthSpeedNormT, gensym("linksLengthSpeedNormT"), A_GIMME, 0); + +/* class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosMean, gensym("linksPosMean"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthMean, gensym("linksLengthMean"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosSpeedMean, gensym("linksPosSpeedMean"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthSpeedMean, gensym("linksLengthSpeedMean"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosStd, gensym("linksPosStd"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthStd, gensym("linksLengthStd"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksPosSpeedStd, gensym("linksPosSpeedStd"), A_GIMME, 0); + class_addmethod(pmpd3d_class, (t_method)pmpd3d_linksLengthSpeedStd0, gensym("linksLengthSpeedStd"), A_GIMME, 0); +*/ + class_addmethod(pmpd3d_class, (t_method)pmpd3d_grabMass, gensym("grabMass"), A_DEFFLOAT, A_DEFFLOAT, A_DEFFLOAT, A_DEFFLOAT, 0); } -- cgit v1.2.1