#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*/ NPdunpack_row(size_t, size_t, double *, double *);
int /*<<< orphan*/ NPdunpack_tril(int, double *, double *, int);
int /*<<< orphan*/ free(double *);
double *malloc(int);
int /*<<< orphan*/ memcpy(double *, double *, int);
void AO2MOrestore_nr8to1(double *eri8, double *eri1, int norb) {
  size_t npair = norb * (norb + 1) / 2;
  size_t i, j, ij;
  size_t d2 = norb * norb;
  size_t d3 = norb * norb * norb;
  double *buf = malloc(sizeof(double) * npair);
  for (ij = 0, i = 0; i < norb; i++) {
    for (j = 0; j < i + 1; j++, ij++) {
      NPdunpack_row(npair, ij, eri8, buf);
      NPdunpack_tril(norb, buf, eri1 + i * d3 + j * d2, 1);
      if (i > j) {
        memcpy(eri1 + j * d3 + i * d2, eri1 + i * d3 + j * d2,
               sizeof(double) * norb * norb);
      }
    }
  }
  free(buf);
}


x86 assembly (clang, -O0):
<AO2MOrestore_nr8to1>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cltq
mov    %rax,-0x20(%rbp)
mov    -0x14(%rbp),%eax
imul   -0x14(%rbp),%eax
cltq
mov    %rax,-0x40(%rbp)
mov    -0x14(%rbp),%eax
imul   -0x14(%rbp),%eax
imul   -0x14(%rbp),%eax
cltq
mov    %rax,-0x48(%rbp)
mov    -0x20(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x50(%rbp)
movq   $0x0,-0x38(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x14(%rbp),%rcx
cmp    %rcx,%rax
jae    12db <AO2MOrestore_nr8to1+0x18b>
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
mov    -0x28(%rbp),%rcx
add    $0x1,%rcx
cmp    %rcx,%rax
jae    12c5 <AO2MOrestore_nr8to1+0x175>
mov    -0x20(%rbp),%rdi
mov    -0x38(%rbp),%rsi
mov    -0x8(%rbp),%rdx
mov    -0x50(%rbp),%rcx
call   1040 <NPdunpack_row@plt>
mov    -0x14(%rbp),%edi
mov    -0x50(%rbp),%rsi
mov    -0x10(%rbp),%rdx
mov    -0x28(%rbp),%rax
imul   -0x48(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdx
mov    -0x30(%rbp),%rax
imul   -0x40(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdx
mov    $0x1,%ecx
call   1050 <NPdunpack_tril@plt>
mov    -0x28(%rbp),%rax
cmp    -0x30(%rbp),%rax
jbe    12a3 <AO2MOrestore_nr8to1+0x153>
mov    -0x10(%rbp),%rdi
mov    -0x30(%rbp),%rax
imul   -0x48(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x28(%rbp),%rax
imul   -0x40(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x10(%rbp),%rsi
mov    -0x28(%rbp),%rax
imul   -0x48(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
mov    -0x30(%rbp),%rax
imul   -0x40(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
movslq -0x14(%rbp),%rax
shl    $0x3,%rax
movslq -0x14(%rbp),%rcx
imul   %rcx,%rax
mov    %eax,%edx
call   1060 <memcpy@plt>
jmp    12a8 <AO2MOrestore_nr8to1+0x158>
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11d7 <AO2MOrestore_nr8to1+0x87>
jmp    12ca <AO2MOrestore_nr8to1+0x17a>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    11be <AO2MOrestore_nr8to1+0x6e>
mov    -0x50(%rbp),%rdi
call   1030 <free@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<AO2MOrestore_nr8to1>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %edx,%r14d
mov    %rsi,%rbx
mov    %rdi,0x50(%rsp)
lea    0x1(%r14),%eax
imul   %edx,%eax
mov    %eax,%ebp
shr    $0x1f,%ebp
add    %eax,%ebp
sar    %ebp
lea    0x0(,%rbp,8),%edi
call   1070 <malloc@plt>
mov    %rbx,%rdx
mov    %rax,%rbx
mov    %r14,0x10(%rsp)
test   %r14d,%r14d
je     1286 <AO2MOrestore_nr8to1+0x136>
movslq %ebp,%rax
mov    %rax,0x48(%rsp)
mov    0x10(%rsp),%rcx
mov    %ecx,%eax
imul   %eax,%eax
mov    %eax,%esi
lea    0x0(,%rax,8),%rdi
mov    %rdi,0x28(%rsp)
imul   %ecx,%eax
cltq
shl    $0x3,%esi
mov    %esi,0xc(%rsp)
movslq %ecx,%rcx
mov    %rcx,0x30(%rsp)
shl    $0x3,%rax
mov    %rax,0x20(%rsp)
mov    $0x1,%eax
mov    %rax,0x18(%rsp)
xor    %r15d,%r15d
mov    %rdx,%rbp
xor    %r12d,%r12d
jmp    1218 <AO2MOrestore_nr8to1+0xc8>
nopl   0x0(%rax)
add    $0x1,%r15
addq   $0x1,0x18(%rsp)
mov    0x38(%rsp),%rbp
add    0x20(%rsp),%rbp
mov    0x40(%rsp),%rdx
add    0x28(%rsp),%rdx
add    %r13,%r12
cmp    0x30(%rsp),%r15
je     1286 <AO2MOrestore_nr8to1+0x136>
mov    %rdx,0x40(%rsp)
mov    %rdx,%r14
mov    %rbp,0x38(%rsp)
xor    %r13d,%r13d
jmp    1245 <AO2MOrestore_nr8to1+0xf5>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r13
add    0x28(%rsp),%rbp
add    0x20(%rsp),%r14
cmp    %r13,0x18(%rsp)
je     11f0 <AO2MOrestore_nr8to1+0xa0>
lea    (%r12,%r13,1),%rsi
mov    0x48(%rsp),%rdi
mov    0x50(%rsp),%rdx
mov    %rbx,%rcx
call   1040 <NPdunpack_row@plt>
mov    0x10(%rsp),%rdi
mov    %rbx,%rsi
mov    %rbp,%rdx
mov    $0x1,%ecx
call   1050 <NPdunpack_tril@plt>
cmp    %r13,%r15
jbe    1230 <AO2MOrestore_nr8to1+0xe0>
mov    %r14,%rdi
mov    %rbp,%rsi
mov    0xc(%rsp),%edx
call   1060 <memcpy@plt>
jmp    1230 <AO2MOrestore_nr8to1+0xe0>
mov    %rbx,%rdi
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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