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


x86 assembly (clang, -O0):
<mult_const_r_2>:
push   %rbp
mov    %rsp,%rbp
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     11a2 <mult_const_r_2+0xa2>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x28(%rbp),%ecx
sub    -0x24(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     118f <mult_const_r_2+0x8f>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
mulsd  (%rax,%rcx,8),%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    113e <mult_const_r_2+0x3e>
jmp    1194 <mult_const_r_2+0x94>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1126 <mult_const_r_2+0x26>
mov    -0x8(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<mult_const_r_2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rdi,%rax
sub    %edx,%ecx
jge    111c <mult_const_r_2+0x1c>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
sub    %r8d,%r9d
lea    0x1(%r9),%r15d
add    $0x2,%r9d
add    $0x1,%ecx
mov    %r9d,%r14d
and    $0xfffffffc,%r14d
lea    -0x4(%r14),%r8
mov    %r8,%r10
shr    $0x2,%r10
add    $0x1,%r10
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
mov    %r10,%rdx
and    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x8(%rsp)
mov    %r9d,%r11d
and    $0x3,%r11d
xor    %r13d,%r13d
jmp    1169 <mult_const_r_2+0x69>
add    $0x1,%r13
cmp    %rcx,%r13
je     1111 <mult_const_r_2+0x11>
test   %r15d,%r15d
js     1160 <mult_const_r_2+0x60>
mov    (%rsi,%r13,8),%rdx
mov    (%rax,%r13,8),%rbp
cmp    $0x4,%r9d
jb     119d <mult_const_r_2+0x9d>
lea    (%rdx,%r9,8),%rdi
cmp    %rdi,%rbp
jae    121d <mult_const_r_2+0x11d>
lea    0x0(,%r9,8),%rdi
add    %rbp,%rdi
cmp    %rdi,%rdx
jae    121d <mult_const_r_2+0x11d>
xor    %ebx,%ebx
mov    %rbx,%rdi
not    %rdi
add    %r9,%rdi
test   %r11,%r11
je     11c9 <mult_const_r_2+0xc9>
mov    %r11,%r12
movsd  (%rdx,%rbx,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x0(%rbp,%rbx,8)
add    $0x1,%rbx
add    $0xffffffffffffffff,%r12
jne    11b0 <mult_const_r_2+0xb0>
cmp    $0x3,%rdi
jb     1160 <mult_const_r_2+0x60>
nop
movsd  (%rdx,%rbx,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x0(%rbp,%rbx,8)
movsd  0x8(%rdx,%rbx,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x8(%rbp,%rbx,8)
movsd  0x10(%rdx,%rbx,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x10(%rbp,%rbx,8)
movsd  0x18(%rdx,%rbx,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x18(%rbp,%rbx,8)
add    $0x4,%rbx
cmp    %rbx,%r9
jne    11d0 <mult_const_r_2+0xd0>
jmp    1160 <mult_const_r_2+0x60>
test   %r8,%r8
je     12af <mult_const_r_2+0x1af>
mov    -0x8(%rsp),%rdi
xor    %ebx,%ebx
nopl   (%rax)
movupd (%rdx,%rbx,8),%xmm2
movupd 0x10(%rdx,%rbx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,0x0(%rbp,%rbx,8)
movupd %xmm3,0x10(%rbp,%rbx,8)
movupd 0x20(%rdx,%rbx,8),%xmm2
movupd 0x30(%rdx,%rbx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,0x20(%rbp,%rbx,8)
movupd %xmm3,0x30(%rbp,%rbx,8)
add    $0x8,%rbx
add    $0xfffffffffffffffe,%rdi
jne    1230 <mult_const_r_2+0x130>
test   $0x1,%r10b
je     129e <mult_const_r_2+0x19e>
movupd (%rdx,%rbx,8),%xmm2
movupd 0x10(%rdx,%rbx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,0x0(%rbp,%rbx,8)
movupd %xmm3,0x10(%rbp,%rbx,8)
mov    %r14,%rbx
cmp    %r9,%r14
je     1160 <mult_const_r_2+0x60>
jmp    119f <mult_const_r_2+0x9f>
xor    %ebx,%ebx
test   $0x1,%r10b
jne    127f <mult_const_r_2+0x17f>
jmp    129e <mult_const_r_2+0x19e>

LLVM IR (-O0):
define double** @mult_const_r_2(double** %0, double %1, double** %2, i32 %3, i32 %4, i32 %5, i32 %6) {
%8 = alloca double**
%9 = alloca double
%10 = alloca double**
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
store double** %0, double*** %8
store double %1, double* %9
store double** %2, double*** %10
store i32 %3, i32* %11
store i32 %4, i32* %12
store i32 %5, i32* %13
store i32 %6, i32* %14
store i32 0, i32* %15
br label %17
17: ; preds = %55, %7
%18 = load i32, i32* %15
%19 = load i32, i32* %12
%20 = load i32, i32* %11
%21 = sub nsw i32 %19, %20
%22 = icmp sle i32 %18, %21
br i1 %22, label %23, label %58
23: ; preds = %17
store i32 0, i32* %16
br label %24
24: ; preds = %51, %23
%25 = load i32, i32* %16
%26 = load i32, i32* %14
%27 = load i32, i32* %13
%28 = sub nsw i32 %26, %27
%29 = add nsw i32 %28, 1
%30 = icmp sle i32 %25, %29
br i1 %30, label %31, label %54
31: ; preds = %24
%32 = load double, double* %9
%33 = load double**, double*** %10
%34 = load i32, i32* %15
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double*, double** %33, i64 %35
%37 = load double*, double** %36
%38 = load i32, i32* %16
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds double, double* %37, i64 %39
%41 = load double, double* %40
%42 = fmul double %32, %41
%43 = load double**, double*** %8
%44 = load i32, i32* %15
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds double*, double** %43, i64 %45
%47 = load double*, double** %46
%48 = load i32, i32* %16
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %47, i64 %49
store double %42, double* %50
br label %51
51: ; preds = %31
%52 = load i32, i32* %16
%53 = add nsw i32 %52, 1
store i32 %53, i32* %16
br label %24
54: ; preds = %24
br label %55
55: ; preds = %54
%56 = load i32, i32* %15
%57 = add nsw i32 %56, 1
store i32 %57, i32* %15
br label %17
58: ; preds = %17
%59 = load double**, double*** %8
ret double** %59
}

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