Algo

Pensez vous quil y ait possiblité de simplifier cette routine???

Recherche :

Mot : Pseudo : Filtrer
Bas de page
Auteur	Sujet : Pensez vous quil y ait possiblité de simplifier cette routine???

red faction

pour un intro vraiment tres petite (256b), j'ai besoin que cette routine prenne le moins de place possible et j'aimerai savoir s'il y a possiblilité d'optimiser (en taille) ce calcul:

Code :

l1=abs(1/dx);
l2=abs(1/dy);
l3=abs(1/dz);
l=min(min(max(l1,l2),max(l2,l3)),max(l1,l3));
ix=dx*l;
iy=dy*l;
iz=dz*l;

Publicité

_darkalt3_

Proctopathe

L'idée c'est pas de voir le code assembleur qui est généré à partir de ce code et de le réduire si c'est possible ?

---------------
Töp of the plöp

red faction

L'idée c surtout deviter de tout coder en assembleur , se rendre compte que on peut optimiser des lignes en simplifiants les calculs et de tout recoder apres (car en asm tout est souvent tres lié).....

dividee

Tu peux déjà éviter de calculer 1/dx,1/dy,1/dz en divisant par l au lieu de multiplier par l:

Code :

l1=abs(dx);
l2=abs(dy);
l3=abs(dz);
l=min(min(max(l1,l2),max(l2,l3)),max(l1,l3));
ix=dx/l;
iy=dy/l;
iz=dz/l;

Message cité 1 fois

oliv5

Pourquoi ? Parce que !

Si tes nombres sont des entiers, tu peux utiliser certains algo particuliers pour calculer abs, min et max.

Jette un coup d'oeil là : http://graphics.stanford.edu/~seander/bithacks.html

Il y a des algo amusants et peu évident pour faire tout ça.
Avec un peu de chance, ca te permettra de faire un code plus concis.

Taz

bisounours-codeur

ouais enfin si c'est des entiers, 1/x c'est vite vu ...

franceso

dividee a écrit :

Tu peux déjà éviter de calculer 1/dx,1/dy,1/dz en divisant par l au lieu de multiplier par l:

Bonne idée, mais dans ce cas là il faut faire attention aux min/max : les valeurs directes dx, dy et dz ne sont pas rangées dans le même ordre que leurs inverses...

Message cité 1 fois

---------------
TriScale innov

red faction

bon pour info

on a effectivement viré les 1/dx 1/dy et pour le moment ca pose pas trop de problemes au minmax

pour linstant on est a 300bytes a peu pres tout compris (routine de base , vsync and co)

dividee

franceso a écrit :

Bonne idée, mais dans ce cas là il faut faire attention aux min/max : les valeurs directes dx, dy et dz ne sont pas rangées dans le même ordre que leurs inverses...

Effectivement, elles sont rangées dans l'ordre inverse, mais le calcul avec min/max est en fait le calcul de la valeur médiane, et l'inverse de la médiane est égale à la médiane des inverses...

[edit:] Si le nombre de valeur est impair (pour que la médiane soit bien définie).

Message cité 1 fois
Message édité par dividee le 16-03-2007 à 13:24:42

franceso

dividee a écrit :

:jap:

effectivement, j'avais pas prêté attention à ce détail...

---------------
TriScale innov

Publicité

tbp

Solution tangente, vectorisons.
Disons que l'on a a { l1,l2,l3 }, et par rotation b { l2,l3,l1 }, c { l3,l1,l2 }.
Prenons le masque de mask_ab = a > b et mask_ac = a > c, alors le masque mask = mask_ab ^ mask_ac désigne la valeur médiane.
Extraction du masque etc...

SSE ftw.

Message cité 1 fois
Message édité par tbp le 16-03-2007 à 13:52:42

red faction

tbp a écrit :

tu peut detailller un peu sil te plait??

sinon on a aussi une autre routine qui prend de la place :

Code :

cx=cos(t*0.6);
cy=cos(t*0.1);
cz=cos(t*0.4);
sx=sin(t*0.6);
sy=sin(t*0.1);
sz=sin(t*0.4);
dx1=dx*cz-dy*sz;
dy1=dx*sz+dy*cz;
dz1=dz;
dx2=dx1;
dy2=dy1*cx-dz1*sx;
dz2=dy1*sx+dz1*cx;
dx3=dx2*cy-dz2*sy;
dy3=dy2;
dz3=dx2*sy+dz2*cy;

pour le moment on la simplifié par :

Code :

c=sin(t);
s=sin(t);
rotate(dx,dy,c,s);
rotate(dy,dz,c,s);
rotate(dx,dz,c,s);
void rotate(&dx,&dy,c,s)
fx=dx*c-dy*s;
fy=dx*s+dy*c;
dx=fx;
dy=fy;

bon ca fait pas tout fait la meme chose (calcul de c et s) mais bon

Message cité 1 fois

tbp

red faction a écrit :

tu peut detailller un peu sil te plait??

Feignasse.
En plus les types sont mal définis, mais je vais supposer que a) ya que des flottants b) SSE2 n'est pas verbotten.

Avec les 3 flottants alignés (à 16 octets) sur la pile

Code :

40f0ff: andnps (%esp),%xmm0
40f103: movaps %xmm0,%xmm2
40f106: movaps %xmm0,%xmm1
40f109: shufps $0xc9,%xmm0,%xmm2
40f10d: shufps $0xd2,%xmm0,%xmm1
40f111: cmpltps %xmm0,%xmm2
40f115: cmpltps %xmm0,%xmm1
40f119: xorps %xmm1,%xmm2
40f11c: movmskps %xmm2,%edx
40f11f:

... ce qui nous fait en gros 32 octets pour trouver une valeur absolue + mediane.
Reste 2 points à élucider.
Premierement d'ou vient le masque (dans xmm0 au début) pour virer les bits de signe: 2 solutions, soit le génerer avec

Code :

pxor %xmm0,%xmm0 // 4 octets
pcmpeqd %xmm0,%xmm0 // 4 octets
pslld $0x1f,%xmm0 // 5 octets

ou un truc plus lent dans le genre, soit le charger.

Ensuite, on a un masque dans un GPR mais il faut triturer les doublons.
En gros

Code :

uint_t bits = %edx;
uint_t index = 0;
if (bits & 2) index = 1;
if (bits & 4) index = 2;
const float l = (esp, index, 4);

mais rien de bien fringuant me vient à l'esprit, alors une bonne vielle LUT fera l'affaire.
Voilà.

Joel F

Real men use unique_ptr

et gcc 4.0 inventa les intrinsics SSE2
Ca fait qquelques temps qu'on a plus besoin de descendre au niveau ASM pour vectoriser quoique ce soit

Message cité 2 fois

tbp

Joel F a écrit :

et gcc 4.0 inventa les intrinsics SSE2
Ca fait qquelques temps qu'on a plus besoin de descendre au niveau ASM pour vectoriser quoique ce soit

Quelle intéressante notion!
Soit dit en passant
a) j'ai posté un désassemblage parce que le critère était la taille du code généré et incidemment mis à contribution g++ & intrinsics en amont.
b) les intrinsics viennent en droite ligne de chez Intel, historiquement gcc/g++ exposaient le tout via builtins.

red faction

Joel F a écrit :

et gcc 4.0 inventa les intrinsics SSE2
Ca fait qquelques temps qu'on a plus besoin de descendre au niveau ASM pour vectoriser quoique ce soit

sauf que la c un com qui doit tenir en 256b donc pas question de compiler quoi que se soit de haut niveau....

FORUM HardWare.fr

Programmation

Algo

Pensez vous quil y ait possiblité de simplifier cette routine???

Sujets relatifs
[Important !] Vous en pensez quoi ?	Un service de cartes virtuelles ; Vous en pensez quoi ??
[PHP] Simplifier la récupération de données d'un formulaire	Design vs en pensez quoi? + aide CSS
Que pensez vous de mon nouveau site ?	Que pensez vous de mon htaccess pour un site oscommerce ?
[Batch] Que pensez vous de mes modifications de XP	routine lien cellule-bouton-image
Que pensez vous du design de mon site...	que pensez vous de l'érgonomie de mon site
Plus de sujets relatifs à : Pensez vous quil y ait possiblité de simplifier cette routine???

Page générée en 0.069 secondes