#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 sin(double);
void grib_trigs(int m, double lambda, double *c, double *s) {
  int i;
  double a, b;
  b = sin(lambda);
  a = 1 - 2 * sin(lambda / 2.0) * sin(lambda / 2.0);
  c[0] = 1.0;
  s[0] = 0.0;
  for (i = 1; i <= m; i++) {
    c[i] = a * c[i - 1] - b * s[i - 1];
    s[i] = a * s[i - 1] + b * c[i - 1];
  }
}


x86 assembly (clang, -O0):
<grib_trigs>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
movsd  -0x10(%rbp),%xmm0
call   1030 <sin@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  0xec0(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1030 <sin@plt>
shl    $0x1,%eax
mov    %eax,-0x3c(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  0xea4(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1030 <sin@plt>
mov    -0x3c(%rbp),%ecx
imul   %eax,%ecx
mov    $0x1,%eax
sub    %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movsd  0xe71(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     125c <grib_trigs+0x14c>
movsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11a5 <grib_trigs+0x95>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<grib_trigs>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x10,%rsp
mov    %rdx,%r14
mov    %rsi,%rbx
movsd  %xmm0,0x8(%rsp)
mov    %edi,%r12d
call   1030 <sin@plt>
mov    %eax,%r15d
movsd  0x8(%rsp),%xmm0
mulsd  0xebf(%rip),%xmm0
movsd  %xmm0,0x8(%rsp)
call   1030 <sin@plt>
mov    %eax,%ebp
movsd  0x8(%rsp),%xmm0
call   1030 <sin@plt>
movabs $0x3ff0000000000000,%rcx
mov    %rcx,(%rbx)
movq   $0x0,(%r14)
test   %r12d,%r12d
jle    12e0 <grib_trigs+0x1d0>
xorps  %xmm0,%xmm0
cvtsi2sd %r15d,%xmm0
imul   %eax,%ebp
add    %ebp,%ebp
mov    $0x1,%eax
sub    %ebp,%eax
cvtsi2sd %eax,%xmm1
mov    %r12d,%ecx
lea    0x1(%rcx),%rax
lea    (%r14,%rcx,8),%rdx
add    $0x8,%rdx
cmp    %rbx,%rdx
jbe    1202 <grib_trigs+0xf2>
lea    (%rbx,%rcx,8),%rcx
add    $0x8,%rcx
cmp    %r14,%rcx
jbe    1202 <grib_trigs+0xf2>
mov    $0x1,%ecx
xorpd  %xmm2,%xmm2
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%rbx,%rcx,8),%xmm3
mulsd  %xmm0,%xmm2
movapd %xmm1,%xmm4
mulsd  %xmm3,%xmm4
subsd  %xmm2,%xmm4
movsd  %xmm4,(%rbx,%rcx,8)
mulsd  %xmm0,%xmm3
movsd  -0x8(%r14,%rcx,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm3,%xmm2
movsd  %xmm2,(%r14,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    11c0 <grib_trigs+0xb0>
jmp    12e0 <grib_trigs+0x1d0>
mov    %eax,%ecx
movsd  (%r14),%xmm2
movhpd (%rbx),%xmm2
unpcklpd %xmm1,%xmm1
unpcklpd %xmm0,%xmm0
add    $0xffffffffffffffff,%rax
cmp    $0x2,%ecx
jne    1234 <grib_trigs+0x124>
mov    $0x1,%ecx
xorpd  %xmm3,%xmm3
test   $0x1,%al
jne    12bc <grib_trigs+0x1ac>
jmp    12e0 <grib_trigs+0x1d0>
mov    %rax,%rdx
and    $0xfffffffffffffffe,%rdx
neg    %rdx
mov    $0x1,%ecx
xorpd  %xmm4,%xmm4
movapd 0xdc1(%rip),%xmm3
nop
xorpd  %xmm3,%xmm4
movapd %xmm1,%xmm5
mulpd  %xmm2,%xmm5
shufpd $0x1,%xmm4,%xmm2
mulpd  %xmm0,%xmm2
addpd  %xmm5,%xmm2
movhpd %xmm2,(%rbx,%rcx,8)
movlpd %xmm2,(%r14,%rcx,8)
movapd %xmm2,%xmm4
xorpd  %xmm3,%xmm4
movapd %xmm1,%xmm5
mulpd  %xmm2,%xmm5
shufpd $0x1,%xmm4,%xmm2
mulpd  %xmm0,%xmm2
addpd  %xmm5,%xmm2
movhpd %xmm2,0x8(%rbx,%rcx,8)
movlpd %xmm2,0x8(%r14,%rcx,8)
lea    (%rdx,%rcx,1),%rsi
add    $0x2,%rsi
add    $0x2,%rcx
movapd %xmm2,%xmm4
cmp    $0x1,%rsi
jne    1250 <grib_trigs+0x140>
movapd %xmm2,%xmm3
test   $0x1,%al
je     12e0 <grib_trigs+0x1d0>
xorpd  0xd4c(%rip),%xmm3
mulpd  %xmm2,%xmm1
shufpd $0x1,%xmm3,%xmm2
mulpd  %xmm0,%xmm2
addpd  %xmm1,%xmm2
movhpd %xmm2,(%rbx,%rcx,8)
movlpd %xmm2,(%r14,%rcx,8)
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @grib_trigs(i32 %0, double %1, double* nocapture %2, double* nocapture %3) {
%5 = tail call i32 @sin(double %1)
%6 = sitofp i32 %5 to double
%7 = fmul double %1, 5.000000e-01
%8 = tail call i32 @sin(double %7)
%9 = mul i32 %8, -2
%10 = tail call i32 @sin(double %7)
%11 = mul i32 %9, %10
%12 = or i32 %11, 1
%13 = sitofp i32 %12 to double
store double 1.000000e+00, double* %2
store double 0.000000e+00, double* %3
%14 = icmp slt i32 %0, 1
br i1 %14, label %101, label %15
15: ; preds = %4
%16 = add nuw i32 %0, 1
%17 = zext i32 %16 to i64
%18 = getelementptr double, double* %2, i64 %17
%19 = getelementptr double, double* %3, i64 %17
%20 = icmp ugt double* %19, %2
%21 = icmp ugt double* %18, %3
%22 = and i1 %20, %21
br i1 %22, label %23, label %40
23: ; preds = %15, %23
%24 = phi double [ %36, %23 ], [ 0.000000e+00, %15 ]
%25 = phi i64 [ %38, %23 ], [ 1, %15 ]
%26 = add nsw i64 %25, -1
%27 = getelementptr inbounds double, double* %2, i64 %26
%28 = load double, double* %27
%29 = getelementptr inbounds double, double* %3, i64 %26
%30 = fneg double %24
%31 = fmul double %6, %30
%32 = tail call double @llvm.fmuladd.f64(double %13, double %28, double %31)
%33 = getelementptr inbounds double, double* %2, i64 %25
store double %32, double* %33
%34 = load double, double* %29
%35 = fmul double %28, %6
%36 = tail call double @llvm.fmuladd.f64(double %13, double %34, double %35)
%37 = getelementptr inbounds double, double* %3, i64 %25
store double %36, double* %37
%38 = add nuw nsw i64 %25, 1
%39 = icmp eq i64 %38, %17
br i1 %39, label %101, label %23
40: ; preds = %15
%41 = load double, double* %2
%42 = load double, double* %3
%43 = insertelement <2 x double> poison, double %42, i64 0
%44 = insertelement <2 x double> %43, double %41, i64 1
%45 = insertelement <2 x double> poison, double %13, i64 0
%46 = shufflevector <2 x double> %45, <2 x double> poison, <2 x i32> zeroinitializer
%47 = insertelement <2 x double> poison, double %6, i64 0
%48 = shufflevector <2 x double> %47, <2 x double> poison, <2 x i32> zeroinitializer
%49 = add nsw i64 %17, -1
%50 = and i64 %49, 1
%51 = icmp eq i32 %16, 2
br i1 %51, label %85, label %52
52: ; preds = %40
%53 = and i64 %49, -2
br label %54
54: ; preds = %54, %52
%55 = phi double [ 0.000000e+00, %52 ], [ %79, %54 ]
%56 = phi i64 [ 1, %52 ], [ %80, %54 ]
%57 = phi <2 x double> [ %44, %52 ], [ %76, %54 ]
%58 = phi i64 [ 0, %52 ], [ %81, %54 ]
%59 = fneg double %55
%60 = getelementptr inbounds double, double* %2, i64 %56
%61 = extractelement <2 x double> %57, i64 1
%62 = insertelement <2 x double> poison, double %61, i64 0
%63 = insertelement <2 x double> %62, double %59, i64 1
%64 = fmul <2 x double> %63, %48
%65 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %46, <2 x double> %57, <2 x double> %64)
%66 = extractelement <2 x double> %65, i64 1
store double %66, double* %60
%67 = getelementptr inbounds double, double* %3, i64 %56
%68 = extractelement <2 x double> %65, i64 0
store double %68, double* %67
%69 = add nuw nsw i64 %56, 1
%70 = fneg double %68
%71 = getelementptr inbounds double, double* %2, i64 %69
%72 = extractelement <2 x double> %65, i64 1
%73 = insertelement <2 x double> poison, double %72, i64 0
%74 = insertelement <2 x double> %73, double %70, i64 1
%75 = fmul <2 x double> %74, %48
%76 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %46, <2 x double> %65, <2 x double> %75)
%77 = extractelement <2 x double> %76, i64 1
store double %77, double* %71
%78 = getelementptr inbounds double, double* %3, i64 %69
%79 = extractelement <2 x double> %76, i64 0
store double %79, double* %78
%80 = add nuw nsw i64 %56, 2
%81 = add i64 %58, 2
%82 = icmp eq i64 %81, %53
br i1 %82, label %83, label %54
83: ; preds = %54
%84 = extractelement <2 x double> %76, i64 0
br label %85
85: ; preds = %83, %40
%86 = phi double [ 0.000000e+00, %40 ], [ %84, %83 ]
%87 = phi i64 [ 1, %40 ], [ %80, %83 ]
%88 = phi <2 x double> [ %44, %40 ], [ %76, %83 ]
%89 = icmp eq i64 %50, 0
br i1 %89, label %101, label %90
90: ; preds = %85
%91 = fneg double %86
%92 = getelementptr inbounds double, double* %2, i64 %87
%93 = extractelement <2 x double> %88, i64 1
%94 = insertelement <2 x double> poison, double %93, i64 0
%95 = insertelement <2 x double> %94, double %91, i64 1
%96 = fmul <2 x double> %95, %48
%97 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %46, <2 x double> %88, <2 x double> %96)
%98 = extractelement <2 x double> %97, i64 1
store double %98, double* %92
%99 = getelementptr inbounds double, double* %3, i64 %87
%100 = extractelement <2 x double> %97, i64 0
store double %100, double* %99
br label %101
101: ; preds = %90, %85, %23, %4
ret void
}