#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*/ fgets(char *, int, int /*<<< orphan*/);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *);
int /*<<< orphan*/ printf(char *, ...);
double sqrt(float);
int sscanf(char *, char *, float *);
int /*<<< orphan*/ stderr;
int /*<<< orphan*/ stdin;
int main(int argc, char *argv[]) {
  float texp;
  float npix;
  float rsrc;
  float rsky;
  float rtherm;
  float rdc;
  float Nsrc;
  float Nsky;
  float Nbkgd;
  float Ndc;
  float sn;
  float sn2;
  float noise;
  float tmp;
  float tmpa;
  float tmpb;
  float tmpc;
  char line[100];
  printf("\nCount rate from source in counts/sec\n");
  printf("  Get this value from wfc3tp.pro in IDL/destiny directory.\n");
  printf("  Enter value here:  ");
  fgets(line, 100, stdin);
  while (sscanf(line, "%f", &rsrc) != 1) {
    fprintf(stderr, "ERROR: bad input\n");
    fprintf(stderr, "Enter count rate again:  ");
    fgets(line, 100, stdin);
  }
  printf("Enter number of pixels covered by source:  ");
  fgets(line, 100, stdin);
  while (sscanf(line, "%f", &npix) != 1) {
    fprintf(stderr, "ERROR: bad input\n");
    fprintf(stderr, "Enter number of pixels again:  ");
    fgets(line, 100, stdin);
  }
  printf("Enter desired SNR:  ");
  fgets(line, 100, stdin);
  while (sscanf(line, "%f", &sn) != 1) {
    fprintf(stderr, "ERROR: bad input\n");
    fprintf(stderr, "Enter SNR again:  ");
    fgets(line, 100, stdin);
  }
  rsky = 0.61 * 0.449 * 0.023 * 5060.0 * npix * 0.13 * 0.13;
  rdc = 0.12 * npix;
  rtherm = 0.449 * 0.1 * npix;
  tmp = rsrc + rsky + rdc + rtherm;
  sn2 = sn * sn;
  tmpa = rsrc * rsrc;
  tmpb = -1.0 * sn2 * tmp;
  tmpc = -1.0 * sn2 * npix * 16.0 * 16.0;
  texp = (-1.0 * tmpb + sqrt(tmpb * tmpb - 4.0 * tmpa * tmpc)) / (2.0 * tmpa);
  noise = sqrt(tmp * texp + npix * 16.0 * 16.0);
  printf("\nFor this observation the input parameters are:\n");
  printf(" Source count rate     = %g e-/sec\n", rsrc);
  printf(" Sky count rate        = %g e-/sec\n", rsky);
  printf("    Sky                = %f phot/sec/Ang/arcsec^2\n", 0.023);
  printf("    QE                 = %f\n", 0.449);
  printf("    Filter bandwidth   = %f Ang\n", 5060.0);
  printf("    Pixel size         = %f arcsec\n", 0.13);
  printf("    N_pix in source    = %f\n", npix);
  printf(" Thermal count rate    = %f e-/sec\n", rtherm);
  printf(" Dark current rate     = %f e-/sec\n", rdc);
  printf(" Read noise            = %f e-\n", 16.0);
  printf("Giving:\n");
  printf("  Source counts     = %g cts\n", rsrc * texp);
  printf("  Sky counts        = %g cts\n", rsky * texp);
  printf("  Thermal counts    = %g cts\n", rtherm * texp);
  printf("  Dark counts       = %g cts\n", rdc * texp);
  printf("Associated noise\n");
  printf("  source            = %g cts\n", sqrt(rsrc * texp));
  printf("  sky               = %g cts\n", sqrt(rsky * texp));
  printf("  thermal           = %g cts\n", sqrt(rtherm * texp));
  printf("  dark current      = %g cts\n", sqrt(rdc * texp));
  printf("  read noise        = %g cts\n", 16.0 * sqrt(npix));
  printf(" Total noise        = %g cts\n\n", noise);
  printf("Which means you need an exposure time of %6.0f sec to get\n", texp);
  printf(" a SNR of %6.2f\n", sn);
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0xd0,%rsp
movl   $0x0,-0x4(%rbp)
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
lea    0xee8(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
lea    0xf01(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
lea    0xf2f(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
lea    -0xc0(%rbp),%rdi
mov    0x2e47(%rip),%rax
mov    (%rax),%edx
mov    $0x64,%esi
call   1040 <fgets@plt>
lea    -0xc0(%rbp),%rdi
lea    0xf16(%rip),%rsi
lea    -0x1c(%rbp),%rdx
call   1060 <sscanf@plt>
cmp    $0x1,%eax
je     1216 <main+0xc6>
mov    0x2e04(%rip),%rax
mov    (%rax),%edi
lea    0xef7(%rip),%rsi
call   1050 <fprintf@plt>
mov    0x2def(%rip),%rax
mov    (%rax),%edi
lea    0xef4(%rip),%rsi
call   1050 <fprintf@plt>
lea    -0xc0(%rbp),%rdi
mov    0x2de3(%rip),%rax
mov    (%rax),%edx
mov    $0x64,%esi
call   1040 <fgets@plt>
jmp    11ad <main+0x5d>
lea    0xee3(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
lea    -0xc0(%rbp),%rdi
mov    0x2db6(%rip),%rax
mov    (%rax),%edx
mov    $0x64,%esi
call   1040 <fgets@plt>
lea    -0xc0(%rbp),%rdi
lea    0xe85(%rip),%rsi
lea    -0x18(%rbp),%rdx
call   1060 <sscanf@plt>
cmp    $0x1,%eax
je     12a7 <main+0x157>
mov    0x2d73(%rip),%rax
mov    (%rax),%edi
lea    0xe66(%rip),%rsi
call   1050 <fprintf@plt>
mov    0x2d5e(%rip),%rax
mov    (%rax),%edi
lea    0xea9(%rip),%rsi
call   1050 <fprintf@plt>
lea    -0xc0(%rbp),%rdi
mov    0x2d52(%rip),%rax
mov    (%rax),%edx
mov    $0x64,%esi
call   1040 <fgets@plt>
jmp    123e <main+0xee>
lea    0xe9e(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
lea    -0xc0(%rbp),%rdi
mov    0x2d25(%rip),%rax
mov    (%rax),%edx
mov    $0x64,%esi
call   1040 <fgets@plt>
lea    -0xc0(%rbp),%rdi
lea    0xdf4(%rip),%rsi
lea    -0x3c(%rbp),%rdx
call   1060 <sscanf@plt>
cmp    $0x1,%eax
je     1338 <main+0x1e8>
mov    0x2ce2(%rip),%rax
mov    (%rax),%edi
lea    0xdd5(%rip),%rsi
call   1050 <fprintf@plt>
mov    0x2ccd(%rip),%rax
mov    (%rax),%edi
lea    0xe4d(%rip),%rsi
call   1050 <fprintf@plt>
lea    -0xc0(%rbp),%rdi
mov    0x2cc1(%rip),%rax
mov    (%rax),%edx
mov    $0x64,%esi
call   1040 <fgets@plt>
jmp    12cf <main+0x17f>
movss  -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd07(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  0xcb3(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xca7(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xcc9(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xca3(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x1c(%rbp),%xmm0
addss  -0x20(%rbp),%xmm0
addss  -0x28(%rbp),%xmm0
addss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x3c(%rbp),%xmm0
mulss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x1c(%rbp),%xmm0
mulss  -0x1c(%rbp),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x40(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xc46(%rip),%xmm0
mulsd  %xmm1,%xmm0
movss  -0x48(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x40(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xc1b(%rip),%xmm0
mulsd  %xmm1,%xmm0
movss  -0x18(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xbca(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xbbe(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x50(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movss  -0x50(%rbp),%xmm0
mulss  -0x50(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movss  -0x4c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm2
movsd  0xbb1(%rip),%xmm0
mulsd  %xmm2,%xmm0
movss  -0x54(%rbp),%xmm2
cvtss2sd %xmm2,%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1070 <sqrt@plt>
movsd  -0xc8(%rbp),%xmm1
ff
subsd  %xmm1,%xmm0
movss  -0x4c(%rbp),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  0xb57(%rip),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x48(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movss  -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xaff(%rip),%xmm2
mulsd  %xmm2,%xmm0
movsd  0xaf3(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1070 <sqrt@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
lea    0xc46(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movss  -0x1c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xc60(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xc6d(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
lea    0xc83(%rip),%rdi
movsd  0xaae(%rip),%xmm0
mov    $0x1,%al
call   1030 <printf@plt>
lea    0xca0(%rip),%rdi
movsd  0xa90(%rip),%xmm0
mov    $0x1,%al
call   1030 <printf@plt>
lea    0xca7(%rip),%rdi
movsd  0xa72(%rip),%xmm0
mov    $0x1,%al
call   1030 <printf@plt>
lea    0xcb2(%rip),%rdi
movsd  0xa54(%rip),%xmm0
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xcb7(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xcbd(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x28(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xcca(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
lea    0xce0(%rip),%rdi
movsd  0x9f1(%rip),%xmm0
mov    $0x1,%al
call   1030 <printf@plt>
lea    0xcea(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movss  -0x1c(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xcd7(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x20(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xcd9(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x24(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xcdb(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x28(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xcdd(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
lea    0xced(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movss  -0x1c(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
call   1070 <sqrt@plt>
lea    0xce2(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x20(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
call   1070 <sqrt@plt>
lea    0xce3(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x24(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
call   1070 <sqrt@plt>
lea    0xce4(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x28(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
call   1070 <sqrt@plt>
lea    0xce5(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x18(%rbp),%xmm0
call   1070 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0x8d5(%rip),%xmm0
mulsd  %xmm1,%xmm0
lea    0xcdc(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x44(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xce3(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xceb(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x3c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xd0f(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
xor    %eax,%eax
add    $0xd0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x98,%rsp
lea    0x127c(%rip),%rdi
call   1030 <puts@plt>
lea    0x1296(%rip),%rdi
call   1030 <puts@plt>
lea    0xed0(%rip),%rdi
xor    %eax,%eax
call   1040 <printf@plt>
mov    0x2e4a(%rip),%rax
mov    (%rax),%edx
lea    0x30(%rsp),%rbx
mov    %rbx,%rdi
mov    $0x64,%esi
call   1050 <fgets@plt>
lea    0xebd(%rip),%rsi
lea    0x4(%rsp),%rdx
mov    %rbx,%rdi
call   1070 <sscanf@plt>
cmp    $0x1,%eax
je     123f <main+0xdf>
mov    0x2e06(%rip),%rbp
lea    0xea0(%rip),%r14
lea    0xeab(%rip),%r15
lea    0x30(%rsp),%rbx
lea    0xe8a(%rip),%r12
lea    0x4(%rsp),%r13
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    0x0(%rbp),%edi
mov    %r14,%rsi
call   1060 <fprintf@plt>
mov    0x0(%rbp),%edi
mov    %r15,%rsi
call   1060 <fprintf@plt>
mov    0x2dcb(%rip),%rax
mov    (%rax),%edx
mov    %rbx,%rdi
mov    $0x64,%esi
call   1050 <fgets@plt>
mov    %rbx,%rdi
mov    %r12,%rsi
mov    %r13,%rdx
call   1070 <sscanf@plt>
cmp    $0x1,%eax
jne    1200 <main+0xa0>
lea    0xe5f(%rip),%rdi
xor    %eax,%eax
call   1040 <printf@plt>
mov    0x2d94(%rip),%rax
mov    (%rax),%edx
lea    0x30(%rsp),%rbx
mov    %rbx,%rdi
mov    $0x64,%esi
call   1050 <fgets@plt>
lea    0xe07(%rip),%rsi
lea    0x8(%rsp),%rdx
mov    %rbx,%rdi
call   1070 <sscanf@plt>
cmp    $0x1,%eax
je     12ef <main+0x18f>
mov    0x2d50(%rip),%rbp
lea    0xdea(%rip),%r14
lea    0xe3b(%rip),%r15
lea    0x30(%rsp),%rbx
lea    0xdd4(%rip),%r12
lea    0x8(%rsp),%r13
nopw   0x0(%rax,%rax,1)
mov    0x0(%rbp),%edi
mov    %r14,%rsi
call   1060 <fprintf@plt>
mov    0x0(%rbp),%edi
mov    %r15,%rsi
call   1060 <fprintf@plt>
mov    0x2d1b(%rip),%rax
mov    (%rax),%edx
mov    %rbx,%rdi
mov    $0x64,%esi
call   1050 <fgets@plt>
mov    %rbx,%rdi
mov    %r12,%rsi
mov    %r13,%rdx
call   1070 <sscanf@plt>
cmp    $0x1,%eax
jne    12b0 <main+0x150>
lea    0xdfb(%rip),%rdi
xor    %eax,%eax
call   1040 <printf@plt>
mov    0x2ce4(%rip),%rax
mov    (%rax),%edx
lea    0x30(%rsp),%rbx
mov    %rbx,%rdi
mov    $0x64,%esi
call   1050 <fgets@plt>
lea    0xd57(%rip),%rsi
lea    0x24(%rsp),%rdx
mov    %rbx,%rdi
call   1070 <sscanf@plt>
cmp    $0x1,%eax
je     139f <main+0x23f>
mov    0x2ca0(%rip),%rbp
lea    0xd3a(%rip),%r14
lea    0xdc0(%rip),%r15
lea    0x30(%rsp),%rbx
lea    0xd24(%rip),%r12
lea    0x24(%rsp),%r13
nopw   0x0(%rax,%rax,1)
mov    0x0(%rbp),%edi
mov    %r14,%rsi
call   1060 <fprintf@plt>
mov    0x0(%rbp),%edi
mov    %r15,%rsi
call   1060 <fprintf@plt>
mov    0x2c6b(%rip),%rax
mov    (%rax),%edx
mov    %rbx,%rdi
mov    $0x64,%esi
call   1050 <fgets@plt>
mov    %rbx,%rdi
mov    %r12,%rsi
mov    %r13,%rdx
call   1070 <sscanf@plt>
cmp    $0x1,%eax
jne    1360 <main+0x200>
movss  0x8(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xc4f(%rip),%xmm1
mulsd  %xmm0,%xmm1
movsd  0xc4b(%rip),%xmm2
mulsd  %xmm2,%xmm1
mulsd  %xmm2,%xmm1
cvtsd2ss %xmm1,%xmm4
movss  %xmm4,0xc(%rsp)
movsd  0xc39(%rip),%xmm1
mulsd  %xmm0,%xmm1
xorps  %xmm2,%xmm2
cvtsd2ss %xmm1,%xmm2
movss  %xmm2,0x14(%rsp)
movsd  0xc28(%rip),%xmm1
mulsd  %xmm0,%xmm1
cvtsd2ss %xmm1,%xmm3
movss  %xmm3,0x10(%rsp)
movss  0x4(%rsp),%xmm1
movaps %xmm1,%xmm5
addss  %xmm4,%xmm5
addss  %xmm2,%xmm5
addss  %xmm3,%xmm5
movss  %xmm5,0x20(%rsp)
movss  0x24(%rsp),%xmm2
mulss  %xmm2,%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm2,%xmm0
xorpd  0xc1d(%rip),%xmm2
xorps  %xmm3,%xmm3
cvtss2sd %xmm5,%xmm3
mulsd  %xmm2,%xmm3
xorps  %xmm2,%xmm2
cvtsd2ss %xmm3,%xmm2
movsd  0xbd3(%rip),%xmm3
mulsd  %xmm3,%xmm0
mulsd  %xmm3,%xmm0
cvtsd2ss %xmm0,%xmm0
mulss  %xmm1,%xmm1
xorps  %xmm3,%xmm3
cvtss2sd %xmm2,%xmm3
movsd  %xmm3,0x18(%rsp)
mulss  %xmm2,%xmm2
cvtss2sd %xmm2,%xmm2
xorps  %xmm3,%xmm3
cvtss2sd %xmm1,%xmm3
movsd  %xmm3,0x28(%rsp)
cvtss2sd %xmm0,%xmm0
movsd  0xb9d(%rip),%xmm1
mulsd  %xmm3,%xmm1
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
call   1080 <sqrt@plt>
subsd  0x18(%rsp),%xmm0
movsd  0x28(%rsp),%xmm1
addsd  %xmm1,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x18(%rsp)
movss  0x20(%rsp),%xmm1
mulss  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
movss  0x8(%rsp),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xb3c(%rip),%xmm2
mulsd  %xmm2,%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
call   1080 <sqrt@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x20(%rsp)
lea    0xf48(%rip),%rdi
call   1030 <puts@plt>
movss  0x4(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xbf6(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movss  0xc(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xc02(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
lea    0xc18(%rip),%rdi
movsd  0xadf(%rip),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
lea    0xc35(%rip),%rdi
movsd  0xad1(%rip),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
lea    0xc3c(%rip),%rdi
movsd  0xac3(%rip),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
lea    0xc47(%rip),%rdi
movsd  0xa75(%rip),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
movss  0x8(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xc4b(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movss  0x10(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xc50(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movss  0x14(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xc5c(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
lea    0xc72(%rip),%rdi
movsd  0xa2f(%rip),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
lea    0xe86(%rip),%rdi
call   1030 <puts@plt>
movss  0x4(%rsp),%xmm0
mulss  0x18(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xc60(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movss  0xc(%rsp),%xmm0
mulss  0x18(%rsp),%xmm0
movss  %xmm0,0xc(%rsp)
cvtss2sd %xmm0,%xmm0
lea    0xc5a(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movss  0x10(%rsp),%xmm0
mulss  0x18(%rsp),%xmm0
movss  %xmm0,0x10(%rsp)
cvtss2sd %xmm0,%xmm0
lea    0xc54(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movss  0x14(%rsp),%xmm0
mulss  0x18(%rsp),%xmm0
movss  %xmm0,0x14(%rsp)
cvtss2sd %xmm0,%xmm0
lea    0xc4e(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
lea    0xdf8(%rip),%rdi
call   1030 <puts@plt>
movss  0x4(%rsp),%xmm0
mulss  0x18(%rsp),%xmm0
call   1080 <sqrt@plt>
lea    0xc41(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movss  0xc(%rsp),%xmm0
call   1080 <sqrt@plt>
lea    0xc46(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movss  0x10(%rsp),%xmm0
call   1080 <sqrt@plt>
lea    0xc4b(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movss  0x14(%rsp),%xmm0
call   1080 <sqrt@plt>
lea    0xc50(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movss  0x8(%rsp),%xmm0
call   1080 <sqrt@plt>
mulsd  0x909(%rip),%xmm0
lea    0xc4d(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movss  0x20(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xc53(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movss  0x18(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xc5a(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movss  0x24(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xc7d(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
xor    %eax,%eax
add    $0x98,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [39 x i8] c"\0ACount rate from source in counts/sec\0A\00"
@.str.1 = [60 x i8] c" Get this value from wfc3tp.pro in IDL/destiny directory.\0A\00"
@.str.2 = [22 x i8] c" Enter value here: \00"
@stdin = global i32 0
@.str.3 = [3 x i8] c"%f\00"
@stderr = global i32 0
@.str.4 = [18 x i8] c"ERROR: bad input\0A\00"
@.str.5 = [26 x i8] c"Enter count rate again: \00"
@.str.6 = [44 x i8] c"Enter number of pixels covered by source: \00"
@.str.7 = [32 x i8] c"Enter number of pixels again: \00"
@.str.8 = [21 x i8] c"Enter desired SNR: \00"
@.str.9 = [19 x i8] c"Enter SNR again: \00"
@.str.10 = [49 x i8] c"\0AFor this observation the input parameters are:\0A\00"
@.str.11 = [36 x i8] c" Source count rate = %g e-/sec\0A\00"
@.str.12 = [36 x i8] c" Sky count rate = %g e-/sec\0A\00"
@.str.13 = [51 x i8] c" Sky = %f phot/sec/Ang/arcsec^2\0A\00"
@.str.14 = [29 x i8] c" QE = %f\0A\00"
@.str.15 = [33 x i8] c" Filter bandwidth = %f Ang\0A\00"
@.str.16 = [36 x i8] c" Pixel size = %f arcsec\0A\00"
@.str.17 = [29 x i8] c" N_pix in source = %f\0A\00"
@.str.18 = [36 x i8] c" Thermal count rate = %f e-/sec\0A\00"
@.str.19 = [36 x i8] c" Dark current rate = %f e-/sec\0A\00"
@.str.20 = [32 x i8] c" Read noise = %f e-\0A\00"
@.str.21 = [9 x i8] c"Giving:\0A\00"
@.str.22 = [30 x i8] c" Source counts = %g cts\0A\00"
@.str.23 = [30 x i8] c" Sky counts = %g cts\0A\00"
@.str.24 = [30 x i8] c" Thermal counts = %g cts\0A\00"
@.str.25 = [30 x i8] c" Dark counts = %g cts\0A\00"
@.str.26 = [18 x i8] c"Associated noise\0A\00"
@.str.27 = [30 x i8] c" source = %g cts\0A\00"
@.str.28 = [30 x i8] c" sky = %g cts\0A\00"
@.str.29 = [30 x i8] c" thermal = %g cts\0A\00"
@.str.30 = [30 x i8] c" dark current = %g cts\0A\00"
@.str.31 = [30 x i8] c" read noise = %g cts\0A\00"
@.str.32 = [31 x i8] c" Total noise = %g cts\0A\0A\00"
@.str.33 = [59 x i8] c"Which means you need an exposure time of %6.0f sec to get\0A\00"
@.str.34 = [17 x i8] c" a SNR of %6.2f\0A\00"
define i32 @main(i32 %0, i8** %1) {
%3 = alloca i32
%4 = alloca i32
%5 = alloca i8**
%6 = alloca float
%7 = alloca float
%8 = alloca float
%9 = alloca float
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca float
%14 = alloca float
%15 = alloca float
%16 = alloca float
%17 = alloca float
%18 = alloca float
%19 = alloca float
%20 = alloca float
%21 = alloca float
%22 = alloca float
%23 = alloca [100 x i8]
store i32 0, i32* %3
store i32 %0, i32* %4
store i8** %1, i8*** %5
%24 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([39 x i8], [39 x i8]* @.str, i64 0, i64 0))
%25 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([60 x i8], [60 x i8]* @.str.1, i64 0, i64 0))
%26 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([22 x i8], [22 x i8]* @.str.2, i64 0, i64 0))
%27 = getelementptr inbounds [100 x i8], [100 x i8]* %23, i64 0, i64 0
%28 = load i32, i32* @stdin
%29 = call i32 @fgets(i8* %27, i32 100, i32 %28)
br label %30
30: ; preds = %34, %2
%31 = getelementptr inbounds [100 x i8], [100 x i8]* %23, i64 0, i64 0
%32 = call i32 @sscanf(i8* %31, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.3, i64 0, i64 0), float* %8)
%33 = icmp ne i32 %32, 1
br i1 %33, label %34, label %42
34: ; preds = %30
%35 = load i32, i32* @stderr
%36 = call i32 @fprintf(i32 %35, i8* getelementptr inbounds ([18 x i8], [18 x i8]* @.str.4, i64 0, i64 0))
%37 = load i32, i32* @stderr
%38 = call i32 @fprintf(i32 %37, i8* getelementptr inbounds ([26 x i8], [26 x i8]* @.str.5, i64 0, i64 0))
%39 = getelementptr inbounds [100 x i8], [100 x i8]* %23, i64 0, i64 0
%40 = load i32, i32* @stdin
%41 = call i32 @fgets(i8* %39, i32 100, i32 %40)
br label %30
42: ; preds = %30
%43 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([44 x i8], [44 x i8]* @.str.6, i64 0, i64 0))
%44 = getelementptr inbounds [100 x i8], [100 x i8]* %23, i64 0, i64 0
%45 = load i32, i32* @stdin
%46 = call i32 @fgets(i8* %44, i32 100, i32 %45)
br label %47
47: ; preds = %51, %42
%48 = getelementptr inbounds [100 x i8], [100 x i8]* %23, i64 0, i64 0
%49 = call i32 @sscanf(i8* %48, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.3, i64 0, i64 0), float* %7)
%50 = icmp ne i32 %49, 1
br i1 %50, label %51, label %59
51: ; preds = %47
%52 = load i32, i32* @stderr
%53 = call i32 @fprintf(i32 %52, i8* getelementptr inbounds ([18 x i8], [18 x i8]* @.str.4, i64 0, i64 0))
%54 = load i32, i32* @stderr
%55 = call i32 @fprintf(i32 %54, i8* getelementptr inbounds ([32 x i8], [32 x i8]* @.str.7, i64 0, i64 0))
%56 = getelementptr inbounds [100 x i8], [100 x i8]* %23, i64 0, i64 0
%57 = load i32, i32* @stdin
%58 = call i32 @fgets(i8* %56, i32 100, i32 %57)
br label %47
59: ; preds = %47
%60 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([21 x i8], [21 x i8]* @.str.8, i64 0, i64 0))
%61 = getelementptr inbounds [100 x i8], [100 x i8]* %23, i64 0, i64 0
%62 = load i32, i32* @stdin
%63 = call i32 @fgets(i8* %61, i32 100, i32 %62)
br label %64
64: ; preds = %68, %59
%65 = getelementptr inbounds [100 x i8], [100 x i8]* %23, i64 0, i64 0
%66 = call i32 @sscanf(i8* %65, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.3, i64 0, i64 0), float* %16)
%67 = icmp ne i32 %66, 1
br i1 %67, label %68, label %76
68: ; preds = %64
%69 = load i32, i32* @stderr
%70 = call i32 @fprintf(i32 %69, i8* getelementptr inbounds ([18 x i8], [18 x i8]* @.str.4, i64 0, i64 0))
%71 = load i32, i32* @stderr
%72 = call i32 @fprintf(i32 %71, i8* getelementptr inbounds ([19 x i8], [19 x i8]* @.str.9, i64 0, i64 0))
%73 = getelementptr inbounds [100 x i8], [100 x i8]* %23, i64 0, i64 0
%74 = load i32, i32* @stdin
%75 = call i32 @fgets(i8* %73, i32 100, i32 %74)
br label %64
76: ; preds = %64
%77 = load float, float* %7
%78 = fpext float %77 to double
%79 = fmul double 0x403FE014DA8297F6, %78
%80 = fmul double %79, 1.300000e-01
%81 = fmul double %80, 1.300000e-01
%82 = fptrunc double %81 to float
store float %82, float* %9
%83 = load float, float* %7
%84 = fpext float %83 to double
%85 = fmul double 1.200000e-01, %84
%86 = fptrunc double %85 to float
store float %86, float* %11
%87 = load float, float* %7
%88 = fpext float %87 to double
%89 = fmul double 4.490000e-02, %88
%90 = fptrunc double %89 to float
store float %90, float* %10
%91 = load float, float* %8
%92 = load float, float* %9
%93 = fadd float %91, %92
%94 = load float, float* %11
%95 = fadd float %93, %94
%96 = load float, float* %10
%97 = fadd float %95, %96
store float %97, float* %19
%98 = load float, float* %16
%99 = load float, float* %16
%100 = fmul float %98, %99
store float %100, float* %17
%101 = load float, float* %8
%102 = load float, float* %8
%103 = fmul float %101, %102
store float %103, float* %20
%104 = load float, float* %17
%105 = fpext float %104 to double
%106 = fmul double -1.000000e+00, %105
%107 = load float, float* %19
%108 = fpext float %107 to double
%109 = fmul double %106, %108
%110 = fptrunc double %109 to float
store float %110, float* %21
%111 = load float, float* %17
%112 = fpext float %111 to double
%113 = fmul double -1.000000e+00, %112
%114 = load float, float* %7
%115 = fpext float %114 to double
%116 = fmul double %113, %115
%117 = fmul double %116, 1.600000e+01
%118 = fmul double %117, 1.600000e+01
%119 = fptrunc double %118 to float
store float %119, float* %22
%120 = load float, float* %21
%121 = fpext float %120 to double
%122 = load float, float* %21
%123 = load float, float* %21
%124 = fmul float %122, %123
%125 = fpext float %124 to double
%126 = load float, float* %20
%127 = fpext float %126 to double
%128 = fmul double 4.000000e+00, %127
%129 = load float, float* %22
%130 = fpext float %129 to double
%131 = fneg double %128
%132 = call double @llvm.fmuladd.f64(double %131, double %130, double %125)
%133 = fptrunc double %132 to float
%134 = call double @sqrt(float %133)
%135 = call double @llvm.fmuladd.f64(double -1.000000e+00, double %121, double %134)
%136 = load float, float* %20
%137 = fpext float %136 to double
%138 = fmul double 2.000000e+00, %137
%139 = fdiv double %135, %138
%140 = fptrunc double %139 to float
store float %140, float* %6
%141 = load float, float* %19
%142 = load float, float* %6
%143 = fmul float %141, %142
%144 = fpext float %143 to double
%145 = load float, float* %7
%146 = fpext float %145 to double
%147 = fmul double %146, 1.600000e+01
%148 = call double @llvm.fmuladd.f64(double %147, double 1.600000e+01, double %144)
%149 = fptrunc double %148 to float
%150 = call double @sqrt(float %149)
%151 = fptrunc double %150 to float
store float %151, float* %18
%152 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([49 x i8], [49 x i8]* @.str.10, i64 0, i64 0))
%153 = load float, float* %8
%154 = fpext float %153 to double
%155 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([36 x i8], [36 x i8]* @.str.11, i64 0, i64 0), double %154)
%156 = load float, float* %9
%157 = fpext float %156 to double
%158 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([36 x i8], [36 x i8]* @.str.12, i64 0, i64 0), double %157)
%159 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([51 x i8], [51 x i8]* @.str.13, i64 0, i64 0), double 2.300000e-02)
%160 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([29 x i8], [29 x i8]* @.str.14, i64 0, i64 0), double 4.490000e-01)
%161 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([33 x i8], [33 x i8]* @.str.15, i64 0, i64 0), double 5.060000e+03)
%162 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([36 x i8], [36 x i8]* @.str.16, i64 0, i64 0), double 1.300000e-01)
%163 = load float, float* %7
%164 = fpext float %163 to double
%165 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([29 x i8], [29 x i8]* @.str.17, i64 0, i64 0), double %164)
%166 = load float, float* %10
%167 = fpext float %166 to double
%168 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([36 x i8], [36 x i8]* @.str.18, i64 0, i64 0), double %167)
%169 = load float, float* %11
%170 = fpext float %169 to double
%171 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([36 x i8], [36 x i8]* @.str.19, i64 0, i64 0), double %170)
%172 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([32 x i8], [32 x i8]* @.str.20, i64 0, i64 0), double 1.600000e+01)
%173 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([9 x i8], [9 x i8]* @.str.21, i64 0, i64 0))
%174 = load float, float* %8
%175 = load float, float* %6
%176 = fmul float %174, %175
%177 = fpext float %176 to double
%178 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str.22, i64 0, i64 0), double %177)
%179 = load float, float* %9
%180 = load float, float* %6
%181 = fmul float %179, %180
%182 = fpext float %181 to double
%183 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str.23, i64 0, i64 0), double %182)
%184 = load float, float* %10
%185 = load float, float* %6
%186 = fmul float %184, %185
%187 = fpext float %186 to double
%188 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str.24, i64 0, i64 0), double %187)
%189 = load float, float* %11
%190 = load float, float* %6
%191 = fmul float %189, %190
%192 = fpext float %191 to double
%193 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str.25, i64 0, i64 0), double %192)
%194 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([18 x i8], [18 x i8]* @.str.26, i64 0, i64 0))
%195 = load float, float* %8
%196 = load float, float* %6
%197 = fmul float %195, %196
%198 = call double @sqrt(float %197)
%199 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str.27, i64 0, i64 0), double %198)
%200 = load float, float* %9
%201 = load float, float* %6
%202 = fmul float %200, %201
%203 = call double @sqrt(float %202)
%204 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str.28, i64 0, i64 0), double %203)
%205 = load float, float* %10
%206 = load float, float* %6
%207 = fmul float %205, %206
%208 = call double @sqrt(float %207)
%209 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str.29, i64 0, i64 0), double %208)
%210 = load float, float* %11
%211 = load float, float* %6
%212 = fmul float %210, %211
%213 = call double @sqrt(float %212)
%214 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str.30, i64 0, i64 0), double %213)
%215 = load float, float* %7
%216 = call double @sqrt(float %215)
%217 = fmul double 1.600000e+01, %216
%218 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str.31, i64 0, i64 0), double %217)
%219 = load float, float* %18
%220 = fpext float %219 to double
%221 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([31 x i8], [31 x i8]* @.str.32, i64 0, i64 0), double %220)
%222 = load float, float* %6
%223 = fpext float %222 to double
%224 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([59 x i8], [59 x i8]* @.str.33, i64 0, i64 0), double %223)
%225 = load float, float* %16
%226 = fpext float %225 to double
%227 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([17 x i8], [17 x i8]* @.str.34, i64 0, i64 0), double %226)
ret i32 0
}

LLVM IR (-O2):
@.str.2 = [22 x i8] c" Enter value here: \00"
@stdin = global i32 0
@.str.3 = [3 x i8] c"%f\00"
@stderr = global i32 0
@.str.4 = [18 x i8] c"ERROR: bad input\0A\00"
@.str.5 = [26 x i8] c"Enter count rate again: \00"
@.str.6 = [44 x i8] c"Enter number of pixels covered by source: \00"
@.str.7 = [32 x i8] c"Enter number of pixels again: \00"
@.str.8 = [21 x i8] c"Enter desired SNR: \00"
@.str.9 = [19 x i8] c"Enter SNR again: \00"
@.str.11 = [36 x i8] c" Source count rate = %g e-/sec\0A\00"
@.str.12 = [36 x i8] c" Sky count rate = %g e-/sec\0A\00"
@.str.13 = [51 x i8] c" Sky = %f phot/sec/Ang/arcsec^2\0A\00"
@.str.14 = [29 x i8] c" QE = %f\0A\00"
@.str.15 = [33 x i8] c" Filter bandwidth = %f Ang\0A\00"
@.str.16 = [36 x i8] c" Pixel size = %f arcsec\0A\00"
@.str.17 = [29 x i8] c" N_pix in source = %f\0A\00"
@.str.18 = [36 x i8] c" Thermal count rate = %f e-/sec\0A\00"
@.str.19 = [36 x i8] c" Dark current rate = %f e-/sec\0A\00"
@.str.20 = [32 x i8] c" Read noise = %f e-\0A\00"
@.str.22 = [30 x i8] c" Source counts = %g cts\0A\00"
@.str.23 = [30 x i8] c" Sky counts = %g cts\0A\00"
@.str.24 = [30 x i8] c" Thermal counts = %g cts\0A\00"
@.str.25 = [30 x i8] c" Dark counts = %g cts\0A\00"
@.str.27 = [30 x i8] c" source = %g cts\0A\00"
@.str.28 = [30 x i8] c" sky = %g cts\0A\00"
@.str.29 = [30 x i8] c" thermal = %g cts\0A\00"
@.str.30 = [30 x i8] c" dark current = %g cts\0A\00"
@.str.31 = [30 x i8] c" read noise = %g cts\0A\00"
@.str.32 = [31 x i8] c" Total noise = %g cts\0A\0A\00"
@.str.33 = [59 x i8] c"Which means you need an exposure time of %6.0f sec to get\0A\00"
@.str.34 = [17 x i8] c" a SNR of %6.2f\0A\00"
@str = [38 x i8] c"\0ACount rate from source in counts/sec\00"
@str.35 = [59 x i8] c" Get this value from wfc3tp.pro in IDL/destiny directory.\00"
@str.36 = [48 x i8] c"\0AFor this observation the input parameters are:\00"
@str.37 = [8 x i8] c"Giving:\00"
@str.38 = [17 x i8] c"Associated noise\00"
define i32 @main(i32 %0, i8** nocapture readnone %1) {
%3 = alloca float
%4 = alloca float
%5 = alloca float
%6 = alloca [100 x i8]
%7 = bitcast float* %3 to i8*
%8 = bitcast float* %4 to i8*
%9 = bitcast float* %5 to i8*
%10 = getelementptr inbounds [100 x i8], [100 x i8]* %6, i64 0, i64 0
%11 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([38 x i8], [38 x i8]* @str, i64 0, i64 0))
%12 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([59 x i8], [59 x i8]* @str.35, i64 0, i64 0))
%13 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([22 x i8], [22 x i8]* @.str.2, i64 0, i64 0))
%14 = load i32, i32* @stdin
%15 = call i32 @fgets(i8* %10, i32 100, i32 %14)
%16 = call i32 @sscanf(i8* %10, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.3, i64 0, i64 0), float* %4)
%17 = icmp eq i32 %16, 1
br i1 %17, label %27, label %18
18: ; preds = %2, %18
%19 = load i32, i32* @stderr
%20 = call i32 @fprintf(i32 %19, i8* getelementptr inbounds ([18 x i8], [18 x i8]* @.str.4, i64 0, i64 0))
%21 = load i32, i32* @stderr
%22 = call i32 @fprintf(i32 %21, i8* getelementptr inbounds ([26 x i8], [26 x i8]* @.str.5, i64 0, i64 0))
%23 = load i32, i32* @stdin
%24 = call i32 @fgets(i8* %10, i32 100, i32 %23)
%25 = call i32 @sscanf(i8* %10, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.3, i64 0, i64 0), float* %4)
%26 = icmp eq i32 %25, 1
br i1 %26, label %27, label %18
27: ; preds = %18, %2
%28 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([44 x i8], [44 x i8]* @.str.6, i64 0, i64 0))
%29 = load i32, i32* @stdin
%30 = call i32 @fgets(i8* %10, i32 100, i32 %29)
%31 = call i32 @sscanf(i8* %10, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.3, i64 0, i64 0), float* %3)
%32 = icmp eq i32 %31, 1
br i1 %32, label %42, label %33
33: ; preds = %27, %33
%34 = load i32, i32* @stderr
%35 = call i32 @fprintf(i32 %34, i8* getelementptr inbounds ([18 x i8], [18 x i8]* @.str.4, i64 0, i64 0))
%36 = load i32, i32* @stderr
%37 = call i32 @fprintf(i32 %36, i8* getelementptr inbounds ([32 x i8], [32 x i8]* @.str.7, i64 0, i64 0))
%38 = load i32, i32* @stdin
%39 = call i32 @fgets(i8* %10, i32 100, i32 %38)
%40 = call i32 @sscanf(i8* %10, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.3, i64 0, i64 0), float* %3)
%41 = icmp eq i32 %40, 1
br i1 %41, label %42, label %33
42: ; preds = %33, %27
%43 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([21 x i8], [21 x i8]* @.str.8, i64 0, i64 0))
%44 = load i32, i32* @stdin
%45 = call i32 @fgets(i8* %10, i32 100, i32 %44)
%46 = call i32 @sscanf(i8* %10, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.3, i64 0, i64 0), float* %5)
%47 = icmp eq i32 %46, 1
br i1 %47, label %57, label %48
48: ; preds = %42, %48
%49 = load i32, i32* @stderr
%50 = call i32 @fprintf(i32 %49, i8* getelementptr inbounds ([18 x i8], [18 x i8]* @.str.4, i64 0, i64 0))
%51 = load i32, i32* @stderr
%52 = call i32 @fprintf(i32 %51, i8* getelementptr inbounds ([19 x i8], [19 x i8]* @.str.9, i64 0, i64 0))
%53 = load i32, i32* @stdin
%54 = call i32 @fgets(i8* %10, i32 100, i32 %53)
%55 = call i32 @sscanf(i8* %10, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.3, i64 0, i64 0), float* %5)
%56 = icmp eq i32 %55, 1
br i1 %56, label %57, label %48
57: ; preds = %48, %42
%58 = load float, float* %3
%59 = fpext float %58 to double
%60 = fmul double %59, 0x403FE014DA8297F6
%61 = fmul double %60, 1.300000e-01
%62 = fmul double %61, 1.300000e-01
%63 = fptrunc double %62 to float
%64 = fmul double %59, 1.200000e-01
%65 = fptrunc double %64 to float
%66 = fmul double %59, 4.490000e-02
%67 = fptrunc double %66 to float
%68 = load float, float* %4
%69 = fadd float %68, %63
%70 = fadd float %69, %65
%71 = fadd float %70, %67
%72 = load float, float* %5
%73 = fmul float %72, %72
%74 = fmul float %68, %68
%75 = fpext float %73 to double
%76 = fneg double %75
%77 = fpext float %71 to double
%78 = fmul double %76, %77
%79 = fptrunc double %78 to float
%80 = fmul double %59, %76
%81 = fmul double %80, 1.600000e+01
%82 = fmul double %81, 1.600000e+01
%83 = fptrunc double %82 to float
%84 = fpext float %79 to double
%85 = fmul float %79, %79
%86 = fpext float %85 to double
%87 = fpext float %74 to double
%88 = fpext float %83 to double
%89 = fmul double %87, -4.000000e+00
%90 = call double @llvm.fmuladd.f64(double %89, double %88, double %86)
%91 = fptrunc double %90 to float
%92 = call double @sqrt(float %91)
%93 = call double @llvm.fmuladd.f64(double %84, double -1.000000e+00, double %92)
%94 = fmul double %87, 2.000000e+00
%95 = fdiv double %93, %94
%96 = fptrunc double %95 to float
%97 = fmul float %71, %96
%98 = fpext float %97 to double
%99 = load float, float* %3
%100 = fpext float %99 to double
%101 = fmul double %100, 1.600000e+01
%102 = call double @llvm.fmuladd.f64(double %101, double 1.600000e+01, double %98)
%103 = fptrunc double %102 to float
%104 = call double @sqrt(float %103)
%105 = fptrunc double %104 to float
%106 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([48 x i8], [48 x i8]* @str.36, i64 0, i64 0))
%107 = load float, float* %4
%108 = fpext float %107 to double
%109 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([36 x i8], [36 x i8]* @.str.11, i64 0, i64 0), double %108)
%110 = fpext float %63 to double
%111 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([36 x i8], [36 x i8]* @.str.12, i64 0, i64 0), double %110)
%112 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([51 x i8], [51 x i8]* @.str.13, i64 0, i64 0), double 2.300000e-02)
%113 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([29 x i8], [29 x i8]* @.str.14, i64 0, i64 0), double 4.490000e-01)
%114 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([33 x i8], [33 x i8]* @.str.15, i64 0, i64 0), double 5.060000e+03)
%115 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([36 x i8], [36 x i8]* @.str.16, i64 0, i64 0), double 1.300000e-01)
%116 = load float, float* %3
%117 = fpext float %116 to double
%118 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([29 x i8], [29 x i8]* @.str.17, i64 0, i64 0), double %117)
%119 = fpext float %67 to double
%120 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([36 x i8], [36 x i8]* @.str.18, i64 0, i64 0), double %119)
%121 = fpext float %65 to double
%122 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([36 x i8], [36 x i8]* @.str.19, i64 0, i64 0), double %121)
%123 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([32 x i8], [32 x i8]* @.str.20, i64 0, i64 0), double 1.600000e+01)
%124 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([8 x i8], [8 x i8]* @str.37, i64 0, i64 0))
%125 = load float, float* %4
%126 = fmul float %125, %96
%127 = fpext float %126 to double
%128 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([30 x i8], [30 x i8]* @.str.22, i64 0, i64 0), double %127)
%129 = fmul float %63, %96
%130 = fpext float %129 to double
%131 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([30 x i8], [30 x i8]* @.str.23, i64 0, i64 0), double %130)
%132 = fmul float %67, %96
%133 = fpext float %132 to double
%134 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([30 x i8], [30 x i8]* @.str.24, i64 0, i64 0), double %133)
%135 = fmul float %65, %96
%136 = fpext float %135 to double
%137 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([30 x i8], [30 x i8]* @.str.25, i64 0, i64 0), double %136)
%138 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([17 x i8], [17 x i8]* @str.38, i64 0, i64 0))
%139 = load float, float* %4
%140 = fmul float %139, %96
%141 = call double @sqrt(float %140)
%142 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([30 x i8], [30 x i8]* @.str.27, i64 0, i64 0), double %141)
%143 = call double @sqrt(float %129)
%144 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([30 x i8], [30 x i8]* @.str.28, i64 0, i64 0), double %143)
%145 = call double @sqrt(float %132)
%146 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([30 x i8], [30 x i8]* @.str.29, i64 0, i64 0), double %145)
%147 = call double @sqrt(float %135)
%148 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([30 x i8], [30 x i8]* @.str.30, i64 0, i64 0), double %147)
%149 = load float, float* %3
%150 = call double @sqrt(float %149)
%151 = fmul double %150, 1.600000e+01
%152 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([30 x i8], [30 x i8]* @.str.31, i64 0, i64 0), double %151)
%153 = fpext float %105 to double
%154 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([31 x i8], [31 x i8]* @.str.32, i64 0, i64 0), double %153)
%155 = fpext float %96 to double
%156 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([59 x i8], [59 x i8]* @.str.33, i64 0, i64 0), double %155)
%157 = load float, float* %5
%158 = fpext float %157 to double
%159 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([17 x i8], [17 x i8]* @.str.34, i64 0, i64 0), double %158)
ret i32 0
}