void get_factor(double *vec, double *fact, int p) {
  int j, m;
  double tmp[2];
  for (j = 0; j < p - 1; ++j) {
    tmp[0] = 0.;
    tmp[1] = 0.;
    for (m = 0; m < (p - (j + 1)); ++m) {
      tmp[0] += vec[(m + j * p)];
      tmp[1] += vec[(m + (1 + j) * p)];
    }
    fact[j] = tmp[1] / tmp[0];
  }
}


x86 assembly (clang, -O0):
<get_factor>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    11ce <get_factor+0xce>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x14(%rbp),%ecx
mov    -0x18(%rbp),%edx
add    $0x1,%edx
sub    %edx,%ecx
cmp    %ecx,%eax
jge    11a9 <get_factor+0xa9>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    -0x18(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    -0x18(%rbp),%edx
add    $0x1,%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    113e <get_factor+0x3e>
movsd  -0x28(%rbp),%xmm0
divsd  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1116 <get_factor+0x16>
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_factor>:
push   %r15
push   %r14
push   %rbx
cmp    $0x2,%edx
jl     11dc <get_factor+0xdc>
mov    %edx,%r10d
add    $0xffffffff,%edx
mov    %edx,%r8d
lea    0x8(%rdi),%rcx
lea    0x0(,%r10,8),%r9
xor    %r14d,%r14d
movsd  0xed2(%rip),%xmm0
jmp    114c <get_factor+0x4c>
divsd  %xmm1,%xmm2
movsd  %xmm2,(%rsi,%r14,8)
add    $0xffffffff,%edx
add    %r9,%rcx
mov    %r11,%r14
cmp    %r8,%r11
je     11dc <get_factor+0xdc>
mov    %edx,%edx
lea    0x1(%r14),%r11
movapd %xmm0,%xmm2
cmp    %r10,%r11
jae    1134 <get_factor+0x34>
cmp    $0x1,%edx
jne    1170 <get_factor+0x70>
xorpd  %xmm1,%xmm1
xor    %ebx,%ebx
xorpd  %xmm2,%xmm2
jmp    11af <get_factor+0xaf>
nopl   0x0(%rax)
mov    %edx,%r15d
and    $0xfffffffe,%r15d
lea    (%rcx,%r9,1),%rax
xorpd  %xmm1,%xmm1
xor    %ebx,%ebx
xorpd  %xmm2,%xmm2
cs nopw 0x0(%rax,%rax,1)
nop
addsd  -0x8(%rcx,%rbx,8),%xmm1
addsd  -0x8(%rax,%rbx,8),%xmm2
addsd  (%rcx,%rbx,8),%xmm1
addsd  (%rax,%rbx,8),%xmm2
add    $0x2,%rbx
cmp    %rbx,%r15
jne    1190 <get_factor+0x90>
test   $0x1,%dl
je     1130 <get_factor+0x30>
mov    %r14,%rax
imul   %r10,%rax
mov    %r11,%r15
add    %rbx,%rax
addsd  (%rdi,%rax,8),%xmm1
imul   %r10,%r15
add    %rbx,%r15
addsd  (%rdi,%r15,8),%xmm2
jmp    1130 <get_factor+0x30>
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @get_factor(double* %0, double* %1, i32 %2) {
%4 = alloca double*
%5 = alloca double*
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca [2 x double]
store double* %0, double** %4
store double* %1, double** %5
store i32 %2, i32* %6
store i32 0, i32* %7
br label %10
10: ; preds = %64, %3
%11 = load i32, i32* %7
%12 = load i32, i32* %6
%13 = sub nsw i32 %12, 1
%14 = icmp slt i32 %11, %13
br i1 %14, label %15, label %67
15: ; preds = %10
%16 = getelementptr inbounds [2 x double], [2 x double]* %9, i64 0, i64 0
store double 0.000000e+00, double* %16
%17 = getelementptr inbounds [2 x double], [2 x double]* %9, i64 0, i64 1
store double 0.000000e+00, double* %17
store i32 0, i32* %8
br label %18
18: ; preds = %51, %15
%19 = load i32, i32* %8
%20 = load i32, i32* %6
%21 = load i32, i32* %7
%22 = add nsw i32 %21, 1
%23 = sub nsw i32 %20, %22
%24 = icmp slt i32 %19, %23
br i1 %24, label %25, label %54
25: ; preds = %18
%26 = load double*, double** %4
%27 = load i32, i32* %8
%28 = load i32, i32* %7
%29 = load i32, i32* %6
%30 = mul nsw i32 %28, %29
%31 = add nsw i32 %27, %30
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %26, i64 %32
%34 = load double, double* %33
%35 = getelementptr inbounds [2 x double], [2 x double]* %9, i64 0, i64 0
%36 = load double, double* %35
%37 = fadd double %36, %34
store double %37, double* %35
%38 = load double*, double** %4
%39 = load i32, i32* %8
%40 = load i32, i32* %7
%41 = add nsw i32 1, %40
%42 = load i32, i32* %6
%43 = mul nsw i32 %41, %42
%44 = add nsw i32 %39, %43
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds double, double* %38, i64 %45
%47 = load double, double* %46
%48 = getelementptr inbounds [2 x double], [2 x double]* %9, i64 0, i64 1
%49 = load double, double* %48
%50 = fadd double %49, %47
store double %50, double* %48
br label %51
51: ; preds = %25
%52 = load i32, i32* %8
%53 = add nsw i32 %52, 1
store i32 %53, i32* %8
br label %18
54: ; preds = %18
%55 = getelementptr inbounds [2 x double], [2 x double]* %9, i64 0, i64 1
%56 = load double, double* %55
%57 = getelementptr inbounds [2 x double], [2 x double]* %9, i64 0, i64 0
%58 = load double, double* %57
%59 = fdiv double %56, %58
%60 = load double*, double** %5
%61 = load i32, i32* %7
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %60, i64 %62
store double %59, double* %63
br label %64
64: ; preds = %54
%65 = load i32, i32* %7
%66 = add nsw i32 %65, 1
store i32 %66, i32* %7
br label %10
67: ; preds = %10
ret void
}

LLVM IR (-O2):
define void @get_factor(double* nocapture %0, double* nocapture writeonly %1, i32 %2) {
%4 = add i32 %2, -1
%5 = icmp sgt i32 %2, 1
br i1 %5, label %6, label %72
6: ; preds = %3
%7 = zext i32 %2 to i64
%8 = zext i32 %4 to i64
br label %9
9: ; preds = %6, %67
%10 = phi i64 [ 0, %6 ], [ %13, %67 ]
%11 = phi i32 [ %4, %6 ], [ %70, %67 ]
%12 = zext i32 %11 to i64
%13 = add nuw nsw i64 %10, 1
%14 = icmp ult i64 %13, %7
br i1 %14, label %15, label %67
15: ; preds = %9
%16 = mul nsw i64 %10, %7
%17 = mul nsw i64 %13, %7
%18 = and i64 %12, 1
%19 = icmp eq i32 %11, 1
br i1 %19, label %47, label %20
20: ; preds = %15
%21 = and i64 %12, 4294967294
br label %22
22: ; preds = %22, %20
%23 = phi i64 [ 0, %20 ], [ %44, %22 ]
%24 = phi double [ 0.000000e+00, %20 ], [ %39, %22 ]
%25 = phi double [ 0.000000e+00, %20 ], [ %43, %22 ]
%26 = phi i64 [ 0, %20 ], [ %45, %22 ]
%27 = add nuw nsw i64 %23, %16
%28 = getelementptr inbounds double, double* %0, i64 %27
%29 = load double, double* %28
%30 = fadd double %24, %29
%31 = add nuw nsw i64 %23, %17
%32 = getelementptr inbounds double, double* %0, i64 %31
%33 = load double, double* %32
%34 = fadd double %25, %33
%35 = or i64 %23, 1
%36 = add nuw nsw i64 %35, %16
%37 = getelementptr inbounds double, double* %0, i64 %36
%38 = load double, double* %37
%39 = fadd double %30, %38
%40 = add nuw nsw i64 %35, %17
%41 = getelementptr inbounds double, double* %0, i64 %40
%42 = load double, double* %41
%43 = fadd double %34, %42
%44 = add nuw nsw i64 %23, 2
%45 = add i64 %26, 2
%46 = icmp eq i64 %45, %21
br i1 %46, label %47, label %22
47: ; preds = %22, %15
%48 = phi double [ undef, %15 ], [ %39, %22 ]
%49 = phi double [ undef, %15 ], [ %43, %22 ]
%50 = phi i64 [ 0, %15 ], [ %44, %22 ]
%51 = phi double [ 0.000000e+00, %15 ], [ %39, %22 ]
%52 = phi double [ 0.000000e+00, %15 ], [ %43, %22 ]
%53 = icmp eq i64 %18, 0
br i1 %53, label %63, label %54
54: ; preds = %47
%55 = add nuw nsw i64 %50, %16
%56 = getelementptr inbounds double, double* %0, i64 %55
%57 = load double, double* %56
%58 = fadd double %51, %57
%59 = add nuw nsw i64 %50, %17
%60 = getelementptr inbounds double, double* %0, i64 %59
%61 = load double, double* %60
%62 = fadd double %52, %61
br label %63
63: ; preds = %47, %54
%64 = phi double [ %48, %47 ], [ %58, %54 ]
%65 = phi double [ %49, %47 ], [ %62, %54 ]
%66 = fdiv double %65, %64
br label %67
67: ; preds = %63, %9
%68 = phi double [ %66, %63 ], [ 0x7FF8000000000000, %9 ]
%69 = getelementptr inbounds double, double* %1, i64 %10
store double %68, double* %69
%70 = add i32 %11, -1
%71 = icmp eq i64 %13, %8
br i1 %71, label %72, label %9
72: ; preds = %67, %3
ret void
}