void inv_q_mult(const long *q1, const long *q2, long *qProd) {
  typedef int invensensePutFunctionCallsHere;
  qProd[0] = (long)(((long long)q1[0] * q2[0] - (long long)q1[1] * q2[1] -
                     (long long)q1[2] * q2[2] - (long long)q1[3] * q2[3]) >>
                    30);
  qProd[1] = (int)(((long long)q1[0] * q2[1] + (long long)q1[1] * q2[0] +
                    (long long)q1[2] * q2[3] - (long long)q1[3] * q2[2]) >>
                   30);
  qProd[2] = (long)(((long long)q1[0] * q2[2] - (long long)q1[1] * q2[3] +
                     (long long)q1[2] * q2[0] + (long long)q1[3] * q2[1]) >>
                    30);
  qProd[3] = (long)(((long long)q1[0] * q2[3] + (long long)q1[1] * q2[2] -
                     (long long)q1[2] * q2[1] + (long long)q1[3] * q2[0]) >>
                    30);
}


x86 assembly (clang, -O0):
<inv_q_mult>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%rcx
mov    -0x10(%rbp),%rax
imul   (%rax),%rcx
mov    -0x8(%rbp),%rax
mov    0x8(%rax),%rax
mov    -0x10(%rbp),%rdx
imul   0x8(%rdx),%rax
sub    %rax,%rcx
mov    -0x8(%rbp),%rax
mov    0x10(%rax),%rax
mov    -0x10(%rbp),%rdx
imul   0x10(%rdx),%rax
sub    %rax,%rcx
mov    -0x8(%rbp),%rax
mov    0x18(%rax),%rax
mov    -0x10(%rbp),%rdx
imul   0x18(%rdx),%rax
sub    %rax,%rcx
sar    $0x1e,%rcx
mov    -0x18(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
mov    -0x10(%rbp),%rcx
imul   0x8(%rcx),%rax
mov    -0x8(%rbp),%rcx
mov    0x8(%rcx),%rcx
mov    -0x10(%rbp),%rdx
imul   (%rdx),%rcx
add    %rcx,%rax
mov    -0x8(%rbp),%rcx
mov    0x10(%rcx),%rcx
mov    -0x10(%rbp),%rdx
imul   0x18(%rdx),%rcx
add    %rcx,%rax
mov    -0x8(%rbp),%rcx
mov    0x18(%rcx),%rcx
mov    -0x10(%rbp),%rdx
imul   0x10(%rdx),%rcx
sub    %rcx,%rax
sar    $0x1e,%rax
movslq %eax,%rcx
mov    -0x18(%rbp),%rax
mov    %rcx,0x8(%rax)
mov    -0x8(%rbp),%rax
mov    (%rax),%rcx
mov    -0x10(%rbp),%rax
imul   0x10(%rax),%rcx
mov    -0x8(%rbp),%rax
mov    0x8(%rax),%rax
mov    -0x10(%rbp),%rdx
imul   0x18(%rdx),%rax
sub    %rax,%rcx
mov    -0x8(%rbp),%rax
mov    0x10(%rax),%rax
mov    -0x10(%rbp),%rdx
imul   (%rdx),%rax
add    %rax,%rcx
mov    -0x8(%rbp),%rax
mov    0x18(%rax),%rax
mov    -0x10(%rbp),%rdx
imul   0x8(%rdx),%rax
add    %rax,%rcx
sar    $0x1e,%rcx
mov    -0x18(%rbp),%rax
mov    %rcx,0x10(%rax)
mov    -0x8(%rbp),%rax
mov    (%rax),%rcx
mov    -0x10(%rbp),%rax
imul   0x18(%rax),%rcx
mov    -0x8(%rbp),%rax
mov    0x8(%rax),%rax
mov    -0x10(%rbp),%rdx
imul   0x10(%rdx),%rax
add    %rax,%rcx
mov    -0x8(%rbp),%rax
mov    0x10(%rax),%rax
mov    -0x10(%rbp),%rdx
imul   0x8(%rdx),%rax
sub    %rax,%rcx
mov    -0x8(%rbp),%rax
mov    0x18(%rax),%rax
mov    -0x10(%rbp),%rdx
imul   (%rdx),%rax
add    %rax,%rcx
sar    $0x1e,%rcx
mov    -0x18(%rbp),%rax
mov    %rcx,0x18(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<inv_q_mult>:
mov    (%rsi),%rax
imul   (%rdi),%rax
mov    0x8(%rsi),%r8
mov    0x10(%rsi),%r9
imul   0x8(%rdi),%r8
mov    0x18(%rsi),%rcx
imul   0x10(%rdi),%r9
add    %r8,%r9
imul   0x18(%rdi),%rcx
add    %r9,%rcx
sub    %rcx,%rax
sar    $0x1e,%rax
mov    %rax,(%rdx)
mov    (%rsi),%rax
mov    0x8(%rsi),%rcx
imul   (%rdi),%rcx
imul   0x8(%rdi),%rax
add    %rcx,%rax
mov    0x18(%rsi),%rcx
imul   0x10(%rdi),%rcx
add    %rax,%rcx
mov    0x10(%rsi),%rax
imul   0x18(%rdi),%rax
sub    %rax,%rcx
shr    $0x1e,%rcx
movslq %ecx,%rax
mov    %rax,0x8(%rdx)
mov    0x10(%rsi),%r8
imul   (%rdi),%r8
mov    0x18(%rsi),%rcx
imul   0x8(%rdi),%rcx
sub    %rcx,%r8
mov    (%rsi),%rcx
mov    0x8(%rsi),%rax
imul   0x10(%rdi),%rcx
add    %r8,%rcx
imul   0x18(%rdi),%rax
add    %rcx,%rax
sar    $0x1e,%rax
mov    %rax,0x10(%rdx)
mov    0x18(%rsi),%rax
imul   (%rdi),%rax
mov    0x10(%rsi),%rcx
imul   0x8(%rdi),%rcx
add    %rax,%rcx
mov    (%rsi),%rax
mov    0x8(%rsi),%rsi
imul   0x10(%rdi),%rsi
imul   0x18(%rdi),%rax
sub    %rsi,%rcx
add    %rcx,%rax
sar    $0x1e,%rax
mov    %rax,0x18(%rdx)
ret

LLVM IR (-O0):
define void @inv_q_mult(i64* %0, i64* %1, i64* %2) {
%4 = alloca i64*
%5 = alloca i64*
%6 = alloca i64*
store i64* %0, i64** %4
store i64* %1, i64** %5
store i64* %2, i64** %6
%7 = load i64*, i64** %4
%8 = getelementptr inbounds i64, i64* %7, i64 0
%9 = load i64, i64* %8
%10 = load i64*, i64** %5
%11 = getelementptr inbounds i64, i64* %10, i64 0
%12 = load i64, i64* %11
%13 = mul nsw i64 %9, %12
%14 = load i64*, i64** %4
%15 = getelementptr inbounds i64, i64* %14, i64 1
%16 = load i64, i64* %15
%17 = load i64*, i64** %5
%18 = getelementptr inbounds i64, i64* %17, i64 1
%19 = load i64, i64* %18
%20 = mul nsw i64 %16, %19
%21 = sub nsw i64 %13, %20
%22 = load i64*, i64** %4
%23 = getelementptr inbounds i64, i64* %22, i64 2
%24 = load i64, i64* %23
%25 = load i64*, i64** %5
%26 = getelementptr inbounds i64, i64* %25, i64 2
%27 = load i64, i64* %26
%28 = mul nsw i64 %24, %27
%29 = sub nsw i64 %21, %28
%30 = load i64*, i64** %4
%31 = getelementptr inbounds i64, i64* %30, i64 3
%32 = load i64, i64* %31
%33 = load i64*, i64** %5
%34 = getelementptr inbounds i64, i64* %33, i64 3
%35 = load i64, i64* %34
%36 = mul nsw i64 %32, %35
%37 = sub nsw i64 %29, %36
%38 = ashr i64 %37, 30
%39 = load i64*, i64** %6
%40 = getelementptr inbounds i64, i64* %39, i64 0
store i64 %38, i64* %40
%41 = load i64*, i64** %4
%42 = getelementptr inbounds i64, i64* %41, i64 0
%43 = load i64, i64* %42
%44 = load i64*, i64** %5
%45 = getelementptr inbounds i64, i64* %44, i64 1
%46 = load i64, i64* %45
%47 = mul nsw i64 %43, %46
%48 = load i64*, i64** %4
%49 = getelementptr inbounds i64, i64* %48, i64 1
%50 = load i64, i64* %49
%51 = load i64*, i64** %5
%52 = getelementptr inbounds i64, i64* %51, i64 0
%53 = load i64, i64* %52
%54 = mul nsw i64 %50, %53
%55 = add nsw i64 %47, %54
%56 = load i64*, i64** %4
%57 = getelementptr inbounds i64, i64* %56, i64 2
%58 = load i64, i64* %57
%59 = load i64*, i64** %5
%60 = getelementptr inbounds i64, i64* %59, i64 3
%61 = load i64, i64* %60
%62 = mul nsw i64 %58, %61
%63 = add nsw i64 %55, %62
%64 = load i64*, i64** %4
%65 = getelementptr inbounds i64, i64* %64, i64 3
%66 = load i64, i64* %65
%67 = load i64*, i64** %5
%68 = getelementptr inbounds i64, i64* %67, i64 2
%69 = load i64, i64* %68
%70 = mul nsw i64 %66, %69
%71 = sub nsw i64 %63, %70
%72 = ashr i64 %71, 30
%73 = trunc i64 %72 to i32
%74 = sext i32 %73 to i64
%75 = load i64*, i64** %6
%76 = getelementptr inbounds i64, i64* %75, i64 1
store i64 %74, i64* %76
%77 = load i64*, i64** %4
%78 = getelementptr inbounds i64, i64* %77, i64 0
%79 = load i64, i64* %78
%80 = load i64*, i64** %5
%81 = getelementptr inbounds i64, i64* %80, i64 2
%82 = load i64, i64* %81
%83 = mul nsw i64 %79, %82
%84 = load i64*, i64** %4
%85 = getelementptr inbounds i64, i64* %84, i64 1
%86 = load i64, i64* %85
%87 = load i64*, i64** %5
%88 = getelementptr inbounds i64, i64* %87, i64 3
%89 = load i64, i64* %88
%90 = mul nsw i64 %86, %89
%91 = sub nsw i64 %83, %90
%92 = load i64*, i64** %4
%93 = getelementptr inbounds i64, i64* %92, i64 2
%94 = load i64, i64* %93
%95 = load i64*, i64** %5
%96 = getelementptr inbounds i64, i64* %95, i64 0
%97 = load i64, i64* %96
%98 = mul nsw i64 %94, %97
%99 = add nsw i64 %91, %98
%100 = load i64*, i64** %4
%101 = getelementptr inbounds i64, i64* %100, i64 3
%102 = load i64, i64* %101
%103 = load i64*, i64** %5
%104 = getelementptr inbounds i64, i64* %103, i64 1
%105 = load i64, i64* %104
%106 = mul nsw i64 %102, %105
%107 = add nsw i64 %99, %106
%108 = ashr i64 %107, 30
%109 = load i64*, i64** %6
%110 = getelementptr inbounds i64, i64* %109, i64 2
store i64 %108, i64* %110
%111 = load i64*, i64** %4
%112 = getelementptr inbounds i64, i64* %111, i64 0
%113 = load i64, i64* %112
%114 = load i64*, i64** %5
%115 = getelementptr inbounds i64, i64* %114, i64 3
%116 = load i64, i64* %115
%117 = mul nsw i64 %113, %116
%118 = load i64*, i64** %4
%119 = getelementptr inbounds i64, i64* %118, i64 1
%120 = load i64, i64* %119
%121 = load i64*, i64** %5
%122 = getelementptr inbounds i64, i64* %121, i64 2
%123 = load i64, i64* %122
%124 = mul nsw i64 %120, %123
%125 = add nsw i64 %117, %124
%126 = load i64*, i64** %4
%127 = getelementptr inbounds i64, i64* %126, i64 2
%128 = load i64, i64* %127
%129 = load i64*, i64** %5
%130 = getelementptr inbounds i64, i64* %129, i64 1
%131 = load i64, i64* %130
%132 = mul nsw i64 %128, %131
%133 = sub nsw i64 %125, %132
%134 = load i64*, i64** %4
%135 = getelementptr inbounds i64, i64* %134, i64 3
%136 = load i64, i64* %135
%137 = load i64*, i64** %5
%138 = getelementptr inbounds i64, i64* %137, i64 0
%139 = load i64, i64* %138
%140 = mul nsw i64 %136, %139
%141 = add nsw i64 %133, %140
%142 = ashr i64 %141, 30
%143 = load i64*, i64** %6
%144 = getelementptr inbounds i64, i64* %143, i64 3
store i64 %142, i64* %144
ret void
}

LLVM IR (-O2):
define void @inv_q_mult(i64* nocapture %0, i64* nocapture %1, i64* nocapture writeonly %2) {
%4 = load i64, i64* %0
%5 = load i64, i64* %1
%6 = mul nsw i64 %5, %4
%7 = getelementptr inbounds i64, i64* %0, i64 1
%8 = load i64, i64* %7
%9 = getelementptr inbounds i64, i64* %1, i64 1
%10 = load i64, i64* %9
%11 = getelementptr inbounds i64, i64* %0, i64 2
%12 = load i64, i64* %11
%13 = getelementptr inbounds i64, i64* %1, i64 2
%14 = load i64, i64* %13
%15 = getelementptr inbounds i64, i64* %0, i64 3
%16 = load i64, i64* %15
%17 = getelementptr inbounds i64, i64* %1, i64 3
%18 = load i64, i64* %17
%19 = mul i64 %10, %8
%20 = mul i64 %14, %12
%21 = mul i64 %18, %16
%22 = add i64 %20, %19
%23 = add i64 %22, %21
%24 = sub i64 %6, %23
%25 = ashr i64 %24, 30
store i64 %25, i64* %2
%26 = load i64, i64* %0
%27 = load i64, i64* %9
%28 = mul nsw i64 %27, %26
%29 = load i64, i64* %7
%30 = load i64, i64* %1
%31 = mul nsw i64 %30, %29
%32 = add nsw i64 %31, %28
%33 = load i64, i64* %11
%34 = load i64, i64* %17
%35 = mul nsw i64 %34, %33
%36 = add nsw i64 %32, %35
%37 = load i64, i64* %15
%38 = load i64, i64* %13
%39 = mul nsw i64 %38, %37
%40 = sub i64 %36, %39
%41 = shl i64 %40, 2
%42 = ashr i64 %41, 32
%43 = getelementptr inbounds i64, i64* %2, i64 1
store i64 %42, i64* %43
%44 = load i64, i64* %0
%45 = load i64, i64* %13
%46 = mul nsw i64 %45, %44
%47 = load i64, i64* %7
%48 = load i64, i64* %17
%49 = mul nsw i64 %48, %47
%50 = sub i64 %46, %49
%51 = load i64, i64* %11
%52 = load i64, i64* %1
%53 = mul nsw i64 %52, %51
%54 = add nsw i64 %50, %53
%55 = load i64, i64* %15
%56 = load i64, i64* %9
%57 = mul nsw i64 %56, %55
%58 = add nsw i64 %54, %57
%59 = ashr i64 %58, 30
%60 = getelementptr inbounds i64, i64* %2, i64 2
store i64 %59, i64* %60
%61 = load i64, i64* %0
%62 = load i64, i64* %17
%63 = mul nsw i64 %62, %61
%64 = load i64, i64* %7
%65 = load i64, i64* %13
%66 = mul nsw i64 %65, %64
%67 = add nsw i64 %66, %63
%68 = load i64, i64* %11
%69 = load i64, i64* %9
%70 = mul nsw i64 %69, %68
%71 = sub i64 %67, %70
%72 = load i64, i64* %15
%73 = load i64, i64* %1
%74 = mul nsw i64 %73, %72
%75 = add nsw i64 %71, %74
%76 = ashr i64 %75, 30
%77 = getelementptr inbounds i64, i64* %2, i64 3
store i64 %76, i64* %77
ret void
}