void vflip1_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 - i][j];
      S[n - i][j] = S[i][j];
      S[i][j] = t;
    }
  }
}


x86 assembly (clang, -O0):
<vflip1_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 <vflip1_dmatrix+0xd6>
mov    -0x20(%rbp),%rax
mov    %eax,-0x30(%rbp)
movslq -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
jg     11c3 <vflip1_dmatrix+0xc3>
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%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    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%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 <vflip1_dmatrix+0x3f>
jmp    11c8 <vflip1_dmatrix+0xc8>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    112a <vflip1_dmatrix+0x2a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<vflip1_dmatrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
movslq %esi,%rax
cmp    %rdx,%rax
jle    111d <vflip1_dmatrix+0x1d>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
add    %edx,%esi
movslq %ecx,%rbx
movslq %esi,%r14
lea    0x1(%r8),%rcx
lea    0x1(%rbx),%r9
cmp    %r9,%rcx
cmovg  %rcx,%r9
mov    %r9,%r10
sub    %rbx,%r10
lea    -0x4(%r10),%rcx
mov    %rcx,-0x8(%rsp)
mov    %rcx,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %r10,%rcx
and    $0xfffffffffffffffc,%rcx
mov    %rcx,-0x18(%rsp)
add    %rbx,%rcx
mov    %rcx,-0x20(%rsp)
mov    %rsi,%rcx
mov    %rsi,-0x10(%rsp)
and    $0xfffffffffffffffe,%rsi
mov    %rsi,-0x28(%rsp)
lea    0x6(%rbx),%rcx
mov    %rcx,-0x30(%rsp)
jmp    1189 <vflip1_dmatrix+0x89>
xchg   %ax,%ax
add    $0x1,%rax
cmp    %rdx,%rax
jg     1112 <vflip1_dmatrix+0x12>
cmp    %r8,%rbx
jg     1180 <vflip1_dmatrix+0x80>
mov    %r14,%rcx
sub    %rax,%rcx
mov    (%rdi,%rcx,8),%rbp
mov    (%rdi,%rax,8),%r15
mov    %rbx,%rcx
cmp    $0x4,%r10
jb     12d0 <vflip1_dmatrix+0x1d0>
lea    0x0(,%rbx,8),%rcx
add    %rbp,%rcx
lea    (%r15,%r9,8),%rsi
cmp    %rsi,%rcx
jae    11d8 <vflip1_dmatrix+0xd8>
lea    0x0(,%r9,8),%r11
add    %rbp,%r11
lea    (%r15,%rbx,8),%rsi
mov    %rbx,%rcx
cmp    %r11,%rsi
jb     12d0 <vflip1_dmatrix+0x1d0>
mov    %r14,%rcx
cmpq   $0x4,-0x8(%rsp)
jae    11eb <vflip1_dmatrix+0xeb>
xor    %r12d,%r12d
jmp    1278 <vflip1_dmatrix+0x178>
mov    -0x30(%rsp),%rsi
lea    (%r15,%rsi,8),%r13
lea    0x0(,%rsi,8),%r11
add    %rbp,%r11
mov    -0x28(%rsp),%r14
xor    %r12d,%r12d
nopw   0x0(%rax,%rax,1)
movups -0x30(%r11,%r12,8),%xmm0
movups -0x20(%r11,%r12,8),%xmm1
movups -0x30(%r13,%r12,8),%xmm2
movups -0x20(%r13,%r12,8),%xmm3
movups %xmm2,-0x30(%r11,%r12,8)
movups %xmm3,-0x20(%r11,%r12,8)
movups %xmm0,-0x30(%r13,%r12,8)
movups %xmm1,-0x20(%r13,%r12,8)
movups -0x10(%r11,%r12,8),%xmm0
movups (%r11,%r12,8),%xmm1
movups -0x10(%r13,%r12,8),%xmm2
movups 0x0(%r13,%r12,8),%xmm3
movups %xmm2,-0x10(%r11,%r12,8)
movups %xmm3,(%r11,%r12,8)
movups %xmm0,-0x10(%r13,%r12,8)
movups %xmm1,0x0(%r13,%r12,8)
add    $0x8,%r12
add    $0xfffffffffffffffe,%r14
jne    1210 <vflip1_dmatrix+0x110>
testb  $0x1,-0x10(%rsp)
mov    %rcx,%r14
je     12b3 <vflip1_dmatrix+0x1b3>
add    %rbx,%r12
movups 0x0(%rbp,%r12,8),%xmm0
movups 0x10(%rbp,%r12,8),%xmm1
movups (%r15,%r12,8),%xmm2
movups 0x10(%r15,%r12,8),%xmm3
movups %xmm2,0x0(%rbp,%r12,8)
movups %xmm3,0x10(%rbp,%r12,8)
movups %xmm0,(%r15,%r12,8)
movups %xmm1,0x10(%r15,%r12,8)
mov    -0x20(%rsp),%rcx
cmp    -0x18(%rsp),%r10
je     1180 <vflip1_dmatrix+0x80>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  0x0(%rbp,%rcx,8),%xmm0
movsd  (%r15,%rcx,8),%xmm1
movsd  %xmm1,0x0(%rbp,%rcx,8)
movsd  %xmm0,(%r15,%rcx,8)
add    $0x1,%rcx
cmp    %r8,%rcx
jle    12d0 <vflip1_dmatrix+0x1d0>
jmp    1180 <vflip1_dmatrix+0x80>

LLVM IR (-O0):
define void @vflip1_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* %11
%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* %12
%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* %11
%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* %12
%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 @vflip1_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 %126, 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
%16 = add nsw i64 %4, 1
%17 = add nsw i64 %12, 1
%18 = call i64 @llvm.smax.i64(i64 %16, i64 %17)
%19 = add nsw i64 %4, 1
%20 = add nsw i64 %12, 1
%21 = call i64 @llvm.smax.i64(i64 %19, i64 %20)
%22 = sub i64 %21, %12
%23 = add i64 %22, -4
%24 = lshr i64 %23, 2
%25 = add nuw nsw i64 %24, 1
%26 = icmp ult i64 %22, 4
%27 = and i64 %22, -4
%28 = add i64 %12, %27
%29 = and i64 %25, 1
%30 = icmp ult i64 %23, 4
%31 = and i64 %25, 9223372036854775806
%32 = icmp eq i64 %29, 0
%33 = icmp eq i64 %22, %27
br label %34
34: ; preds = %9, %123
%35 = phi i64 [ %7, %9 ], [ %124, %123 ]
br i1 %13, label %123, label %36
36: ; preds = %34
%37 = sub nsw i64 %15, %35
%38 = getelementptr inbounds double*, double** %0, i64 %37
%39 = load double*, double** %38
%40 = getelementptr inbounds double*, double** %0, i64 %35
%41 = load double*, double** %40
br i1 %26, label %113, label %42
42: ; preds = %36
%43 = getelementptr double, double* %39, i64 %12
%44 = getelementptr double, double* %39, i64 %18
%45 = getelementptr double, double* %41, i64 %12
%46 = getelementptr double, double* %41, i64 %18
%47 = icmp ult double* %43, %46
%48 = icmp ult double* %45, %44
%49 = and i1 %47, %48
br i1 %49, label %113, label %50
50: ; preds = %42
br i1 %30, label %92, label %51
51: ; preds = %50, %51
%52 = phi i64 [ %89, %51 ], [ 0, %50 ]
%53 = phi i64 [ %90, %51 ], [ 0, %50 ]
%54 = add i64 %12, %52
%55 = getelementptr inbounds double, double* %39, i64 %54
%56 = bitcast double* %55 to <2 x double>*
%57 = load <2 x double>, <2 x double>* %56
%58 = getelementptr inbounds double, double* %55, i64 2
%59 = bitcast double* %58 to <2 x double>*
%60 = load <2 x double>, <2 x double>* %59
%61 = getelementptr inbounds double, double* %41, i64 %54
%62 = bitcast double* %61 to <2 x double>*
%63 = load <2 x double>, <2 x double>* %62
%64 = getelementptr inbounds double, double* %61, i64 2
%65 = bitcast double* %64 to <2 x double>*
%66 = load <2 x double>, <2 x double>* %65
%67 = bitcast double* %55 to <2 x double>*
store <2 x double> %63, <2 x double>* %67
%68 = bitcast double* %58 to <2 x double>*
store <2 x double> %66, <2 x double>* %68
%69 = bitcast double* %61 to <2 x double>*
store <2 x double> %57, <2 x double>* %69
%70 = bitcast double* %64 to <2 x double>*
store <2 x double> %60, <2 x double>* %70
%71 = or i64 %52, 4
%72 = add i64 %12, %71
%73 = getelementptr inbounds double, double* %39, i64 %72
%74 = bitcast double* %73 to <2 x double>*
%75 = load <2 x double>, <2 x double>* %74
%76 = getelementptr inbounds double, double* %73, i64 2
%77 = bitcast double* %76 to <2 x double>*
%78 = load <2 x double>, <2 x double>* %77
%79 = getelementptr inbounds double, double* %41, i64 %72
%80 = bitcast double* %79 to <2 x double>*
%81 = load <2 x double>, <2 x double>* %80
%82 = getelementptr inbounds double, double* %79, i64 2
%83 = bitcast double* %82 to <2 x double>*
%84 = load <2 x double>, <2 x double>* %83
%85 = bitcast double* %73 to <2 x double>*
store <2 x double> %81, <2 x double>* %85
%86 = bitcast double* %76 to <2 x double>*
store <2 x double> %84, <2 x double>* %86
%87 = bitcast double* %79 to <2 x double>*
store <2 x double> %75, <2 x double>* %87
%88 = bitcast double* %82 to <2 x double>*
store <2 x double> %78, <2 x double>* %88
%89 = add nuw i64 %52, 8
%90 = add i64 %53, 2
%91 = icmp eq i64 %90, %31
br i1 %91, label %92, label %51
92: ; preds = %51, %50
%93 = phi i64 [ 0, %50 ], [ %89, %51 ]
br i1 %32, label %112, label %94
94: ; preds = %92
%95 = add i64 %12, %93
%96 = getelementptr inbounds double, double* %39, i64 %95
%97 = bitcast double* %96 to <2 x double>*
%98 = load <2 x double>, <2 x double>* %97
%99 = getelementptr inbounds double, double* %96, i64 2
%100 = bitcast double* %99 to <2 x double>*
%101 = load <2 x double>, <2 x double>* %100
%102 = getelementptr inbounds double, double* %41, i64 %95
%103 = bitcast double* %102 to <2 x double>*
%104 = load <2 x double>, <2 x double>* %103
%105 = getelementptr inbounds double, double* %102, i64 2
%106 = bitcast double* %105 to <2 x double>*
%107 = load <2 x double>, <2 x double>* %106
%108 = bitcast double* %96 to <2 x double>*
store <2 x double> %104, <2 x double>* %108
%109 = bitcast double* %99 to <2 x double>*
store <2 x double> %107, <2 x double>* %109
%110 = bitcast double* %102 to <2 x double>*
store <2 x double> %98, <2 x double>* %110
%111 = bitcast double* %105 to <2 x double>*
store <2 x double> %101, <2 x double>* %111
br label %112
112: ; preds = %92, %94
br i1 %33, label %123, label %113
113: ; preds = %42, %36, %112
%114 = phi i64 [ %12, %42 ], [ %12, %36 ], [ %28, %112 ]
br label %115
115: ; preds = %113, %115
%116 = phi i64 [ %121, %115 ], [ %114, %113 ]
%117 = getelementptr inbounds double, double* %39, i64 %116
%118 = load double, double* %117
%119 = getelementptr inbounds double, double* %41, i64 %116
%120 = load double, double* %119
store double %120, double* %117
store double %118, double* %119
%121 = add i64 %116, 1
%122 = icmp sgt i64 %121, %4
br i1 %122, label %123, label %115
123: ; preds = %115, %112, %34
%124 = add i64 %35, 1
%125 = icmp sgt i64 %124, %2
br i1 %125, label %126, label %34
126: ; preds = %123, %5
ret void
}