Discussion :

[ERCO503]

Bonjour,

Existe-t-il des algorithmes ou astuces pour déterminer quelles faces seront atteintes par un rayon partant d'un point p à l'aide d'un mesh (dessiné à l'aide de la fonction patch) contenant une liste de faces et vertices ?

En ce moment, le mieux que j'aie pu faire comme estimation est de boucler itérativement chaque face en y projetant un rayon en son centre et déterminer s'il n'y a pas une autre face qui est atteinte par ce rayon avant d'arriver à cette face, mais cela prend environ 50 secondes pour 18 000 faces en vectorisant une partie. Tenter de tout vectoriser serait trop lourd en mémoire malheureusement.

Merci,

Éric

[Jerome Briot]

Si tu nous montres ton code, nous pourrons essayer de l'optimiser ensemble.

Si tu es sûr d'avoir optimisé au maximum, il faut utiliser le profiler de MATLAB pour déterminer les lignes gourmandes :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

doc profile

A partir de la, tu peux éventuellement essayer d'utiliser des fichiers MEX pour accélérer certaine partie du code.

Pour finir, ce type de calcul peut généralement être parallélisé. Possèdes-tu la Parallel Procesing Toolbox ?

[ERCO503]

Malheureusement, je ne peux pas fournir ce code.

Je voulais plutôt savoir s'il n'existait pas déjà des astuces, des fonctions, algorithmes ou autres qui permettaient d'effectuer cela.

Autrement, en ce moment, mon bout de code le plus long dans ma boucle est celui-ci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

mask = ((s>=0) & (t>=0) & ((s+t)<=1));

Où s et t sont des vecteurs d'environ 18 000 éléments. Cela prend environ 10 secondes au total pour environ 14 000 exécutions.

Merci,

Éric