double ***sub_const_r_3(double ***res, double cn, double ***vec, int lo1,
                        int hi1, int lo2, int hi2, int lo3, int hi3) {
  int i, j, k;
  for (i = 0; i <= hi1 - lo1; i++)
    for (j = 0; j <= hi2 - lo2; j++) {
      for (k = 0; k <= hi3 - lo3; k++)
        res[i][j][k] = cn - vec[i][j][k];
      res[i][j][hi3 - lo3 + 1] -= cn * (hi3 - lo3 + 1);
    }
  return res;
}


x86 assembly (clang, -O0):
<sub_const_r_3>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %ecx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9d,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x20(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jg     1236 <sub_const_r_3+0x136>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x28(%rbp),%ecx
sub    -0x24(%rbp),%ecx
cmp    %ecx,%eax
jg     1223 <sub_const_r_3+0x123>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    0x18(%rbp),%ecx
sub    0x10(%rbp),%ecx
cmp    %ecx,%eax
jg     11ba <sub_const_r_3+0xba>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    115c <sub_const_r_3+0x5c>
movsd  -0x10(%rbp),%xmm0
mov    0x18(%rbp),%eax
sub    0x10(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm2
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    0x18(%rbp),%ecx
sub    0x10(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1144 <sub_const_r_3+0x44>
jmp    1228 <sub_const_r_3+0x128>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    112c <sub_const_r_3+0x2c>
mov    -0x8(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<sub_const_r_3>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r9d,-0x2c(%rsp)
mov    %r8d,-0x30(%rsp)
mov    %rdi,-0x20(%rsp)
sub    %edx,%ecx
mov    %rcx,-0x28(%rsp)
jge    1132 <sub_const_r_3+0x32>
mov    -0x20(%rsp),%rax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    0x40(%rsp),%r10d
sub    0x38(%rsp),%r10d
lea    0x1(%r10),%edi
cvtsi2sd %edi,%xmm3
movslq %edi,%rdx
movapd 0xeb1(%rip),%xmm1
xorpd  %xmm0,%xmm1
mov    -0x2c(%rsp),%eax
mov    %eax,%r11d
sub    -0x30(%rsp),%r11d
add    $0x1,%r11d
mov    -0x28(%rsp),%rax
add    $0x1,%eax
mov    %rax,-0x28(%rsp)
mov    %edi,%r14d
and    $0xfffffffc,%r14d
lea    -0x4(%r14),%rax
mov    %rax,-0x8(%rsp)
mov    %rax,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
movapd %xmm0,%xmm2
unpcklpd %xmm0,%xmm2
mov    %rcx,%rax
mov    %rcx,-0x10(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,-0x18(%rsp)
mov    %edi,%r15d
and    $0x3,%r15d
xor    %ebp,%ebp
mulsd  %xmm3,%xmm1
jmp    11cf <sub_const_r_3+0xcf>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rbp
cmp    -0x28(%rsp),%rbp
je     1122 <sub_const_r_3+0x22>
mov    -0x30(%rsp),%eax
cmp    %eax,-0x2c(%rsp)
jl     11c0 <sub_const_r_3+0xc0>
mov    -0x20(%rsp),%rax
mov    (%rax,%rbp,8),%r12
xor    %eax,%eax
jmp    120d <sub_const_r_3+0x10d>
cs nopw 0x0(%rax,%rax,1)
mov    (%r12,%rax,8),%r9
movsd  (%r9,%rdx,8),%xmm3
addsd  %xmm1,%xmm3
movsd  %xmm3,(%r9,%rdx,8)
add    $0x1,%rax
cmp    %r11,%rax
je     11c0 <sub_const_r_3+0xc0>
test   %r10d,%r10d
js     11f0 <sub_const_r_3+0xf0>
mov    (%rsi,%rbp,8),%rcx
mov    (%rcx,%rax,8),%r8
mov    (%r12,%rax,8),%r9
cmp    $0x4,%edi
jb     123d <sub_const_r_3+0x13d>
lea    (%r8,%rdi,8),%rcx
cmp    %rcx,%r9
jae    12e4 <sub_const_r_3+0x1e4>
lea    (%r9,%rdi,8),%rcx
cmp    %rcx,%r8
jae    12e4 <sub_const_r_3+0x1e4>
xor    %r13d,%r13d
mov    %r13,%rcx
not    %rcx
add    %rdi,%rcx
test   %r15,%r15
je     127a <sub_const_r_3+0x17a>
mov    %r15,%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movapd %xmm0,%xmm3
subsd  (%r8,%r13,8),%xmm3
movsd  %xmm3,(%r9,%r13,8)
add    $0x1,%r13
add    $0xffffffffffffffff,%rbx
jne    1260 <sub_const_r_3+0x160>
cmp    $0x3,%rcx
jb     11f4 <sub_const_r_3+0xf4>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movapd %xmm0,%xmm3
subsd  (%r8,%r13,8),%xmm3
movsd  %xmm3,(%r9,%r13,8)
movapd %xmm0,%xmm3
subsd  0x8(%r8,%r13,8),%xmm3
movsd  %xmm3,0x8(%r9,%r13,8)
movapd %xmm0,%xmm3
subsd  0x10(%r8,%r13,8),%xmm3
movsd  %xmm3,0x10(%r9,%r13,8)
movapd %xmm0,%xmm3
subsd  0x18(%r8,%r13,8),%xmm3
movsd  %xmm3,0x18(%r9,%r13,8)
add    $0x4,%r13
cmp    %r13,%rdi
jne    1290 <sub_const_r_3+0x190>
jmp    11f4 <sub_const_r_3+0xf4>
cmpq   $0x0,-0x8(%rsp)
je     13a2 <sub_const_r_3+0x2a2>
mov    -0x18(%rsp),%rbx
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movupd (%r8,%rcx,8),%xmm3
movupd 0x10(%r8,%rcx,8),%xmm4
movapd %xmm2,%xmm5
subpd  %xmm3,%xmm5
movapd %xmm2,%xmm3
subpd  %xmm4,%xmm3
movupd %xmm5,(%r9,%rcx,8)
movupd %xmm3,0x10(%r9,%rcx,8)
movupd 0x20(%r8,%rcx,8),%xmm3
movupd 0x30(%r8,%rcx,8),%xmm4
movapd %xmm2,%xmm5
subpd  %xmm3,%xmm5
movapd %xmm2,%xmm3
subpd  %xmm4,%xmm3
movupd %xmm5,0x20(%r9,%rcx,8)
movupd %xmm3,0x30(%r9,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rbx
jne    1300 <sub_const_r_3+0x200>
testb  $0x1,-0x10(%rsp)
je     1391 <sub_const_r_3+0x291>
movupd (%r8,%rcx,8),%xmm3
movupd 0x10(%r8,%rcx,8),%xmm4
movapd %xmm2,%xmm5
subpd  %xmm3,%xmm5
movapd %xmm2,%xmm3
subpd  %xmm4,%xmm3
movupd %xmm5,(%r9,%rcx,8)
movupd %xmm3,0x10(%r9,%rcx,8)
mov    %r14,%r13
cmp    %rdi,%r14
je     11f4 <sub_const_r_3+0xf4>
jmp    1240 <sub_const_r_3+0x140>
xor    %ecx,%ecx
testb  $0x1,-0x10(%rsp)
jne    1367 <sub_const_r_3+0x267>
jmp    1391 <sub_const_r_3+0x291>

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

LLVM IR (-O2):
define double*** @sub_const_r_3(double*** returned %0, double %1, double*** nocapture %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7, i32 %8) {
%10 = icmp slt i32 %4, %3
br i1 %10, label %166, label %11
11: ; preds = %9
%12 = icmp slt i32 %6, %5
%13 = sub nsw i32 %8, %7
%14 = icmp slt i32 %13, 0
%15 = add nsw i32 %13, 1
%16 = sitofp i32 %15 to double
%17 = sext i32 %15 to i64
%18 = fneg double %1
%19 = add i32 %8, 1
%20 = sub i32 %19, %7
%21 = add i32 %6, 1
%22 = sub i32 %21, %5
%23 = add i32 %4, 1
%24 = sub i32 %23, %3
%25 = zext i32 %24 to i64
%26 = zext i32 %22 to i64
%27 = zext i32 %20 to i64
%28 = and i64 %27, 4294967292
%29 = add nsw i64 %28, -4
%30 = lshr exact i64 %29, 2
%31 = add nuw nsw i64 %30, 1
%32 = icmp ult i32 %20, 4
%33 = and i64 %27, 4294967292
%34 = insertelement <2 x double> poison, double %1, i64 0
%35 = shufflevector <2 x double> %34, <2 x double> poison, <2 x i32> zeroinitializer
%36 = insertelement <2 x double> poison, double %1, i64 0
%37 = shufflevector <2 x double> %36, <2 x double> poison, <2 x i32> zeroinitializer
%38 = and i64 %31, 1
%39 = icmp eq i64 %29, 0
%40 = and i64 %31, 9223372036854775806
%41 = icmp eq i64 %38, 0
%42 = icmp eq i64 %33, %27
%43 = and i64 %27, 3
%44 = icmp eq i64 %43, 0
br label %45
45: ; preds = %11, %163
%46 = phi i64 [ 0, %11 ], [ %164, %163 ]
br i1 %12, label %163, label %47
47: ; preds = %45
%48 = getelementptr inbounds double**, double*** %2, i64 %46
%49 = getelementptr inbounds double**, double*** %0, i64 %46
%50 = load double**, double*** %49
br label %51
51: ; preds = %47, %156
%52 = phi i64 [ 0, %47 ], [ %161, %156 ]
br i1 %14, label %53, label %56
53: ; preds = %51
%54 = getelementptr inbounds double*, double** %50, i64 %52
%55 = load double*, double** %54
br label %156
56: ; preds = %51
%57 = load double**, double*** %48
%58 = getelementptr inbounds double*, double** %57, i64 %52
%59 = load double*, double** %58
%60 = getelementptr inbounds double*, double** %50, i64 %52
%61 = load double*, double** %60
br i1 %32, label %116, label %62
62: ; preds = %56
%63 = getelementptr double, double* %61, i64 %27
%64 = getelementptr double, double* %59, i64 %27
%65 = icmp ult double* %61, %64
%66 = icmp ult double* %59, %63
%67 = and i1 %65, %66
br i1 %67, label %116, label %68
68: ; preds = %62
br i1 %39, label %100, label %69
69: ; preds = %68, %69
%70 = phi i64 [ %97, %69 ], [ 0, %68 ]
%71 = phi i64 [ %98, %69 ], [ 0, %68 ]
%72 = getelementptr inbounds double, double* %59, i64 %70
%73 = bitcast double* %72 to <2 x double>*
%74 = load <2 x double>, <2 x double>* %73
%75 = getelementptr inbounds double, double* %72, i64 2
%76 = bitcast double* %75 to <2 x double>*
%77 = load <2 x double>, <2 x double>* %76
%78 = fsub <2 x double> %35, %74
%79 = fsub <2 x double> %37, %77
%80 = getelementptr inbounds double, double* %61, i64 %70
%81 = bitcast double* %80 to <2 x double>*
store <2 x double> %78, <2 x double>* %81
%82 = getelementptr inbounds double, double* %80, i64 2
%83 = bitcast double* %82 to <2 x double>*
store <2 x double> %79, <2 x double>* %83
%84 = or i64 %70, 4
%85 = getelementptr inbounds double, double* %59, i64 %84
%86 = bitcast double* %85 to <2 x double>*
%87 = load <2 x double>, <2 x double>* %86
%88 = getelementptr inbounds double, double* %85, i64 2
%89 = bitcast double* %88 to <2 x double>*
%90 = load <2 x double>, <2 x double>* %89
%91 = fsub <2 x double> %35, %87
%92 = fsub <2 x double> %37, %90
%93 = getelementptr inbounds double, double* %61, i64 %84
%94 = bitcast double* %93 to <2 x double>*
store <2 x double> %91, <2 x double>* %94
%95 = getelementptr inbounds double, double* %93, i64 2
%96 = bitcast double* %95 to <2 x double>*
store <2 x double> %92, <2 x double>* %96
%97 = add nuw i64 %70, 8
%98 = add i64 %71, 2
%99 = icmp eq i64 %98, %40
br i1 %99, label %100, label %69
100: ; preds = %69, %68
%101 = phi i64 [ 0, %68 ], [ %97, %69 ]
br i1 %41, label %115, label %102
102: ; preds = %100
%103 = getelementptr inbounds double, double* %59, i64 %101
%104 = bitcast double* %103 to <2 x double>*
%105 = load <2 x double>, <2 x double>* %104
%106 = getelementptr inbounds double, double* %103, i64 2
%107 = bitcast double* %106 to <2 x double>*
%108 = load <2 x double>, <2 x double>* %107
%109 = fsub <2 x double> %35, %105
%110 = fsub <2 x double> %37, %108
%111 = getelementptr inbounds double, double* %61, i64 %101
%112 = bitcast double* %111 to <2 x double>*
store <2 x double> %109, <2 x double>* %112
%113 = getelementptr inbounds double, double* %111, i64 2
%114 = bitcast double* %113 to <2 x double>*
store <2 x double> %110, <2 x double>* %114
br label %115
115: ; preds = %100, %102
br i1 %42, label %156, label %116
116: ; preds = %62, %56, %115
%117 = phi i64 [ 0, %62 ], [ 0, %56 ], [ %33, %115 ]
%118 = xor i64 %117, -1
%119 = add nsw i64 %118, %27
br i1 %44, label %130, label %120
120: ; preds = %116, %120
%121 = phi i64 [ %127, %120 ], [ %117, %116 ]
%122 = phi i64 [ %128, %120 ], [ 0, %116 ]
%123 = getelementptr inbounds double, double* %59, i64 %121
%124 = load double, double* %123
%125 = fsub double %1, %124
%126 = getelementptr inbounds double, double* %61, i64 %121
store double %125, double* %126
%127 = add nuw nsw i64 %121, 1
%128 = add i64 %122, 1
%129 = icmp eq i64 %128, %43
br i1 %129, label %130, label %120
130: ; preds = %120, %116
%131 = phi i64 [ %117, %116 ], [ %127, %120 ]
%132 = icmp ult i64 %119, 3
br i1 %132, label %156, label %133
133: ; preds = %130, %133
%134 = phi i64 [ %154, %133 ], [ %131, %130 ]
%135 = getelementptr inbounds double, double* %59, i64 %134
%136 = load double, double* %135
%137 = fsub double %1, %136
%138 = getelementptr inbounds double, double* %61, i64 %134
store double %137, double* %138
%139 = add nuw nsw i64 %134, 1
%140 = getelementptr inbounds double, double* %59, i64 %139
%141 = load double, double* %140
%142 = fsub double %1, %141
%143 = getelementptr inbounds double, double* %61, i64 %139
store double %142, double* %143
%144 = add nuw nsw i64 %134, 2
%145 = getelementptr inbounds double, double* %59, i64 %144
%146 = load double, double* %145
%147 = fsub double %1, %146
%148 = getelementptr inbounds double, double* %61, i64 %144
store double %147, double* %148
%149 = add nuw nsw i64 %134, 3
%150 = getelementptr inbounds double, double* %59, i64 %149
%151 = load double, double* %150
%152 = fsub double %1, %151
%153 = getelementptr inbounds double, double* %61, i64 %149
store double %152, double* %153
%154 = add nuw nsw i64 %134, 4
%155 = icmp eq i64 %154, %27
br i1 %155, label %156, label %133
156: ; preds = %130, %133, %115, %53
%157 = phi double* [ %55, %53 ], [ %61, %115 ], [ %61, %133 ], [ %61, %130 ]
%158 = getelementptr inbounds double, double* %157, i64 %17
%159 = load double, double* %158
%160 = tail call double @llvm.fmuladd.f64(double %18, double %16, double %159)
store double %160, double* %158
%161 = add nuw nsw i64 %52, 1
%162 = icmp eq i64 %161, %26
br i1 %162, label %163, label %51
163: ; preds = %156, %45
%164 = add nuw nsw i64 %46, 1
%165 = icmp eq i64 %164, %25
br i1 %165, label %166, label %45
166: ; preds = %163, %9
ret double*** %0
}