double DEVpnjlim(vnew, vold, vt, vcrit, icheck)
double vnew;
double vold;
double vt;
double vcrit;
int *icheck;
{
  double arg;
  if ((vnew > vcrit) &&
      ((((vnew - vold) < 0) ? -(vnew - vold) : (vnew - vold)) > (vt + vt))) {
    if (vold > 0) {
      arg = 1 + (vnew - vold) / vt;
      if (arg > 0) {
        vnew = vold + vt * log(arg);
      } else {
        vnew = vcrit;
      }
    } else {
      vnew = vt * log(vnew / vt);
    }
    *icheck = 1;
  } else {
    *icheck = 0;
  }
  return (vnew);
}


x86 assembly (clang, -O0):
<DEVpnjlim>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
movsd  %xmm3,-0x20(%rbp)
mov    %rdi,-0x28(%rbp)
movsd  -0x8(%rbp),%xmm0
ucomisd -0x20(%rbp),%xmm0
jbe    1273 <DEVpnjlim+0x163>
movsd  -0x8(%rbp),%xmm1
subsd  -0x10(%rbp),%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    1182 <DEVpnjlim+0x72>
movsd  -0x8(%rbp),%xmm0
subsd  -0x10(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    1191 <DEVpnjlim+0x81>
movsd  -0x8(%rbp),%xmm0
subsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
addsd  -0x18(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jbe    1273 <DEVpnjlim+0x163>
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    123a <DEVpnjlim+0x12a>
movsd  -0x8(%rbp),%xmm1
subsd  -0x10(%rbp),%xmm1
divsd  -0x18(%rbp),%xmm1
movsd  0xe2d(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    122b <DEVpnjlim+0x11b>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x30(%rbp),%xmm0
call   1030 <log@plt>
movsd  -0x48(%rbp),%xmm1
movaps %xmm0,%xmm2
movsd  -0x40(%rbp),%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1235 <DEVpnjlim+0x125>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1264 <DEVpnjlim+0x154>
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x8(%rbp),%xmm0
divsd  -0x18(%rbp),%xmm0
call   1030 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x50(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
mov    -0x28(%rbp),%rax
movl   $0x1,(%rax)
jmp    127d <DEVpnjlim+0x16d>
mov    -0x28(%rbp),%rax
movl   $0x0,(%rax)
movsd  -0x8(%rbp),%xmm0
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<DEVpnjlim>:
push   %rbp
push   %rbx
sub    $0x18,%rsp
mov    %rdi,%rbx
xor    %ebp,%ebp
ucomisd %xmm3,%xmm0
jbe    11cc <DEVpnjlim+0xbc>
movapd %xmm0,%xmm7
subsd  %xmm1,%xmm7
movapd 0xecb(%rip),%xmm4
xorpd  %xmm7,%xmm4
xorpd  %xmm6,%xmm6
movapd %xmm7,%xmm5
cmpltsd %xmm6,%xmm5
andpd  %xmm5,%xmm4
andnpd %xmm7,%xmm5
orpd   %xmm4,%xmm5
movapd %xmm2,%xmm4
addsd  %xmm2,%xmm4
ucomisd %xmm4,%xmm5
jbe    11cc <DEVpnjlim+0xbc>
ucomisd %xmm6,%xmm1
jbe    11ac <DEVpnjlim+0x9c>
divsd  %xmm2,%xmm7
addsd  0xe9e(%rip),%xmm7
mov    $0x1,%ebp
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm7
jbe    11c8 <DEVpnjlim+0xb8>
movapd %xmm7,%xmm0
movsd  %xmm2,0x8(%rsp)
movsd  %xmm1,0x10(%rsp)
call   1030 <log@plt>
movsd  0x8(%rsp),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x10(%rsp),%xmm1
movapd %xmm1,%xmm0
jmp    11cc <DEVpnjlim+0xbc>
divsd  %xmm2,%xmm0
movsd  %xmm2,0x8(%rsp)
call   1030 <log@plt>
mulsd  0x8(%rsp),%xmm0
mov    $0x1,%ebp
jmp    11cc <DEVpnjlim+0xbc>
movapd %xmm3,%xmm0
mov    %ebp,(%rbx)
add    $0x18,%rsp
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
define double @DEVpnjlim(double %0, double %1, double %2, double %3, i32* %4) {
%6 = alloca double
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca i32*
%11 = alloca double
store double %0, double* %6
store double %1, double* %7
store double %2, double* %8
store double %3, double* %9
store i32* %4, i32** %10
%12 = load double, double* %6
%13 = load double, double* %9
%14 = fcmp ogt double %12, %13
br i1 %14, label %15, label %65
15: ; preds = %5
%16 = load double, double* %6
%17 = load double, double* %7
%18 = fsub double %16, %17
%19 = fcmp olt double %18, 0.000000e+00
br i1 %19, label %20, label %25
20: ; preds = %15
%21 = load double, double* %6
%22 = load double, double* %7
%23 = fsub double %21, %22
%24 = fneg double %23
br label %29
25: ; preds = %15
%26 = load double, double* %6
%27 = load double, double* %7
%28 = fsub double %26, %27
br label %29
29: ; preds = %25, %20
%30 = phi double [ %24, %20 ], [ %28, %25 ]
%31 = load double, double* %8
%32 = load double, double* %8
%33 = fadd double %31, %32
%34 = fcmp ogt double %30, %33
br i1 %34, label %35, label %65
35: ; preds = %29
%36 = load double, double* %7
%37 = fcmp ogt double %36, 0.000000e+00
br i1 %37, label %38, label %56
38: ; preds = %35
%39 = load double, double* %6
%40 = load double, double* %7
%41 = fsub double %39, %40
%42 = load double, double* %8
%43 = fdiv double %41, %42
%44 = fadd double 1.000000e+00, %43
store double %44, double* %11
%45 = load double, double* %11
%46 = fcmp ogt double %45, 0.000000e+00
br i1 %46, label %47, label %53
47: ; preds = %38
%48 = load double, double* %7
%49 = load double, double* %8
%50 = load double, double* %11
%51 = call double @log(double %50)
%52 = call double @llvm.fmuladd.f64(double %49, double %51, double %48)
store double %52, double* %6
br label %55
53: ; preds = %38
%54 = load double, double* %9
store double %54, double* %6
br label %55
55: ; preds = %53, %47
br label %63
56: ; preds = %35
%57 = load double, double* %8
%58 = load double, double* %6
%59 = load double, double* %8
%60 = fdiv double %58, %59
%61 = call double @log(double %60)
%62 = fmul double %57, %61
store double %62, double* %6
br label %63
63: ; preds = %56, %55
%64 = load i32*, i32** %10
store i32 1, i32* %64
br label %67
65: ; preds = %29, %5
%66 = load i32*, i32** %10
store i32 0, i32* %66
br label %67
67: ; preds = %65, %63
%68 = load double, double* %6
ret double %68
}

LLVM IR (-O2):
define double @DEVpnjlim(double %0, double %1, double %2, double %3, i32* nocapture writeonly %4) {
%6 = fcmp ogt double %0, %3
br i1 %6, label %7, label %27
7: ; preds = %5
%8 = fsub double %0, %1
%9 = fcmp olt double %8, 0.000000e+00
%10 = fneg double %8
%11 = select i1 %9, double %10, double %8
%12 = fadd double %2, %2
%13 = fcmp ogt double %11, %12
br i1 %13, label %14, label %27
14: ; preds = %7
%15 = fcmp ogt double %1, 0.000000e+00
br i1 %15, label %16, label %23
16: ; preds = %14
%17 = fdiv double %8, %2
%18 = fadd double %17, 1.000000e+00
%19 = fcmp ogt double %18, 0.000000e+00
br i1 %19, label %20, label %27
20: ; preds = %16
%21 = tail call double @log(double %18)
%22 = tail call double @llvm.fmuladd.f64(double %2, double %21, double %1)
br label %27
23: ; preds = %14
%24 = fdiv double %0, %2
%25 = tail call double @log(double %24)
%26 = fmul double %25, %2
br label %27
27: ; preds = %5, %7, %23, %20, %16
%28 = phi i32 [ 1, %16 ], [ 1, %20 ], [ 1, %23 ], [ 0, %7 ], [ 0, %5 ]
%29 = phi double [ %3, %16 ], [ %22, %20 ], [ %26, %23 ], [ %0, %7 ], [ %0, %5 ]
store i32 %28, i32* %4
ret double %29
}