#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 sqrt(double);
int tpd1(int inverse, const int i[], const double p[], int Nhat,
         const double rawcrd[], double *discrd) {
  if (i[3 + inverse] != 4 || 2 < Nhat) {
    return 1;
  }
  double r, s;
  double u = rawcrd[0];
  double v = rawcrd[1];
  if (i[5]) {
    r = p[0] + p[1] * u + p[2] * v;
    v = p[3] + p[4] * u + p[5] * v;
    u = r;
    p += 6;
  }
  if (inverse)
    p += i[3];
  *discrd = p[0] + u * p[1];
  if (Nhat == 1)
    return 0;
  *discrd += v * p[2];
  if (i[6]) {
    s = u * u + v * v;
    r = sqrt(s);
    *discrd += r * p[3];
  }
  return 0;
}


x86 assembly (clang, -O0):
<tpd1>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x8(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
cmpl   $0x4,(%rax,%rcx,4)
jne    1153 <tpd1+0x43>
mov    $0x2,%eax
cmp    -0x1c(%rbp),%eax
jge    115f <tpd1+0x4f>
movl   $0x1,-0x4(%rbp)
jmp    12f9 <tpd1+0x1e9>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x28(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
cmpl   $0x0,0x14(%rax)
je     1211 <tpd1+0x101>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm1
movsd  -0x48(%rbp),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  -0x50(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x20(%rax),%xmm1
movsd  -0x48(%rbp),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movsd  0x28(%rax),%xmm0
movsd  -0x50(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%rax
add    $0x30,%rax
mov    %rax,-0x18(%rbp)
cmpl   $0x0,-0x8(%rbp)
je     1234 <tpd1+0x124>
mov    -0x10(%rbp),%rax
mov    0xc(%rax),%ecx
mov    -0x18(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
cmpl   $0x1,-0x1c(%rbp)
jne    1270 <tpd1+0x160>
movl   $0x0,-0x4(%rbp)
jmp    12f9 <tpd1+0x1e9>
movsd  -0x50(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm2
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
cmpl   $0x0,0x18(%rax)
je     12f2 <tpd1+0x1e2>
movsd  -0x48(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm2
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<tpd1>:
push   %rbp
push   %r14
push   %rbx
mov    $0x1,%ebp
cmp    $0x2,%ecx
jg     11f3 <tpd1+0xe3>
movslq %edi,%rax
cmpl   $0x4,0xc(%rsi,%rax,4)
jne    11f3 <tpd1+0xe3>
mov    %r9,%r14
mov    %rdx,%rbx
movupd (%r8),%xmm0
cmpl   $0x0,0x14(%rsi)
je     117e <tpd1+0x6e>
movsd  (%rbx),%xmm1
movsd  0x8(%rbx),%xmm2
movhpd 0x20(%rbx),%xmm2
movsd  0x10(%rbx),%xmm3
movapd %xmm0,%xmm4
unpcklpd %xmm0,%xmm4
movhpd 0x18(%rbx),%xmm1
mulpd  %xmm2,%xmm4
addpd  %xmm4,%xmm1
movhpd 0x28(%rbx),%xmm3
unpckhpd %xmm0,%xmm0
mulpd  %xmm3,%xmm0
addpd  %xmm1,%xmm0
add    $0x30,%rbx
test   %edi,%edi
je     118a <tpd1+0x7a>
movslq 0xc(%rsi),%rax
lea    (%rbx,%rax,8),%rbx
movsd  0x8(%rbx),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rbx),%xmm1
movsd  %xmm1,(%r14)
xor    %ebp,%ebp
cmp    $0x1,%ecx
je     11f3 <tpd1+0xe3>
movapd %xmm0,%xmm2
unpckhpd %xmm0,%xmm2
movsd  0x10(%rbx),%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm1,%xmm3
movsd  %xmm3,(%r14)
cmpl   $0x0,0x18(%rsi)
je     11f3 <tpd1+0xe3>
mulsd  %xmm2,%xmm2
mulsd  %xmm0,%xmm0
addsd  %xmm2,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     11df <tpd1+0xcf>
sqrtsd %xmm0,%xmm0
jmp    11e4 <tpd1+0xd4>
call   1030 <sqrt@plt>
mulsd  0x18(%rbx),%xmm0
addsd  (%r14),%xmm0
movsd  %xmm0,(%r14)
mov    %ebp,%eax
pop    %rbx
pop    %r14
pop    %rbp
ret

LLVM IR (-O0):
define i32 @tpd1(i32 %0, i32* %1, double* %2, i32 %3, double* %4, double* %5) {
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32*
%10 = alloca double*
%11 = alloca i32
%12 = alloca double*
%13 = alloca double*
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
store i32 %0, i32* %8
store i32* %1, i32** %9
store double* %2, double** %10
store i32 %3, i32* %11
store double* %4, double** %12
store double* %5, double** %13
%18 = load i32*, i32** %9
%19 = load i32, i32* %8
%20 = add nsw i32 3, %19
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds i32, i32* %18, i64 %21
%23 = load i32, i32* %22
%24 = icmp ne i32 %23, 4
br i1 %24, label %28, label %25
25: ; preds = %6
%26 = load i32, i32* %11
%27 = icmp slt i32 2, %26
br i1 %27, label %28, label %29
28: ; preds = %25, %6
store i32 1, i32* %7
br label %122
29: ; preds = %25
%30 = load double*, double** %12
%31 = getelementptr inbounds double, double* %30, i64 0
%32 = load double, double* %31
store double %32, double* %16
%33 = load double*, double** %12
%34 = getelementptr inbounds double, double* %33, i64 1
%35 = load double, double* %34
store double %35, double* %17
%36 = load i32*, i32** %9
%37 = getelementptr inbounds i32, i32* %36, i64 5
%38 = load i32, i32* %37
%39 = icmp ne i32 %38, 0
br i1 %39, label %40, label %70
40: ; preds = %29
%41 = load double*, double** %10
%42 = getelementptr inbounds double, double* %41, i64 0
%43 = load double, double* %42
%44 = load double*, double** %10
%45 = getelementptr inbounds double, double* %44, i64 1
%46 = load double, double* %45
%47 = load double, double* %16
%48 = call double @llvm.fmuladd.f64(double %46, double %47, double %43)
%49 = load double*, double** %10
%50 = getelementptr inbounds double, double* %49, i64 2
%51 = load double, double* %50
%52 = load double, double* %17
%53 = call double @llvm.fmuladd.f64(double %51, double %52, double %48)
store double %53, double* %14
%54 = load double*, double** %10
%55 = getelementptr inbounds double, double* %54, i64 3
%56 = load double, double* %55
%57 = load double*, double** %10
%58 = getelementptr inbounds double, double* %57, i64 4
%59 = load double, double* %58
%60 = load double, double* %16
%61 = call double @llvm.fmuladd.f64(double %59, double %60, double %56)
%62 = load double*, double** %10
%63 = getelementptr inbounds double, double* %62, i64 5
%64 = load double, double* %63
%65 = load double, double* %17
%66 = call double @llvm.fmuladd.f64(double %64, double %65, double %61)
store double %66, double* %17
%67 = load double, double* %14
store double %67, double* %16
%68 = load double*, double** %10
%69 = getelementptr inbounds double, double* %68, i64 6
store double* %69, double** %10
br label %70
70: ; preds = %40, %29
%71 = load i32, i32* %8
%72 = icmp ne i32 %71, 0
br i1 %72, label %73, label %80
73: ; preds = %70
%74 = load i32*, i32** %9
%75 = getelementptr inbounds i32, i32* %74, i64 3
%76 = load i32, i32* %75
%77 = load double*, double** %10
%78 = sext i32 %76 to i64
%79 = getelementptr inbounds double, double* %77, i64 %78
store double* %79, double** %10
br label %80
80: ; preds = %73, %70
%81 = load double*, double** %10
%82 = getelementptr inbounds double, double* %81, i64 0
%83 = load double, double* %82
%84 = load double, double* %16
%85 = load double*, double** %10
%86 = getelementptr inbounds double, double* %85, i64 1
%87 = load double, double* %86
%88 = call double @llvm.fmuladd.f64(double %84, double %87, double %83)
%89 = load double*, double** %13
store double %88, double* %89
%90 = load i32, i32* %11
%91 = icmp eq i32 %90, 1
br i1 %91, label %92, label %93
92: ; preds = %80
store i32 0, i32* %7
br label %122
93: ; preds = %80
%94 = load double, double* %17
%95 = load double*, double** %10
%96 = getelementptr inbounds double, double* %95, i64 2
%97 = load double, double* %96
%98 = load double*, double** %13
%99 = load double, double* %98
%100 = call double @llvm.fmuladd.f64(double %94, double %97, double %99)
store double %100, double* %98
%101 = load i32*, i32** %9
%102 = getelementptr inbounds i32, i32* %101, i64 6
%103 = load i32, i32* %102
%104 = icmp ne i32 %103, 0
br i1 %104, label %105, label %121
105: ; preds = %93
%106 = load double, double* %16
%107 = load double, double* %16
%108 = load double, double* %17
%109 = load double, double* %17
%110 = fmul double %108, %109
%111 = call double @llvm.fmuladd.f64(double %106, double %107, double %110)
store double %111, double* %15
%112 = load double, double* %15
%113 = call double @sqrt(double %112)
store double %113, double* %14
%114 = load double, double* %14
%115 = load double*, double** %10
%116 = getelementptr inbounds double, double* %115, i64 3
%117 = load double, double* %116
%118 = load double*, double** %13
%119 = load double, double* %118
%120 = call double @llvm.fmuladd.f64(double %114, double %117, double %119)
store double %120, double* %118
br label %121
121: ; preds = %105, %93
store i32 0, i32* %7
br label %122
122: ; preds = %121, %92, %28
%123 = load i32, i32* %7
ret i32 %123
}

LLVM IR (-O2):
define i32 @tpd1(i32 %0, i32* nocapture %1, double* nocapture %2, i32 %3, double* nocapture %4, double* nocapture %5) {
%7 = add nsw i32 %0, 3
%8 = sext i32 %7 to i64
%9 = getelementptr inbounds i32, i32* %1, i64 %8
%10 = load i32, i32* %9
%11 = icmp ne i32 %10, 4
%12 = icmp sgt i32 %3, 2
%13 = or i1 %12, %11
br i1 %13, label %76, label %14
14: ; preds = %6
%15 = bitcast double* %4 to <2 x double>*
%16 = load <2 x double>, <2 x double>* %15
%17 = getelementptr inbounds i32, i32* %1, i64 5
%18 = load i32, i32* %17
%19 = icmp eq i32 %18, 0
br i1 %19, label %43, label %20
20: ; preds = %14
%21 = load double, double* %2
%22 = getelementptr inbounds double, double* %2, i64 1
%23 = load double, double* %22
%24 = getelementptr inbounds double, double* %2, i64 2
%25 = load double, double* %24
%26 = getelementptr inbounds double, double* %2, i64 3
%27 = load double, double* %26
%28 = getelementptr inbounds double, double* %2, i64 4
%29 = load double, double* %28
%30 = insertelement <2 x double> poison, double %23, i64 0
%31 = insertelement <2 x double> %30, double %29, i64 1
%32 = shufflevector <2 x double> %16, <2 x double> undef, <2 x i32> zeroinitializer
%33 = insertelement <2 x double> poison, double %21, i64 0
%34 = insertelement <2 x double> %33, double %27, i64 1
%35 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %31, <2 x double> %32, <2 x double> %34)
%36 = getelementptr inbounds double, double* %2, i64 5
%37 = load double, double* %36
%38 = insertelement <2 x double> poison, double %25, i64 0
%39 = insertelement <2 x double> %38, double %37, i64 1
%40 = shufflevector <2 x double> %16, <2 x double> undef, <2 x i32> <i32 1, i32 1>
%41 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %39, <2 x double> %40, <2 x double> %35)
%42 = getelementptr inbounds double, double* %2, i64 6
br label %43
43: ; preds = %20, %14
%44 = phi double* [ %42, %20 ], [ %2, %14 ]
%45 = phi <2 x double> [ %41, %20 ], [ %16, %14 ]
%46 = icmp eq i32 %0, 0
br i1 %46, label %52, label %47
47: ; preds = %43
%48 = getelementptr inbounds i32, i32* %1, i64 3
%49 = load i32, i32* %48
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double, double* %44, i64 %50
br label %52
52: ; preds = %47, %43
%53 = phi double* [ %51, %47 ], [ %44, %43 ]
%54 = load double, double* %53
%55 = getelementptr inbounds double, double* %53, i64 1
%56 = load double, double* %55
%57 = extractelement <2 x double> %45, i64 0
%58 = tail call double @llvm.fmuladd.f64(double %57, double %56, double %54)
store double %58, double* %5
%59 = icmp eq i32 %3, 1
br i1 %59, label %76, label %60
60: ; preds = %52
%61 = getelementptr inbounds double, double* %53, i64 2
%62 = load double, double* %61
%63 = extractelement <2 x double> %45, i64 1
%64 = tail call double @llvm.fmuladd.f64(double %63, double %62, double %58)
store double %64, double* %5
%65 = getelementptr inbounds i32, i32* %1, i64 6
%66 = load i32, i32* %65
%67 = icmp eq i32 %66, 0
br i1 %67, label %76, label %68
68: ; preds = %60
%69 = fmul double %63, %63
%70 = tail call double @llvm.fmuladd.f64(double %57, double %57, double %69)
%71 = tail call double @sqrt(double %70)
%72 = getelementptr inbounds double, double* %53, i64 3
%73 = load double, double* %72
%74 = load double, double* %5
%75 = tail call double @llvm.fmuladd.f64(double %71, double %73, double %74)
store double %75, double* %5
br label %76
76: ; preds = %52, %68, %60, %6
%77 = phi i32 [ 1, %6 ], [ 0, %60 ], [ 0, %68 ], [ 0, %52 ]
ret i32 %77
}