void c_model2(int *nin, int *nsin, double *x, double *x1, double *titre,
              double *titrepred, double *dd, double *mu) {
  int n = nin[0];
  int nsamp = nsin[0];
  int k;
  int i;
  double xx2;
  for (k = 0; k < nsamp; k++) {
    xx2 = 0;
    for (i = 0; i < n; i++) {
      x1[i] = mu[0] * dd[k * n + i] * x[i];
    }
    for (i = 0; i < n; i++) {
      xx2 = xx2 + x1[i];
    }
    titrepred[k] = xx2;
  }
}


x86 assembly (clang, -O0):
<c_model2>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x34(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    1202 <c_model2+0x102>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    11aa <c_model2+0xaa>
mov    0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x34(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1158 <c_model2+0x58>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    11e2 <c_model2+0xe2>
movsd  -0x48(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11b1 <c_model2+0xb1>
movsd  -0x48(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    113d <c_model2+0x3d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<c_model2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r9,-0x8(%rsp)
mov    (%rsi),%eax
mov    %rax,-0x10(%rsp)
test   %eax,%eax
jle    1373 <c_model2+0x273>
mov    0x40(%rsp),%rsi
mov    0x38(%rsp),%r8
movslq (%rdi),%r9
mov    %r9d,%eax
lea    (%rcx,%rax,8),%rbp
lea    0x8(%rsi),%rdi
lea    (%rdx,%rax,8),%r10
lea    -0x1(%rax),%rbx
mov    %rbx,-0x30(%rsp)
cmp    %rcx,%rdi
seta   %dil
cmp    %rsi,%rbp
seta   %bl
and    %dil,%bl
mov    %bl,-0x31(%rsp)
cmp    %rcx,%r10
seta   %dil
mov    %rbp,-0x20(%rsp)
cmp    %rdx,%rbp
seta   %bl
and    %dil,%bl
mov    %bl,-0x32(%rsp)
mov    %eax,%edi
and    $0xfffffffc,%edi
mov    %rdi,-0x28(%rsp)
mov    %eax,%ebp
and    $0x7,%ebp
mov    %eax,%r10d
and    $0xfffffff8,%r10d
lea    0x10(%r8),%r11
lea    0x0(,%r9,8),%r15
lea    0x0(,%rax,8),%r14
and    $0xffffffffffffffe0,%r14
lea    0x8(%r8),%r12
xor    %r8d,%r8d
mov    %r9,-0x18(%rsp)
jmp    11d0 <c_model2+0xd0>
nopl   0x0(%rax)
mov    -0x8(%rsp),%rdi
movsd  %xmm0,(%rdi,%r8,8)
add    $0x1,%r8
add    %r15,%r11
add    %r15,%r12
cmp    -0x10(%rsp),%r8
je     1373 <c_model2+0x273>
xorpd  %xmm0,%xmm0
test   %r9d,%r9d
jle    11b0 <c_model2+0xb0>
mov    %r8,%r13
imul   %r9,%r13
cmp    $0x4,%eax
jb     1212 <c_model2+0x112>
lea    (%rax,%r13,1),%rdi
mov    0x38(%rsp),%r9
lea    (%r9,%r13,8),%rbx
lea    (%r9,%rdi,8),%rdi
cmp    %rcx,%rdi
seta   %dil
cmp    -0x20(%rsp),%rbx
setb   %bl
and    %dil,%bl
or     -0x31(%rsp),%bl
or     -0x32(%rsp),%bl
je     1286 <c_model2+0x186>
xor    %r9d,%r9d
mov    %r9,%rdi
test   $0x1,%al
je     1241 <c_model2+0x141>
movsd  (%rsi),%xmm1
add    %r9,%r13
mov    0x38(%rsp),%rdi
mulsd  (%rdi,%r13,8),%xmm1
mulsd  (%rdx,%r9,8),%xmm1
movsd  %xmm1,(%rcx,%r9,8)
mov    %r9,%rdi
or     $0x1,%rdi
cmp    %r9,-0x30(%rsp)
je     12e0 <c_model2+0x1e0>
nopl   0x0(%rax)
movsd  (%rsi),%xmm1
mulsd  -0x8(%r12,%rdi,8),%xmm1
mulsd  (%rdx,%rdi,8),%xmm1
movsd  %xmm1,(%rcx,%rdi,8)
movsd  (%rsi),%xmm1
mulsd  (%r12,%rdi,8),%xmm1
mulsd  0x8(%rdx,%rdi,8),%xmm1
movsd  %xmm1,0x8(%rcx,%rdi,8)
add    $0x2,%rdi
cmp    %rdi,%rax
jne    1250 <c_model2+0x150>
jmp    12e0 <c_model2+0x1e0>
movsd  (%rsi),%xmm1
unpcklpd %xmm1,%xmm1
xor    %edi,%edi
movupd -0x10(%r11,%rdi,1),%xmm2
movupd (%r11,%rdi,1),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd (%rdx,%rdi,1),%xmm4
mulpd  %xmm2,%xmm4
movupd 0x10(%rdx,%rdi,1),%xmm2
mulpd  %xmm3,%xmm2
movupd %xmm4,(%rcx,%rdi,1)
movupd %xmm2,0x10(%rcx,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%r14
jne    1290 <c_model2+0x190>
mov    -0x28(%rsp),%rdi
mov    %rdi,%r9
cmp    %rax,%rdi
jne    1215 <c_model2+0x115>
nopl   (%rax)
test   %eax,%eax
mov    -0x18(%rsp),%r9
jle    11b0 <c_model2+0xb0>
cmpq   $0x7,-0x30(%rsp)
jae    12fd <c_model2+0x1fd>
xorpd  %xmm0,%xmm0
xor    %edi,%edi
jmp    1348 <c_model2+0x248>
xorpd  %xmm0,%xmm0
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rcx,%rdi,8),%xmm0
addsd  0x8(%rcx,%rdi,8),%xmm0
addsd  0x10(%rcx,%rdi,8),%xmm0
addsd  0x18(%rcx,%rdi,8),%xmm0
addsd  0x20(%rcx,%rdi,8),%xmm0
addsd  0x28(%rcx,%rdi,8),%xmm0
addsd  0x30(%rcx,%rdi,8),%xmm0
addsd  0x38(%rcx,%rdi,8),%xmm0
add    $0x8,%rdi
cmp    %rdi,%r10
jne    1310 <c_model2+0x210>
test   %rbp,%rbp
je     11b0 <c_model2+0xb0>
lea    (%rcx,%rdi,8),%rdi
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
addsd  (%rdi,%rbx,8),%xmm0
add    $0x1,%rbx
cmp    %rbx,%rbp
jne    1360 <c_model2+0x260>
jmp    11b0 <c_model2+0xb0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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