Comment puis-je multiplier des opérandes 64 bits et obtenir un résultat 128 bits de manière portable?

Pour x64, je peux utiliser ceci:

{ uint64_t hi, lo; // hi,lo = 64bit x 64bit multiply of c[0] and b[0] __asm__("mulq %3\n\t" : "=d" (hi), "=a" (lo) : "%a" (c[0]), "rm" (b[0]) : "cc" ); a[0] += hi; a[1] += lo; } 

Mais j’aimerais effectuer le même calcul de manière portable. Par exemple pour travailler sur x86.

Si je comprends bien la question, vous voulez une implémentation en C pur et portable de la multiplication à 64 bits, avec une sortie de 128 bits stockée dans deux valeurs de 64 bits. Dans ce cas, cet article prétend avoir ce dont vous avez besoin. Ce code est écrit pour C ++. Il ne faut pas grand chose pour le transformer en code C:

 void mult64to128(uint64_t op1, uint64_t op2, uint64_t *hi, uint64_t *lo) { uint64_t u1 = (op1 & 0xffffffff); uint64_t v1 = (op2 & 0xffffffff); uint64_t t = (u1 * v1); uint64_t w3 = (t & 0xffffffff); uint64_t k = (t >> 32); op1 >>= 32; t = (op1 * v1) + k; k = (t & 0xffffffff); uint64_t w1 = (t >> 32); op2 >>= 32; t = (u1 * op2) + k; k = (t >> 32); *hi = (op1 * op2) + w1 + k; *lo = (t << 32) + w3; } 

Puisque vous avez gcc comme balise, notez que vous pouvez simplement utiliser le type entier 128 bits de gcc :

 typedef unsigned __int128 uint128_t; // ... uint64_t x, y; // ... uint128_t result = (uint128_t)x * y; uint64_t lo = result; uint64_t hi = result >> 64;