#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int quatMultiply(double *, double *, double *, int);
int quatMultLoop(double *q1, double *q2, double *q1q2, unsigned long long nq1,
                 unsigned long long nq2, int p) {
  unsigned long long i, q1step, q2step, q1q2step, nq;
  int trash;
  q1q2step = 4;
  q1step = (nq1 == 1) ? 0 : 4;
  q2step = (nq2 == 1) ? 0 : 4;
  if ((nq1 > 1) && (nq2 > 1)) {
    nq = (nq1 <= nq2) ? nq1 : nq2;
  } else {
    nq = (nq1 >= nq2) ? nq1 : nq2;
  }
#pragma omp parallel shared(q1, q2, q1q2) private(i, trash)                    \
    firstprivate(q1step, q2step, q1q2step, nq, p)
  {
#pragma omp for schedule(static) nowait
    for (i = 0; i < nq; i++) {
      trash = quatMultiply(&(q1[i * q1step]), &(q2[i * q2step]),
                           &(q1q2[i * q1q2step]), p);
    }
  }
  return 1;
}


x86 assembly (clang, -O0):
<quatMultLoop>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
movq   $0x4,-0x50(%rbp)
mov    -0x20(%rbp),%rdx
mov    $0x4,%eax
xor    %ecx,%ecx
cmp    $0x1,%rdx
cmove  %ecx,%eax
cltq
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%rdx
mov    $0x4,%eax
xor    %ecx,%ecx
cmp    $0x1,%rdx
cmove  %ecx,%eax
cltq
mov    %rax,-0x48(%rbp)
cmpq   $0x1,-0x20(%rbp)
jbe    11ae <quatMultLoop+0x9e>
cmpq   $0x1,-0x28(%rbp)
jbe    11ae <quatMultLoop+0x9e>
mov    -0x20(%rbp),%rax
cmp    -0x28(%rbp),%rax
ja     1199 <quatMultLoop+0x89>
mov    -0x20(%rbp),%rax
mov    %rax,-0x68(%rbp)
jmp    11a1 <quatMultLoop+0x91>
mov    -0x28(%rbp),%rax
mov    %rax,-0x68(%rbp)
mov    -0x68(%rbp),%rax
mov    %rax,-0x58(%rbp)
jmp    11d9 <quatMultLoop+0xc9>
mov    -0x20(%rbp),%rax
cmp    -0x28(%rbp),%rax
jb     11c9 <quatMultLoop+0xb9>
mov    -0x20(%rbp),%rax
mov    %rax,-0x70(%rbp)
jmp    11d1 <quatMultLoop+0xc1>
mov    -0x28(%rbp),%rax
mov    %rax,-0x70(%rbp)
mov    -0x70(%rbp),%rax
mov    %rax,-0x58(%rbp)
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x58(%rbp),%rax
jae    1247 <quatMultLoop+0x137>
mov    -0x8(%rbp),%rdi
mov    -0x38(%rbp),%rax
imul   -0x40(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x10(%rbp),%rsi
mov    -0x38(%rbp),%rax
imul   -0x48(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
mov    -0x18(%rbp),%rdx
mov    -0x38(%rbp),%rax
imul   -0x50(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdx
mov    -0x2c(%rbp),%ecx
call   1030 <quatMultiply@plt>
mov    %eax,-0x5c(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11e1 <quatMultLoop+0xd1>
mov    $0x1,%eax
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<quatMultLoop>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
cmp    $0x1,%r8
setne  %r10b
cmp    $0x1,%rcx
setne  %r11b
cmp    %r8,%rcx
mov    %rcx,%r14
cmova  %r8,%r14
mov    %rcx,%rax
cmovb  %r8,%rax
cmp    $0x2,%r8
cmovb  %rax,%r14
cmp    $0x2,%rcx
cmovb  %rax,%r14
test   %r14,%r14
je     11a2 <quatMultLoop+0x92>
mov    %r9d,%r12d
mov    %rdx,%r15
mov    %rsi,%rbp
mov    %rdi,%rbx
xor    %r13d,%r13d
mov    %r11b,%r13b
shl    $0x2,%r13
xor    %eax,%eax
mov    %r10b,%al
shl    $0x2,%rax
shl    $0x3,%rax
mov    %rax,(%rsp)
shl    $0x3,%r13
nopl   0x0(%rax)
mov    %rbx,%rdi
mov    %rbp,%rsi
mov    %r15,%rdx
mov    %r12d,%ecx
call   1030 <quatMultiply@plt>
add    (%rsp),%rbp
add    %r13,%rbx
add    $0x20,%r15
add    $0xffffffffffffffff,%r14
jne    1180 <quatMultLoop+0x70>
mov    $0x1,%eax
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @quatMultLoop(double* %0, double* %1, double* %2, i64 %3, i64 %4, i32 %5) {
%7 = alloca double*
%8 = alloca double*
%9 = alloca double*
%10 = alloca i64
%11 = alloca i64
%12 = alloca i32
%13 = alloca i64
%14 = alloca i64
%15 = alloca i64
%16 = alloca i64
%17 = alloca i64
%18 = alloca i32
store double* %0, double** %7
store double* %1, double** %8
store double* %2, double** %9
store i64 %3, i64* %10
store i64 %4, i64* %11
store i32 %5, i32* %12
store i64 4, i64* %16
%19 = load i64, i64* %10
%20 = icmp eq i64 %19, 1
%21 = zext i1 %20 to i64
%22 = select i1 %20, i32 0, i32 4
%23 = sext i32 %22 to i64
store i64 %23, i64* %14
%24 = load i64, i64* %11
%25 = icmp eq i64 %24, 1
%26 = zext i1 %25 to i64
%27 = select i1 %25, i32 0, i32 4
%28 = sext i32 %27 to i64
store i64 %28, i64* %15
%29 = load i64, i64* %10
%30 = icmp ugt i64 %29, 1
br i1 %30, label %31, label %44
31: ; preds = %6
%32 = load i64, i64* %11
%33 = icmp ugt i64 %32, 1
br i1 %33, label %34, label %44
34: ; preds = %31
%35 = load i64, i64* %10
%36 = load i64, i64* %11
%37 = icmp ule i64 %35, %36
br i1 %37, label %38, label %40
38: ; preds = %34
%39 = load i64, i64* %10
br label %42
40: ; preds = %34
%41 = load i64, i64* %11
br label %42
42: ; preds = %40, %38
%43 = phi i64 [ %39, %38 ], [ %41, %40 ]
store i64 %43, i64* %17
br label %54
44: ; preds = %31, %6
%45 = load i64, i64* %10
%46 = load i64, i64* %11
%47 = icmp uge i64 %45, %46
br i1 %47, label %48, label %50
48: ; preds = %44
%49 = load i64, i64* %10
br label %52
50: ; preds = %44
%51 = load i64, i64* %11
br label %52
52: ; preds = %50, %48
%53 = phi i64 [ %49, %48 ], [ %51, %50 ]
store i64 %53, i64* %17
br label %54
54: ; preds = %52, %42
store i64 0, i64* %13
br label %55
55: ; preds = %77, %54
%56 = load i64, i64* %13
%57 = load i64, i64* %17
%58 = icmp ult i64 %56, %57
br i1 %58, label %59, label %80
59: ; preds = %55
%60 = load double*, double** %7
%61 = load i64, i64* %13
%62 = load i64, i64* %14
%63 = mul i64 %61, %62
%64 = getelementptr inbounds double, double* %60, i64 %63
%65 = load double*, double** %8
%66 = load i64, i64* %13
%67 = load i64, i64* %15
%68 = mul i64 %66, %67
%69 = getelementptr inbounds double, double* %65, i64 %68
%70 = load double*, double** %9
%71 = load i64, i64* %13
%72 = load i64, i64* %16
%73 = mul i64 %71, %72
%74 = getelementptr inbounds double, double* %70, i64 %73
%75 = load i32, i32* %12
%76 = call i32 @quatMultiply(double* %64, double* %69, double* %74, i32 %75)
store i32 %76, i32* %18
br label %77
77: ; preds = %59
%78 = load i64, i64* %13
%79 = add i64 %78, 1
store i64 %79, i64* %13
br label %55
80: ; preds = %55
ret i32 1
}

LLVM IR (-O2):
define i32 @quatMultLoop(double* %0, double* %1, double* %2, i64 %3, i64 %4, i32 %5) {
%7 = icmp eq i64 %3, 1
%8 = select i1 %7, i64 0, i64 4
%9 = icmp eq i64 %4, 1
%10 = select i1 %9, i64 0, i64 4
%11 = icmp ugt i64 %3, 1
%12 = icmp ugt i64 %4, 1
%13 = and i1 %11, %12
%14 = icmp ugt i64 %3, %4
%15 = select i1 %14, i64 %4, i64 %3
%16 = icmp ult i64 %3, %4
%17 = select i1 %16, i64 %4, i64 %3
%18 = select i1 %13, i64 %15, i64 %17
%19 = icmp eq i64 %18, 0
br i1 %19, label %31, label %20
20: ; preds = %6, %20
%21 = phi i64 [ %29, %20 ], [ 0, %6 ]
%22 = mul i64 %21, %8
%23 = getelementptr inbounds double, double* %0, i64 %22
%24 = mul i64 %21, %10
%25 = getelementptr inbounds double, double* %1, i64 %24
%26 = shl i64 %21, 2
%27 = getelementptr inbounds double, double* %2, i64 %26
%28 = tail call i32 @quatMultiply(double* %23, double* %25, double* %27, i32 %5)
%29 = add nuw i64 %21, 1
%30 = icmp eq i64 %29, %18
br i1 %30, label %31, label %20
31: ; preds = %20, %6
ret i32 1
}