rot2(a, cose, sine, c1, c2) double a[4][4], cose, sine;
int c1, c2;
{
  register int i;
  double t1, t2;
  for (i = 0; i < 4; i++) {
    t1 = a[i][c1] * cose;
    t2 = a[i][c1] * sine;
    a[i][c1] = t1 - a[i][c2] * sine;
    a[i][c2] = t2 + a[i][c2] * cose;
  }
}


x86 assembly (clang, -O0):
<rot2>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    %esi,-0x24(%rbp)
mov    %edx,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
cmpl   $0x4,-0x2c(%rbp)
jge    1216 <rot2+0x116>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x20(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x18(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    111f <rot2+0x1f>
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<rot2>:
movslq %esi,%rcx
movslq %edx,%rax
movsd  (%rdi,%rcx,8),%xmm2
movapd %xmm2,%xmm3
mulsd  %xmm0,%xmm3
mulsd  %xmm1,%xmm2
movsd  (%rdi,%rax,8),%xmm4
mulsd  %xmm1,%xmm4
subsd  %xmm4,%xmm3
movsd  %xmm3,(%rdi,%rcx,8)
movsd  (%rdi,%rax,8),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%rdi,%rax,8)
movsd  0x20(%rdi,%rcx,8),%xmm2
movapd %xmm2,%xmm3
mulsd  %xmm0,%xmm3
mulsd  %xmm1,%xmm2
movsd  0x20(%rdi,%rax,8),%xmm4
mulsd  %xmm1,%xmm4
subsd  %xmm4,%xmm3
movsd  %xmm3,0x20(%rdi,%rcx,8)
movsd  0x20(%rdi,%rax,8),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,0x20(%rdi,%rax,8)
movsd  0x40(%rdi,%rcx,8),%xmm2
movapd %xmm2,%xmm3
mulsd  %xmm0,%xmm3
mulsd  %xmm1,%xmm2
movsd  0x40(%rdi,%rax,8),%xmm4
mulsd  %xmm1,%xmm4
subsd  %xmm4,%xmm3
movsd  %xmm3,0x40(%rdi,%rcx,8)
movsd  0x40(%rdi,%rax,8),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,0x40(%rdi,%rax,8)
movsd  0x60(%rdi,%rcx,8),%xmm2
movapd %xmm2,%xmm3
mulsd  %xmm0,%xmm3
mulsd  %xmm1,%xmm2
mulsd  0x60(%rdi,%rax,8),%xmm1
subsd  %xmm1,%xmm3
movsd  %xmm3,0x60(%rdi,%rcx,8)
mulsd  0x60(%rdi,%rax,8),%xmm0
addsd  %xmm2,%xmm0
movsd  %xmm0,0x60(%rdi,%rax,8)
ret

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

LLVM IR (-O2):
define i32 @rot2([4 x double]* nocapture %0, double %1, double %2, i32 %3, i32 %4) {
%6 = sext i32 %3 to i64
%7 = sext i32 %4 to i64
%8 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 0, i64 %6
%9 = load double, double* %8
%10 = fmul double %9, %1
%11 = fmul double %9, %2
%12 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 0, i64 %7
%13 = load double, double* %12
%14 = fneg double %13
%15 = tail call double @llvm.fmuladd.f64(double %14, double %2, double %10)
store double %15, double* %8
%16 = load double, double* %12
%17 = tail call double @llvm.fmuladd.f64(double %16, double %1, double %11)
store double %17, double* %12
%18 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 1, i64 %6
%19 = load double, double* %18
%20 = fmul double %19, %1
%21 = fmul double %19, %2
%22 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 1, i64 %7
%23 = load double, double* %22
%24 = fneg double %23
%25 = tail call double @llvm.fmuladd.f64(double %24, double %2, double %20)
store double %25, double* %18
%26 = load double, double* %22
%27 = tail call double @llvm.fmuladd.f64(double %26, double %1, double %21)
store double %27, double* %22
%28 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 2, i64 %6
%29 = load double, double* %28
%30 = fmul double %29, %1
%31 = fmul double %29, %2
%32 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 2, i64 %7
%33 = load double, double* %32
%34 = fneg double %33
%35 = tail call double @llvm.fmuladd.f64(double %34, double %2, double %30)
store double %35, double* %28
%36 = load double, double* %32
%37 = tail call double @llvm.fmuladd.f64(double %36, double %1, double %31)
store double %37, double* %32
%38 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 3, i64 %6
%39 = load double, double* %38
%40 = fmul double %39, %1
%41 = fmul double %39, %2
%42 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 3, i64 %7
%43 = load double, double* %42
%44 = fneg double %43
%45 = tail call double @llvm.fmuladd.f64(double %44, double %2, double %40)
store double %45, double* %38
%46 = load double, double* %42
%47 = tail call double @llvm.fmuladd.f64(double %46, double %1, double %41)
store double %47, double* %42
ret i32 undef
}