#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 */
double legendre(int, double);
int /*<<< orphan*/ lobatto_nodes(double *, int);
void lobatto_quad(double *restrict z, double *restrict w, int n) {
  int i, j;
  lobatto_nodes(z, n);
  for (i = 0; i <= (n - 1) / 2; ++i) {
    double d = legendre(n - 1, z[i]);
    w[i] = 2 / ((n - 1) * n * d * d);
  }
  for (j = (n + 1) / 2, i = n / 2 - 1; j < n; ++j, --i)
    w[j] = w[i];
}


x86 assembly (clang, -O0):
<lobatto_quad>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x14(%rbp),%esi
call   1030 <lobatto_nodes@plt>
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x2c(%rbp),%eax
cmp    %ecx,%eax
jg     11c3 <lobatto_quad+0xa3>
mov    -0x14(%rbp),%edi
sub    $0x1,%edi
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1040 <legendre@plt>
movsd  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
imul   -0x14(%rbp),%eax
cvtsi2sd %eax,%xmm1
mulsd  -0x28(%rbp),%xmm1
mulsd  -0x28(%rbp),%xmm1
movsd  0xe5c(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1146 <lobatto_quad+0x26>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x1c(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
sub    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1222 <lobatto_quad+0x102>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x18(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x18(%rbp)
jmp    11e5 <lobatto_quad+0xc5>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<lobatto_quad>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %edx,%ebx
mov    %rsi,%r14
mov    %rdi,%r15
mov    %edx,%esi
call   1030 <lobatto_nodes@plt>
test   %ebx,%ebx
js     11a6 <lobatto_quad+0x86>
lea    -0x1(%rbx),%r12d
mov    %r12d,%r13d
shr    $0x1f,%r13d
add    %r12d,%r13d
mov    %r12d,%eax
imul   %ebx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%rsp)
sar    %r13d
add    $0x1,%r13d
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%r15,%rbp,8),%xmm0
mov    %r12d,%edi
call   1040 <legendre@plt>
movapd %xmm0,%xmm1
mulsd  (%rsp),%xmm1
mulsd  %xmm0,%xmm1
movsd  0xe6d(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%r14,%rbp,8)
add    $0x1,%rbp
cmp    %rbp,%r13
jne    1170 <lobatto_quad+0x50>
lea    0x1(%rbx),%eax
shr    $0x1f,%eax
lea    (%rbx,%rax,1),%esi
add    $0x1,%esi
sar    %esi
cmp    %ebx,%esi
jge    1257 <lobatto_quad+0x137>
mov    %ebx,%eax
shr    $0x1f,%eax
add    %ebx,%eax
sar    %eax
movslq %eax,%rcx
movslq %esi,%rax
movslq %ebx,%rdx
sub    %esi,%ebx
mov    %rax,%rsi
not    %rsi
add    %rdx,%rsi
and    $0x3,%rbx
je     11fb <lobatto_quad+0xdb>
nop
movsd  -0x8(%r14,%rcx,8),%xmm0
add    $0xffffffffffffffff,%rcx
movsd  %xmm0,(%r14,%rax,8)
add    $0x1,%rax
add    $0xffffffffffffffff,%rbx
jne    11e0 <lobatto_quad+0xc0>
cmp    $0x3,%rsi
jb     1257 <lobatto_quad+0x137>
lea    (%r14,%rax,8),%rsi
add    $0x18,%rsi
sub    %rdx,%rax
lea    (%r14,%rcx,8),%rcx
add    $0xfffffffffffffff8,%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rcx,%rdx,8),%xmm0
movsd  %xmm0,-0x18(%rsi)
movsd  -0x8(%rcx,%rdx,8),%xmm0
movsd  %xmm0,-0x10(%rsi)
movsd  -0x10(%rcx,%rdx,8),%xmm0
movsd  %xmm0,-0x8(%rsi)
movsd  -0x18(%rcx,%rdx,8),%xmm0
movsd  %xmm0,(%rsi)
add    $0xfffffffffffffffc,%rdx
add    $0x20,%rsi
cmp    %rdx,%rax
jne    1220 <lobatto_quad+0x100>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @lobatto_quad(double* noalias %0, double* noalias nocapture %1, i32 %2) {
%4 = tail call i32 @lobatto_nodes(double* %0, i32 %2)
%5 = add nsw i32 %2, -1
%6 = icmp slt i32 %2, 0
br i1 %6, label %24, label %7
7: ; preds = %3
%8 = sdiv i32 %5, 2
%9 = mul nsw i32 %5, %2
%10 = sitofp i32 %9 to double
%11 = add nuw nsw i32 %8, 1
%12 = zext i32 %11 to i64
br label %13
13: ; preds = %7, %13
%14 = phi i64 [ 0, %7 ], [ %22, %13 ]
%15 = getelementptr inbounds double, double* %0, i64 %14
%16 = load double, double* %15
%17 = tail call double @legendre(i32 %5, double %16)
%18 = fmul double %17, %10
%19 = fmul double %17, %18
%20 = fdiv double 2.000000e+00, %19
%21 = getelementptr inbounds double, double* %1, i64 %14
store double %20, double* %21
%22 = add nuw nsw i64 %14, 1
%23 = icmp eq i64 %22, %12
br i1 %23, label %24, label %13
24: ; preds = %13, %3
%25 = add nsw i32 %2, 1
%26 = sdiv i32 %25, 2
%27 = icmp slt i32 %26, %2
br i1 %27, label %28, label %77
28: ; preds = %24
%29 = sdiv i32 %2, 2
%30 = sext i32 %29 to i64
%31 = sext i32 %26 to i64
%32 = sext i32 %2 to i64
%33 = sub nsw i64 %32, %31
%34 = xor i64 %31, -1
%35 = add nsw i64 %34, %32
%36 = and i64 %33, 3
%37 = icmp eq i64 %36, 0
br i1 %37, label %49, label %38
38: ; preds = %28, %38
%39 = phi i64 [ %46, %38 ], [ %31, %28 ]
%40 = phi i64 [ %42, %38 ], [ %30, %28 ]
%41 = phi i64 [ %47, %38 ], [ 0, %28 ]
%42 = add nsw i64 %40, -1
%43 = getelementptr inbounds double, double* %1, i64 %42
%44 = load double, double* %43
%45 = getelementptr inbounds double, double* %1, i64 %39
store double %44, double* %45
%46 = add nsw i64 %39, 1
%47 = add i64 %41, 1
%48 = icmp eq i64 %47, %36
br i1 %48, label %49, label %38
49: ; preds = %38, %28
%50 = phi i64 [ %31, %28 ], [ %46, %38 ]
%51 = phi i64 [ %30, %28 ], [ %42, %38 ]
%52 = icmp ult i64 %35, 3
br i1 %52, label %77, label %53
53: ; preds = %49, %53
%54 = phi i64 [ %75, %53 ], [ %50, %49 ]
%55 = phi i64 [ %71, %53 ], [ %51, %49 ]
%56 = add nsw i64 %55, -1
%57 = getelementptr inbounds double, double* %1, i64 %56
%58 = load double, double* %57
%59 = getelementptr inbounds double, double* %1, i64 %54
store double %58, double* %59
%60 = add nsw i64 %54, 1
%61 = add nsw i64 %55, -2
%62 = getelementptr inbounds double, double* %1, i64 %61
%63 = load double, double* %62
%64 = getelementptr inbounds double, double* %1, i64 %60
store double %63, double* %64
%65 = add nsw i64 %54, 2
%66 = add nsw i64 %55, -3
%67 = getelementptr inbounds double, double* %1, i64 %66
%68 = load double, double* %67
%69 = getelementptr inbounds double, double* %1, i64 %65
store double %68, double* %69
%70 = add nsw i64 %54, 3
%71 = add nsw i64 %55, -4
%72 = getelementptr inbounds double, double* %1, i64 %71
%73 = load double, double* %72
%74 = getelementptr inbounds double, double* %1, i64 %70
store double %73, double* %74
%75 = add nsw i64 %54, 4
%76 = icmp eq i64 %75, %32
br i1 %76, label %77, label %53
77: ; preds = %49, %53, %24
ret void
}