void lsx_ordered_convolve(int n, void *not_used, double *a, const double *b) {
  int i;
  a[0] *= b[0];
  a[1] *= b[1];
  for (i = 2; i < n; i += 2) {
    double tmp = a[i];
    a[i] = b[i] * tmp - b[i + 1] * a[i + 1];
    a[i + 1] = b[i + 1] * tmp + b[i] * a[i + 1];
  }
  (void)not_used;
}


x86 assembly (clang, -O0):
<lsx_ordered_convolve>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulsd  0x8(%rax),%xmm0
movsd  %xmm0,0x8(%rax)
movl   $0x2,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1223 <lsx_ordered_convolve+0x123>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x30(%rbp),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x30(%rbp),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x24(%rbp)
jmp    1145 <lsx_ordered_convolve+0x45>
pop    %rbp
ret

x86 assembly (clang, -O2):
<lsx_ordered_convolve>:
movsd  (%rcx),%xmm0
mulsd  (%rdx),%xmm0
movsd  %xmm0,(%rdx)
movsd  0x8(%rcx),%xmm0
mulsd  0x8(%rdx),%xmm0
movsd  %xmm0,0x8(%rdx)
cmp    $0x3,%edi
jl     125d <lsx_ordered_convolve+0x15d>
mov    %edi,%r10d
cmp    $0x5,%r10
mov    $0x4,%r9d
cmovae %r10,%r9
mov    $0x4,%eax
add    $0xfffffffffffffffd,%r9
mov    $0x2,%esi
cmp    $0x2,%r9
jb     1210 <lsx_ordered_convolve+0x110>
lea    0x10(%rdx),%r8
cmp    $0x5,%r10
cmovae %r10,%rax
add    $0x1,%rax
and    $0xfffffffffffffffe,%rax
lea    (%rcx,%rax,8),%rdi
cmp    %rdi,%r8
jae    117b <lsx_ordered_convolve+0x7b>
lea    (%rdx,%rax,8),%rax
lea    0x10(%rcx),%rdi
cmp    %rax,%rdi
jb     1210 <lsx_ordered_convolve+0x110>
shr    %r9
add    $0x1,%r9
movabs $0x7ffffffffffffffe,%r8
ff ff 7f
and    %r9,%r8
lea    (%r9,%r9,1),%rsi
or     $0x2,%rsi
mov    $0x10,%edi
mov    %r8,%rax
nop
movupd (%rdx,%rdi,1),%xmm0
movupd 0x10(%rdx,%rdi,1),%xmm1
movapd %xmm0,%xmm2
unpcklpd %xmm1,%xmm2
unpckhpd %xmm1,%xmm0
movupd (%rcx,%rdi,1),%xmm1
movupd 0x10(%rcx,%rdi,1),%xmm3
movapd %xmm1,%xmm4
unpcklpd %xmm3,%xmm4
unpckhpd %xmm3,%xmm1
movapd %xmm0,%xmm3
mulpd  %xmm1,%xmm3
mulpd  %xmm4,%xmm0
mulpd  %xmm2,%xmm4
subpd  %xmm3,%xmm4
mulpd  %xmm2,%xmm1
addpd  %xmm0,%xmm1
movapd %xmm4,%xmm0
unpckhpd %xmm1,%xmm0
unpcklpd %xmm1,%xmm4
movupd %xmm4,(%rdx,%rdi,1)
movupd %xmm0,0x10(%rdx,%rdi,1)
add    $0x20,%rdi
add    $0xfffffffffffffffe,%rax
jne    11a0 <lsx_ordered_convolve+0xa0>
cmp    %r8,%r9
je     125d <lsx_ordered_convolve+0x15d>
lea    0x0(,%rsi,8),%rax
movsd  (%rdx,%rsi,8),%xmm0
or     $0x8,%rax
movsd  (%rdx,%rax,1),%xmm1
movsd  (%rcx,%rax,1),%xmm2
mulsd  %xmm1,%xmm2
movsd  (%rcx,%rsi,8),%xmm3
mulsd  %xmm0,%xmm3
subsd  %xmm2,%xmm3
movsd  %xmm3,(%rdx,%rsi,8)
mulsd  (%rcx,%rsi,8),%xmm1
mulsd  (%rcx,%rax,1),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rdx,%rax,1)
add    $0x2,%rsi
cmp    %r10,%rsi
jb     1210 <lsx_ordered_convolve+0x110>
ret

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

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