#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 */
double signx(double);
double sqrt(double);
void jrotate(double *A, int N, double a, double b, int i) {
  int j, r, r1;
  double c, s, den, y, w;
  r = i * N;
  r1 = r + N;
  if (a == 0.0) {
    c = 0.0;
    s = signx(b);
  } else {
    den = sqrt(a * a + b * b);
    c = a / den;
    s = b / den;
  }
  for (j = i; j < N; ++j) {
    y = A[r + j];
    w = A[r1 + j];
    A[r + j] = c * y - s * w;
    A[r1 + j] = s * y + c * w;
  }
}


x86 assembly (clang, -O0):
<jrotate>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    %edx,-0x24(%rbp)
mov    -0x24(%rbp),%eax
imul   -0xc(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
add    -0xc(%rbp),%eax
mov    %eax,-0x30(%rbp)
movsd  -0x18(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1183 <jrotate+0x63>
jp     1183 <jrotate+0x63>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x20(%rbp),%xmm0
call   1030 <signx@plt>
movsd  %xmm0,-0x40(%rbp)
jmp    11c7 <jrotate+0xa7>
movsd  -0x18(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
movsd  -0x20(%rbp),%xmm1
mulsd  -0x20(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movsd  %xmm0,-0x48(%rbp)
movsd  -0x18(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x20(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1288 <jrotate+0x168>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  -0x58(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11cd <jrotate+0xad>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<jrotate>:
push   %r15
push   %r14
push   %rbx
sub    $0x20,%rsp
mov    %edx,%r15d
mov    %esi,%r14d
mov    %rdi,%rbx
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm0
jne    115b <jrotate+0x3b>
jp     115b <jrotate+0x3b>
movapd %xmm1,%xmm0
call   1030 <signx@plt>
xorpd  %xmm9,%xmm9
unpcklpd %xmm0,%xmm9
cmp    %r14d,%r15d
jl     11b1 <jrotate+0x91>
jmp    1317 <jrotate+0x1f7>
movapd %xmm0,%xmm9
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
mulsd  %xmm9,%xmm0
addsd  %xmm2,%xmm0
ucomisd %xmm3,%xmm0
jb     117d <jrotate+0x5d>
sqrtsd %xmm0,%xmm0
jmp    119a <jrotate+0x7a>
movapd %xmm9,0x10(%rsp)
movapd %xmm1,(%rsp)
call   1040 <sqrt@plt>
movapd (%rsp),%xmm1
movapd 0x10(%rsp),%xmm9
unpcklpd %xmm1,%xmm9
unpcklpd %xmm0,%xmm0
divpd  %xmm0,%xmm9
cmp    %r14d,%r15d
jge    1317 <jrotate+0x1f7>
mov    %r15d,%ecx
imul   %r14d,%ecx
movapd %xmm9,%xmm0
unpckhpd %xmm9,%xmm0
movapd 0xe35(%rip),%xmm8
xorpd  %xmm0,%xmm8
movslq %r15d,%rax
movslq %ecx,%r8
movslq %r14d,%rcx
mov    %rcx,%r9
sub    %rax,%r9
cmp    $0x2,%r9
jb     12b9 <jrotate+0x199>
mov    %rax,%rdx
not    %rdx
add    %rcx,%rdx
lea    0x1(%r15),%esi
imul   %r14d,%esi
add    %r15d,%esi
lea    (%rdx,%rsi,1),%edi
cmp    %esi,%edi
jl     12b9 <jrotate+0x199>
shr    $0x20,%rdx
jne    12b9 <jrotate+0x199>
lea    (%rax,%r8,1),%rdx
lea    (%rbx,%rdx,8),%rdi
movslq %esi,%rdx
lea    (%rcx,%rdx,1),%rsi
sub    %rax,%rsi
lea    (%rbx,%rsi,8),%rsi
cmp    %rsi,%rdi
jae    123e <jrotate+0x11e>
lea    (%rcx,%r8,1),%rsi
lea    (%rbx,%rsi,8),%rsi
lea    (%rbx,%rdx,8),%rdx
cmp    %rsi,%rdx
jb     12b9 <jrotate+0x199>
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
add    %r10,%rax
movapd %xmm8,%xmm2
unpcklpd %xmm8,%xmm2
movapd %xmm9,%xmm3
unpcklpd %xmm9,%xmm3
lea    0x1(%r14),%r11d
imul   %r15d,%r11d
add    %r14d,%r11d
xor    %esi,%esi
nopl   0x0(%rax)
movupd (%rdi,%rsi,8),%xmm4
lea    (%r11,%rsi,1),%edx
movslq %edx,%rdx
movupd (%rbx,%rdx,8),%xmm5
movapd %xmm5,%xmm6
mulpd  %xmm2,%xmm6
movapd %xmm3,%xmm7
mulpd  %xmm4,%xmm7
addpd  %xmm6,%xmm7
movupd %xmm7,(%rdi,%rsi,8)
mulpd  %xmm3,%xmm5
mulpd  %xmm0,%xmm4
addpd  %xmm5,%xmm4
movupd %xmm4,(%rbx,%rdx,8)
add    $0x2,%rsi
cmp    %rsi,%r10
jne    1270 <jrotate+0x150>
cmp    %r10,%r9
je     1317 <jrotate+0x1f7>
add    $0x1,%r15d
imul   %r15d,%r14d
lea    (%rbx,%r8,8),%rdx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rdx,%rax,8),%xmm2
lea    (%r14,%rax,1),%esi
movslq %esi,%rsi
movsd  (%rbx,%rsi,8),%xmm3
movapd %xmm3,%xmm4
mulsd  %xmm8,%xmm4
movapd %xmm9,%xmm5
mulsd  %xmm2,%xmm5
addsd  %xmm4,%xmm5
movsd  %xmm5,(%rdx,%rax,8)
mulsd  %xmm9,%xmm3
mulsd  %xmm0,%xmm2
addsd  %xmm3,%xmm2
movsd  %xmm2,(%rbx,%rsi,8)
add    $0x1,%rax
cmp    %rax,%rcx
jne    12d0 <jrotate+0x1b0>
add    $0x20,%rsp
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @jrotate(double* %0, i32 %1, double %2, double %3, i32 %4) {
%6 = alloca double*
%7 = alloca i32
%8 = alloca double
%9 = alloca double
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
store double* %0, double** %6
store i32 %1, i32* %7
store double %2, double* %8
store double %3, double* %9
store i32 %4, i32* %10
%19 = load i32, i32* %10
%20 = load i32, i32* %7
%21 = mul nsw i32 %19, %20
store i32 %21, i32* %12
%22 = load i32, i32* %12
%23 = load i32, i32* %7
%24 = add nsw i32 %22, %23
store i32 %24, i32* %13
%25 = load double, double* %8
%26 = fcmp oeq double %25, 0.000000e+00
br i1 %26, label %27, label %30
27: ; preds = %5
store double 0.000000e+00, double* %14
%28 = load double, double* %9
%29 = call double @signx(double %28)
store double %29, double* %15
br label %44
30: ; preds = %5
%31 = load double, double* %8
%32 = load double, double* %8
%33 = load double, double* %9
%34 = load double, double* %9
%35 = fmul double %33, %34
%36 = call double @llvm.fmuladd.f64(double %31, double %32, double %35)
%37 = call double @sqrt(double %36)
store double %37, double* %16
%38 = load double, double* %8
%39 = load double, double* %16
%40 = fdiv double %38, %39
store double %40, double* %14
%41 = load double, double* %9
%42 = load double, double* %16
%43 = fdiv double %41, %42
store double %43, double* %15
br label %44
44: ; preds = %30, %27
%45 = load i32, i32* %10
store i32 %45, i32* %11
br label %46
46: ; preds = %90, %44
%47 = load i32, i32* %11
%48 = load i32, i32* %7
%49 = icmp slt i32 %47, %48
br i1 %49, label %50, label %93
50: ; preds = %46
%51 = load double*, double** %6
%52 = load i32, i32* %12
%53 = load i32, i32* %11
%54 = add nsw i32 %52, %53
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %51, i64 %55
%57 = load double, double* %56
store double %57, double* %17
%58 = load double*, double** %6
%59 = load i32, i32* %13
%60 = load i32, i32* %11
%61 = add nsw i32 %59, %60
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %58, i64 %62
%64 = load double, double* %63
store double %64, double* %18
%65 = load double, double* %14
%66 = load double, double* %17
%67 = load double, double* %15
%68 = load double, double* %18
%69 = fmul double %67, %68
%70 = fneg double %69
%71 = call double @llvm.fmuladd.f64(double %65, double %66, double %70)
%72 = load double*, double** %6
%73 = load i32, i32* %12
%74 = load i32, i32* %11
%75 = add nsw i32 %73, %74
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %72, i64 %76
store double %71, double* %77
%78 = load double, double* %15
%79 = load double, double* %17
%80 = load double, double* %14
%81 = load double, double* %18
%82 = fmul double %80, %81
%83 = call double @llvm.fmuladd.f64(double %78, double %79, double %82)
%84 = load double*, double** %6
%85 = load i32, i32* %13
%86 = load i32, i32* %11
%87 = add nsw i32 %85, %86
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds double, double* %84, i64 %88
store double %83, double* %89
br label %90
90: ; preds = %50
%91 = load i32, i32* %11
%92 = add nsw i32 %91, 1
store i32 %92, i32* %11
br label %46
93: ; preds = %46
ret void
}

LLVM IR (-O2):
define void @jrotate(double* nocapture %0, i32 %1, double %2, double %3, i32 %4) {
%6 = mul nsw i32 %4, %1
%7 = fcmp oeq double %2, 0.000000e+00
br i1 %7, label %8, label %11
8: ; preds = %5
%9 = tail call double @signx(double %3)
%10 = insertelement <2 x double> <double 0.000000e+00, double poison>, double %9, i64 1
br label %20
11: ; preds = %5
%12 = fmul double %3, %3
%13 = tail call double @llvm.fmuladd.f64(double %2, double %2, double %12)
%14 = tail call double @sqrt(double %13)
%15 = insertelement <2 x double> poison, double %2, i64 0
%16 = insertelement <2 x double> %15, double %3, i64 1
%17 = insertelement <2 x double> poison, double %14, i64 0
%18 = shufflevector <2 x double> %17, <2 x double> poison, <2 x i32> zeroinitializer
%19 = fdiv <2 x double> %16, %18
br label %20
20: ; preds = %11, %8
%21 = phi <2 x double> [ %10, %8 ], [ %19, %11 ]
%22 = icmp slt i32 %4, %1
br i1 %22, label %23, label %110
23: ; preds = %20
%24 = extractelement <2 x double> %21, i64 1
%25 = fneg double %24
%26 = sext i32 %4 to i64
%27 = sext i32 %6 to i64
%28 = sext i32 %1 to i64
%29 = sub nsw i64 %28, %26
%30 = icmp ult i64 %29, 2
br i1 %30, label %90, label %31
31: ; preds = %23
%32 = xor i64 %26, -1
%33 = add nsw i64 %32, %28
%34 = add i32 %4, 1
%35 = mul i32 %34, %1
%36 = add i32 %35, %4
%37 = trunc i64 %33 to i32
%38 = add i32 %36, %37
%39 = icmp slt i32 %38, %36
%40 = icmp ugt i64 %33, 4294967295
%41 = or i1 %39, %40
br i1 %41, label %90, label %42
42: ; preds = %31
%43 = add nsw i64 %26, %27
%44 = getelementptr double, double* %0, i64 %43
%45 = add nsw i64 %28, %27
%46 = getelementptr double, double* %0, i64 %45
%47 = add i32 %4, 1
%48 = mul i32 %47, %1
%49 = add i32 %48, %4
%50 = sext i32 %49 to i64
%51 = getelementptr double, double* %0, i64 %50
%52 = add nsw i64 %28, %50
%53 = sub nsw i64 %52, %26
%54 = getelementptr double, double* %0, i64 %53
%55 = icmp ult double* %44, %54
%56 = icmp ult double* %51, %46
%57 = and i1 %55, %56
br i1 %57, label %90, label %58
58: ; preds = %42
%59 = and i64 %29, -2
%60 = add nsw i64 %59, %26
%61 = insertelement <2 x double> poison, double %25, i64 0
%62 = shufflevector <2 x double> %61, <2 x double> poison, <2 x i32> zeroinitializer
%63 = shufflevector <2 x double> %21, <2 x double> poison, <2 x i32> zeroinitializer
%64 = shufflevector <2 x double> %21, <2 x double> undef, <2 x i32> <i32 1, i32 1>
br label %65
65: ; preds = %65, %58
%66 = phi i64 [ 0, %58 ], [ %86, %65 ]
%67 = add i64 %66, %26
%68 = trunc i64 %66 to i32
%69 = add i32 %68, %4
%70 = add nsw i64 %67, %27
%71 = getelementptr inbounds double, double* %0, i64 %70
%72 = bitcast double* %71 to <2 x double>*
%73 = load <2 x double>, <2 x double>* %72
%74 = add i32 %6, %69
%75 = add i32 %74, %1
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %0, i64 %76
%78 = bitcast double* %77 to <2 x double>*
%79 = load <2 x double>, <2 x double>* %78
%80 = fmul <2 x double> %79, %62
%81 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %63, <2 x double> %73, <2 x double> %80)
%82 = bitcast double* %71 to <2 x double>*
store <2 x double> %81, <2 x double>* %82
%83 = fmul <2 x double> %63, %79
%84 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %64, <2 x double> %73, <2 x double> %83)
%85 = bitcast double* %77 to <2 x double>*
store <2 x double> %84, <2 x double>* %85
%86 = add nuw i64 %66, 2
%87 = icmp eq i64 %86, %59
br i1 %87, label %88, label %65
88: ; preds = %65
%89 = icmp eq i64 %29, %59
br i1 %89, label %110, label %90
90: ; preds = %42, %31, %23, %88
%91 = phi i64 [ %26, %42 ], [ %26, %31 ], [ %26, %23 ], [ %60, %88 ]
%92 = extractelement <2 x double> %21, i64 0
br label %93
93: ; preds = %90, %93
%94 = phi i64 [ %108, %93 ], [ %91, %90 ]
%95 = trunc i64 %94 to i32
%96 = add nsw i64 %94, %27
%97 = getelementptr inbounds double, double* %0, i64 %96
%98 = load double, double* %97
%99 = add i32 %6, %95
%100 = add i32 %99, %1
%101 = sext i32 %100 to i64
%102 = getelementptr inbounds double, double* %0, i64 %101
%103 = load double, double* %102
%104 = fmul double %103, %25
%105 = tail call double @llvm.fmuladd.f64(double %92, double %98, double %104)
store double %105, double* %97
%106 = fmul double %92, %103
%107 = tail call double @llvm.fmuladd.f64(double %24, double %98, double %106)
store double %107, double* %102
%108 = add nsw i64 %94, 1
%109 = icmp eq i64 %108, %28
br i1 %109, label %110, label %93
110: ; preds = %93, %88, %20
ret void
}