#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*/ assert(int);
void cl_fill(int l, int nfunc, int ntheta, int lmax, const double *wigd,
             double *out_cl, const double *in_cf) {
  int f, i;
  double cl;
  assert(l <= lmax);
  for (f = 0; f < nfunc; f++) {
    cl = 0.0;
    for (i = 0; i < ntheta; i++)
      cl += wigd[i] * in_cf[f * ntheta + i];
    out_cl[f * (lmax + 1) + l] = cl;
  }
}


x86 assembly (clang, -O0):
<cl_fill>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
mov    -0x4(%rbp),%eax
cmp    -0x10(%rbp),%eax
setle  %al
and    $0x1,%al
movzbl %al,%edi
call   1030 <assert@plt>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11e2 <cl_fill+0xd2>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11b4 <cl_fill+0xa4>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1165 <cl_fill+0x55>
movsd  -0x30(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x10(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
add    -0x4(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    114a <cl_fill+0x3a>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cl_fill>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %r9,%r14
mov    %r8,%rbx
mov    %ecx,%ebp
mov    %edx,%r15d
mov    %esi,%r12d
mov    %edi,%r13d
xor    %edi,%edi
cmp    %ecx,%r13d
setle  %dil
call   1030 <assert@plt>
test   %r12d,%r12d
jle    124e <cl_fill+0x13e>
mov    %r12d,%eax
mov    0x40(%rsp),%r12
add    $0x1,%ebp
mov    %r13d,%ecx
movslq %r15d,%r13
movslq %ebp,%r9
movslq %ecx,%r10
mov    %eax,%r11d
mov    %r13d,%edx
lea    -0x1(%rdx),%rax
mov    %rax,(%rsp)
mov    %edx,%ebp
and    $0x3,%ebp
and    $0xfffffffc,%edx
lea    0x18(%r12),%rsi
shl    $0x3,%r13
xor    %edi,%edi
jmp    11a7 <cl_fill+0x97>
nopl   (%rax)
xorpd  %xmm0,%xmm0
mov    %rdi,%rax
imul   %r9,%rax
add    %r10,%rax
movsd  %xmm0,(%r14,%rax,8)
add    $0x1,%rdi
add    %r13,%rsi
add    %r13,%r12
cmp    %r11,%rdi
je     124e <cl_fill+0x13e>
test   %r15d,%r15d
jle    1180 <cl_fill+0x70>
cmpq   $0x3,(%rsp)
jae    11c0 <cl_fill+0xb0>
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
jmp    1217 <cl_fill+0x107>
nopl   0x0(%rax,%rax,1)
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rbx,%rcx,8),%xmm1
movsd  0x8(%rbx,%rcx,8),%xmm2
mulsd  -0x18(%rsi,%rcx,8),%xmm1
mulsd  -0x10(%rsi,%rcx,8),%xmm2
addsd  %xmm0,%xmm1
movsd  0x10(%rbx,%rcx,8),%xmm3
mulsd  -0x8(%rsi,%rcx,8),%xmm3
addsd  %xmm1,%xmm2
movsd  0x18(%rbx,%rcx,8),%xmm0
mulsd  (%rsi,%rcx,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    11d0 <cl_fill+0xc0>
test   %rbp,%rbp
je     1184 <cl_fill+0x74>
lea    (%r12,%rcx,8),%rax
lea    (%rbx,%rcx,8),%rcx
xor    %r8d,%r8d
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%r8,8),%xmm1
mulsd  (%rax,%r8,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%r8
cmp    %r8,%rbp
jne    1230 <cl_fill+0x120>
jmp    1184 <cl_fill+0x74>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @cl_fill(i32 %0, i32 %1, i32 %2, i32 %3, double* nocapture %4, double* nocapture writeonly %5, double* nocapture %6) {
%8 = icmp sle i32 %0, %3
%9 = zext i1 %8 to i32
%10 = tail call i32 @assert(i32 %9)
%11 = icmp sgt i32 %1, 0
br i1 %11, label %12, label %87
12: ; preds = %7
%13 = icmp sgt i32 %2, 0
%14 = add nsw i32 %3, 1
%15 = sext i32 %2 to i64
%16 = sext i32 %14 to i64
%17 = sext i32 %0 to i64
%18 = zext i32 %1 to i64
%19 = zext i32 %2 to i64
%20 = add nsw i64 %19, -1
%21 = and i64 %19, 3
%22 = icmp ult i64 %20, 3
%23 = and i64 %19, 4294967292
%24 = icmp eq i64 %21, 0
br label %25
25: ; preds = %12, %80
%26 = phi i64 [ 0, %12 ], [ %85, %80 ]
br i1 %13, label %27, label %80
27: ; preds = %25
%28 = mul nsw i64 %26, %15
br i1 %22, label %63, label %29
29: ; preds = %27, %29
%30 = phi i64 [ %60, %29 ], [ 0, %27 ]
%31 = phi double [ %59, %29 ], [ 0.000000e+00, %27 ]
%32 = phi i64 [ %61, %29 ], [ 0, %27 ]
%33 = getelementptr inbounds double, double* %4, i64 %30
%34 = load double, double* %33
%35 = add nsw i64 %30, %28
%36 = getelementptr inbounds double, double* %6, i64 %35
%37 = load double, double* %36
%38 = tail call double @llvm.fmuladd.f64(double %34, double %37, double %31)
%39 = or i64 %30, 1
%40 = getelementptr inbounds double, double* %4, i64 %39
%41 = load double, double* %40
%42 = add nsw i64 %39, %28
%43 = getelementptr inbounds double, double* %6, i64 %42
%44 = load double, double* %43
%45 = tail call double @llvm.fmuladd.f64(double %41, double %44, double %38)
%46 = or i64 %30, 2
%47 = getelementptr inbounds double, double* %4, i64 %46
%48 = load double, double* %47
%49 = add nsw i64 %46, %28
%50 = getelementptr inbounds double, double* %6, i64 %49
%51 = load double, double* %50
%52 = tail call double @llvm.fmuladd.f64(double %48, double %51, double %45)
%53 = or i64 %30, 3
%54 = getelementptr inbounds double, double* %4, i64 %53
%55 = load double, double* %54
%56 = add nsw i64 %53, %28
%57 = getelementptr inbounds double, double* %6, i64 %56
%58 = load double, double* %57
%59 = tail call double @llvm.fmuladd.f64(double %55, double %58, double %52)
%60 = add nuw nsw i64 %30, 4
%61 = add i64 %32, 4
%62 = icmp eq i64 %61, %23
br i1 %62, label %63, label %29
63: ; preds = %29, %27
%64 = phi double [ undef, %27 ], [ %59, %29 ]
%65 = phi i64 [ 0, %27 ], [ %60, %29 ]
%66 = phi double [ 0.000000e+00, %27 ], [ %59, %29 ]
br i1 %24, label %80, label %67
67: ; preds = %63, %67
%68 = phi i64 [ %77, %67 ], [ %65, %63 ]
%69 = phi double [ %76, %67 ], [ %66, %63 ]
%70 = phi i64 [ %78, %67 ], [ 0, %63 ]
%71 = getelementptr inbounds double, double* %4, i64 %68
%72 = load double, double* %71
%73 = add nsw i64 %68, %28
%74 = getelementptr inbounds double, double* %6, i64 %73
%75 = load double, double* %74
%76 = tail call double @llvm.fmuladd.f64(double %72, double %75, double %69)
%77 = add nuw nsw i64 %68, 1
%78 = add i64 %70, 1
%79 = icmp eq i64 %78, %21
br i1 %79, label %80, label %67
80: ; preds = %63, %67, %25
%81 = phi double [ 0.000000e+00, %25 ], [ %64, %63 ], [ %76, %67 ]
%82 = mul nsw i64 %26, %16
%83 = add nsw i64 %82, %17
%84 = getelementptr inbounds double, double* %5, i64 %83
store double %81, double* %84
%85 = add nuw nsw i64 %26, 1
%86 = icmp eq i64 %85, %18
br i1 %86, label %87, label %25
87: ; preds = %80, %7
ret void
}