Accueil > > > COLLISIONS EN 3D AVEC VISUAL PYTHON
COLLISIONS EN 3D AVEC VISUAL PYTHON
Information sur la source
Description
Ce code fait suite à "COLLISIONS en cascades" déposé fin 2005 sur le site. Il met en oeuvre l'extension VISUAL de PYTHON, très intéressante pour les représentation graphiques en 3D.
Source
- #! /usr/bin/env python
- # -*- coding: Latin-1 -*-
-
- "COLLISIONS_VPYTHON"
- "Ce script est démarqué de l'exemple fourni dans 'GAS' de l'extension VISUAL de PYTHON"
-
- from visual import *
- from time import clock
- from visual.graph import *
- from random import random
- from visual.controls import*
-
- #################################################################################
-
- win=500
- # Valeurs typiques
- L=1.0
- gray=(0.7,0.7,0.7) # couleur des bordures
- Raxes=0.01 # rayon des bordures
- scene=display(title="COLLISIONS_VPYTHON",width=win,height=win,background=(0,1,1),x=0,y=0,range=L/4,center=(L/2.,L/2.,L/2.))
- scene.lights=[vector(1,1,1)]
- scene.ambient=0
- xaxis=curve(pos=[(0,0,0),(L,0,0)],color=gray,radius=Raxes)
- yaxis=curve(pos=[(0,0,0),(0,L,0)],color=gray,radius=Raxes)
- xaxis2=curve(pos=[(0,L,0),(L,L,0)],color=gray,radius=Raxes)
- yaxis2=curve(pos=[(L,0,0),(L,L,0)],color=gray,radius=Raxes)
- dt=0.002
-
- #################################################################################
-
- def boules():
- ""
- global pos,v,radius,boule,N,R,Nsteps,t,dt
- N =10 # valeur à modifier pour augmenter ou diminuer le nombre de boules
- R=0.05
- boule=[]
- colors=[color.red,color.green,color.blue,color.yellow,color.cyan,color.magenta,color.white]
- poslist,vlist,rlist=[],[],[]
- for i in range(N):
- Lmin=R
- Lmax=L-Lmin
- x,y,z,r=R+(Lmax-Lmin)*random(),R+(Lmax-Lmin)*random(),0,R
- boule=boule+[sphere(pos=(x,y,z),radius=r,color=colors[i % 7])]
- theta,phi=pi*random(),2*pi*random()
- vx,vy,vz=sin(theta)*cos(phi),sin(theta)*sin(phi),0
- poslist.append((x,y,z))
- vlist.append((vx,vy,vz))
- rlist.append(r)
- pos=array(poslist)
- v=array(vlist)
- radius=array(rlist)
- t=0.0
- Nsteps=0
- pos=pos+v*(dt/2.)
- time=clock()
-
- ################################################################################
-
- def move_boules():
- ""
- global pos,v,radius,boule,Nsteps,t,dt
- while 1:
- rate(100)
- pos =pos+ v*dt
- r = pos-pos[:,NewAxis]
- rmag = sqrt(add.reduce(r*r,-1))
- hit = less_equal(rmag,radius+radius[:,NewAxis])-identity(N)
- hitlist = sort(nonzero(hit.flat)).tolist()
- # si collisions:
- for ij in hitlist:
- i,j= divmod(ij,N)
- hitlist.remove(j*N+i)
- vtot=v[i]+v[j]
- vi,vj=v[i],v[j]
- ri,rj=boule[i].radius,boule[j].radius
- a=mag(vj-vi)**2
- if a==0: continue
- b=2*dot(pos[i]-pos[j],vj-vi)
- c=mag(pos[i]-pos[j])**2-(ri+rj)**2
- d=b**2-4.*a*c
- if d<0: continue
- deltat=(-b+sqrt(d))/(2.*a) # t-deltat lorsque il y a contact
- pos[i],pos[j]=pos[i]-v[i]*deltat,pos[j]-v[j]*deltat # retour à la configuration de contact
- vcmi,vcmj=v[i]-vtot/2,v[j]-vtot/2
- rrel=norm(pos[j]-pos[i])
- vcmi,vcmj=vcmi-2*dot(vcmi,rrel)*rrel,vcmj-2*dot(vcmj,rrel)*rrel
- v[i],v[j]=vcmi+vtot/2,vcmj+vtot/2
- pos[i],pos[j]=pos[i]+v[i]*deltat,pos[j]+v[j]*deltat
- # rebonds sur les bordures
- outside=less_equal(pos,R+Raxes)
- v1=v*outside
- v=v-v1+abs(v1)
- outside=greater_equal(pos,L-R-Raxes)
- v1=v*outside
- v=v-v1-abs(v1)
- # mise à jour des positions
- for i in range(N):
- boule[i].pos=pos[i]
- Nsteps=Nsteps+1
- t=t+dt
-
- def lancer():
- ""
- boules()
- move_boules()
- lancer()
-
#! /usr/bin/env python
# -*- coding: Latin-1 -*-
"COLLISIONS_VPYTHON"
"Ce script est démarqué de l'exemple fourni dans 'GAS' de l'extension VISUAL de PYTHON"
from visual import *
from time import clock
from visual.graph import *
from random import random
from visual.controls import*
#################################################################################
win=500
# Valeurs typiques
L=1.0
gray=(0.7,0.7,0.7) # couleur des bordures
Raxes=0.01 # rayon des bordures
scene=display(title="COLLISIONS_VPYTHON",width=win,height=win,background=(0,1,1),x=0,y=0,range=L/4,center=(L/2.,L/2.,L/2.))
scene.lights=[vector(1,1,1)]
scene.ambient=0
xaxis=curve(pos=[(0,0,0),(L,0,0)],color=gray,radius=Raxes)
yaxis=curve(pos=[(0,0,0),(0,L,0)],color=gray,radius=Raxes)
xaxis2=curve(pos=[(0,L,0),(L,L,0)],color=gray,radius=Raxes)
yaxis2=curve(pos=[(L,0,0),(L,L,0)],color=gray,radius=Raxes)
dt=0.002
#################################################################################
def boules():
""
global pos,v,radius,boule,N,R,Nsteps,t,dt
N =10 # valeur à modifier pour augmenter ou diminuer le nombre de boules
R=0.05
boule=[]
colors=[color.red,color.green,color.blue,color.yellow,color.cyan,color.magenta,color.white]
poslist,vlist,rlist=[],[],[]
for i in range(N):
Lmin=R
Lmax=L-Lmin
x,y,z,r=R+(Lmax-Lmin)*random(),R+(Lmax-Lmin)*random(),0,R
boule=boule+[sphere(pos=(x,y,z),radius=r,color=colors[i % 7])]
theta,phi=pi*random(),2*pi*random()
vx,vy,vz=sin(theta)*cos(phi),sin(theta)*sin(phi),0
poslist.append((x,y,z))
vlist.append((vx,vy,vz))
rlist.append(r)
pos=array(poslist)
v=array(vlist)
radius=array(rlist)
t=0.0
Nsteps=0
pos=pos+v*(dt/2.)
time=clock()
################################################################################
def move_boules():
""
global pos,v,radius,boule,Nsteps,t,dt
while 1:
rate(100)
pos =pos+ v*dt
r = pos-pos[:,NewAxis]
rmag = sqrt(add.reduce(r*r,-1))
hit = less_equal(rmag,radius+radius[:,NewAxis])-identity(N)
hitlist = sort(nonzero(hit.flat)).tolist()
# si collisions:
for ij in hitlist:
i,j= divmod(ij,N)
hitlist.remove(j*N+i)
vtot=v[i]+v[j]
vi,vj=v[i],v[j]
ri,rj=boule[i].radius,boule[j].radius
a=mag(vj-vi)**2
if a==0: continue
b=2*dot(pos[i]-pos[j],vj-vi)
c=mag(pos[i]-pos[j])**2-(ri+rj)**2
d=b**2-4.*a*c
if d<0: continue
deltat=(-b+sqrt(d))/(2.*a) # t-deltat lorsque il y a contact
pos[i],pos[j]=pos[i]-v[i]*deltat,pos[j]-v[j]*deltat # retour à la configuration de contact
vcmi,vcmj=v[i]-vtot/2,v[j]-vtot/2
rrel=norm(pos[j]-pos[i])
vcmi,vcmj=vcmi-2*dot(vcmi,rrel)*rrel,vcmj-2*dot(vcmj,rrel)*rrel
v[i],v[j]=vcmi+vtot/2,vcmj+vtot/2
pos[i],pos[j]=pos[i]+v[i]*deltat,pos[j]+v[j]*deltat
# rebonds sur les bordures
outside=less_equal(pos,R+Raxes)
v1=v*outside
v=v-v1+abs(v1)
outside=greater_equal(pos,L-R-Raxes)
v1=v*outside
v=v-v1-abs(v1)
# mise à jour des positions
for i in range(N):
boule[i].pos=pos[i]
Nsteps=Nsteps+1
t=t+dt
def lancer():
""
boules()
move_boules()
lancer()
Sources du même auteur
Sources de la même categorie
Commentaires et avis
|
Derniers Blogs
ENUMERABLECOLLECTIONENUMERABLECOLLECTION par Matthieu MEZIL
Prenons le scénario suivant. On utilise MVVM. On a les deux classes suivantes dans le model : public class Child { } public class Parent { private ObservableCollection < Child > _children; public ObservableCollection < Child > Children { get {...
Cliquez pour lire la suite de l'article par Matthieu MEZIL [HS] CHROME 6 + MOI = COUP DE GUEULE ![HS] CHROME 6 + MOI = COUP DE GUEULE ! par JeremyJeanson
Attention, le poste qui suit n'est pas la complainte d'une personne : Qui n'aime pas Chrome. D'un anti Google. D'un développeur qui a un poil énorme dans la main. Ceux qui me fréquentent savent que je change de navigateur favori tous les 2 ou 3 mois afin ...
Cliquez pour lire la suite de l'article par JeremyJeanson [WP7] UTILISER UN WRAPPANEL DANS UNE APPLICATION WINDOWS PHONE 7[WP7] UTILISER UN WRAPPANEL DANS UNE APPLICATION WINDOWS PHONE 7 par Audrey
Lors de la réalisation de ma 2ème application Windows Phone 7, j'ai souhaité utiliser un WrapPanel pour afficher plusieurs photos. Mais le contrôle WrapPanel ne fait pas parti de la liste des contrôles inclus dans le SDK de la version Beta des outils pour...
Cliquez pour lire la suite de l'article par Audrey [WP7] BESOIN D'AVOIR DES DONNéES EN CACHE[WP7] BESOIN D'AVOIR DES DONNéES EN CACHE par Nicolas
Les développeurs ASP.NET ont l'habitude de mettre des données en cache pour éviter de requêter a chaque fois la base de données. Et il est toujours utilie de penser que vos utilisateurs mobiles n'ont pas troujours une super connexion 3G/WIFI et un for...
Cliquez pour lire la suite de l'article par Nicolas [TFS] COMMENT FORCER LA SAISIE D'UN AREA OU ITERATION[TFS] COMMENT FORCER LA SAISIE D'UN AREA OU ITERATION par cyril
Lorsque l'on créé un Work Item dans TFS, il est possible de le classer dans un "area" et dans une "iteration". Dans la plupart des types de projet, un "area" correspond à une catégorie, une "iteration" à un numéro de version. Il est possible de cré...
Cliquez pour lire la suite de l'article par cyril
Forum
RE : PYTHON 3.0RE : PYTHON 3.0 par aera group
Cliquez pour lire la suite par aera group RE : PYTHON 3.0RE : PYTHON 3.0 par xeolin
Cliquez pour lire la suite par xeolin RE : PYTHON 3.0RE : PYTHON 3.0 par aera group
Cliquez pour lire la suite par aera group
Logiciels
uTorrent (2.0.4)UTORRENT (2.0.4)C'est un client BitTorrent très puissant et très performant. Comme son nom l'indique, uTorrent (m... Cliquez pour télécharger uTorrent Bureau de Gestion - ERP Devis Facturation (2.02)BUREAU DE GESTION - ERP DEVIS FACTURATION (2.02)- Version gratuite du 10/06/2010
Le Bureau de Gestion est un logiciel dédié à la gestion de l'en... Cliquez pour télécharger Bureau de Gestion - ERP Devis Facturation 4Videosoft Transfert iPod Mac (3.2.08)4VIDEOSOFT TRANSFERT IPOD MAC (3.2.08)4Videosoft Transfert iPod-Mac caractérise principalement à transférer les fichiers iPod vers Mac.... Cliquez pour télécharger 4Videosoft Transfert iPod Mac 4Videosoft HD Convertisseur (3.3.08)4VIDEOSOFT HD CONVERTISSEUR (3.3.08)Etant le meilleur HD Vidéo Convertisseur, 4Videosoft HD Convertisseur, vous pouvez regarder la vi... Cliquez pour télécharger 4Videosoft HD Convertisseur 4Videosoft Transfert iPad Mac (3.2.08)4VIDEOSOFT TRANSFERT IPAD MAC (3.2.08)4Videosoft Transfert iPad-Mac est un logiciel managérial iPad professionnel qui a des fonctions i... Cliquez pour télécharger 4Videosoft Transfert iPad Mac
|