#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*/ exit(int);
double const exp(double);
double for_large_negative(double, double);
double log(double);
double pow(double, double);
int /*<<< orphan*/ printf(char *);
double sqrt(double);
double th_Gamma(double);
double th_Hyper1F1GG(double, double, double, double);
scalar_t__ th_Uax9Strip(double, double, double, double);
double ParCylU(double a, double x, double eps) {
  double aPh, aD2, Dw, integral, epsI, wStart, xS, xSD2, xSD4, wMax, DwStart,
      Uax, epsK;
  const double RRt2 = 0.707106781187, RtPi = 1.77245385091;
  int kMax, k;
  aPh = a + 0.5;
  aD2 = 0.5 * a;
  xS = x * x;
  xSD2 = 0.5 * xS;
  xSD4 = 0.0;
  if ((x > 3.0) && (a > 1.0)) {
    wMax = 0.5 * x * (sqrt(x * x + 4.0 * a - 2.0) - x);
    wMax = (wMax > 2.5) ? 8.0 * wMax : 25.0;
    kMax = 90;
    DwStart = wMax / ((double)kMax);
    Dw = DwStart / 27.0;
    epsI = 0.05;
    integral =
        pow(epsI, aPh) * (1.0 / aPh - epsI / (aPh + 1.0) +
                          0.5 * epsI * epsI *
                              ((1.0 - 1.0 / xS) / (aPh + 2.0) +
                               epsI * (1.0 / xS - 1.0 / 3.0) / (aPh + 3.0)));
    wStart = epsI;
    for (k = 1; k <= kMax; k++) {
      integral += th_Uax9Strip(a, x, wStart, Dw);
      wStart += DwStart;
    }
    Uax = integral / (exp(xSD4 + aPh * log(x)) * th_Gamma(aPh));
    return Uax;
  }
  if (a <= 1.0) {
    printf("The custom made parametric cylinder function is\n");
    printf("is not garantied to work form a <= 1.0\n\n");
    exit(1);
  }
  if (x < -10) {
    return for_large_negative(a, x) + x * x / 4.0;
  }
  epsK = 0.001 * eps;
  Uax = RtPi * exp(-xSD4) *
        (pow(2.0, -0.25 - aD2) * th_Hyper1F1GG(aD2 + 0.25, 0.5, xSD2, epsK) /
             th_Gamma(aD2 + 0.75) -
         x * pow(2.0, 0.25 - aD2) * th_Hyper1F1GG(aD2 + 0.75, 1.5, xSD2, epsK) /
             th_Gamma(aD2 + 0.25));
  return log(Uax);
}


x86 assembly (clang, -O0):
<ParCylU>:
push   %rbp
mov    %rsp,%rbp
sub    $0x110,%rsp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
movsd  0xe56(%rip),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  0xe3e(%rip),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0xe26(%rip),%xmm0
addsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xe14(%rip),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  0xdf3(%rip),%xmm0
mulsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xdcc(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1506 <ParCylU+0x366>
movsd  -0x10(%rbp),%xmm0
movsd  0xdd5(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1506 <ParCylU+0x366>
movsd  0xdab(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
movsd  0xde4(%rip),%xmm1
mulsd  -0x10(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xdb7(%rip),%xmm1
subsd  %xmm1,%xmm0
call   10b0 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0xa8(%rbp),%xmm0
ff
subsd  -0x18(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  0xda8(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    12e4 <ParCylU+0x144>
movsd  0xda6(%rip),%xmm0
mulsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
jmp    12f9 <ParCylU+0x159>
movsd  0xd84(%rip),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
jmp    12f9 <ParCylU+0x159>
movsd  -0xb0(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x5a,-0x9c(%rbp)
movsd  -0x70(%rbp),%xmm0
cvtsi2sdl -0x9c(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  0xd5d(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xd44(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm1
call   10a0 <pow@plt>
movsd  0xcc0(%rip),%xmm2
divsd  -0x28(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
movsd  0xcae(%rip),%xmm3
addsd  -0x28(%rbp),%xmm3
divsd  %xmm3,%xmm1
subsd  %xmm1,%xmm2
movsd  0xc81(%rip),%xmm1
mulsd  -0x48(%rbp),%xmm1
mulsd  -0x48(%rbp),%xmm1
movsd  0xc87(%rip),%xmm4
divsd  -0x58(%rbp),%xmm4
movsd  0xc7a(%rip),%xmm3
subsd  %xmm4,%xmm3
movsd  0xc96(%rip),%xmm4
addsd  -0x28(%rbp),%xmm4
divsd  %xmm4,%xmm3
movsd  -0x48(%rbp),%xmm4
movsd  0xc58(%rip),%xmm5
divsd  -0x58(%rbp),%xmm5
movsd  0xcab(%rip),%xmm6
subsd  %xmm6,%xmm5
mulsd  %xmm5,%xmm4
movsd  0xc1b(%rip),%xmm5
addsd  -0x28(%rbp),%xmm5
divsd  %xmm5,%xmm4
addsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x1,-0xa0(%rbp)
mov    -0xa0(%rbp),%eax
cmp    -0x9c(%rbp),%eax
jg     1474 <ParCylU+0x2d4>
movsd  -0x10(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
movsd  -0x50(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm3
call   1060 <th_Uax9Strip@plt>
cvtsi2sd %rax,%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x78(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0xa0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa0(%rbp)
jmp    1417 <ParCylU+0x277>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
call   1090 <log@plt>
movsd  -0xd0(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xc8(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1080 <exp@plt>
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
call   1070 <th_Gamma@plt>
movsd  -0xc0(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xb8(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1720 <ParCylU+0x580>
movsd  0xb12(%rip),%xmm0
ucomisd -0x10(%rbp),%xmm0
jb     153b <ParCylU+0x39b>
lea    0xb78(%rip),%rdi
call   1030 <printf@plt>
lea    0xb9d(%rip),%rdi
call   1030 <printf@plt>
mov    $0x1,%edi
call   10c0 <exit@plt>
movsd  0xae5(%rip),%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    1581 <ParCylU+0x3e1>
movsd  -0x10(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
call   1040 <for_large_negative@plt>
movsd  -0x18(%rbp),%xmm1
mulsd  -0x18(%rbp),%xmm1
movsd  0xaf1(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1720 <ParCylU+0x580>
movsd  0xacf(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
call   1080 <exp@plt>
movaps %xmm0,%xmm1
movsd  0xa4e(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movsd  0xa7a(%rip),%xmm1
subsd  -0x30(%rbp),%xmm1
movsd  0xa65(%rip),%xmm0
call   10a0 <pow@plt>
movsd  %xmm0,-0x108(%rbp)
ff
movsd  0xa38(%rip),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm2
movsd  -0x88(%rbp),%xmm3
ff
movsd  0x9f6(%rip),%xmm1
call   1050 <th_Hyper1F1GG@plt>
movaps %xmm0,%xmm1
movsd  -0x108(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
movsd  0xa0a(%rip),%xmm0
addsd  -0x30(%rbp),%xmm0
call   1070 <th_Gamma@plt>
movaps %xmm0,%xmm1
movsd  -0x100(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
movsd  0x9c4(%rip),%xmm1
subsd  -0x30(%rbp),%xmm1
movsd  0x9cf(%rip),%xmm0
call   10a0 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xf8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  0x9a3(%rip),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm2
movsd  -0x88(%rbp),%xmm3
ff
movsd  0x981(%rip),%xmm1
call   1050 <th_Hyper1F1GG@plt>
movaps %xmm0,%xmm1
movsd  -0xf0(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
movsd  0x955(%rip),%xmm0
addsd  -0x30(%rbp),%xmm0
call   1070 <th_Gamma@plt>
movsd  -0xe8(%rbp),%xmm2
ff
movsd  -0xe0(%rbp),%xmm1
ff
movaps %xmm0,%xmm3
movsd  -0xd8(%rbp),%xmm0
ff
divsd  %xmm3,%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
call   1090 <log@plt>
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x110,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ParCylU>:
push   %rbx
sub    $0x80,%rsp
movapd %xmm0,%xmm4
movsd  0xe3c(%rip),%xmm5
movapd %xmm1,%xmm6
mulsd  %xmm1,%xmm6
ucomisd 0xe34(%rip),%xmm0
jbe    1226 <ParCylU+0x76>
ucomisd 0xe32(%rip),%xmm1
jbe    1226 <ParCylU+0x76>
movsd  0xe18(%rip),%xmm2
movsd  0xe60(%rip),%xmm0
mulsd  %xmm4,%xmm0
addsd  %xmm6,%xmm0
addsd  0xe58(%rip),%xmm0
mulsd  %xmm1,%xmm2
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm0
movapd %xmm6,0x10(%rsp)
movsd  %xmm1,0x30(%rsp)
movsd  %xmm4,0x8(%rsp)
jb     1270 <ParCylU+0xc0>
sqrtsd %xmm0,%xmm0
jmp    1293 <ParCylU+0xe3>
movsd  0xdda(%rip),%xmm0
ucomisd %xmm4,%xmm0
jae    1561 <ParCylU+0x3b1>
movsd  0xdd8(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1437 <ParCylU+0x287>
movapd %xmm4,%xmm0
movapd %xmm6,0x10(%rsp)
call   1050 <for_large_negative@plt>
movapd 0x10(%rsp),%xmm1
mulsd  0xdc9(%rip),%xmm1
addsd  %xmm0,%xmm1
jmp    142a <ParCylU+0x27a>
movsd  %xmm2,(%rsp)
call   10c0 <sqrt@plt>
movsd  (%rsp),%xmm2
movsd  0xd79(%rip),%xmm5
movsd  0x8(%rsp),%xmm4
movsd  0x30(%rsp),%xmm1
addsd  %xmm4,%xmm5
subsd  %xmm1,%xmm0
mulsd  %xmm0,%xmm2
ucomisd 0xdb9(%rip),%xmm2
jbe    12bb <ParCylU+0x10b>
mulsd  0xdb7(%rip),%xmm2
divsd  0xdbf(%rip),%xmm2
jmp    12c3 <ParCylU+0x113>
movsd  0xdad(%rip),%xmm2
movsd  %xmm2,(%rsp)
movapd %xmm2,%xmm0
divsd  0xdac(%rip),%xmm0
movsd  %xmm0,0x48(%rsp)
movsd  0xda6(%rip),%xmm0
movsd  %xmm0,0x20(%rsp)
movapd %xmm5,%xmm1
movapd %xmm5,0x50(%rsp)
call   10b0 <pow@plt>
movsd  0xd09(%rip),%xmm1
movapd 0x50(%rsp),%xmm4
movapd %xmm4,%xmm2
addsd  %xmm1,%xmm2
movapd %xmm4,%xmm3
movapd %xmm4,%xmm5
unpcklpd %xmm2,%xmm3
movapd 0xd7f(%rip),%xmm2
divpd  %xmm3,%xmm2
movapd %xmm1,%xmm3
divsd  0x10(%rsp),%xmm3
movapd %xmm2,%xmm4
unpckhpd %xmm2,%xmm4
subsd  %xmm3,%xmm1
addsd  0xd4d(%rip),%xmm3
unpcklpd %xmm3,%xmm3
mulpd  0xd61(%rip),%xmm3
movsd  %xmm1,%xmm3
movapd %xmm5,%xmm1
unpcklpd %xmm5,%xmm1
addpd  0xd5d(%rip),%xmm1
divpd  %xmm1,%xmm3
movapd %xmm3,%xmm1
unpckhpd %xmm3,%xmm1
addsd  %xmm3,%xmm1
mulsd  0xd1d(%rip),%xmm1
subsd  %xmm4,%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
mov    $0x5a,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  0x8(%rsp),%xmm0
movsd  0x30(%rsp),%xmm1
movsd  0x20(%rsp),%xmm2
movsd  0x48(%rsp),%xmm3
call   1070 <th_Uax9Strip@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movapd 0x10(%rsp),%xmm1
addsd  %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
movsd  0x20(%rsp),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,0x20(%rsp)
add    $0xffffffff,%ebx
jne    13a0 <ParCylU+0x1f0>
movsd  0x30(%rsp),%xmm0
call   10a0 <log@plt>
mulsd  0x50(%rsp),%xmm0
xorpd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
call   1090 <exp@plt>
movsd  %xmm0,0x20(%rsp)
movapd 0x50(%rsp),%xmm0
call   1080 <th_Gamma@plt>
mulsd  0x20(%rsp),%xmm0
movapd 0x10(%rsp),%xmm1
divsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
add    $0x80,%rsp
pop    %rbx
ret
mulsd  %xmm5,%xmm4
movsd  %xmm4,0x8(%rsp)
mulsd  %xmm5,%xmm6
movapd %xmm6,0x10(%rsp)
mulsd  0xbcd(%rip),%xmm2
movsd  %xmm2,(%rsp)
movsd  0xbc8(%rip),%xmm0
subsd  %xmm4,%xmm0
movsd  %xmm1,0x30(%rsp)
call   1040 <exp2@plt>
movapd %xmm0,0x20(%rsp)
movsd  0x8(%rsp),%xmm0
addsd  0xbad(%rip),%xmm0
movsd  %xmm0,0x48(%rsp)
movsd  0xb6f(%rip),%xmm1
movaps 0x10(%rsp),%xmm2
movsd  (%rsp),%xmm3
call   1060 <th_Hyper1F1GG@plt>
movapd %xmm0,0x70(%rsp)
movsd  0xb8a(%rip),%xmm0
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x68(%rsp)
call   1080 <th_Gamma@plt>
movapd %xmm0,0x50(%rsp)
movsd  0xb63(%rip),%xmm0
subsd  0x8(%rsp),%xmm0
call   1040 <exp2@plt>
mulsd  0x30(%rsp),%xmm0
movapd %xmm0,0x30(%rsp)
movsd  0xb54(%rip),%xmm1
movsd  0x68(%rsp),%xmm0
movaps 0x10(%rsp),%xmm2
movsd  (%rsp),%xmm3
call   1060 <th_Hyper1F1GG@plt>
movapd 0x70(%rsp),%xmm2
unpcklpd 0x30(%rsp),%xmm2
movapd 0x20(%rsp),%xmm1
unpcklpd %xmm0,%xmm1
mulpd  %xmm2,%xmm1
movapd %xmm1,0x20(%rsp)
movsd  0x48(%rsp),%xmm0
call   1080 <th_Gamma@plt>
movapd 0x50(%rsp),%xmm1
unpcklpd %xmm0,%xmm1
movapd 0x20(%rsp),%xmm0
divpd  %xmm1,%xmm0
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
subsd  %xmm1,%xmm0
mulsd  0xaf4(%rip),%xmm0
add    $0x80,%rsp
pop    %rbx
jmp    10a0 <log@plt>
lea    0xb68(%rip),%rdi
call   1030 <puts@plt>
lea    0xb8c(%rip),%rdi
call   1030 <puts@plt>
mov    $0x1,%edi
call   10d0 <exit@plt>

LLVM IR (-O0):
@.str = [49 x i8] c"The custom made parametric cylinder function is\0A\00"
@.str.1 = [41 x i8] c"is not garantied to work form a <= 1.0\0A\0A\00"
define double @ParCylU(double %0, double %1, double %2) {
%4 = alloca double
%5 = alloca double
%6 = alloca double
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
%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 i32
%24 = alloca i32
store double %0, double* %5
store double %1, double* %6
store double %2, double* %7
store double 0x3FE6A09E667F4BB8, double* %21
store double 0x3FFC5BF891B53E4D, double* %22
%25 = load double, double* %5
%26 = fadd double %25, 5.000000e-01
store double %26, double* %8
%27 = load double, double* %5
%28 = fmul double 5.000000e-01, %27
store double %28, double* %9
%29 = load double, double* %6
%30 = load double, double* %6
%31 = fmul double %29, %30
store double %31, double* %14
%32 = load double, double* %14
%33 = fmul double 5.000000e-01, %32
store double %33, double* %15
store double 0.000000e+00, double* %16
%34 = load double, double* %6
%35 = fcmp ogt double %34, 3.000000e+00
br i1 %35, label %36, label %130
36: ; preds = %3
%37 = load double, double* %5
%38 = fcmp ogt double %37, 1.000000e+00
br i1 %38, label %39, label %130
39: ; preds = %36
%40 = load double, double* %6
%41 = fmul double 5.000000e-01, %40
%42 = load double, double* %6
%43 = load double, double* %6
%44 = load double, double* %5
%45 = fmul double 4.000000e+00, %44
%46 = call double @llvm.fmuladd.f64(double %42, double %43, double %45)
%47 = fsub double %46, 2.000000e+00
%48 = call double @sqrt(double %47)
%49 = load double, double* %6
%50 = fsub double %48, %49
%51 = fmul double %41, %50
store double %51, double* %17
%52 = load double, double* %17
%53 = fcmp ogt double %52, 2.500000e+00
br i1 %53, label %54, label %57
54: ; preds = %39
%55 = load double, double* %17
%56 = fmul double 8.000000e+00, %55
br label %58
57: ; preds = %39
br label %58
58: ; preds = %57, %54
%59 = phi double [ %56, %54 ], [ 2.500000e+01, %57 ]
store double %59, double* %17
store i32 90, i32* %23
%60 = load double, double* %17
%61 = load i32, i32* %23
%62 = sitofp i32 %61 to double
%63 = fdiv double %60, %62
store double %63, double* %18
%64 = load double, double* %18
%65 = fdiv double %64, 2.700000e+01
store double %65, double* %10
store double 5.000000e-02, double* %12
%66 = load double, double* %12
%67 = load double, double* %8
%68 = call double @pow(double %66, double %67)
%69 = load double, double* %8
%70 = fdiv double 1.000000e+00, %69
%71 = load double, double* %12
%72 = load double, double* %8
%73 = fadd double %72, 1.000000e+00
%74 = fdiv double %71, %73
%75 = fsub double %70, %74
%76 = load double, double* %12
%77 = fmul double 5.000000e-01, %76
%78 = load double, double* %12
%79 = fmul double %77, %78
%80 = load double, double* %14
%81 = fdiv double 1.000000e+00, %80
%82 = fsub double 1.000000e+00, %81
%83 = load double, double* %8
%84 = fadd double %83, 2.000000e+00
%85 = fdiv double %82, %84
%86 = load double, double* %12
%87 = load double, double* %14
%88 = fdiv double 1.000000e+00, %87
%89 = fsub double %88, 0x3FD5555555555555
%90 = fmul double %86, %89
%91 = load double, double* %8
%92 = fadd double %91, 3.000000e+00
%93 = fdiv double %90, %92
%94 = fadd double %85, %93
%95 = call double @llvm.fmuladd.f64(double %79, double %94, double %75)
%96 = fmul double %68, %95
store double %96, double* %11
%97 = load double, double* %12
store double %97, double* %13
store i32 1, i32* %24
br label %98
98: ; preds = %114, %58
%99 = load i32, i32* %24
%100 = load i32, i32* %23
%101 = icmp sle i32 %99, %100
br i1 %101, label %102, label %117
102: ; preds = %98
%103 = load double, double* %5
%104 = load double, double* %6
%105 = load double, double* %13
%106 = load double, double* %10
%107 = call i64 @th_Uax9Strip(double %103, double %104, double %105, double %106)
%108 = sitofp i64 %107 to double
%109 = load double, double* %11
%110 = fadd double %109, %108
store double %110, double* %11
%111 = load double, double* %18
%112 = load double, double* %13
%113 = fadd double %112, %111
store double %113, double* %13
br label %114
114: ; preds = %102
%115 = load i32, i32* %24
%116 = add nsw i32 %115, 1
store i32 %116, i32* %24
br label %98
117: ; preds = %98
%118 = load double, double* %11
%119 = load double, double* %16
%120 = load double, double* %8
%121 = load double, double* %6
%122 = call double @log(double %121)
%123 = call double @llvm.fmuladd.f64(double %120, double %122, double %119)
%124 = call double @exp(double %123)
%125 = load double, double* %8
%126 = call double @th_Gamma(double %125)
%127 = fmul double %124, %126
%128 = fdiv double %118, %127
store double %128, double* %19
%129 = load double, double* %19
store double %129, double* %4
br label %188
130: ; preds = %36, %3
%131 = load double, double* %5
%132 = fcmp ole double %131, 1.000000e+00
br i1 %132, label %133, label %137
133: ; preds = %130
%134 = call i32 @printf(i8* getelementptr inbounds ([49 x i8], [49 x i8]* @.str, i64 0, i64 0))
%135 = call i32 @printf(i8* getelementptr inbounds ([41 x i8], [41 x i8]* @.str.1, i64 0, i64 0))
%136 = call i32 @exit(i32 1)
unreachable
137: ; preds = %130
%138 = load double, double* %6
%139 = fcmp olt double %138, -1.000000e+01
br i1 %139, label %140, label %149
140: ; preds = %137
%141 = load double, double* %5
%142 = load double, double* %6
%143 = call double @for_large_negative(double %141, double %142)
%144 = load double, double* %6
%145 = load double, double* %6
%146 = fmul double %144, %145
%147 = fdiv double %146, 4.000000e+00
%148 = fadd double %143, %147
store double %148, double* %4
br label %188
149: ; preds = %137
%150 = load double, double* %7
%151 = fmul double 1.000000e-03, %150
store double %151, double* %20
%152 = load double, double* %16
%153 = fneg double %152
%154 = call double @exp(double %153)
%155 = fmul double 0x3FFC5BF891B53E4D, %154
%156 = load double, double* %9
%157 = fsub double -2.500000e-01, %156
%158 = call double @pow(double 2.000000e+00, double %157)
%159 = load double, double* %9
%160 = fadd double %159, 2.500000e-01
%161 = load double, double* %15
%162 = load double, double* %20
%163 = call double @th_Hyper1F1GG(double %160, double 5.000000e-01, double %161, double %162)
%164 = fmul double %158, %163
%165 = load double, double* %9
%166 = fadd double %165, 7.500000e-01
%167 = call double @th_Gamma(double %166)
%168 = fdiv double %164, %167
%169 = load double, double* %6
%170 = load double, double* %9
%171 = fsub double 2.500000e-01, %170
%172 = call double @pow(double 2.000000e+00, double %171)
%173 = fmul double %169, %172
%174 = load double, double* %9
%175 = fadd double %174, 7.500000e-01
%176 = load double, double* %15
%177 = load double, double* %20
%178 = call double @th_Hyper1F1GG(double %175, double 1.500000e+00, double %176, double %177)
%179 = fmul double %173, %178
%180 = load double, double* %9
%181 = fadd double %180, 2.500000e-01
%182 = call double @th_Gamma(double %181)
%183 = fdiv double %179, %182
%184 = fsub double %168, %183
%185 = fmul double %155, %184
store double %185, double* %19
%186 = load double, double* %19
%187 = call double @log(double %186)
store double %187, double* %4
br label %188
188: ; preds = %149, %140, %117
%189 = load double, double* %4
ret double %189
}

LLVM IR (-O2):
@str = [48 x i8] c"The custom made parametric cylinder function is\00"
@str.2 = [40 x i8] c"is not garantied to work form a <= 1.0\0A\00"
define double @ParCylU(double %0, double %1, double %2) {
%4 = fadd double %0, 5.000000e-01
%5 = fmul double %0, 5.000000e-01
%6 = fmul double %1, %1
%7 = fmul double %6, 5.000000e-01
%8 = fcmp ogt double %1, 3.000000e+00
%9 = fcmp ogt double %0, 1.000000e+00
%10 = and i1 %9, %8
br i1 %10, label %11, label %64
11: ; preds = %3
%12 = fmul double %1, 5.000000e-01
%13 = fmul double %0, 4.000000e+00
%14 = tail call double @llvm.fmuladd.f64(double %1, double %1, double %13)
%15 = fadd double %14, -2.000000e+00
%16 = tail call double @sqrt(double %15)
%17 = fsub double %16, %1
%18 = fmul double %12, %17
%19 = fcmp ogt double %18, 2.500000e+00
%20 = fmul double %18, 8.000000e+00
%21 = fdiv double %20, 9.000000e+01
%22 = select i1 %19, double %21, double 0x3FD1C71C71C71C72
%23 = fdiv double %22, 2.700000e+01
%24 = tail call double @pow(double 5.000000e-02, double %4)
%25 = fadd double %4, 1.000000e+00
%26 = insertelement <2 x double> poison, double %4, i64 0
%27 = insertelement <2 x double> %26, double %25, i64 1
%28 = fdiv <2 x double> <double 1.000000e+00, double 5.000000e-02>, %27
%29 = shufflevector <2 x double> %28, <2 x double> poison, <2 x i32> <i32 1, i32 undef>
%30 = fsub <2 x double> %28, %29
%31 = extractelement <2 x double> %30, i64 0
%32 = fdiv double 1.000000e+00, %6
%33 = fadd double %32, 0xBFD5555555555555
%34 = insertelement <2 x double> poison, double %32, i64 0
%35 = insertelement <2 x double> %34, double %33, i64 1
%36 = fsub <2 x double> <double 1.000000e+00, double poison>, %35
%37 = fmul <2 x double> %35, <double poison, double 5.000000e-02>
%38 = shufflevector <2 x double> %36, <2 x double> %37, <2 x i32> <i32 0, i32 3>
%39 = shufflevector <2 x double> %26, <2 x double> poison, <2 x i32> zeroinitializer
%40 = fadd <2 x double> %39, <double 2.000000e+00, double 3.000000e+00>
%41 = fdiv <2 x double> %38, %40
%42 = shufflevector <2 x double> %41, <2 x double> poison, <2 x i32> <i32 1, i32 undef>
%43 = fadd <2 x double> %41, %42
%44 = extractelement <2 x double> %43, i64 0
%45 = tail call double @llvm.fmuladd.f64(double %44, double 0x3F547AE147AE147C, double %31)
%46 = fmul double %45, %24
br label %47
47: ; preds = %11, %47
%48 = phi i32 [ 1, %11 ], [ %55, %47 ]
%49 = phi double [ %46, %11 ], [ %53, %47 ]
%50 = phi double [ 5.000000e-02, %11 ], [ %54, %47 ]
%51 = tail call i64 @th_Uax9Strip(double %0, double %1, double %50, double %23)
%52 = sitofp i64 %51 to double
%53 = fadd double %49, %52
%54 = fadd double %22, %50
%55 = add nuw nsw i32 %48, 1
%56 = icmp eq i32 %55, 91
br i1 %56, label %57, label %47
57: ; preds = %47
%58 = tail call double @log(double %1)
%59 = tail call double @llvm.fmuladd.f64(double %4, double %58, double 0.000000e+00)
%60 = tail call double @exp(double %59)
%61 = tail call double @th_Gamma(double %4)
%62 = fmul double %60, %61
%63 = fdiv double %53, %62
br label %102
64: ; preds = %3
%65 = fcmp ugt double %0, 1.000000e+00
br i1 %65, label %70, label %66
66: ; preds = %64
%67 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([48 x i8], [48 x i8]* @str, i64 0, i64 0))
%68 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([40 x i8], [40 x i8]* @str.2, i64 0, i64 0))
%69 = tail call i32 @exit(i32 1)
unreachable
70: ; preds = %64
%71 = fcmp olt double %1, -1.000000e+01
br i1 %71, label %72, label %76
72: ; preds = %70
%73 = tail call double @for_large_negative(double %0, double %1)
%74 = fmul double %6, 2.500000e-01
%75 = fadd double %74, %73
br label %102
76: ; preds = %70
%77 = fmul double %2, 1.000000e-03
%78 = fsub double -2.500000e-01, %5
%79 = tail call double @exp2(double %78)
%80 = fadd double %5, 2.500000e-01
%81 = tail call double @th_Hyper1F1GG(double %80, double 5.000000e-01, double %7, double %77)
%82 = fadd double %5, 7.500000e-01
%83 = tail call double @th_Gamma(double %82)
%84 = fsub double 2.500000e-01, %5
%85 = tail call double @exp2(double %84)
%86 = fmul double %85, %1
%87 = tail call double @th_Hyper1F1GG(double %82, double 1.500000e+00, double %7, double %77)
%88 = insertelement <2 x double> poison, double %81, i64 0
%89 = insertelement <2 x double> %88, double %86, i64 1
%90 = insertelement <2 x double> poison, double %79, i64 0
%91 = insertelement <2 x double> %90, double %87, i64 1
%92 = fmul <2 x double> %89, %91
%93 = tail call double @th_Gamma(double %80)
%94 = insertelement <2 x double> poison, double %83, i64 0
%95 = insertelement <2 x double> %94, double %93, i64 1
%96 = fdiv <2 x double> %92, %95
%97 = shufflevector <2 x double> %96, <2 x double> poison, <2 x i32> <i32 1, i32 undef>
%98 = fsub <2 x double> %96, %97
%99 = extractelement <2 x double> %98, i64 0
%100 = fmul double %99, 0x3FFC5BF891B53E4D
%101 = tail call double @log(double %100)
br label %102
102: ; preds = %76, %72, %57
%103 = phi double [ %63, %57 ], [ %75, %72 ], [ %101, %76 ]
ret double %103
}