void permuteInv(int nPoint, int nFact, int fact[], int remain[], double xRe[],
                double xIm[], double yRe[], double yIm[]) {
  int i, j, k;
  int count[20];
  for (i = 1; i <= nFact; i++)
    count[i] = 0;
  k = 0;
  for (i = 0; i <= nPoint - 2; i++) {
    yRe[i] = xRe[k];
    yIm[i] = -xIm[k];
    j = 1;
    k = k + remain[j];
    count[1] = count[1] + 1;
    while (count[j] >= fact[j]) {
      count[j] = 0;
      k = k - remain[j - 1] + remain[j + 1];
      j = j + 1;
      count[j] = count[j] + 1;
    }
  }
  yRe[nPoint - 1] = xRe[nPoint - 1];
  yIm[nPoint - 1] = -xIm[nPoint - 1];
}


x86 assembly (clang, -O0):
<permuteInv>:
push   %rbp
mov    %rsp,%rbp
sub    $0x10,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jg     1156 <permuteInv+0x56>
movslq -0x2c(%rbp),%rax
movl   $0x0,-0x90(%rbp,%rax,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    112d <permuteInv+0x2d>
movl   $0x0,-0x34(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jg     1272 <permuteInv+0x172>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x30(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x18(%rbp),%rcx
movslq -0x30(%rbp),%rdx
add    (%rcx,%rdx,4),%eax
mov    %eax,-0x34(%rbp)
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
movslq -0x30(%rbp),%rax
mov    -0x90(%rbp,%rax,4),%eax
mov    -0x10(%rbp),%rcx
movslq -0x30(%rbp),%rdx
cmp    (%rcx,%rdx,4),%eax
jl     125f <permuteInv+0x15f>
movslq -0x30(%rbp),%rax
movl   $0x0,-0x90(%rbp,%rax,4)
mov    -0x34(%rbp),%eax
mov    -0x18(%rbp),%rcx
mov    -0x30(%rbp),%edx
sub    $0x1,%edx
movslq %edx,%rdx
sub    (%rcx,%rdx,4),%eax
mov    -0x18(%rbp),%rcx
mov    -0x30(%rbp),%edx
add    $0x1,%edx
movslq %edx,%rdx
add    (%rcx,%rdx,4),%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
movslq -0x30(%rbp),%rax
mov    -0x90(%rbp,%rax,4),%ecx
add    $0x1,%ecx
movslq -0x30(%rbp),%rax
mov    %ecx,-0x90(%rbp,%rax,4)
jmp    11e7 <permuteInv+0xe7>
jmp    1264 <permuteInv+0x164>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1164 <permuteInv+0x64>
mov    -0x20(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
add    $0x10,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<permuteInv>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x50,%rsp
mov    %r9,%r14
mov    %r8,%r15
mov    %rcx,%rbx
mov    %rdx,%r13
mov    %edi,%r12d
test   %esi,%esi
jle    1142 <permuteInv+0x32>
lea    0x4(%rsp),%rdi
mov    %esi,%edx
shl    $0x2,%rdx
xor    %esi,%esi
call   1030 <memset@plt>
mov    0x88(%rsp),%r8
mov    0x80(%rsp),%r9
cmp    $0x1,%r12d
jg     118b <permuteInv+0x7b>
add    $0xffffffff,%r12d
movslq %r12d,%rax
movsd  (%r15,%rax,8),%xmm0
movsd  %xmm0,(%r9,%rax,8)
movsd  (%r14,%rax,8),%xmm0
xorps  0xe88(%rip),%xmm0
movlps %xmm0,(%r8,%rax,8)
add    $0x50,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret
mov    0x4(%rbx),%r10d
mov    0x4(%r13),%r11d
add    $0xffffffff,%r12d
xor    %edi,%edi
movaps 0xe60(%rip),%xmm0
xor    %eax,%eax
jmp    11b9 <permuteInv+0xa9>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%rdi
cmp    %r12,%rdi
je     115c <permuteInv+0x4c>
cltq
movsd  (%r15,%rax,8),%xmm1
movsd  %xmm1,(%r9,%rdi,8)
movsd  (%r14,%rax,8),%xmm1
xorps  %xmm0,%xmm1
movlps %xmm1,(%r8,%rdi,8)
add    %r10d,%eax
mov    0x4(%rsp),%ecx
add    $0x1,%ecx
mov    %ecx,0x4(%rsp)
cmp    %r11d,%ecx
jl     11b0 <permuteInv+0xa0>
mov    $0x8,%ecx
nopl   (%rax)
movl   $0x0,-0x4(%rsp,%rcx,1)
sub    -0x8(%rbx,%rcx,1),%eax
add    (%rbx,%rcx,1),%eax
mov    (%rsp,%rcx,1),%edx
add    $0x1,%edx
mov    %edx,(%rsp,%rcx,1)
lea    0x4(%rcx),%rsi
cmp    0x0(%r13,%rcx,1),%edx
mov    %rsi,%rcx
jge    11f0 <permuteInv+0xe0>
jmp    11b0 <permuteInv+0xa0>

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

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