¿Cómo puedo multiplicar los operandos de 64 bits y obtener el resultado de 128 bits de forma portátil?

Para x64 puedo usar esto:

{ 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; } 

Pero me gustaría realizar el mismo cálculo de forma portátil. Por ejemplo, para trabajar en x86.

Según entiendo la pregunta, desea una implementación C portátil pura de multiplicación de 64 bits, con salida a un valor de 128 bits, almacenada en dos valores de 64 bits. En cuyo caso, este artículo pretende tener lo que necesita. Ese código está escrito para C ++. No hace falta mucho para convertirlo en código 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; } 

Como tiene gcc como etiqueta, tenga en cuenta que solo puede usar el tipo de entero de 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;