#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 sqrt(double);
double getPeriod(const double *m, double **k, long n) {
  if (m == NULL || n <= 0)
    return 0;
  double max = m[n - 1] / k[n - 1][n - 1];
  for (n -= 2; n >= 0; n--) {
    if (m[n] / k[n][n] > max) {
      max = m[n] / k[n][n];
    }
  }
  return 2 * 3.14159265358979323846 * sqrt(max);
}


x86 assembly (clang, -O0):
<getPeriod>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
cmpq   $0x0,-0x10(%rbp)
je     113a <getPeriod+0x2a>
cmpq   $0x0,-0x20(%rbp)
jg     1147 <getPeriod+0x37>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    121a <getPeriod+0x10a>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    $0x1,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%rcx
sub    $0x1,%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x20(%rbp),%rax
sub    $0x2,%rax
mov    %rax,-0x20(%rbp)
cmpq   $0x0,-0x20(%rbp)
jl     11fb <getPeriod+0xeb>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    11e5 <getPeriod+0xd5>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    11ea <getPeriod+0xda>
mov    -0x20(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x20(%rbp)
jmp    1186 <getPeriod+0x76>
movsd  -0x28(%rbp),%xmm0
call   1030 <sqrt@plt>
cvtsi2sd %eax,%xmm1
movsd  0xdef(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<getPeriod>:
xorpd  %xmm0,%xmm0
test   %rdi,%rdi
je     11cd <getPeriod+0xbd>
test   %rdx,%rdx
jle    11cd <getPeriod+0xbd>
push   %rax
movsd  -0x8(%rdi,%rdx,8),%xmm1
mov    -0x8(%rsi,%rdx,8),%rax
divsd  -0x8(%rax,%rdx,8),%xmm1
cmp    $0x2,%rdx
jb     1155 <getPeriod+0x45>
lea    -0x1(%rdx),%rcx
lea    -0x2(%rdx),%rax
test   $0x1,%cl
jne    115b <getPeriod+0x4b>
mov    %rax,%rdx
test   %rax,%rax
jne    1180 <getPeriod+0x70>
jmp    11b5 <getPeriod+0xa5>
movapd %xmm1,%xmm0
jmp    11b5 <getPeriod+0xa5>
movsd  -0x10(%rdi,%rdx,8),%xmm0
mov    -0x10(%rsi,%rdx,8),%rcx
divsd  -0x10(%rcx,%rdx,8),%xmm0
maxsd  %xmm1,%xmm0
add    $0xfffffffffffffffd,%rdx
movapd %xmm0,%xmm1
test   %rax,%rax
je     11b5 <getPeriod+0xa5>
nopl   (%rax)
movsd  -0x8(%rdi,%rdx,8),%xmm0
movsd  (%rdi,%rdx,8),%xmm2
mov    -0x8(%rsi,%rdx,8),%rax
mov    (%rsi,%rdx,8),%rcx
divsd  (%rcx,%rdx,8),%xmm2
divsd  -0x8(%rax,%rdx,8),%xmm0
maxsd  %xmm1,%xmm2
maxsd  %xmm2,%xmm0
add    $0xfffffffffffffffe,%rdx
movapd %xmm0,%xmm1
cmp    $0xffffffffffffffff,%rdx
jne    1180 <getPeriod+0x70>
call   1030 <sqrt@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0xe37(%rip),%xmm0
add    $0x8,%rsp
ret

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

LLVM IR (-O2):
define double @getPeriod(double* %0, double** nocapture %1, i64 %2) {
%4 = icmp eq double* %0, null
%5 = icmp slt i64 %2, 1
%6 = or i1 %4, %5
br i1 %6, label %66, label %7
7: ; preds = %3
%8 = add nsw i64 %2, -1
%9 = getelementptr inbounds double, double* %0, i64 %8
%10 = load double, double* %9
%11 = getelementptr inbounds double*, double** %1, i64 %8
%12 = load double*, double** %11
%13 = getelementptr inbounds double, double* %12, i64 %8
%14 = load double, double* %13
%15 = fdiv double %10, %14
%16 = icmp ugt i64 %2, 1
br i1 %16, label %17, label %61
17: ; preds = %7
%18 = add i64 %2, -2
%19 = and i64 %8, 1
%20 = icmp eq i64 %19, 0
br i1 %20, label %32, label %21
21: ; preds = %17
%22 = getelementptr inbounds double, double* %0, i64 %18
%23 = load double, double* %22
%24 = getelementptr inbounds double*, double** %1, i64 %18
%25 = load double*, double** %24
%26 = getelementptr inbounds double, double* %25, i64 %18
%27 = load double, double* %26
%28 = fdiv double %23, %27
%29 = fcmp ogt double %28, %15
%30 = select i1 %29, double %28, double %15
%31 = add i64 %2, -3
br label %32
32: ; preds = %21, %17
%33 = phi double [ undef, %17 ], [ %30, %21 ]
%34 = phi double [ %15, %17 ], [ %30, %21 ]
%35 = phi i64 [ %18, %17 ], [ %31, %21 ]
%36 = icmp eq i64 %18, 0
br i1 %36, label %61, label %37
37: ; preds = %32, %37
%38 = phi double [ %58, %37 ], [ %34, %32 ]
%39 = phi i64 [ %59, %37 ], [ %35, %32 ]
%40 = getelementptr inbounds double, double* %0, i64 %39
%41 = load double, double* %40
%42 = getelementptr inbounds double*, double** %1, i64 %39
%43 = load double*, double** %42
%44 = getelementptr inbounds double, double* %43, i64 %39
%45 = load double, double* %44
%46 = fdiv double %41, %45
%47 = fcmp ogt double %46, %38
%48 = select i1 %47, double %46, double %38
%49 = add nsw i64 %39, -1
%50 = getelementptr inbounds double, double* %0, i64 %49
%51 = load double, double* %50
%52 = getelementptr inbounds double*, double** %1, i64 %49
%53 = load double*, double** %52
%54 = getelementptr inbounds double, double* %53, i64 %49
%55 = load double, double* %54
%56 = fdiv double %51, %55
%57 = fcmp ogt double %56, %48
%58 = select i1 %57, double %56, double %48
%59 = add nsw i64 %39, -2
%60 = icmp eq i64 %49, 0
br i1 %60, label %61, label %37
61: ; preds = %32, %37, %7
%62 = phi double [ %15, %7 ], [ %33, %32 ], [ %58, %37 ]
%63 = tail call i32 @sqrt(double %62)
%64 = sitofp i32 %63 to double
%65 = fmul double %64, 0x401921FB54442D18
br label %66
66: ; preds = %3, %61
%67 = phi double [ %65, %61 ], [ 0.000000e+00, %3 ]
ret double %67
}