#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ free(double *);
double *malloc(int);
int /*<<< orphan*/ memcpy(double *, double *, int);
int reorder_dmat_c(double *dmat, int *permutation, int num_dof) {
  int result, size, i, j, ip, jp, k, kp;
  double *tmp;
  result = 0;
  size = ((num_dof + 1) * num_dof) / 2;
  tmp = malloc(size * sizeof(double));
  if (tmp == NULL) {
    result = -1;
    goto EXIT;
  }
  k = 0;
  for (i = 0; i < num_dof; i++) {
    ip = permutation[i];
    for (j = 0; j <= i; j++) {
      jp = permutation[j];
      if (ip >= jp) {
        kp = (ip * (ip + 1)) / 2 + jp;
      } else {
        kp = (jp * (jp + 1)) / 2 + ip;
      }
      tmp[k] = dmat[kp];
      k++;
    }
  }
  memcpy(dmat, tmp, size * sizeof(double));
EXIT:
  free(tmp);
  return result;
}


x86 assembly (clang, -O0):
<reorder_dmat_c>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
imul   -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x1c(%rbp)
movslq -0x1c(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x40(%rbp)
cmpq   $0x0,-0x40(%rbp)
jne    1189 <reorder_dmat_c+0x59>
movl   $0xffffffff,-0x18(%rbp)
jmp    1272 <reorder_dmat_c+0x142>
movl   $0x0,-0x30(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    125b <reorder_dmat_c+0x12b>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x28(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     1248 <reorder_dmat_c+0x118>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jl     11fd <reorder_dmat_c+0xcd>
mov    -0x28(%rbp),%eax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
jmp    1217 <reorder_dmat_c+0xe7>
mov    -0x2c(%rbp),%eax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    -0x28(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11b8 <reorder_dmat_c+0x88>
jmp    124d <reorder_dmat_c+0x11d>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1197 <reorder_dmat_c+0x67>
mov    -0x8(%rbp),%rdi
mov    -0x40(%rbp),%rsi
movslq -0x1c(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edx
call   1040 <memcpy@plt>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
mov    -0x18(%rbp),%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<reorder_dmat_c>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %edx,%ebx
mov    %rsi,%r12
mov    %rdi,%r14
lea    0x1(%rbx),%eax
imul   %edx,%eax
mov    %eax,%r15d
shr    $0x1f,%r15d
add    %eax,%r15d
shl    $0x2,%r15d
and    $0xfffffff8,%r15d
mov    %r15d,%edi
call   1050 <malloc@plt>
mov    %rax,%r13
test   %rax,%rax
je     1216 <reorder_dmat_c+0xe6>
test   %ebx,%ebx
jle    1204 <reorder_dmat_c+0xd4>
mov    %ebx,%r8d
mov    $0x1,%ecx
xor    %r9d,%r9d
xor    %r10d,%r10d
jmp    11a0 <reorder_dmat_c+0x70>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r9
add    $0x1,%rcx
add    %rbp,%r10
cmp    %r8,%r9
je     1204 <reorder_dmat_c+0xd4>
mov    (%r12,%r9,4),%edi
lea    0x1(%rdi),%edx
imul   %edi,%edx
mov    %edx,%eax
shr    $0x1f,%eax
add    %edx,%eax
sar    %eax
movslq %r10d,%r10
lea    0x0(,%r10,8),%rdx
add    %r13,%rdx
xor    %ebp,%ebp
jmp    11f8 <reorder_dmat_c+0xc8>
cs nopw 0x0(%rax,%rax,1)
nop
lea    0x1(%rsi),%ebx
imul   %esi,%ebx
mov    %ebx,%esi
shr    $0x1f,%esi
add    %ebx,%esi
sar    %esi
add    %edi,%esi
movslq %esi,%rsi
movsd  (%r14,%rsi,8),%xmm0
movsd  %xmm0,(%rdx,%rbp,8)
add    $0x1,%rbp
cmp    %rbp,%rcx
je     1190 <reorder_dmat_c+0x60>
mov    (%r12,%rbp,4),%esi
cmp    %esi,%edi
jl     11d0 <reorder_dmat_c+0xa0>
add    %eax,%esi
jmp    11e1 <reorder_dmat_c+0xb1>
mov    %r14,%rdi
mov    %r13,%rsi
mov    %r15d,%edx
call   1040 <memcpy@plt>
xor    %ebx,%ebx
jmp    121b <reorder_dmat_c+0xeb>
mov    $0xffffffff,%ebx
mov    %r13,%rdi
call   1030 <free@plt>
mov    %ebx,%eax
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @reorder_dmat_c(double* %0, i32* nocapture %1, i32 %2) {
%4 = add nsw i32 %2, 1
%5 = mul nsw i32 %4, %2
%6 = sdiv i32 %5, 2
%7 = shl i32 %6, 3
%8 = tail call double* @malloc(i32 %7)
%9 = icmp eq double* %8, null
br i1 %9, label %53, label %10
10: ; preds = %3
%11 = icmp sgt i32 %2, 0
br i1 %11, label %12, label %51
12: ; preds = %10
%13 = zext i32 %2 to i64
br label %14
14: ; preds = %12, %47
%15 = phi i64 [ 0, %12 ], [ %48, %47 ]
%16 = phi i64 [ 1, %12 ], [ %49, %47 ]
%17 = phi i64 [ 0, %12 ], [ %44, %47 ]
%18 = getelementptr inbounds i32, i32* %1, i64 %15
%19 = load i32, i32* %18
%20 = add nsw i32 %19, 1
%21 = mul nsw i32 %20, %19
%22 = sdiv i32 %21, 2
%23 = shl i64 %17, 32
%24 = ashr exact i64 %23, 32
br label %25
25: ; preds = %14, %38
%26 = phi i64 [ 0, %14 ], [ %45, %38 ]
%27 = phi i64 [ %24, %14 ], [ %44, %38 ]
%28 = getelementptr inbounds i32, i32* %1, i64 %26
%29 = load i32, i32* %28
%30 = icmp slt i32 %19, %29
br i1 %30, label %33, label %31
31: ; preds = %25
%32 = add nsw i32 %29, %22
br label %38
33: ; preds = %25
%34 = add nsw i32 %29, 1
%35 = mul nsw i32 %34, %29
%36 = sdiv i32 %35, 2
%37 = add nsw i32 %36, %19
br label %38
38: ; preds = %33, %31
%39 = phi i32 [ %32, %31 ], [ %37, %33 ]
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %0, i64 %40
%42 = load double, double* %41
%43 = getelementptr inbounds double, double* %8, i64 %27
store double %42, double* %43
%44 = add nsw i64 %27, 1
%45 = add nuw nsw i64 %26, 1
%46 = icmp eq i64 %45, %16
br i1 %46, label %47, label %25
47: ; preds = %38
%48 = add nuw nsw i64 %15, 1
%49 = add nuw nsw i64 %16, 1
%50 = icmp eq i64 %48, %13
br i1 %50, label %51, label %14
51: ; preds = %47, %10
%52 = tail call i32 @memcpy(double* %0, double* %8, i32 %7)
br label %53
53: ; preds = %3, %51
%54 = phi i32 [ 0, %51 ], [ -1, %3 ]
%55 = tail call i32 @free(double* %8)
ret i32 %54
}