#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 */
typedef int /*<<< orphan*/ FILE;

/* Variables and functions */
int BUFSIZ;
double M_PI;
double atof(char *);
int cos(double);
int /*<<< orphan*/ exit(int);
int /*<<< orphan*/ feof(int /*<<< orphan*/ *);
int /*<<< orphan*/ *fopen(char *, char *);
int fread(char *, int, int, int /*<<< orphan*/ *);
int /*<<< orphan*/ printf(char *, ...);
double sin(double);
int /*<<< orphan*/ sscanf(char *, char *, ...);
scalar_t__ strncmp(char *, char *, int);
int main(int argc, char *argv[]) {
  char key[10] = "";
  FILE *fp;
  char buffer[81] = "";
  char val[70] = "";
  int i, j, flag = 0;
  char instname[BUFSIZ];
  double a1 = 1.7330e-4, a3 = -2.09e-7, a5 = 7.95e-10;
  double ZD0 = -1., ZD1 = -1.;
  double INR0 = -1000.0, INR1 = -1000.0;
  double sec2z, tan2z;
  double a, b, z, t, PA, cP, sP, ct, st, b11, b12, b21, b22;
  if (argc < 2) {
    printf("Usage: airmass2 FITS [dtheta]\n");
    exit(-1);
  }
  if (argc == 2) {
    t = 0;
  } else {
    t = atof(argv[2]);
  }
  fp = fopen(argv[1], "r");
  if (NULL == fp) {
    printf("Cannot find file \"%s\"\n", argv[1]);
    exit(-1);
  }
  while (!feof(fp)) {
    if (80 != fread(buffer, 1, 80, fp)) {
      exit(-1);
    }
    if (strncmp("END     ", buffer, 8) == 0) {
      break;
    } else if (strncmp("ZD-STR  =", buffer, 9) == 0) {
      sscanf(buffer + 9, "%lf", &ZD0);
    } else if (strncmp("ZD-END  =", buffer, 9) == 0) {
      sscanf(buffer + 9, "%lf", &ZD1);
    } else if (strncmp("INR-STR =", buffer, 9) == 0) {
      sscanf(buffer + 9, "%lf", &INR0);
    } else if (strncmp("INR-END =", buffer, 9) == 0) {
      sscanf(buffer + 9, "%lf", &INR1);
    } else if (strncmp("INSTRUME=", buffer, 9) == 0) {
      sscanf(buffer + 9, "%s", instname);
      if (strncmp(instname, "'SuprimeCam'", 12) != 0) {
        printf(" Error: INSTRUME=%s is not supported. quit.\n", instname);
        exit(-1);
      }
    }
  }
  if (INR0 > -360 && INR1 > -360) {
    PA = (0.5 * (INR0 + INR1) - 90.) / 180. * M_PI;
  } else {
    exit(1);
  }
  if (ZD0 > 0 && ZD1 > 0) {
    z = 0.5 * (ZD0 + ZD1) / 180. * M_PI;
  } else {
    exit(1);
  }
  sec2z = 1 / cos(z) / cos(z);
  tan2z = sin(z) * sin(z) * sec2z;
  b = (a1 + tan2z * (3. * a3 + 5. * a5 * tan2z)) * sec2z;
  a = 1. - b;
  cP = cos(PA);
  sP = sin(PA);
  ct = cos(t);
  st = sin(t);
  b11 = ct * (cP * cP + a * sP * sP) - st * b * sP * cP;
  b12 = ct * b * sP * cP - st * (a * cP * cP + sP * sP);
  b21 = st * (cP * cP + a * sP * sP) + ct * b * sP * cP;
  b22 = st * b * sP * cP + ct * (a * cP * cP + sP * sP);
  printf("-b11=%g -b12=%g -b21=%g -b22=%g\n", b11, b12, b21, b22);
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x1d0,%rsp
movl   $0x0,-0x4(%rbp)
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movw   $0x0,-0x18(%rbp)
movq   $0x0,-0x20(%rbp)
xorps  %xmm0,%xmm0
movaps %xmm0,-0x1c0(%rbp)
movaps %xmm0,-0x40(%rbp)
movaps %xmm0,-0x50(%rbp)
movaps %xmm0,-0x60(%rbp)
movaps %xmm0,-0x70(%rbp)
movaps %xmm0,-0x80(%rbp)
movb   $0x0,-0x30(%rbp)
movaps %xmm0,-0xa0(%rbp)
movaps %xmm0,-0xb0(%rbp)
movaps %xmm0,-0xc0(%rbp)
movaps %xmm0,-0xd0(%rbp)
movq   $0x0,-0x92(%rbp)
movl   $0x0,-0xdc(%rbp)
mov    0x2db7(%rip),%rax
mov    (%rax),%eax
mov    %rsp,%rcx
mov    %rcx,-0xe8(%rbp)
mov    %rax,%rdx
add    $0xf,%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x1a8(%rbp)
mov    %rcx,%rsp
mov    %rax,-0xf0(%rbp)
movsd  0xdc9(%rip),%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
movsd  0xdb1(%rip),%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
movsd  0xd99(%rip),%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
movsd  0xd81(%rip),%xmm0
movsd  %xmm0,-0x110(%rbp)
ff
movsd  0xd71(%rip),%xmm0
movsd  %xmm0,-0x118(%rbp)
ff
movsd  0xd59(%rip),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
movsd  0xd49(%rip),%xmm0
movsd  %xmm0,-0x128(%rbp)
ff
cmpl   $0x2,-0x8(%rbp)
jge    12e1 <main+0x141>
lea    0xd90(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
mov    $0xffffffff,%edi
call   10c0 <exit@plt>
cmpl   $0x2,-0x8(%rbp)
jne    12fb <main+0x15b>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x158(%rbp)
ff
jmp    1310 <main+0x170>
mov    -0x10(%rbp),%rax
mov    0x10(%rax),%rdi
call   1050 <atof@plt>
movsd  %xmm0,-0x158(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rdi
lea    0xd60(%rip),%rsi
call   10b0 <fopen@plt>
mov    %rax,-0x28(%rbp)
xor    %eax,%eax
cmp    -0x28(%rbp),%rax
jne    1354 <main+0x1b4>
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rsi
lea    0xd3e(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
mov    $0xffffffff,%edi
call   10c0 <exit@plt>
jmp    1359 <main+0x1b9>
mov    -0x28(%rbp),%rdi
call   1070 <feof@plt>
cmp    $0x0,%eax
setne  %al
xor    $0xff,%al
test   $0x1,%al
jne    1377 <main+0x1d7>
jmp    156f <main+0x3cf>
lea    -0x80(%rbp),%rdi
mov    -0x28(%rbp),%rcx
mov    $0x1,%esi
mov    $0x50,%edx
call   1040 <fread@plt>
mov    %eax,%ecx
mov    $0x50,%eax
cmp    %ecx,%eax
je     13a7 <main+0x207>
mov    $0xffffffff,%edi
call   10c0 <exit@plt>
lea    -0x80(%rbp),%rsi
lea    0xce6(%rip),%rdi
mov    $0x8,%edx
call   1030 <strncmp@plt>
cmp    $0x0,%rax
jne    13cb <main+0x22b>
jmp    156f <main+0x3cf>
lea    -0x80(%rbp),%rsi
lea    0xccb(%rip),%rdi
mov    $0x9,%edx
call   1030 <strncmp@plt>
cmp    $0x0,%rax
jne    140c <main+0x26c>
lea    -0x80(%rbp),%rdi
add    $0x9,%rdi
lea    0xcb2(%rip),%rsi
lea    -0x110(%rbp),%rdx
mov    $0x0,%al
call   1090 <sscanf@plt>
jmp    1565 <main+0x3c5>
lea    -0x80(%rbp),%rsi
lea    0xc98(%rip),%rdi
mov    $0x9,%edx
call   1030 <strncmp@plt>
cmp    $0x0,%rax
jne    144d <main+0x2ad>
lea    -0x80(%rbp),%rdi
add    $0x9,%rdi
lea    0xc71(%rip),%rsi
lea    -0x118(%rbp),%rdx
mov    $0x0,%al
call   1090 <sscanf@plt>
jmp    1560 <main+0x3c0>
lea    -0x80(%rbp),%rsi
lea    0xc61(%rip),%rdi
mov    $0x9,%edx
call   1030 <strncmp@plt>
cmp    $0x0,%rax
jne    148e <main+0x2ee>
lea    -0x80(%rbp),%rdi
add    $0x9,%rdi
lea    0xc30(%rip),%rsi
lea    -0x120(%rbp),%rdx
mov    $0x0,%al
call   1090 <sscanf@plt>
jmp    155b <main+0x3bb>
lea    -0x80(%rbp),%rsi
lea    0xc2a(%rip),%rdi
mov    $0x9,%edx
call   1030 <strncmp@plt>
cmp    $0x0,%rax
jne    14cf <main+0x32f>
lea    -0x80(%rbp),%rdi
add    $0x9,%rdi
lea    0xbef(%rip),%rsi
lea    -0x128(%rbp),%rdx
mov    $0x0,%al
call   1090 <sscanf@plt>
jmp    1556 <main+0x3b6>
lea    -0x80(%rbp),%rsi
lea    0xbf3(%rip),%rdi
mov    $0x9,%edx
call   1030 <strncmp@plt>
cmp    $0x0,%rax
jne    1551 <main+0x3b1>
mov    -0x1a8(%rbp),%rdx
lea    -0x80(%rbp),%rdi
add    $0x9,%rdi
lea    0xbd3(%rip),%rsi
mov    $0x0,%al
call   1090 <sscanf@plt>
mov    -0x1a8(%rbp),%rdi
lea    0xbc1(%rip),%rsi
mov    $0xc,%edx
call   1030 <strncmp@plt>
cmp    $0x0,%rax
je     154c <main+0x3ac>
mov    -0x1a8(%rbp),%rsi
lea    0xbac(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
mov    $0xffffffff,%edi
call   10c0 <exit@plt>
jmp    1551 <main+0x3b1>
jmp    1556 <main+0x3b6>
jmp    155b <main+0x3bb>
jmp    1560 <main+0x3c0>
jmp    1565 <main+0x3c5>
jmp    156a <main+0x3ca>
jmp    1359 <main+0x1b9>
movsd  -0x120(%rbp),%xmm0
ff
movsd  0xaa9(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    15ef <main+0x44f>
movsd  -0x128(%rbp),%xmm0
ff
movsd  0xa8f(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    15ef <main+0x44f>
movsd  -0x120(%rbp),%xmm0
ff
addsd  -0x128(%rbp),%xmm0
ff
movsd  0xa7d(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xa79(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  0xa5d(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    0x2a0a(%rip),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x160(%rbp)
ff
jmp    15f9 <main+0x459>
mov    $0x1,%edi
call   10c0 <exit@plt>
movsd  -0x110(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1663 <main+0x4c3>
movsd  -0x118(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1663 <main+0x4c3>
movsd  -0x110(%rbp),%xmm1
ff
addsd  -0x118(%rbp),%xmm1
ff
movsd  0x9fd(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  0x9e9(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    0x2996(%rip),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x150(%rbp)
ff
jmp    166d <main+0x4cd>
mov    $0x1,%edi
call   10c0 <exit@plt>
movsd  -0x150(%rbp),%xmm0
ff
call   10a0 <cos@plt>
mov    %eax,%ecx
mov    $0x1,%eax
cltd
idiv   %ecx
mov    %eax,-0x1cc(%rbp)
movsd  -0x150(%rbp),%xmm0
ff
call   10a0 <cos@plt>
mov    %eax,%ecx
mov    -0x1cc(%rbp),%eax
cltd
idiv   %ecx
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
movsd  -0x150(%rbp),%xmm0
ff
call   1080 <sin@plt>
movsd  %xmm0,-0x1c8(%rbp)
ff
movsd  -0x150(%rbp),%xmm0
ff
call   1080 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x1c8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
mulsd  -0x130(%rbp),%xmm0
ff
movsd  %xmm0,-0x138(%rbp)
ff
movsd  -0xf8(%rbp),%xmm1
ff
movsd  -0x138(%rbp),%xmm0
ff
movsd  -0x100(%rbp),%xmm2
ff
movsd  0x949(%rip),%xmm3
mulsd  -0x108(%rbp),%xmm3
ff
mulsd  -0x138(%rbp),%xmm3
ff
movsd  0x929(%rip),%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mulsd  -0x130(%rbp),%xmm0
ff
movsd  %xmm0,-0x148(%rbp)
ff
movsd  0x8f9(%rip),%xmm0
subsd  -0x148(%rbp),%xmm0
ff
movsd  %xmm0,-0x140(%rbp)
ff
movsd  -0x160(%rbp),%xmm0
ff
call   10a0 <cos@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x168(%rbp)
ff
movsd  -0x160(%rbp),%xmm0
ff
call   1080 <sin@plt>
movsd  %xmm0,-0x170(%rbp)
ff
movsd  -0x158(%rbp),%xmm0
ff
call   10a0 <cos@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x178(%rbp)
ff
movsd  -0x158(%rbp),%xmm0
ff
call   1080 <sin@plt>
movsd  %xmm0,-0x180(%rbp)
ff
movsd  -0x178(%rbp),%xmm0
ff
movsd  -0x168(%rbp),%xmm2
ff
movsd  -0x168(%rbp),%xmm3
ff
movsd  -0x140(%rbp),%xmm1
ff
mulsd  -0x170(%rbp),%xmm1
ff
mulsd  -0x170(%rbp),%xmm1
ff
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
movsd  -0x180(%rbp),%xmm1
ff
mulsd  -0x148(%rbp),%xmm1
ff
mulsd  -0x170(%rbp),%xmm1
ff
mulsd  -0x168(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x188(%rbp)
ff
movsd  -0x178(%rbp),%xmm0
ff
mulsd  -0x148(%rbp),%xmm0
ff
mulsd  -0x170(%rbp),%xmm0
ff
movsd  -0x168(%rbp),%xmm2
ff
movsd  -0x180(%rbp),%xmm1
ff
movsd  -0x140(%rbp),%xmm3
ff
mulsd  -0x168(%rbp),%xmm3
ff
movsd  -0x168(%rbp),%xmm5
ff
movsd  -0x170(%rbp),%xmm4
ff
mulsd  -0x170(%rbp),%xmm4
ff
mulsd  %xmm5,%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x190(%rbp)
ff
movsd  -0x180(%rbp),%xmm0
ff
movsd  -0x168(%rbp),%xmm2
ff
movsd  -0x168(%rbp),%xmm3
ff
movsd  -0x140(%rbp),%xmm1
ff
mulsd  -0x170(%rbp),%xmm1
ff
mulsd  -0x170(%rbp),%xmm1
ff
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
movsd  -0x178(%rbp),%xmm1
ff
mulsd  -0x148(%rbp),%xmm1
ff
mulsd  -0x170(%rbp),%xmm1
ff
mulsd  -0x168(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x198(%rbp)
ff
movsd  -0x180(%rbp),%xmm0
ff
mulsd  -0x148(%rbp),%xmm0
ff
mulsd  -0x170(%rbp),%xmm0
ff
movsd  -0x168(%rbp),%xmm2
ff
movsd  -0x178(%rbp),%xmm1
ff
movsd  -0x140(%rbp),%xmm3
ff
mulsd  -0x168(%rbp),%xmm3
ff
movsd  -0x168(%rbp),%xmm5
ff
movsd  -0x170(%rbp),%xmm4
ff
mulsd  -0x170(%rbp),%xmm4
ff
mulsd  %xmm5,%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1a0(%rbp)
ff
movsd  -0x188(%rbp),%xmm0
ff
movsd  -0x190(%rbp),%xmm1
ff
movsd  -0x198(%rbp),%xmm2
ff
movsd  -0x1a0(%rbp),%xmm3
ff
lea    0x75c(%rip),%rdi
mov    $0x4,%al
call   1060 <printf@plt>
movl   $0x0,-0x4(%rbp)
mov    -0xe8(%rbp),%rax
mov    %rax,%rsp
mov    -0x4(%rbp),%eax
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xb8,%rsp
mov    %rsp,%r14
xorpd  %xmm0,%xmm0
movapd %xmm0,-0xa0(%rbp)
ff
movapd %xmm0,-0xb0(%rbp)
ff
movapd %xmm0,-0xc0(%rbp)
ff
movapd %xmm0,-0xd0(%rbp)
ff
movapd %xmm0,-0xe0(%rbp)
ff
movb   $0x0,-0x90(%rbp)
mov    0x2dd7(%rip),%rax
mov    (%rax),%eax
mov    %rsp,%rcx
add    $0xf,%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%rcx
mov    %rcx,-0x30(%rbp)
mov    %rcx,%rsp
movabs $0xbff0000000000000,%rax
mov    %rax,-0x70(%rbp)
mov    %rax,-0x68(%rbp)
movabs $0xc08f400000000000,%rax
40 8f c0
mov    %rax,-0x60(%rbp)
mov    %rax,-0x58(%rbp)
cmp    $0x1,%edi
jle    163c <main+0x48c>
mov    %rsi,%rbx
xorpd  %xmm0,%xmm0
cmp    $0x2,%edi
je     125a <main+0xaa>
mov    0x10(%rbx),%rdi
call   1060 <atof@plt>
movsd  %xmm0,-0x50(%rbp)
mov    0x8(%rbx),%rdi
lea    0xdd6(%rip),%rsi
call   10c0 <fopen@plt>
test   %rax,%rax
je     1652 <main+0x4a2>
mov    %rax,%r12
mov    %r14,-0x80(%rbp)
mov    %rax,%rdi
call   1080 <feof@plt>
test   %eax,%eax
jne    1632 <main+0x482>
lea    -0xd7(%rbp),%r14
lea    -0xe0(%rbp),%rbx
lea    0xdb5(%rip),%r13
lea    0xdb7(%rip),%r15
jmp    12d1 <main+0x121>
nopl   (%rax)
mov    %r14,%rdi
lea    0xdb2(%rip),%rsi
xor    %eax,%eax
call   10a0 <sscanf@plt>
mov    %r12,%rdi
call   1080 <feof@plt>
test   %eax,%eax
jne    13f9 <main+0x249>
mov    %rbx,%rdi
mov    $0x1,%esi
mov    $0x50,%edx
mov    %r12,%rcx
call   1050 <fread@plt>
cmp    $0x50,%eax
jne    1628 <main+0x478>
mov    %r13,%rdi
mov    %rbx,%rsi
mov    $0x8,%edx
call   1030 <strncmp@plt>
test   %rax,%rax
je     13f9 <main+0x249>
mov    %r15,%rdi
mov    %rbx,%rsi
mov    $0x9,%edx
call   1030 <strncmp@plt>
lea    -0x70(%rbp),%rdx
test   %rax,%rax
je     12b0 <main+0x100>
lea    0xd48(%rip),%rdi
mov    %rbx,%rsi
mov    $0x9,%edx
call   1030 <strncmp@plt>
lea    -0x68(%rbp),%rdx
test   %rax,%rax
je     12b0 <main+0x100>
lea    0xd31(%rip),%rdi
mov    %rbx,%rsi
mov    $0x9,%edx
call   1030 <strncmp@plt>
lea    -0x60(%rbp),%rdx
test   %rax,%rax
je     12b0 <main+0x100>
lea    0xd1a(%rip),%rdi
mov    %rbx,%rsi
mov    $0x9,%edx
call   1030 <strncmp@plt>
lea    -0x58(%rbp),%rdx
test   %rax,%rax
je     12b0 <main+0x100>
lea    0xd03(%rip),%rdi
mov    %rbx,%rsi
mov    $0x9,%edx
call   1030 <strncmp@plt>
test   %rax,%rax
jne    12c1 <main+0x111>
mov    %r14,%rdi
lea    0xced(%rip),%rsi
mov    %r15,%r13
mov    %r14,%r15
mov    -0x30(%rbp),%r14
mov    %r14,%rdx
xor    %eax,%eax
call   10a0 <sscanf@plt>
mov    %r14,%rdi
mov    %r15,%r14
mov    %r13,%r15
lea    0xc8a(%rip),%r13
lea    0xcc5(%rip),%rsi
mov    $0xc,%edx
call   1030 <strncmp@plt>
test   %rax,%rax
je     12c1 <main+0x111>
lea    0xcb8(%rip),%rdi
mov    -0x30(%rbp),%rsi
jmp    165d <main+0x4ad>
movsd  -0x60(%rbp),%xmm3
ucomisd 0xbfa(%rip),%xmm3
jbe    1632 <main+0x482>
movsd  -0x58(%rbp),%xmm1
ucomisd 0xbe7(%rip),%xmm1
jbe    1632 <main+0x482>
movsd  -0x70(%rbp),%xmm0
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm0
mov    -0x80(%rbp),%r14
jbe    1632 <main+0x482>
movsd  -0x68(%rbp),%xmm4
ucomisd %xmm2,%xmm4
jbe    1632 <main+0x482>
addsd  %xmm1,%xmm3
mulsd  0xbb7(%rip),%xmm3
addsd  0xbb7(%rip),%xmm3
divsd  0xbb7(%rip),%xmm3
mov    0x2b80(%rip),%rax
movsd  (%rax),%xmm1
addsd  %xmm4,%xmm0
mulsd  0xb90(%rip),%xmm0
mulsd  %xmm1,%xmm3
movsd  %xmm3,-0x38(%rbp)
divsd  0xb8f(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
call   10b0 <cos@plt>
mov    %eax,%ebx
lea    0x1(%rbx),%eax
xor    %ecx,%ecx
cmp    $0x3,%eax
cmovae %ecx,%ebx
movsd  -0x30(%rbp),%xmm0
call   10b0 <cos@plt>
mov    %eax,%ecx
mov    %ebx,%eax
cltd
idiv   %ecx
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x30(%rbp),%xmm0
call   1090 <sin@plt>
movsd  %xmm0,-0x40(%rbp)
movsd  -0x30(%rbp),%xmm0
call   1090 <sin@plt>
movsd  -0x40(%rbp),%xmm2
mulsd  %xmm0,%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm1,%xmm2
movsd  0xb2c(%rip),%xmm0
mulsd  %xmm2,%xmm0
addsd  0xb28(%rip),%xmm0
mulsd  %xmm2,%xmm0
addsd  0xb24(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xb1b(%rip),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
call   10b0 <cos@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1090 <sin@plt>
movsd  %xmm0,-0x78(%rbp)
movsd  -0x50(%rbp),%xmm0
call   10b0 <cos@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x50(%rbp),%xmm0
call   1090 <sin@plt>
movsd  -0x78(%rbp),%xmm6
movapd %xmm6,%xmm1
movsd  -0x48(%rbp),%xmm5
mulsd  %xmm5,%xmm1
mulsd  %xmm6,%xmm1
movsd  -0x40(%rbp),%xmm4
movapd %xmm4,%xmm2
mulsd  %xmm4,%xmm2
addsd  %xmm1,%xmm2
movapd %xmm0,%xmm3
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm1,%xmm3
mulsd  %xmm6,%xmm3
movsd  -0x38(%rbp),%xmm7
movapd %xmm7,%xmm9
mulsd  %xmm2,%xmm9
mulsd  %xmm4,%xmm3
subsd  %xmm3,%xmm9
mulsd  %xmm7,%xmm1
movapd %xmm7,%xmm8
mulsd  %xmm6,%xmm1
mulsd  %xmm4,%xmm5
mulsd  %xmm6,%xmm6
mulsd  %xmm4,%xmm5
movapd %xmm4,%xmm7
addsd  %xmm6,%xmm5
movapd %xmm5,%xmm4
mulsd  %xmm0,%xmm4
mulsd  %xmm7,%xmm1
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
subsd  %xmm4,%xmm1
mulsd  %xmm8,%xmm5
addsd  %xmm3,%xmm5
movapd %xmm5,%xmm3
lea    0xacd(%rip),%rdi
movapd %xmm9,%xmm0
mov    $0x4,%al
call   1070 <printf@plt>
mov    %r14,%rsp
xor    %eax,%eax
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    $0xffffffff,%edi
call   10d0 <exit@plt>
mov    $0x1,%edi
call   10d0 <exit@plt>
lea    0xab3(%rip),%rdi
call   1040 <puts@plt>
mov    $0xffffffff,%edi
call   10d0 <exit@plt>
mov    0x8(%rbx),%rsi
lea    0x9e5(%rip),%rdi
xor    %eax,%eax
call   1070 <printf@plt>
mov    $0xffffffff,%edi
call   10d0 <exit@plt>

LLVM IR (-O0):
@BUFSIZ = global i32 0
@.str = [31 x i8] c"Usage: airmass2 FITS [dtheta]\0A\00"
@.str.1 = [2 x i8] c"r\00"
@.str.2 = [23 x i8] c"Cannot find file \22%s\22\0A\00"
@.str.3 = [9 x i8] c"END \00"
@.str.4 = [10 x i8] c"ZD-STR =\00"
@.str.5 = [4 x i8] c"%lf\00"
@.str.6 = [10 x i8] c"ZD-END =\00"
@.str.7 = [10 x i8] c"INR-STR =\00"
@.str.8 = [10 x i8] c"INR-END =\00"
@.str.9 = [10 x i8] c"INSTRUME=\00"
@.str.10 = [3 x i8] c"%s\00"
@.str.11 = [13 x i8] c"'SuprimeCam'\00"
@.str.12 = [45 x i8] c" Error: INSTRUME=%s is not supported. quit.\0A\00"
@M_PI = global double 0.000000e+00
@.str.13 = [33 x i8] c"-b11=%g -b12=%g -b21=%g -b22=%g\0A\00"
define i32 @main(i32 %0, i8** %1) {
%3 = alloca i32
%4 = alloca i32
%5 = alloca i8**
%6 = alloca [10 x i8]
%7 = alloca i32*
%8 = alloca [81 x i8]
%9 = alloca [70 x i8]
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i8*
%14 = alloca i64
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double
%26 = alloca double
%27 = alloca double
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca double
%32 = alloca double
%33 = alloca double
%34 = alloca double
%35 = alloca double
%36 = alloca double
store i32 0, i32* %3
store i32 %0, i32* %4
store i8** %1, i8*** %5
%37 = bitcast [10 x i8]* %6 to i8*
%38 = bitcast [81 x i8]* %8 to i8*
%39 = bitcast [70 x i8]* %9 to i8*
store i32 0, i32* %12
%40 = load i32, i32* @BUFSIZ
%41 = zext i32 %40 to i64
%42 = call i8* @llvm.stacksave()
store i8* %42, i8** %13
%43 = alloca i8, i64 %41
store i64 %41, i64* %14
store double 1.733000e-04, double* %15
store double -2.090000e-07, double* %16
store double 7.950000e-10, double* %17
store double -1.000000e+00, double* %18
store double -1.000000e+00, double* %19
store double -1.000000e+03, double* %20
store double -1.000000e+03, double* %21
%44 = load i32, i32* %4
%45 = icmp slt i32 %44, 2
br i1 %45, label %46, label %49
46: ; preds = %2
%47 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([31 x i8], [31 x i8]* @.str, i64 0, i64 0))
%48 = call i32 @exit(i32 -1)
unreachable
49: ; preds = %2
%50 = load i32, i32* %4
%51 = icmp eq i32 %50, 2
br i1 %51, label %52, label %53
52: ; preds = %49
store double 0.000000e+00, double* %27
br label %58
53: ; preds = %49
%54 = load i8**, i8*** %5
%55 = getelementptr inbounds i8*, i8** %54, i64 2
%56 = load i8*, i8** %55
%57 = call double @atof(i8* %56)
store double %57, double* %27
br label %58
58: ; preds = %53, %52
%59 = load i8**, i8*** %5
%60 = getelementptr inbounds i8*, i8** %59, i64 1
%61 = load i8*, i8** %60
%62 = call i32* @fopen(i8* %61, i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.1, i64 0, i64 0))
store i32* %62, i32** %7
%63 = load i32*, i32** %7
%64 = icmp eq i32* null, %63
br i1 %64, label %65, label %71
65: ; preds = %58
%66 = load i8**, i8*** %5
%67 = getelementptr inbounds i8*, i8** %66, i64 1
%68 = load i8*, i8** %67
%69 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([23 x i8], [23 x i8]* @.str.2, i64 0, i64 0), i8* %68)
%70 = call i32 @exit(i32 -1)
unreachable
71: ; preds = %58
br label %72
72: ; preds = %140, %71
%73 = load i32*, i32** %7
%74 = call i32 @feof(i32* %73)
%75 = icmp ne i32 %74, 0
%76 = xor i1 %75, true
br i1 %76, label %77, label %141
77: ; preds = %72
%78 = getelementptr inbounds [81 x i8], [81 x i8]* %8, i64 0, i64 0
%79 = load i32*, i32** %7
%80 = call i32 @fread(i8* %78, i32 1, i32 80, i32* %79)
%81 = icmp ne i32 80, %80
br i1 %81, label %82, label %84
82: ; preds = %77
%83 = call i32 @exit(i32 -1)
unreachable
84: ; preds = %77
%85 = getelementptr inbounds [81 x i8], [81 x i8]* %8, i64 0, i64 0
%86 = call i64 @strncmp(i8* getelementptr inbounds ([9 x i8], [9 x i8]* @.str.3, i64 0, i64 0), i8* %85, i32 8)
%87 = icmp eq i64 %86, 0
br i1 %87, label %88, label %89
88: ; preds = %84
br label %141
89: ; preds = %84
%90 = getelementptr inbounds [81 x i8], [81 x i8]* %8, i64 0, i64 0
%91 = call i64 @strncmp(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str.4, i64 0, i64 0), i8* %90, i32 9)
%92 = icmp eq i64 %91, 0
br i1 %92, label %93, label %97
93: ; preds = %89
%94 = getelementptr inbounds [81 x i8], [81 x i8]* %8, i64 0, i64 0
%95 = getelementptr inbounds i8, i8* %94, i64 9
%96 = call i32 (i8*, i8*, ...) @sscanf(i8* %95, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.5, i64 0, i64 0), double* %18)
br label %139
97: ; preds = %89
%98 = getelementptr inbounds [81 x i8], [81 x i8]* %8, i64 0, i64 0
%99 = call i64 @strncmp(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str.6, i64 0, i64 0), i8* %98, i32 9)
%100 = icmp eq i64 %99, 0
br i1 %100, label %101, label %105
101: ; preds = %97
%102 = getelementptr inbounds [81 x i8], [81 x i8]* %8, i64 0, i64 0
%103 = getelementptr inbounds i8, i8* %102, i64 9
%104 = call i32 (i8*, i8*, ...) @sscanf(i8* %103, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.5, i64 0, i64 0), double* %19)
br label %138
105: ; preds = %97
%106 = getelementptr inbounds [81 x i8], [81 x i8]* %8, i64 0, i64 0
%107 = call i64 @strncmp(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str.7, i64 0, i64 0), i8* %106, i32 9)
%108 = icmp eq i64 %107, 0
br i1 %108, label %109, label %113
109: ; preds = %105
%110 = getelementptr inbounds [81 x i8], [81 x i8]* %8, i64 0, i64 0
%111 = getelementptr inbounds i8, i8* %110, i64 9
%112 = call i32 (i8*, i8*, ...) @sscanf(i8* %111, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.5, i64 0, i64 0), double* %20)
br label %137
113: ; preds = %105
%114 = getelementptr inbounds [81 x i8], [81 x i8]* %8, i64 0, i64 0
%115 = call i64 @strncmp(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str.8, i64 0, i64 0), i8* %114, i32 9)
%116 = icmp eq i64 %115, 0
br i1 %116, label %117, label %121
117: ; preds = %113
%118 = getelementptr inbounds [81 x i8], [81 x i8]* %8, i64 0, i64 0
%119 = getelementptr inbounds i8, i8* %118, i64 9
%120 = call i32 (i8*, i8*, ...) @sscanf(i8* %119, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.5, i64 0, i64 0), double* %21)
br label %136
121: ; preds = %113
%122 = getelementptr inbounds [81 x i8], [81 x i8]* %8, i64 0, i64 0
%123 = call i64 @strncmp(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str.9, i64 0, i64 0), i8* %122, i32 9)
%124 = icmp eq i64 %123, 0
br i1 %124, label %125, label %135
125: ; preds = %121
%126 = getelementptr inbounds [81 x i8], [81 x i8]* %8, i64 0, i64 0
%127 = getelementptr inbounds i8, i8* %126, i64 9
%128 = call i32 (i8*, i8*, ...) @sscanf(i8* %127, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.10, i64 0, i64 0), i8* %43)
%129 = call i64 @strncmp(i8* %43, i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str.11, i64 0, i64 0), i32 12)
%130 = icmp ne i64 %129, 0
br i1 %130, label %131, label %134
131: ; preds = %125
%132 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([45 x i8], [45 x i8]* @.str.12, i64 0, i64 0), i8* %43)
%133 = call i32 @exit(i32 -1)
unreachable
134: ; preds = %125
br label %135
135: ; preds = %134, %121
br label %136
136: ; preds = %135, %117
br label %137
137: ; preds = %136, %109
br label %138
138: ; preds = %137, %101
br label %139
139: ; preds = %138, %93
br label %140
140: ; preds = %139
br label %72
141: ; preds = %88, %72
%142 = load double, double* %20
%143 = fcmp ogt double %142, -3.600000e+02
br i1 %143, label %144, label %155
144: ; preds = %141
%145 = load double, double* %21
%146 = fcmp ogt double %145, -3.600000e+02
br i1 %146, label %147, label %155
147: ; preds = %144
%148 = load double, double* %20
%149 = load double, double* %21
%150 = fadd double %148, %149
%151 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %150, double -9.000000e+01)
%152 = fdiv double %151, 1.800000e+02
%153 = load double, double* @M_PI
%154 = fmul double %152, %153
store double %154, double* %28
br label %157
155: ; preds = %144, %141
%156 = call i32 @exit(i32 1)
unreachable
157: ; preds = %147
%158 = load double, double* %18
%159 = fcmp ogt double %158, 0.000000e+00
br i1 %159, label %160, label %171
160: ; preds = %157
%161 = load double, double* %19
%162 = fcmp ogt double %161, 0.000000e+00
br i1 %162, label %163, label %171
163: ; preds = %160
%164 = load double, double* %18
%165 = load double, double* %19
%166 = fadd double %164, %165
%167 = fmul double 5.000000e-01, %166
%168 = fdiv double %167, 1.800000e+02
%169 = load double, double* @M_PI
%170 = fmul double %168, %169
store double %170, double* %26
br label %173
171: ; preds = %160, %157
%172 = call i32 @exit(i32 1)
unreachable
173: ; preds = %163
%174 = load double, double* %26
%175 = call i32 @cos(double %174)
%176 = sdiv i32 1, %175
%177 = load double, double* %26
%178 = call i32 @cos(double %177)
%179 = sdiv i32 %176, %178
%180 = sitofp i32 %179 to double
store double %180, double* %22
%181 = load double, double* %26
%182 = call double @sin(double %181)
%183 = load double, double* %26
%184 = call double @sin(double %183)
%185 = fmul double %182, %184
%186 = load double, double* %22
%187 = fmul double %185, %186
store double %187, double* %23
%188 = load double, double* %15
%189 = load double, double* %23
%190 = load double, double* %16
%191 = load double, double* %17
%192 = fmul double 5.000000e+00, %191
%193 = load double, double* %23
%194 = fmul double %192, %193
%195 = call double @llvm.fmuladd.f64(double 3.000000e+00, double %190, double %194)
%196 = call double @llvm.fmuladd.f64(double %189, double %195, double %188)
%197 = load double, double* %22
%198 = fmul double %196, %197
store double %198, double* %25
%199 = load double, double* %25
%200 = fsub double 1.000000e+00, %199
store double %200, double* %24
%201 = load double, double* %28
%202 = call i32 @cos(double %201)
%203 = sitofp i32 %202 to double
store double %203, double* %29
%204 = load double, double* %28
%205 = call double @sin(double %204)
store double %205, double* %30
%206 = load double, double* %27
%207 = call i32 @cos(double %206)
%208 = sitofp i32 %207 to double
store double %208, double* %31
%209 = load double, double* %27
%210 = call double @sin(double %209)
store double %210, double* %32
%211 = load double, double* %31
%212 = load double, double* %29
%213 = load double, double* %29
%214 = load double, double* %24
%215 = load double, double* %30
%216 = fmul double %214, %215
%217 = load double, double* %30
%218 = fmul double %216, %217
%219 = call double @llvm.fmuladd.f64(double %212, double %213, double %218)
%220 = load double, double* %32
%221 = load double, double* %25
%222 = fmul double %220, %221
%223 = load double, double* %30
%224 = fmul double %222, %223
%225 = load double, double* %29
%226 = fmul double %224, %225
%227 = fneg double %226
%228 = call double @llvm.fmuladd.f64(double %211, double %219, double %227)
store double %228, double* %33
%229 = load double, double* %31
%230 = load double, double* %25
%231 = fmul double %229, %230
%232 = load double, double* %30
%233 = fmul double %231, %232
%234 = load double, double* %29
%235 = load double, double* %32
%236 = load double, double* %24
%237 = load double, double* %29
%238 = fmul double %236, %237
%239 = load double, double* %29
%240 = load double, double* %30
%241 = load double, double* %30
%242 = fmul double %240, %241
%243 = call double @llvm.fmuladd.f64(double %238, double %239, double %242)
%244 = fmul double %235, %243
%245 = fneg double %244
%246 = call double @llvm.fmuladd.f64(double %233, double %234, double %245)
store double %246, double* %34
%247 = load double, double* %32
%248 = load double, double* %29
%249 = load double, double* %29
%250 = load double, double* %24
%251 = load double, double* %30
%252 = fmul double %250, %251
%253 = load double, double* %30
%254 = fmul double %252, %253
%255 = call double @llvm.fmuladd.f64(double %248, double %249, double %254)
%256 = load double, double* %31
%257 = load double, double* %25
%258 = fmul double %256, %257
%259 = load double, double* %30
%260 = fmul double %258, %259
%261 = load double, double* %29
%262 = fmul double %260, %261
%263 = call double @llvm.fmuladd.f64(double %247, double %255, double %262)
store double %263, double* %35
%264 = load double, double* %32
%265 = load double, double* %25
%266 = fmul double %264, %265
%267 = load double, double* %30
%268 = fmul double %266, %267
%269 = load double, double* %29
%270 = load double, double* %31
%271 = load double, double* %24
%272 = load double, double* %29
%273 = fmul double %271, %272
%274 = load double, double* %29
%275 = load double, double* %30
%276 = load double, double* %30
%277 = fmul double %275, %276
%278 = call double @llvm.fmuladd.f64(double %273, double %274, double %277)
%279 = fmul double %270, %278
%280 = call double @llvm.fmuladd.f64(double %268, double %269, double %279)
store double %280, double* %36
%281 = load double, double* %33
%282 = load double, double* %34
%283 = load double, double* %35
%284 = load double, double* %36
%285 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([33 x i8], [33 x i8]* @.str.13, i64 0, i64 0), double %281, double %282, double %283, double %284)
store i32 0, i32* %3
%286 = load i8*, i8** %13
%287 = load i32, i32* %3
ret i32 %287
}

LLVM IR (-O2):
@BUFSIZ = global i32 0
@.str.1 = [2 x i8] c"r\00"
@.str.2 = [23 x i8] c"Cannot find file \22%s\22\0A\00"
@.str.3 = [9 x i8] c"END \00"
@.str.4 = [10 x i8] c"ZD-STR =\00"
@.str.5 = [4 x i8] c"%lf\00"
@.str.6 = [10 x i8] c"ZD-END =\00"
@.str.7 = [10 x i8] c"INR-STR =\00"
@.str.8 = [10 x i8] c"INR-END =\00"
@.str.9 = [10 x i8] c"INSTRUME=\00"
@.str.10 = [3 x i8] c"%s\00"
@.str.11 = [13 x i8] c"'SuprimeCam'\00"
@.str.12 = [45 x i8] c" Error: INSTRUME=%s is not supported. quit.\0A\00"
@M_PI = global double 0.000000e+00
@.str.13 = [33 x i8] c"-b11=%g -b12=%g -b21=%g -b22=%g\0A\00"
@str = [30 x i8] c"Usage: airmass2 FITS [dtheta]\00"
define i32 @main(i32 %0, i8** nocapture %1) {
%3 = alloca [81 x i8]
%4 = alloca double
%5 = alloca double
%6 = alloca double
%7 = alloca double
%8 = getelementptr inbounds [81 x i8], [81 x i8]* %3, i64 0, i64 0
%9 = load i32, i32* @BUFSIZ
%10 = zext i32 %9 to i64
%11 = tail call i8* @llvm.stacksave()
%12 = alloca i8, i64 %10
%13 = bitcast double* %4 to i8*
store double -1.000000e+00, double* %4
%14 = bitcast double* %5 to i8*
store double -1.000000e+00, double* %5
%15 = bitcast double* %6 to i8*
store double -1.000000e+03, double* %6
%16 = bitcast double* %7 to i8*
store double -1.000000e+03, double* %7
%17 = icmp slt i32 %0, 2
br i1 %17, label %18, label %21
18: ; preds = %2
%19 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([30 x i8], [30 x i8]* @str, i64 0, i64 0))
%20 = tail call i32 @exit(i32 -1)
unreachable
21: ; preds = %2
%22 = icmp eq i32 %0, 2
br i1 %22, label %27, label %23
23: ; preds = %21
%24 = getelementptr inbounds i8*, i8** %1, i64 2
%25 = load i8*, i8** %24
%26 = tail call double @atof(i8* %25)
br label %27
27: ; preds = %21, %23
%28 = phi double [ %26, %23 ], [ 0.000000e+00, %21 ]
%29 = getelementptr inbounds i8*, i8** %1, i64 1
%30 = load i8*, i8** %29
%31 = tail call i32* @fopen(i8* %30, i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.1, i64 0, i64 0))
%32 = icmp eq i32* %31, null
br i1 %32, label %38, label %33
33: ; preds = %27
%34 = call i32 @feof(i32* %31)
%35 = icmp eq i32 %34, 0
br i1 %35, label %36, label %95
36: ; preds = %33
%37 = getelementptr inbounds [81 x i8], [81 x i8]* %3, i64 0, i64 9
br label %42
38: ; preds = %27
%39 = load i8*, i8** %29
%40 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([23 x i8], [23 x i8]* @.str.2, i64 0, i64 0), i8* %39)
%41 = tail call i32 @exit(i32 -1)
unreachable
42: ; preds = %36, %75
%43 = call i32 @fread(i8* %8, i32 1, i32 80, i32* %31)
%44 = icmp eq i32 %43, 80
br i1 %44, label %47, label %45
45: ; preds = %42
%46 = call i32 @exit(i32 -1)
unreachable
47: ; preds = %42
%48 = call i64 @strncmp(i8* getelementptr inbounds ([9 x i8], [9 x i8]* @.str.3, i64 0, i64 0), i8* %8, i32 8)
%49 = icmp eq i64 %48, 0
br i1 %49, label %78, label %50
50: ; preds = %47
%51 = call i64 @strncmp(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str.4, i64 0, i64 0), i8* %8, i32 9)
%52 = icmp eq i64 %51, 0
br i1 %52, label %72, label %53
53: ; preds = %50
%54 = call i64 @strncmp(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str.6, i64 0, i64 0), i8* %8, i32 9)
%55 = icmp eq i64 %54, 0
br i1 %55, label %72, label %56
56: ; preds = %53
%57 = call i64 @strncmp(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str.7, i64 0, i64 0), i8* %8, i32 9)
%58 = icmp eq i64 %57, 0
br i1 %58, label %72, label %59
59: ; preds = %56
%60 = call i64 @strncmp(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str.8, i64 0, i64 0), i8* %8, i32 9)
%61 = icmp eq i64 %60, 0
br i1 %61, label %72, label %62
62: ; preds = %59
%63 = call i64 @strncmp(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str.9, i64 0, i64 0), i8* %8, i32 9)
%64 = icmp eq i64 %63, 0
br i1 %64, label %65, label %75
65: ; preds = %62
%66 = call i32 (i8*, i8*, ...) @sscanf(i8* %37, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.10, i64 0, i64 0), i8* %12)
%67 = call i64 @strncmp(i8* %12, i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str.11, i64 0, i64 0), i32 12)
%68 = icmp eq i64 %67, 0
br i1 %68, label %75, label %69
69: ; preds = %65
%70 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([45 x i8], [45 x i8]* @.str.12, i64 0, i64 0), i8* %12)
%71 = call i32 @exit(i32 -1)
unreachable
72: ; preds = %59, %56, %53, %50
%73 = phi double* [ %4, %50 ], [ %5, %53 ], [ %6, %56 ], [ %7, %59 ]
%74 = call i32 (i8*, i8*, ...) @sscanf(i8* %37, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.5, i64 0, i64 0), double* %73)
br label %75
75: ; preds = %72, %62, %65
%76 = call i32 @feof(i32* %31)
%77 = icmp eq i32 %76, 0
br i1 %77, label %42, label %78
78: ; preds = %75, %47
%79 = load double, double* %6
%80 = load double, double* %7
%81 = fcmp ogt double %79, -3.600000e+02
%82 = fcmp ogt double %80, -3.600000e+02
%83 = select i1 %81, i1 %82, i1 false
br i1 %83, label %84, label %95
84: ; preds = %78
%85 = fadd double %79, %80
%86 = call double @llvm.fmuladd.f64(double %85, double 5.000000e-01, double -9.000000e+01)
%87 = fdiv double %86, 1.800000e+02
%88 = load double, double* @M_PI
%89 = fmul double %87, %88
%90 = load double, double* %4
%91 = fcmp ogt double %90, 0.000000e+00
%92 = load double, double* %5
%93 = fcmp ogt double %92, 0.000000e+00
%94 = select i1 %91, i1 %93, i1 false
br i1 %94, label %97, label %145
95: ; preds = %33, %78
%96 = call i32 @exit(i32 1)
unreachable
97: ; preds = %84
%98 = fadd double %90, %92
%99 = fmul double %98, 5.000000e-01
%100 = fdiv double %99, 1.800000e+02
%101 = fmul double %88, %100
%102 = call i32 @cos(double %101)
%103 = add i32 %102, 1
%104 = icmp ult i32 %103, 3
%105 = select i1 %104, i32 %102, i32 0
%106 = call i32 @cos(double %101)
%107 = sdiv i32 %105, %106
%108 = sitofp i32 %107 to double
%109 = call double @sin(double %101)
%110 = call double @sin(double %101)
%111 = fmul double %109, %110
%112 = fmul double %111, %108
%113 = fmul double %112, 3.975000e-09
%114 = fadd double %113, -6.270000e-07
%115 = call double @llvm.fmuladd.f64(double %112, double %114, double 1.733000e-04)
%116 = fmul double %115, %108
%117 = fsub double 1.000000e+00, %116
%118 = call i32 @cos(double %89)
%119 = sitofp i32 %118 to double
%120 = call double @sin(double %89)
%121 = call i32 @cos(double %28)
%122 = sitofp i32 %121 to double
%123 = call double @sin(double %28)
%124 = fmul double %120, %117
%125 = fmul double %120, %124
%126 = call double @llvm.fmuladd.f64(double %119, double %119, double %125)
%127 = fmul double %123, %116
%128 = fmul double %120, %127
%129 = fneg double %128
%130 = fmul double %119, %129
%131 = call double @llvm.fmuladd.f64(double %122, double %126, double %130)
%132 = fmul double %116, %122
%133 = fmul double %120, %132
%134 = fmul double %117, %119
%135 = fmul double %120, %120
%136 = call double @llvm.fmuladd.f64(double %134, double %119, double %135)
%137 = fneg double %123
%138 = fmul double %136, %137
%139 = call double @llvm.fmuladd.f64(double %133, double %119, double %138)
%140 = fmul double %133, %119
%141 = call double @llvm.fmuladd.f64(double %123, double %126, double %140)
%142 = fmul double %136, %122
%143 = call double @llvm.fmuladd.f64(double %128, double %119, double %142)
%144 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([33 x i8], [33 x i8]* @.str.13, i64 0, i64 0), double %131, double %139, double %141, double %143)
ret i32 0
145: ; preds = %84
%146 = call i32 @exit(i32 1)
unreachable
}