#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*/ floor(double);
scalar_t__ malloc(int);
size_t *Mat_CalcSudirripts2(int rank, size_t *dims, size_t index) {
  int i;
  size_t *subs;
  double l;
  subs = (size_t *)malloc(rank * sizeof(size_t));
  l = (double)index;
  for (i = rank; i--;) {
    int j;
    size_t k = 1;
    for (j = i; j--;)
      k *= dims[j];
    subs[i] = (size_t)floor(l / (double)k);
    l -= subs[i] * k;
    subs[i]++;
  }
  return subs;
}


x86 assembly (clang, -O0):
<Mat_CalcSudirripts2>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x4(%rbp),%edi
shl    $0x3,%edi
call   1030 <malloc@plt>
mov    %rax,-0x28(%rbp)
movaps 0xeb7(%rip),%xmm1
movq   -0x18(%rbp),%xmm0
punpckldq %xmm1,%xmm0
movapd 0xeb6(%rip),%xmm1
subpd  %xmm1,%xmm0
movaps %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x4(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x1c(%rbp)
cmp    $0x0,%eax
je     126b <Mat_CalcSudirripts2+0x14b>
movq   $0x1,-0x40(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x34(%rbp)
cmp    $0x0,%eax
je     11c4 <Mat_CalcSudirripts2+0xa4>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
imul   -0x40(%rbp),%rax
mov    %rax,-0x40(%rbp)
jmp    1196 <Mat_CalcSudirripts2+0x76>
movsd  -0x30(%rbp),%xmm0
movq   -0x40(%rbp),%xmm1
movaps 0xe2b(%rip),%xmm2
movaps %xmm2,-0x60(%rbp)
punpckldq %xmm2,%xmm1
movapd 0xe2b(%rip),%xmm2
movaps %xmm2,-0x50(%rbp)
subpd  %xmm2,%xmm1
movaps %xmm1,%xmm2
unpckhpd %xmm1,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
call   1040 <floor@plt>
movaps -0x60(%rbp),%xmm2
movaps -0x50(%rbp),%xmm0
movslq %eax,%rdx
mov    -0x28(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x40(%rbp),%rcx
imul   %rcx,%rax
movq   %rax,%xmm1
punpckldq %xmm2,%xmm1
subpd  %xmm0,%xmm1
movaps %xmm1,%xmm0
unpckhpd %xmm1,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x30(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
add    $0x1,%rdx
mov    %rdx,(%rax,%rcx,8)
jmp    1174 <Mat_CalcSudirripts2+0x54>
mov    -0x28(%rbp),%rax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Mat_CalcSudirripts2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rdx,%r12
mov    %rsi,%rbx
mov    %edi,%r14d
lea    0x0(,%r14,8),%edi
call   1030 <malloc@plt>
mov    %rax,0x8(%rsp)
test   %r14d,%r14d
je     12d3 <Mat_CalcSudirripts2+0x1b3>
movq   %r12,%xmm0
punpckldq 0xea1(%rip),%xmm0
subpd  0xea9(%rip),%xmm0
movapd %xmm0,%xmm2
unpckhpd %xmm0,%xmm2
addsd  %xmm0,%xmm2
lea    -0x2(%r14),%r12d
movslq %r14d,%r13
add    $0x7,%r14b
xor    %ebp,%ebp
mov    %r12d,0x4(%rsp)
jmp    1224 <Mat_CalcSudirripts2+0x104>
nopw   0x0(%rax,%rax,1)
mov    $0x1,%r15d
movq   %r15,%xmm0
movdqa 0xe5d(%rip),%xmm1
punpckldq %xmm1,%xmm0
movapd 0xe61(%rip),%xmm1
subpd  %xmm1,%xmm0
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
addsd  %xmm0,%xmm1
movapd %xmm2,%xmm0
divsd  %xmm1,%xmm0
movapd %xmm2,0x10(%rsp)
call   1040 <floor@plt>
movapd 0x10(%rsp),%xmm2
cltq
imul   %rax,%r15
movq   %r15,%xmm0
punpckldq 0xe15(%rip),%xmm0
subpd  0xe1d(%rip),%xmm0
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
addsd  %xmm0,%xmm1
subsd  %xmm1,%xmm2
add    $0x1,%rax
mov    0x8(%rsp),%rcx
mov    %rax,(%rcx,%r13,8)
add    $0xffffffff,%r12d
add    $0x1,%ebp
add    $0x7,%r14b
test   %r13d,%r13d
je     12d3 <Mat_CalcSudirripts2+0x1b3>
add    $0xffffffffffffffff,%r13
test   %r13d,%r13d
je     1190 <Mat_CalcSudirripts2+0x70>
movslq %r12d,%rax
mov    0x4(%rsp),%ecx
sub    %ebp,%ecx
movslq %ecx,%rdx
lea    0x1(%rdx),%rcx
mov    $0x1,%r15d
cmp    $0x7,%edx
jb     1299 <Mat_CalcSudirripts2+0x179>
lea    0x1(%rax),%rdx
and    $0xfffffffffffffff8,%rdx
mov    $0x1,%r15d
nopw   0x0(%rax,%rax,1)
imul   (%rbx,%rax,8),%r15
imul   -0x8(%rbx,%rax,8),%r15
imul   -0x10(%rbx,%rax,8),%r15
imul   -0x18(%rbx,%rax,8),%r15
imul   -0x20(%rbx,%rax,8),%r15
imul   -0x28(%rbx,%rax,8),%r15
imul   -0x30(%rbx,%rax,8),%r15
imul   -0x38(%rbx,%rax,8),%r15
add    $0xfffffffffffffff8,%rax
add    $0xfffffffffffffff8,%rdx
jne    1260 <Mat_CalcSudirripts2+0x140>
test   $0x7,%cl
je     1196 <Mat_CalcSudirripts2+0x76>
movzbl %r14b,%ecx
and    $0x7,%ecx
neg    %rcx
lea    (%rbx,%rax,8),%rax
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
imul   (%rax,%rdx,8),%r15
add    $0xffffffffffffffff,%rdx
cmp    %rdx,%rcx
jne    12c0 <Mat_CalcSudirripts2+0x1a0>
jmp    1196 <Mat_CalcSudirripts2+0x76>
mov    0x8(%rsp),%rax
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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