From: Gabriel Peyr <nik...@wa...> - 2001-08-10 08:08:11
|
j'ai pas bien compris ta remarque (c'est pas ironique ! j'ai pas compris) = : oui, je peux calculer le time step pour chaque particule, et prendre la = moyenne, comme =E7a, elle sont toutes updat=E9es de la m=EAme mani=E8re. mais bon, je pense que tu as raison, et qu'on doit y gagner pas mal =E0 a= dapter le h pour chaque particule bein, pour mettre =E0 jour le time step, y'a des algo super efficace pour= =E7a, je pense pas que =E7a soit la peine de r=E9inventer la roue, y'a de= s matheux qui ont pass=E9 leur vie l=E0 dessus. lis donc les numerical recipies, tout est expliqu=E9 ! www.nr.com grosso modo, pour le RK4, y'a un truc de base : c'est faire des =E9valuat= ion tous les h et tous les 2h, et calculer la deviation (tr=E9s grosso mod= o) : c'est pas terrible. la meilleur technique, c'est d'utiliser un RK plus =E9lev=E9 (RK5), pour = calculer une erreur locale. ensuite, il suffit de faire les modifs en fonction de cette erreur et d'un certain treshold gabriel >Messsage du 09/08/2001 21:41 >De : <ori...@li...> >A : <ori...@li...> >Copie =E0 : >Objet : [Orion3d-dev] adaptive > > > >sinon, je vais essayer de faire un time step adaptatif, mais l=E0 enco= re, =E7a > risque d'=EAtre assez > >chaud, au niveau de la synchro avec le temps r=E9el, ainsi qu'=E0 cau= se du > nombre de particules =E0 > >prendre en compte (je veux dire : le choix du timestep doit =EAtre fai= t pour > toutes les particule : > >je choisis une moyenne, ou alors je me calle sur le plus mauvais =3F=3F= ) > > > tu le fais pour chaque particule, pasque si tu fais la moyenne ou si tu= te > calle sur la plus > mauvaise c que tu connais le h de toutes les particules.... :) > > perso mon pb est qd mettre a jour le timestep=3F je suis a la recherche= d'une > methode > efficace ne bouffant pas ce qu'on gagne d'un cot=E9, si qq1 a des conse= ils.... > > + > Nith > > > > > > _______________________________________________ > Orion3d-dev mailing list > Ori...@li... > http://lists.sourceforge.net/lists/listinfo/orion3d-dev > |