#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*/ NPdcopy(double *, double *, int);
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);
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) {
        NPdcopy(eri1 + j * d3 + i * d2, eri1 + i * d3 + j * d2, 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   1060 <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    12d0 <AO2MOrestore_nr8to1+0x180>
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
mov    -0x28(%rbp),%rcx
add    $0x1,%rcx
cmp    %rcx,%rax
jae    12ba <AO2MOrestore_nr8to1+0x16a>
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    1298 <AO2MOrestore_nr8to1+0x148>
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
mov    -0x14(%rbp),%edx
imul   -0x14(%rbp),%edx
call   1070 <NPdcopy@plt>
jmp    129d <AO2MOrestore_nr8to1+0x14d>
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    12bf <AO2MOrestore_nr8to1+0x16f>
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   1060 <malloc@plt>
mov    %rbx,%rdx
mov    %rax,%rbx
mov    %r14,0x8(%rsp)
test   %r14d,%r14d
je     1287 <AO2MOrestore_nr8to1+0x137>
movslq %ebp,%rax
mov    %rax,0x48(%rsp)
mov    0x8(%rsp),%rcx
mov    %ecx,%esi
imul   %esi,%esi
mov    %esi,%eax
imul   %ecx,%eax
cltq
movslq %ecx,%rcx
mov    %rcx,0x28(%rsp)
shl    $0x3,%rax
mov    %rax,0x20(%rsp)
mov    %rsi,0x40(%rsp)
lea    0x0(,%rsi,8),%rax
mov    %rax,0x18(%rsp)
mov    $0x1,%eax
mov    %rax,0x10(%rsp)
xor    %r15d,%r15d
mov    %rdx,%rbp
xor    %r12d,%r12d
jmp    1218 <AO2MOrestore_nr8to1+0xc8>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r15
addq   $0x1,0x10(%rsp)
mov    0x30(%rsp),%rbp
add    0x20(%rsp),%rbp
mov    0x38(%rsp),%rdx
add    0x18(%rsp),%rdx
add    %r13,%r12
cmp    0x28(%rsp),%r15
je     1287 <AO2MOrestore_nr8to1+0x137>
mov    %rdx,0x38(%rsp)
mov    %rdx,%r14
mov    %rbp,0x30(%rsp)
xor    %r13d,%r13d
jmp    1245 <AO2MOrestore_nr8to1+0xf5>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r13
add    0x18(%rsp),%rbp
add    0x20(%rsp),%r14
cmp    %r13,0x10(%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    0x8(%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    0x40(%rsp),%rdx
call   1070 <NPdcopy@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 = %96, %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 %99
39: ; preds = %34
store i64 0, i64* %9
br label %40
40: ; preds = %90, %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 %95
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 %89
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 = load i32, i32* %6
%87 = mul nsw i32 %85, %86
%88 = call i32 @NPdcopy(double* %75, double* %84, i32 %87)
br label %89
89: ; preds = %66, %45
br label %90
90: ; preds = %89
%91 = load i64, i64* %9
%92 = add i64 %91, 1
store i64 %92, i64* %9
%93 = load i64, i64* %10
%94 = add i64 %93, 1
store i64 %94, i64* %10
br label %40
95: ; preds = %40
br label %96
96: ; preds = %95
%97 = load i64, i64* %8
%98 = add i64 %97, 1
store i64 %98, i64* %8
br label %34
99: ; preds = %34
%100 = load double*, double** %13
%101 = call i32 @free(double* %100)
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 %45, label %15
15: ; preds = %3
%16 = sext i32 %2 to i64
br label %21
17: ; preds = %41
%18 = add nuw i64 %23, 1
%19 = add i64 %22, 1
%20 = icmp eq i64 %18, %16
br i1 %20, label %45, label %21
21: ; preds = %15, %17
%22 = phi i64 [ 1, %15 ], [ %19, %17 ]
%23 = phi i64 [ 0, %15 ], [ %18, %17 ]
%24 = phi i64 [ 0, %15 ], [ %43, %17 ]
%25 = mul i64 %23, %11
%26 = getelementptr inbounds double, double* %1, i64 %25
%27 = mul i64 %23, %9
%28 = getelementptr inbounds double, double* %1, i64 %27
br label %29
29: ; preds = %21, %41
%30 = phi i64 [ 0, %21 ], [ %42, %41 ]
%31 = phi i64 [ %24, %21 ], [ %43, %41 ]
%32 = tail call i32 @NPdunpack_row(i64 %7, i64 %31, double* %0, double* %13)
%33 = mul i64 %30, %9
%34 = getelementptr inbounds double, double* %26, i64 %33
%35 = tail call i32 @NPdunpack_tril(i32 %2, double* %13, double* %34, i32 1)
%36 = icmp ugt i64 %23, %30
br i1 %36, label %37, label %41
37: ; preds = %29
%38 = mul i64 %30, %11
%39 = getelementptr inbounds double, double* %28, i64 %38
%40 = tail call i32 @NPdcopy(double* %39, double* %34, i32 %8)
br label %41
41: ; preds = %29, %37
%42 = add nuw i64 %30, 1
%43 = add i64 %31, 1
%44 = icmp eq i64 %42, %22
br i1 %44, label %17, label %29
45: ; preds = %17, %3
%46 = tail call i32 @free(double* %13)
ret void
}