double get_bij(int i, int j, double *adirissa, double *weight) {
  double ret;
  if (j <= i)
    ret = weight[j] * adirissa[j] * (1.0 - adirissa[i]);
  else
    ret = weight[j] * adirissa[i] * (1.0 - adirissa[j]);
  return ret;
}


x86 assembly (clang, -O0):
<get_bij>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    -0x8(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     115b <get_bij+0x5b>
mov    -0x18(%rbp),%rax
movslq -0x8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x8(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
movsd  0xeb8(%rip),%xmm1
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    1193 <get_bij+0x93>
mov    -0x18(%rbp),%rax
movslq -0x8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x8(%rbp),%rcx
movsd  0xe7b(%rip),%xmm1
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_bij>:
movslq %esi,%rax
movsd  (%rcx,%rax,8),%xmm1
cmp    %edi,%eax
mov    %edi,%ecx
cmovg  %eax,%ecx
cmovl  %eax,%edi
movslq %ecx,%rax
movslq %edi,%rcx
mulsd  (%rdx,%rcx,8),%xmm1
movsd  0xedb(%rip),%xmm0
subsd  (%rdx,%rax,8),%xmm0
mulsd  %xmm1,%xmm0
ret

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

LLVM IR (-O2):
define double @get_bij(i32 %0, i32 %1, double* nocapture %2, double* nocapture %3) {
%5 = sext i32 %1 to i64
%6 = getelementptr inbounds double, double* %3, i64 %5
%7 = load double, double* %6
%8 = icmp sgt i32 %1, %0
%9 = select i1 %8, i32 %1, i32 %0
%10 = sext i32 %9 to i64
%11 = icmp slt i32 %1, %0
%12 = select i1 %11, i32 %1, i32 %0
%13 = sext i32 %12 to i64
%14 = getelementptr inbounds double, double* %2, i64 %13
%15 = load double, double* %14
%16 = fmul double %7, %15
%17 = getelementptr inbounds double, double* %2, i64 %10
%18 = load double, double* %17
%19 = fsub double 1.000000e+00, %18
%20 = fmul double %16, %19
ret double %20
}