1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
|
/*
* Test R5900-specific three-operand MADD and MADD1.
*/
#include <stdio.h>
#include <inttypes.h>
#include <assert.h>
int64_t madd(int64_t a, int32_t rs, int32_t rt)
{
int32_t lo = a;
int32_t hi = a >> 32;
int32_t rd;
int64_t r;
__asm__ __volatile__ (
" mtlo %5\n"
" mthi %6\n"
" madd %0, %3, %4\n"
" mflo %1\n"
" mfhi %2\n"
: "=r" (rd), "=r" (lo), "=r" (hi)
: "r" (rs), "r" (rt), "r" (lo), "r" (hi));
r = ((int64_t)hi << 32) | (uint32_t)lo;
assert(a + (int64_t)rs * rt == r);
assert(rd == lo);
return r;
}
int64_t madd1(int64_t a, int32_t rs, int32_t rt)
{
int32_t lo = a;
int32_t hi = a >> 32;
int32_t rd;
int64_t r;
__asm__ __volatile__ (
" mtlo1 %5\n"
" mthi1 %6\n"
" madd1 %0, %3, %4\n"
" mflo1 %1\n"
" mfhi1 %2\n"
: "=r" (rd), "=r" (lo), "=r" (hi)
: "r" (rs), "r" (rt), "r" (lo), "r" (hi));
r = ((int64_t)hi << 32) | (uint32_t)lo;
assert(a + (int64_t)rs * rt == r);
assert(rd == lo);
return r;
}
static int64_t madd_variants(int64_t a, int32_t rs, int32_t rt)
{
int64_t rd = madd(a, rs, rt);
int64_t rd1 = madd1(a, rs, rt);
assert(rd == rd1);
return rd;
}
static void verify_madd(int64_t a, int32_t rs, int32_t rt, int64_t expected)
{
assert(madd_variants(a, rs, rt) == expected);
assert(madd_variants(a, -rs, rt) == a + a - expected);
assert(madd_variants(a, rs, -rt) == a + a - expected);
assert(madd_variants(a, -rs, -rt) == expected);
}
int main()
{
verify_madd(13, 17, 19, 336);
return 0;
}
|
