void sortsig(double *sigma, int ma, int ia[], int mfit) {
  int i, k;
  double swap;
  for (i = mfit + 1; i <= ma; i++)
    sigma[i] = 0.;
  k = mfit;
  for (i = ma; i >= 1; i--) {
    if (ia[i]) {
      {
        swap = (sigma[i]);
        (sigma[i]) = (sigma[k]);
        (sigma[k]) = swap;
      };
      k--;
    };
  };
}


x86 assembly (clang, -O0):
<sortsig>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1145 <sortsig+0x45>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    111b <sortsig+0x1b>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x1,-0x20(%rbp)
jl     11c7 <sortsig+0xc7>
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     11b4 <sortsig+0xb4>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    11b9 <sortsig+0xb9>
mov    -0x20(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
jmp    1151 <sortsig+0x51>
pop    %rbp
ret

x86 assembly (clang, -O2):
<sortsig>:
push   %rbp
push   %r15
push   %r14
push   %rbx
push   %rax
mov    %ecx,%ebx
mov    %rdx,%r14
mov    %esi,%ebp
mov    %rdi,%r15
cmp    %esi,%ecx
jge    1145 <sortsig+0x35>
movslq %ebx,%rax
lea    (%r15,%rax,8),%rdi
add    $0x8,%rdi
mov    %ebx,%eax
not    %eax
add    %ebp,%eax
lea    0x8(,%rax,8),%rdx
xor    %esi,%esi
call   1030 <memset@plt>
test   %ebp,%ebp
jle    117d <sortsig+0x6d>
mov    %ebp,%eax
test   $0x1,%al
je     1178 <sortsig+0x68>
cmpl   $0x0,(%r14,%rax,4)
je     1174 <sortsig+0x64>
movsd  (%r15,%rax,8),%xmm0
movslq %ebx,%rcx
movsd  (%r15,%rcx,8),%xmm1
movsd  %xmm1,(%r15,%rax,8)
movsd  %xmm0,(%r15,%rcx,8)
add    $0xffffffff,%ebx
add    $0xffffffffffffffff,%rax
cmp    $0x1,%ebp
jne    1188 <sortsig+0x78>
add    $0x8,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret
add    $0x2,%rax
jmp    119a <sortsig+0x8a>
xchg   %ax,%ax
add    $0xfffffffffffffffe,%rax
cmp    $0x2,%rax
jle    117d <sortsig+0x6d>
cmpl   $0x0,-0x8(%r14,%rax,4)
je     11c2 <sortsig+0xb2>
movsd  -0x10(%r15,%rax,8),%xmm0
movslq %ebx,%rbx
movsd  (%r15,%rbx,8),%xmm1
movsd  %xmm1,-0x10(%r15,%rax,8)
movsd  %xmm0,(%r15,%rbx,8)
add    $0xffffffff,%ebx
cmpl   $0x0,-0xc(%r14,%rax,4)
je     1190 <sortsig+0x80>
movsd  -0x18(%r15,%rax,8),%xmm0
movslq %ebx,%rbx
movsd  (%r15,%rbx,8),%xmm1
movsd  %xmm1,-0x18(%r15,%rax,8)
movsd  %xmm0,(%r15,%rbx,8)
add    $0xffffffff,%ebx
jmp    1190 <sortsig+0x80>

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

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