Discussion :

[uriotcea]

Bonjour
J'essaye de tester l'implémentation d'instruction ss3 avec l'example ci-dessous qui fonctionne bien

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
 
void testSSE2(float* v1, float* v2)
{
  __m128 v1m, v2m;
  memcpy(&v1m, v1, 4* sizeof(float));
  memcpy(&v2m, v2, 4* sizeof(float));
  v1m = _mm_add_ps(v1m, v2m);
  v1m = _mm_sub_ps(v1m, v2m);
  v1m = _mm_mul_ps(v1m, v2m);
  v1m = _mm_div_ps(v1m, v2m);
  v1m = _mm_sqrt_ps(v1m);
  memcpy(v1, &v1m, 4* sizeof(float));
}
 
float v1[4] = { 1, 2, 3, 4};
float v2[4] = { 1, 2, 3, 4};
for(int i = 0; i < 1000; i++)  testSSE2(v1,v2);
std::cout << v1[0] << " " << v1[1] << " " << v1[2] << " " << v1[3] << std::endl;

Maintenant je voudrais passer tout ca en double mais je n'y arrive pas. Je n'arrive pas non plus à savoir si les instruction ss3 acceptent des doubles.

Merci de votre aide

[Goten]

remplace les intrisics *_ps par *_pd (precision double)

(et pour voir une bonne abstraction bien comme il faut avec un DSEL regarde nt2 :€ )

[uriotcea]

Ok merci, ca fonctionne en remplacant aussi __m128 par __m128d
Au passage, il y a t-il moyen de savoir si un processeur accepte ou non les instructions sse3

[Goten]

Ok merci, ca fonctionne en remplacant aussi __m128 par __m128d
Au passage, il y a t-il moyen de savoir si un processeur accepte ou non les instructions sse3

C'est /chiant/, ie c'est os-dependant. Du genre sous *BSD j'irais lire le dmesg ou sysctl, sur linux via /proc/cpuinfo et la solution *générique* via cpuid (l'instruction asm).

Nt2 fait un peu de tout, notamment compiler un programme qui utilise cpuid.

Si ça t'intéresse j'ai un bout de code (en C) qui fait ça.

[uriotcea]

Oui ca m'interresse, pour ma part j'ai trouvé ca, je ne sais pas si c'est compatible

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 
int cpuid;
    asm (    "movl $1, %%eax;"
             "cpuid;"
        "movl %%edx, %0;"
        :"=&r"(cpuid)
        :
        :"%ebx");
    if(cpuid & 0x800000)
        printf("Found MMX capabilities\n");
    if(cpuid & 0x2000000)
        printf("Found SSE capabilities\n");
    if(cpuid & 0x4000000)
        printf("Found SSE2 capabilities\n");
    if(cpuid & 0x1)
        printf("Found SSE3 capabilities\n");
    if(cpuid & 0x200)
        printf("Found SSE3+ capabilities\n");
    if(cpuid & 0x80000)
        printf("Found SSE4.1 capabilities\n");

Au passage (un de plus) je suis à la recheche d'example simple d'utilisation d'instruction sse3

[Joel F]

pas de memcpy pour charger un registre SSE, on utilise _mm_load_ps/pd.
ta memoire doit etre allouée alignée aussi , sinon seg fault.

EDIT: et bien sur http://github.com/MetaScale/nt2 precisement le module (mal/peu documenté) sdk/simd