void ComputeT(double *h0, double *h1, int M, double *alfa, double sigma,
              double *t0, double *tl, double *tu) {
  int i = 0;
  for (i = 0; i < M; i++) {
    t0[i] = (h0[i] * alfa[i] * sigma) + 1;
    if (i < M - 1) {
      tu[i] = h1[i] * alfa[i + 1] * sigma;
      tl[i] = h1[i] * alfa[i] * sigma;
    }
  }
}


x86 assembly (clang, -O0):
<ComputeT>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11fb <ComputeT+0xfb>
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xe97(%rip),%xmm1
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    11e8 <ComputeT+0xe8>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x28(%rbp),%xmm0
mov    0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x28(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    11ed <ComputeT+0xed>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1132 <ComputeT+0x32>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ComputeT>:
test   %edx,%edx
jle    1178 <ComputeT+0x78>
mov    0x8(%rsp),%r10
lea    -0x1(%rdx),%r11d
mov    %edx,%edx
xor    %eax,%eax
movsd  0xee7(%rip),%xmm1
jmp    1152 <ComputeT+0x52>
nopl   0x0(%rax,%rax,1)
movsd  (%rsi,%rax,8),%xmm2
mulsd  0x8(%rcx,%rax,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,(%r10,%rax,8)
movsd  (%rsi,%rax,8),%xmm2
mulsd  (%rcx,%rax,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,(%r9,%rax,8)
add    $0x1,%rax
cmp    %rdx,%rax
je     1178 <ComputeT+0x78>
movsd  (%rdi,%rax,8),%xmm2
mulsd  (%rcx,%rax,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%r8,%rax,8)
cmp    %r11,%rax
jl     1120 <ComputeT+0x20>
add    $0x1,%rax
cmp    %rdx,%rax
jne    1152 <ComputeT+0x52>
ret

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

LLVM IR (-O2):
define void @ComputeT(double* nocapture %0, double* nocapture %1, i32 %2, double* nocapture %3, double %4, double* nocapture writeonly %5, double* nocapture writeonly %6, double* nocapture writeonly %7) {
%9 = icmp sgt i32 %2, 0
br i1 %9, label %10, label %43
10: ; preds = %8
%11 = add nsw i32 %2, -1
%12 = zext i32 %11 to i64
%13 = zext i32 %2 to i64
br label %14
14: ; preds = %10, %40
%15 = phi i64 [ 0, %10 ], [ %41, %40 ]
%16 = getelementptr inbounds double, double* %0, i64 %15
%17 = load double, double* %16
%18 = getelementptr inbounds double, double* %3, i64 %15
%19 = load double, double* %18
%20 = fmul double %17, %19
%21 = tail call double @llvm.fmuladd.f64(double %20, double %4, double 1.000000e+00)
%22 = getelementptr inbounds double, double* %5, i64 %15
store double %21, double* %22
%23 = icmp slt i64 %15, %12
br i1 %23, label %26, label %24
24: ; preds = %14
%25 = add nuw nsw i64 %15, 1
br label %40
26: ; preds = %14
%27 = getelementptr inbounds double, double* %1, i64 %15
%28 = load double, double* %27
%29 = add nuw nsw i64 %15, 1
%30 = getelementptr inbounds double, double* %3, i64 %29
%31 = load double, double* %30
%32 = fmul double %28, %31
%33 = fmul double %32, %4
%34 = getelementptr inbounds double, double* %7, i64 %15
store double %33, double* %34
%35 = load double, double* %27
%36 = load double, double* %18
%37 = fmul double %35, %36
%38 = fmul double %37, %4
%39 = getelementptr inbounds double, double* %6, i64 %15
store double %38, double* %39
br label %40
40: ; preds = %24, %26
%41 = phi i64 [ %25, %24 ], [ %29, %26 ]
%42 = icmp eq i64 %41, %13
br i1 %42, label %43, label %14
43: ; preds = %40, %8
ret void
}