void dst(x, y, m, tab) double *x, *y, *tab;
int m;
{
  int n, n2, i;
  double *xp1, *xp2;
  double temp;
  n = 1 << m;
  n2 = n >> 1;
  xp1 = x + 1;
  for (i = 0; i < n2; i++) {
    *xp1 = -*xp1;
    xp1 += 2;
  }
  dct(x, y, m, tab);
  xp1 = x;
  xp2 = x + n - 1;
  for (i = 0; i < n2; i++) {
    temp = *xp1;
    *xp1 = *xp2;
    *xp2 = temp;
    xp1++;
    xp2--;
  }
}


x86 assembly (clang, -O0):
<dst>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x14(%rbp),%ecx
mov    $0x1,%eax
shl    %cl,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x8(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    119d <dst+0x8d>
mov    -0x38(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x38(%rbp),%rax
add    $0x10,%rax
mov    %rax,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1150 <dst+0x40>
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x14(%rbp),%edx
mov    -0x20(%rbp),%rcx
mov    $0x0,%al
call   1030 <dct@plt>
mov    -0x8(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    1235 <dst+0x125>
mov    -0x38(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x48(%rbp),%xmm0
mov    -0x40(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x38(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x38(%rbp)
mov    -0x40(%rbp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11d9 <dst+0xc9>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dst>:
push   %rbp
push   %r15
push   %r14
push   %rbx
push   %rax
mov    %rcx,%r8
mov    %rdi,%rbx
mov    $0x1,%r14d
mov    $0x1,%ebp
mov    %edx,%ecx
shl    %cl,%ebp
cmp    $0x1,%ebp
jle    128a <dst+0x17a>
mov    %ebp,%eax
sar    %eax
lea    0x8(%rbx),%rcx
cmp    $0x2,%eax
cmovge %eax,%r14d
lea    -0x1(%r14),%r15d
cmp    $0x3,%r15d
jb     1197 <dst+0x87>
mov    %r14d,%edi
and    $0x7ffffffc,%edi
movaps 0xea2(%rip),%xmm0
xchg   %ax,%ax
movsd  (%rcx),%xmm1
movsd  0x10(%rcx),%xmm2
movsd  0x20(%rcx),%xmm3
movsd  0x30(%rcx),%xmm4
xorps  %xmm0,%xmm1
movlps %xmm1,(%rcx)
xorps  %xmm0,%xmm2
movlps %xmm2,0x10(%rcx)
xorps  %xmm0,%xmm3
movlps %xmm3,0x20(%rcx)
xorps  %xmm0,%xmm4
movlps %xmm4,0x30(%rcx)
add    $0x40,%rcx
add    $0xfffffffc,%edi
jne    1160 <dst+0x50>
test   $0x3,%r14b
je     11c4 <dst+0xb4>
mov    %r14d,%edi
and    $0x3,%edi
shl    $0x4,%rdi
xor    %eax,%eax
movaps 0xe50(%rip),%xmm0
movsd  (%rcx,%rax,1),%xmm1
xorps  %xmm0,%xmm1
movlps %xmm1,(%rcx,%rax,1)
add    $0x10,%rax
cmp    %eax,%edi
jne    11b0 <dst+0xa0>
mov    %rbx,%rdi
mov    %r8,%rcx
xor    %eax,%eax
call   1030 <dct@plt>
cmp    $0x2,%ebp
jl     127f <dst+0x16f>
mov    %ebp,%eax
lea    (%rbx,%rax,8),%rax
cmp    $0x3,%r15d
jb     124b <dst+0x13b>
mov    %r14d,%ecx
and    $0x7ffffffc,%ecx
nop
movsd  (%rbx),%xmm0
movsd  -0x8(%rax),%xmm1
movsd  %xmm1,(%rbx)
movsd  %xmm0,-0x8(%rax)
movsd  0x8(%rbx),%xmm0
movsd  -0x10(%rax),%xmm1
movsd  %xmm1,0x8(%rbx)
movsd  %xmm0,-0x10(%rax)
movsd  0x10(%rbx),%xmm0
movsd  -0x18(%rax),%xmm1
movsd  %xmm1,0x10(%rbx)
movsd  %xmm0,-0x18(%rax)
movsd  0x18(%rbx),%xmm0
movsd  -0x20(%rax),%xmm1
movsd  %xmm1,0x18(%rbx)
movsd  %xmm0,-0x20(%rax)
add    $0xffffffffffffffe0,%rax
add    $0x20,%rbx
add    $0xfffffffc,%ecx
jne    11f0 <dst+0xe0>
test   $0x3,%r14b
je     127f <dst+0x16f>
add    $0xfffffffffffffff8,%rax
and    $0x3,%r14d
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movsd  (%rbx,%rcx,8),%xmm0
movsd  (%rax),%xmm1
movsd  %xmm1,(%rbx,%rcx,8)
movsd  %xmm0,(%rax)
add    $0xfffffffffffffff8,%rax
add    $0x1,%rcx
cmp    %ecx,%r14d
jne    1260 <dst+0x150>
add    $0x8,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %rbx,%rdi
mov    %r8,%rcx
xor    %eax,%eax
add    $0x8,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <dct@plt>

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

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