void lrotate1_dmatrix(double **S, long nrl, long nrh, long ncl, long nch) {
  int i, j, n = nrl + nrh;
  double t;
  for (i = nrl; i <= nrh; i++) {
    for (j = ncl; j <= nch; j++) {
      t = S[n - j][i];
      S[n - j][i] = S[i][j];
      S[i][j] = t;
    }
  }
}


x86 assembly (clang, -O0):
<lrotate1_dmatrix>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x10(%rbp),%rax
add    -0x18(%rbp),%rax
mov    %eax,-0x34(%rbp)
mov    -0x10(%rbp),%rax
mov    %eax,-0x2c(%rbp)
movslq -0x2c(%rbp),%rax
cmp    -0x18(%rbp),%rax
jg     11d6 <lrotate1_dmatrix+0xd6>
mov    -0x20(%rbp),%rax
mov    %eax,-0x30(%rbp)
movslq -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
jg     11c3 <lrotate1_dmatrix+0xc3>
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    113f <lrotate1_dmatrix+0x3f>
jmp    11c8 <lrotate1_dmatrix+0xc8>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    112a <lrotate1_dmatrix+0x2a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<lrotate1_dmatrix>:
push   %rbx
movslq %esi,%r11
cmp    %rdx,%r11
jle    110b <lrotate1_dmatrix+0xb>
pop    %rbx
ret
add    %edx,%esi
movslq %ecx,%r10
movslq %esi,%rax
sub    %r10,%rax
lea    (%rdi,%rax,8),%r9
jmp    1129 <lrotate1_dmatrix+0x29>
nopl   0x0(%rax)
add    $0x1,%r11
cmp    %rdx,%r11
jg     1109 <lrotate1_dmatrix+0x9>
cmp    %r8,%r10
jg     1120 <lrotate1_dmatrix+0x20>
mov    (%rdi,%r11,8),%rsi
mov    %r9,%rcx
mov    %r10,%rax
nopl   0x0(%rax,%rax,1)
mov    (%rcx),%rbx
movsd  (%rbx,%r11,8),%xmm0
movsd  (%rsi,%rax,8),%xmm1
movsd  %xmm1,(%rbx,%r11,8)
movsd  %xmm0,(%rsi,%rax,8)
add    $0x1,%rax
add    $0xfffffffffffffff8,%rcx
cmp    %r8,%rax
jle    1140 <lrotate1_dmatrix+0x40>
jmp    1120 <lrotate1_dmatrix+0x20>

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

LLVM IR (-O2):
define void @lrotate1_dmatrix(double** nocapture %0, i64 %1, i64 %2, i64 %3, i64 %4) {
%6 = shl i64 %1, 32
%7 = ashr exact i64 %6, 32
%8 = icmp sgt i64 %7, %2
br i1 %8, label %35, label %9
9: ; preds = %5
%10 = add nsw i64 %2, %1
%11 = shl i64 %3, 32
%12 = ashr exact i64 %11, 32
%13 = icmp sgt i64 %12, %4
%14 = shl i64 %10, 32
%15 = ashr exact i64 %14, 32
br label %16
16: ; preds = %9, %32
%17 = phi i64 [ %7, %9 ], [ %33, %32 ]
br i1 %13, label %32, label %18
18: ; preds = %16
%19 = getelementptr inbounds double*, double** %0, i64 %17
%20 = load double*, double** %19
br label %21
21: ; preds = %18, %21
%22 = phi i64 [ %12, %18 ], [ %30, %21 ]
%23 = sub nsw i64 %15, %22
%24 = getelementptr inbounds double*, double** %0, i64 %23
%25 = load double*, double** %24
%26 = getelementptr inbounds double, double* %25, i64 %17
%27 = load double, double* %26
%28 = getelementptr inbounds double, double* %20, i64 %22
%29 = load double, double* %28
store double %29, double* %26
store double %27, double* %28
%30 = add i64 %22, 1
%31 = icmp sgt i64 %30, %4
br i1 %31, label %32, label %21
32: ; preds = %21, %16
%33 = add i64 %17, 1
%34 = icmp sgt i64 %33, %2
br i1 %34, label %35, label %16
35: ; preds = %32, %5
ret void
}