void qtrans(double *q, int m, double a) {
  int i;
  double d, o;
  m += m;
  i = 1;
  d = q[i];
  for (i++; i <= m; i++) {
    o = q[i] + a * d;
    d = q[i];
    q[i] = o;
  }
  return;
}


x86 assembly (clang, -O0):
<qtrans>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    -0xc(%rbp),%eax
add    -0xc(%rbp),%eax
mov    %eax,-0xc(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     119d <qtrans+0x9d>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x18(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    113b <qtrans+0x3b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<qtrans>:
test   %esi,%esi
jle    124f <qtrans+0x14f>
movsd  0x8(%rdi),%xmm1
add    %esi,%esi
cmp    $0x3,%esi
mov    $0x2,%eax
cmovge %esi,%eax
or     $0x1,%eax
lea    -0x2(%rax),%rcx
mov    $0x2,%edx
cmp    $0x4,%rcx
jb     1230 <qtrans+0x130>
and    $0xfffffffffffffffc,%rcx
unpcklpd %xmm1,%xmm1
movapd %xmm0,%xmm2
unpcklpd %xmm0,%xmm2
lea    -0x4(%rcx),%rdx
mov    %rdx,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rdx,%rdx
je     1250 <qtrans+0x150>
mov    %r8,%rsi
and    $0xfffffffffffffffe,%rsi
mov    $0x2,%edx
movapd %xmm1,%xmm3
nopl   0x0(%rax,%rax,1)
movupd (%rdi,%rdx,8),%xmm4
movupd 0x10(%rdi,%rdx,8),%xmm5
movupd 0x20(%rdi,%rdx,8),%xmm6
movupd 0x30(%rdi,%rdx,8),%xmm1
shufpd $0x1,%xmm4,%xmm3
mulpd  %xmm2,%xmm3
addpd  %xmm4,%xmm3
shufpd $0x1,%xmm5,%xmm4
mulpd  %xmm2,%xmm4
addpd  %xmm5,%xmm4
movupd %xmm3,(%rdi,%rdx,8)
movupd %xmm4,0x10(%rdi,%rdx,8)
shufpd $0x1,%xmm6,%xmm5
mulpd  %xmm2,%xmm5
addpd  %xmm6,%xmm5
shufpd $0x1,%xmm1,%xmm6
mulpd  %xmm2,%xmm6
addpd  %xmm1,%xmm6
movupd %xmm5,0x20(%rdi,%rdx,8)
movupd %xmm6,0x30(%rdi,%rdx,8)
add    $0x8,%rdx
movapd %xmm1,%xmm3
add    $0xfffffffffffffffe,%rsi
jne    1170 <qtrans+0x70>
test   $0x1,%r8b
je     121a <qtrans+0x11a>
movupd (%rdi,%rdx,8),%xmm3
movupd 0x10(%rdi,%rdx,8),%xmm4
shufpd $0x1,%xmm3,%xmm1
mulpd  %xmm2,%xmm1
addpd  %xmm3,%xmm1
shufpd $0x1,%xmm4,%xmm3
mulpd  %xmm2,%xmm3
addpd  %xmm4,%xmm3
movupd %xmm1,(%rdi,%rdx,8)
movupd %xmm3,0x10(%rdi,%rdx,8)
movapd %xmm4,%xmm1
or     $0x2,%rcx
unpckhpd %xmm1,%xmm1
mov    %rcx,%rdx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rdi,%rdx,8),%xmm2
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,(%rdi,%rdx,8)
add    $0x1,%rdx
movapd %xmm2,%xmm1
cmp    %rdx,%rax
jne    1230 <qtrans+0x130>
ret
mov    $0x2,%edx
test   $0x1,%r8b
jne    11e6 <qtrans+0xe6>
jmp    121a <qtrans+0x11a>

LLVM IR (-O0):
define void @qtrans(double* %0, i32 %1, double %2) {
%4 = alloca double*
%5 = alloca i32
%6 = alloca double
%7 = alloca i32
%8 = alloca double
%9 = alloca double
store double* %0, double** %4
store i32 %1, i32* %5
store double %2, double* %6
%10 = load i32, i32* %5
%11 = load i32, i32* %5
%12 = add nsw i32 %11, %10
store i32 %12, i32* %5
store i32 1, i32* %7
%13 = load double*, double** %4
%14 = load i32, i32* %7
%15 = sext i32 %14 to i64
%16 = getelementptr inbounds double, double* %13, i64 %15
%17 = load double, double* %16
store double %17, double* %8
%18 = load i32, i32* %7
%19 = add nsw i32 %18, 1
store i32 %19, i32* %7
br label %20
20: ; preds = %43, %3
%21 = load i32, i32* %7
%22 = load i32, i32* %5
%23 = icmp sle i32 %21, %22
br i1 %23, label %24, label %46
24: ; preds = %20
%25 = load double*, double** %4
%26 = load i32, i32* %7
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
%29 = load double, double* %28
%30 = load double, double* %6
%31 = load double, double* %8
%32 = call double @llvm.fmuladd.f64(double %30, double %31, double %29)
store double %32, double* %9
%33 = load double*, double** %4
%34 = load i32, i32* %7
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double, double* %33, i64 %35
%37 = load double, double* %36
store double %37, double* %8
%38 = load double, double* %9
%39 = load double*, double** %4
%40 = load i32, i32* %7
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds double, double* %39, i64 %41
store double %38, double* %42
br label %43
43: ; preds = %24
%44 = load i32, i32* %7
%45 = add nsw i32 %44, 1
store i32 %45, i32* %7
br label %20
46: ; preds = %20
ret void
}

LLVM IR (-O2):
define void @qtrans(double* nocapture %0, i32 %1, double %2) {
%4 = icmp slt i32 %1, 1
br i1 %4, label %96, label %5
5: ; preds = %3
%6 = getelementptr inbounds double, double* %0, i64 1
%7 = load double, double* %6
%8 = shl nuw i32 %1, 1
%9 = call i32 @llvm.smax.i32(i32 %8, i32 2)
%10 = or i32 %9, 1
%11 = zext i32 %10 to i64
%12 = add nsw i64 %11, -2
%13 = icmp ult i64 %12, 4
br i1 %13, label %85, label %14
14: ; preds = %5
%15 = and i64 %12, -4
%16 = or i64 %15, 2
%17 = insertelement <2 x double> poison, double %7, i64 1
%18 = insertelement <2 x double> poison, double %2, i64 0
%19 = shufflevector <2 x double> %18, <2 x double> poison, <2 x i32> zeroinitializer
%20 = insertelement <2 x double> poison, double %2, i64 0
%21 = shufflevector <2 x double> %20, <2 x double> poison, <2 x i32> zeroinitializer
%22 = add nsw i64 %15, -4
%23 = lshr exact i64 %22, 2
%24 = add nuw nsw i64 %23, 1
%25 = and i64 %24, 1
%26 = icmp eq i64 %22, 0
br i1 %26, label %64, label %27
27: ; preds = %14
%28 = and i64 %24, 9223372036854775806
br label %29
29: ; preds = %29, %27
%30 = phi i64 [ 0, %27 ], [ %59, %29 ]
%31 = phi <2 x double> [ %17, %27 ], [ %52, %29 ]
%32 = phi i64 [ 0, %27 ], [ %60, %29 ]
%33 = or i64 %30, 2
%34 = getelementptr inbounds double, double* %0, i64 %33
%35 = bitcast double* %34 to <2 x double>*
%36 = load <2 x double>, <2 x double>* %35
%37 = getelementptr inbounds double, double* %34, i64 2
%38 = bitcast double* %37 to <2 x double>*
%39 = load <2 x double>, <2 x double>* %38
%40 = shufflevector <2 x double> %31, <2 x double> %36, <2 x i32> <i32 1, i32 2>
%41 = shufflevector <2 x double> %36, <2 x double> %39, <2 x i32> <i32 1, i32 2>
%42 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %19, <2 x double> %40, <2 x double> %36)
%43 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %21, <2 x double> %41, <2 x double> %39)
%44 = bitcast double* %34 to <2 x double>*
store <2 x double> %42, <2 x double>* %44
%45 = bitcast double* %37 to <2 x double>*
store <2 x double> %43, <2 x double>* %45
%46 = or i64 %30, 6
%47 = getelementptr inbounds double, double* %0, i64 %46
%48 = bitcast double* %47 to <2 x double>*
%49 = load <2 x double>, <2 x double>* %48
%50 = getelementptr inbounds double, double* %47, i64 2
%51 = bitcast double* %50 to <2 x double>*
%52 = load <2 x double>, <2 x double>* %51
%53 = shufflevector <2 x double> %39, <2 x double> %49, <2 x i32> <i32 1, i32 2>
%54 = shufflevector <2 x double> %49, <2 x double> %52, <2 x i32> <i32 1, i32 2>
%55 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %19, <2 x double> %53, <2 x double> %49)
%56 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %21, <2 x double> %54, <2 x double> %52)
%57 = bitcast double* %47 to <2 x double>*
store <2 x double> %55, <2 x double>* %57
%58 = bitcast double* %50 to <2 x double>*
store <2 x double> %56, <2 x double>* %58
%59 = add nuw i64 %30, 8
%60 = add i64 %32, 2
%61 = icmp eq i64 %60, %28
br i1 %61, label %62, label %29
62: ; preds = %29
%63 = or i64 %59, 2
br label %64
64: ; preds = %62, %14
%65 = phi <2 x double> [ undef, %14 ], [ %52, %62 ]
%66 = phi i64 [ 2, %14 ], [ %63, %62 ]
%67 = phi <2 x double> [ %17, %14 ], [ %52, %62 ]
%68 = icmp eq i64 %25, 0
br i1 %68, label %82, label %69
69: ; preds = %64
%70 = getelementptr inbounds double, double* %0, i64 %66
%71 = bitcast double* %70 to <2 x double>*
%72 = load <2 x double>, <2 x double>* %71
%73 = getelementptr inbounds double, double* %70, i64 2
%74 = bitcast double* %73 to <2 x double>*
%75 = load <2 x double>, <2 x double>* %74
%76 = shufflevector <2 x double> %67, <2 x double> %72, <2 x i32> <i32 1, i32 2>
%77 = shufflevector <2 x double> %72, <2 x double> %75, <2 x i32> <i32 1, i32 2>
%78 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %19, <2 x double> %76, <2 x double> %72)
%79 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %21, <2 x double> %77, <2 x double> %75)
%80 = bitcast double* %70 to <2 x double>*
store <2 x double> %78, <2 x double>* %80
%81 = bitcast double* %73 to <2 x double>*
store <2 x double> %79, <2 x double>* %81
br label %82
82: ; preds = %64, %69
%83 = phi <2 x double> [ %65, %64 ], [ %75, %69 ]
%84 = extractelement <2 x double> %83, i64 1
br label %85
85: ; preds = %82, %5
%86 = phi i64 [ %16, %82 ], [ 2, %5 ]
%87 = phi double [ %84, %82 ], [ %7, %5 ]
br label %88
88: ; preds = %85, %88
%89 = phi i64 [ %94, %88 ], [ %86, %85 ]
%90 = phi double [ %92, %88 ], [ %87, %85 ]
%91 = getelementptr inbounds double, double* %0, i64 %89
%92 = load double, double* %91
%93 = tail call double @llvm.fmuladd.f64(double %2, double %90, double %92)
store double %93, double* %91
%94 = add nuw nsw i64 %89, 1
%95 = icmp eq i64 %94, %11
br i1 %95, label %96, label %88
96: ; preds = %88, %3
ret void
}