void diag_contr(int num_SD, int factor, double *sd_diag, double *init,
                double *final) {
  if (num_SD == 0)
    return;
  do {
    *(final++) = factor * (*(init++)) * (*(sd_diag++));
  } while (--num_SD);
  return;
}


x86 assembly (clang, -O0):
<diag_contr>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
cmpl   $0x0,-0x4(%rbp)
jne    1125 <diag_contr+0x25>
jmp    117f <diag_contr+0x7f>
jmp    112a <diag_contr+0x2a>
cvtsi2sdl -0x8(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x18(%rbp)
mulsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x10(%rbp)
mulsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x20(%rbp)
movsd  %xmm0,(%rax)
mov    -0x4(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x4(%rbp)
cmp    $0x0,%eax
jne    112a <diag_contr+0x2a>
jmp    117f <diag_contr+0x7f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<diag_contr>:
test   %edi,%edi
je     1301 <diag_contr+0x201>
cvtsi2sd %esi,%xmm0
lea    -0x1(%rdi),%r11d
cmp    $0x3,%r11d
jae    11d6 <diag_contr+0xd6>
lea    -0x1(%rdi),%r9d
mov    %edi,%esi
and    $0x3,%esi
je     1156 <diag_contr+0x56>
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
movsd  (%rcx),%xmm1
mulsd  %xmm0,%xmm1
add    $0x8,%rcx
mulsd  (%rdx),%xmm1
add    $0x8,%rdx
movsd  %xmm1,(%r8)
add    $0x8,%r8
add    $0x1,%eax
cmp    %eax,%esi
jne    1130 <diag_contr+0x30>
sub    %eax,%edi
cmp    $0x3,%r9d
jb     1301 <diag_contr+0x201>
mov    %edi,%eax
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rcx,%rsi,8),%xmm1
mulsd  %xmm0,%xmm1
mulsd  (%rdx,%rsi,8),%xmm1
movsd  %xmm1,(%r8,%rsi,8)
movsd  0x8(%rcx,%rsi,8),%xmm1
mulsd  %xmm0,%xmm1
mulsd  0x8(%rdx,%rsi,8),%xmm1
movsd  %xmm1,0x8(%r8,%rsi,8)
movsd  0x10(%rcx,%rsi,8),%xmm1
mulsd  %xmm0,%xmm1
mulsd  0x10(%rdx,%rsi,8),%xmm1
movsd  %xmm1,0x10(%r8,%rsi,8)
movsd  0x18(%rcx,%rsi,8),%xmm1
mulsd  %xmm0,%xmm1
mulsd  0x18(%rdx,%rsi,8),%xmm1
movsd  %xmm1,0x18(%r8,%rsi,8)
add    $0x4,%rsi
cmp    %esi,%eax
jne    1170 <diag_contr+0x70>
jmp    1301 <diag_contr+0x201>
lea    (%rcx,%r11,8),%rax
add    $0x8,%rax
cmp    %r8,%rax
seta   %r9b
lea    (%r8,%r11,8),%rsi
add    $0x8,%rsi
cmp    %rcx,%rsi
seta   %r10b
lea    (%rdx,%r11,8),%rax
add    $0x8,%rax
cmp    %r8,%rax
seta   %al
cmp    %rdx,%rsi
seta   %sil
test   %r10b,%r9b
jne    111a <diag_contr+0x1a>
and    %sil,%al
jne    111a <diag_contr+0x1a>
add    $0x1,%r11
mov    %r11,%r9
and    $0xfffffffffffffffc,%r9
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%r9),%rax
mov    %rax,%r10
shr    $0x2,%r10
add    $0x1,%r10
test   %rax,%rax
je     1316 <diag_contr+0x216>
mov    %r10,%rax
and    $0xfffffffffffffffe,%rax
xor    %esi,%esi
nop
movupd (%rcx,%rsi,8),%xmm2
movupd 0x10(%rcx,%rsi,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd (%rdx,%rsi,8),%xmm4
mulpd  %xmm2,%xmm4
movupd 0x10(%rdx,%rsi,8),%xmm2
mulpd  %xmm3,%xmm2
movupd %xmm4,(%r8,%rsi,8)
movupd %xmm2,0x10(%r8,%rsi,8)
movupd 0x20(%rcx,%rsi,8),%xmm2
movupd 0x30(%rcx,%rsi,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd 0x20(%rdx,%rsi,8),%xmm4
mulpd  %xmm2,%xmm4
movupd 0x30(%rdx,%rsi,8),%xmm2
mulpd  %xmm3,%xmm2
movupd %xmm4,0x20(%r8,%rsi,8)
movupd %xmm2,0x30(%r8,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rax
jne    1250 <diag_contr+0x150>
test   $0x1,%r10b
je     12fc <diag_contr+0x1fc>
movupd (%rcx,%rsi,8),%xmm2
movupd 0x10(%rcx,%rsi,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd (%rdx,%rsi,8),%xmm1
mulpd  %xmm2,%xmm1
movupd 0x10(%rdx,%rsi,8),%xmm2
mulpd  %xmm3,%xmm2
movupd %xmm1,(%r8,%rsi,8)
movupd %xmm2,0x10(%r8,%rsi,8)
cmp    %r9,%r11
jne    1302 <diag_contr+0x202>
ret
sub    %r9d,%edi
lea    (%rdx,%r9,8),%rdx
lea    (%rcx,%r9,8),%rcx
lea    (%r8,%r9,8),%r8
jmp    111a <diag_contr+0x1a>
xor    %esi,%esi
test   $0x1,%r10b
jne    12c9 <diag_contr+0x1c9>
jmp    12fc <diag_contr+0x1fc>

LLVM IR (-O0):
define void @diag_contr(i32 %0, i32 %1, double* %2, double* %3, double* %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca double*
%9 = alloca double*
%10 = alloca double*
store i32 %0, i32* %6
store i32 %1, i32* %7
store double* %2, double** %8
store double* %3, double** %9
store double* %4, double** %10
%11 = load i32, i32* %6
%12 = icmp eq i32 %11, 0
br i1 %12, label %13, label %14
13: ; preds = %5
br label %33
14: ; preds = %5
br label %15
15: ; preds = %28, %14
%16 = load i32, i32* %7
%17 = sitofp i32 %16 to double
%18 = load double*, double** %9
%19 = getelementptr inbounds double, double* %18, i32 1
store double* %19, double** %9
%20 = load double, double* %18
%21 = fmul double %17, %20
%22 = load double*, double** %8
%23 = getelementptr inbounds double, double* %22, i32 1
store double* %23, double** %8
%24 = load double, double* %22
%25 = fmul double %21, %24
%26 = load double*, double** %10
%27 = getelementptr inbounds double, double* %26, i32 1
store double* %27, double** %10
store double %25, double* %26
br label %28
28: ; preds = %15
%29 = load i32, i32* %6
%30 = add nsw i32 %29, -1
store i32 %30, i32* %6
%31 = icmp ne i32 %30, 0
br i1 %31, label %15, label %32
32: ; preds = %28
br label %33
33: ; preds = %32, %13
ret void
}

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