void poly_sum(double a[], int Na, double b[], int Nb, double ans[], int Nans) {
  int i;
  Nans = Na > Nb ? Na : Nb;
  for (i = 0; i < Nans; i++) {
    if (Na > Nb) {
      ans[i] = a[i];
      if (i > Na - Nb - 1) {
        ans[i] = ans[i] + b[i - Na + Nb];
      }
    } else {
      ans[i] = b[i];
      if (i > Nb - Na - 1) {
        ans[i] = ans[i] + a[i - Nb + Na];
      }
    }
  }
}


x86 assembly (clang, -O0):
<poly_sum>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jle    1131 <poly_sum+0x31>
mov    -0xc(%rbp),%eax
mov    %eax,-0x34(%rbp)
jmp    1137 <poly_sum+0x37>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    1233 <poly_sum+0x133>
mov    -0xc(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jle    11be <poly_sum+0xbe>
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jle    11b9 <poly_sum+0xb9>
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    -0xc(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1220 <poly_sum+0x120>
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jle    121b <poly_sum+0x11b>
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1220 <poly_sum+0x120>
jmp    1225 <poly_sum+0x125>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1144 <poly_sum+0x44>
pop    %rbp
ret

x86 assembly (clang, -O2):
<poly_sum>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %esi,%r11d
sub    %ecx,%r11d
mov    %ecx,%eax
cmovg  %esi,%eax
test   %eax,%eax
jle    1199 <poly_sum+0x99>
mov    %esi,%ebx
not    %ebx
add    %ecx,%ebx
mov    %ecx,%ebp
not    %ebp
add    %esi,%ebp
mov    %ecx,%r9d
sub    %esi,%r9d
movslq %ebp,%r10
movslq %ebx,%r14
mov    %eax,%r15d
xor    %eax,%eax
mov    %r8,%rbx
jmp    115f <poly_sum+0x5f>
nop
add    %eax,%ebp
movslq %ebp,%rbp
addsd  0x0(%r13,%rbp,8),%xmm0
movsd  %xmm0,(%r12)
add    $0x1,%rax
add    $0x8,%rbx
cmp    %rax,%r15
je     1199 <poly_sum+0x99>
cmp    %ecx,%esi
jle    1180 <poly_sum+0x80>
movsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,(%rbx)
cmp    %r10,%rax
jle    1152 <poly_sum+0x52>
lea    (%r8,%rax,8),%r12
mov    %r9d,%ebp
mov    %rdx,%r13
jmp    1140 <poly_sum+0x40>
nopl   (%rax)
movsd  (%rdx,%rax,8),%xmm0
movsd  %xmm0,(%rbx)
mov    %r11d,%ebp
mov    %rdi,%r13
mov    %rbx,%r12
cmp    %r14,%rax
jg     1140 <poly_sum+0x40>
jmp    1152 <poly_sum+0x52>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @poly_sum(double* %0, i32 %1, double* %2, i32 %3, double* %4, i32 %5) {
%7 = alloca double*
%8 = alloca i32
%9 = alloca double*
%10 = alloca i32
%11 = alloca double*
%12 = alloca i32
%13 = alloca i32
store double* %0, double** %7
store i32 %1, i32* %8
store double* %2, double** %9
store i32 %3, i32* %10
store double* %4, double** %11
store i32 %5, i32* %12
%14 = load i32, i32* %8
%15 = load i32, i32* %10
%16 = icmp sgt i32 %14, %15
br i1 %16, label %17, label %19
17: ; preds = %6
%18 = load i32, i32* %8
br label %21
19: ; preds = %6
%20 = load i32, i32* %10
br label %21
21: ; preds = %19, %17
%22 = phi i32 [ %18, %17 ], [ %20, %19 ]
store i32 %22, i32* %12
store i32 0, i32* %13
br label %23
23: ; preds = %106, %21
%24 = load i32, i32* %13
%25 = load i32, i32* %12
%26 = icmp slt i32 %24, %25
br i1 %26, label %27, label %109
27: ; preds = %23
%28 = load i32, i32* %8
%29 = load i32, i32* %10
%30 = icmp sgt i32 %28, %29
br i1 %30, label %31, label %68
31: ; preds = %27
%32 = load double*, double** %7
%33 = load i32, i32* %13
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds double, double* %32, i64 %34
%36 = load double, double* %35
%37 = load double*, double** %11
%38 = load i32, i32* %13
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds double, double* %37, i64 %39
store double %36, double* %40
%41 = load i32, i32* %13
%42 = load i32, i32* %8
%43 = load i32, i32* %10
%44 = sub nsw i32 %42, %43
%45 = sub nsw i32 %44, 1
%46 = icmp sgt i32 %41, %45
br i1 %46, label %47, label %67
47: ; preds = %31
%48 = load double*, double** %11
%49 = load i32, i32* %13
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double, double* %48, i64 %50
%52 = load double, double* %51
%53 = load double*, double** %9
%54 = load i32, i32* %13
%55 = load i32, i32* %8
%56 = sub nsw i32 %54, %55
%57 = load i32, i32* %10
%58 = add nsw i32 %56, %57
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double, double* %53, i64 %59
%61 = load double, double* %60
%62 = fadd double %52, %61
%63 = load double*, double** %11
%64 = load i32, i32* %13
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
store double %62, double* %66
br label %67
67: ; preds = %47, %31
br label %105
68: ; preds = %27
%69 = load double*, double** %9
%70 = load i32, i32* %13
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds double, double* %69, i64 %71
%73 = load double, double* %72
%74 = load double*, double** %11
%75 = load i32, i32* %13
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %74, i64 %76
store double %73, double* %77
%78 = load i32, i32* %13
%79 = load i32, i32* %10
%80 = load i32, i32* %8
%81 = sub nsw i32 %79, %80
%82 = sub nsw i32 %81, 1
%83 = icmp sgt i32 %78, %82
br i1 %83, label %84, label %104
84: ; preds = %68
%85 = load double*, double** %11
%86 = load i32, i32* %13
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double, double* %85, i64 %87
%89 = load double, double* %88
%90 = load double*, double** %7
%91 = load i32, i32* %13
%92 = load i32, i32* %10
%93 = sub nsw i32 %91, %92
%94 = load i32, i32* %8
%95 = add nsw i32 %93, %94
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds double, double* %90, i64 %96
%98 = load double, double* %97
%99 = fadd double %89, %98
%100 = load double*, double** %11
%101 = load i32, i32* %13
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds double, double* %100, i64 %102
store double %99, double* %103
br label %104
104: ; preds = %84, %68
br label %105
105: ; preds = %104, %67
br label %106
106: ; preds = %105
%107 = load i32, i32* %13
%108 = add nsw i32 %107, 1
store i32 %108, i32* %13
br label %23
109: ; preds = %23
ret void
}

LLVM IR (-O2):
define void @poly_sum(double* nocapture %0, i32 %1, double* nocapture %2, i32 %3, double* nocapture writeonly %4, i32 %5) {
%7 = icmp sgt i32 %1, %3
%8 = select i1 %7, i32 %1, i32 %3
%9 = icmp sgt i32 %8, 0
br i1 %9, label %10, label %46
10: ; preds = %6
%11 = xor i32 %1, -1
%12 = add i32 %11, %3
%13 = sub i32 %1, %3
%14 = xor i32 %3, -1
%15 = add i32 %14, %1
%16 = sub i32 %3, %1
%17 = sext i32 %15 to i64
%18 = sext i32 %12 to i64
%19 = zext i32 %8 to i64
br label %20
20: ; preds = %10, %43
%21 = phi i64 [ 0, %10 ], [ %44, %43 ]
br i1 %7, label %22, label %27
22: ; preds = %20
%23 = getelementptr inbounds double, double* %0, i64 %21
%24 = load double, double* %23
%25 = getelementptr inbounds double, double* %4, i64 %21
store double %24, double* %25
%26 = icmp sgt i64 %21, %17
br i1 %26, label %32, label %43
27: ; preds = %20
%28 = getelementptr inbounds double, double* %2, i64 %21
%29 = load double, double* %28
%30 = getelementptr inbounds double, double* %4, i64 %21
store double %29, double* %30
%31 = icmp sgt i64 %21, %18
br i1 %31, label %32, label %43
32: ; preds = %27, %22
%33 = phi i32 [ %16, %22 ], [ %13, %27 ]
%34 = phi double* [ %2, %22 ], [ %0, %27 ]
%35 = phi double [ %24, %22 ], [ %29, %27 ]
%36 = phi double* [ %25, %22 ], [ %30, %27 ]
%37 = trunc i64 %21 to i32
%38 = add i32 %33, %37
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds double, double* %34, i64 %39
%41 = load double, double* %40
%42 = fadd double %35, %41
store double %42, double* %36
br label %43
43: ; preds = %32, %22, %27
%44 = add nuw nsw i64 %21, 1
%45 = icmp eq i64 %44, %19
br i1 %45, label %46, label %20
46: ; preds = %43, %6
ret void
}