double midpnt(double (*func)(double), double a, double b, int n) {
  double x, tnm, sum, del, ddel;
  static double s;
  int it, j;
  if (n == 1) {
    return (s = (b - a) * ((*func)(0.5 * (a + b))));
  } else {
    for (it = 1, j = 1; j < n - 1; j++)
      it *= 3;
    tnm = it;
    del = (b - a) / (3.0 * tnm);
    ddel = del + del;
    x = a + 0.5 * del;
    sum = 0.0;
    for (j = 1; j <= it; j++) {
      sum += ((*func)(x));
      x += ddel;
      sum += ((*func)(x));
      x += del;
    }
    s = (s + (b - a) * sum / tnm) / 3.0;
    return s;
  }
}


x86 assembly (clang, -O0):
<midpnt>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    %esi,-0x24(%rbp)
cmpl   $0x1,-0x24(%rbp)
jne    116c <midpnt+0x6c>
movsd  -0x20(%rbp),%xmm0
subsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x10(%rbp),%rax
movsd  -0x18(%rbp),%xmm1
addsd  -0x20(%rbp),%xmm1
movsd  0xeb8(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   *%rax
movaps %xmm0,%xmm1
movsd  -0x60(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x2ec6(%rip)
movsd  %xmm0,-0x8(%rbp)
jmp    12aa <midpnt+0x1aa>
movl   $0x1,-0x54(%rbp)
movl   $0x1,-0x58(%rbp)
mov    -0x58(%rbp),%eax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    11a0 <midpnt+0xa0>
imul   $0x3,-0x54(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    117a <midpnt+0x7a>
cvtsi2sdl -0x54(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x20(%rbp),%xmm0
subsd  -0x18(%rbp),%xmm0
movsd  0xe4c(%rip),%xmm1
mulsd  -0x38(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x18(%rbp),%xmm1
movsd  -0x48(%rbp),%xmm0
movsd  0xe15(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x1,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x54(%rbp),%eax
jg     1269 <midpnt+0x169>
mov    -0x10(%rbp),%rax
movsd  -0x30(%rbp),%xmm0
call   *%rax
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x50(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movsd  -0x30(%rbp),%xmm0
call   *%rax
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1207 <midpnt+0x107>
movsd  0x2db7(%rip),%xmm0
movsd  -0x20(%rbp),%xmm1
subsd  -0x18(%rbp),%xmm1
mulsd  -0x40(%rbp),%xmm1
divsd  -0x38(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  0xd77(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x2d8b(%rip)
movsd  0x2d83(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<midpnt>:
push   %rbp
push   %rbx
sub    $0x38,%rsp
mov    %rdi,%rbx
cmp    $0x1,%esi
jne    1135 <midpnt+0x35>
movapd %xmm1,%xmm2
subsd  %xmm0,%xmm2
movsd  %xmm2,0x8(%rsp)
addsd  %xmm1,%xmm0
mulsd  0xf00(%rip),%xmm0
call   *%rbx
mulsd  0x8(%rsp),%xmm0
jmp    1313 <midpnt+0x213>
mov    $0x1,%ebp
mov    $0x1,%edx
cmp    $0x3,%esi
jl     1248 <midpnt+0x148>
lea    -0x2(%rsi),%r8d
mov    $0x1,%eax
mov    $0x1,%edx
cmp    $0x8,%r8d
jb     1231 <midpnt+0x131>
mov    %r8d,%eax
and    $0xfffffff8,%eax
lea    -0x8(%rax),%ecx
mov    %ecx,%edi
shr    $0x3,%edi
add    $0x1,%edi
mov    %edi,%edx
and    $0x7,%edx
cmp    $0x38,%ecx
jae    1189 <midpnt+0x89>
movdqa 0xe7d(%rip),%xmm5
movdqa %xmm5,%xmm2
jmp    11db <midpnt+0xdb>
and    $0xfffffff8,%edi
movdqa 0xe6c(%rip),%xmm5
movdqa 0xe74(%rip),%xmm3
movdqa %xmm5,%xmm2
pshufd $0xf5,%xmm5,%xmm4
pmuludq %xmm3,%xmm5
pshufd $0xe8,%xmm5,%xmm5
pmuludq %xmm3,%xmm4
pshufd $0xe8,%xmm4,%xmm4
punpckldq %xmm4,%xmm5
pshufd $0xf5,%xmm2,%xmm4
pmuludq %xmm3,%xmm2
pshufd $0xe8,%xmm2,%xmm2
pmuludq %xmm3,%xmm4
pshufd $0xe8,%xmm4,%xmm4
punpckldq %xmm4,%xmm2
add    $0xfffffff8,%edi
jne    11a0 <midpnt+0xa0>
test   %edx,%edx
je     11fd <midpnt+0xfd>
nop
movdqa %xmm5,%xmm3
paddd  %xmm5,%xmm3
paddd  %xmm3,%xmm5
movdqa %xmm2,%xmm3
paddd  %xmm2,%xmm3
paddd  %xmm3,%xmm2
add    $0xffffffff,%edx
jne    11e0 <midpnt+0xe0>
pshufd $0xf5,%xmm5,%xmm3
pshufd $0xf5,%xmm2,%xmm4
pmuludq %xmm3,%xmm4
pmuludq %xmm5,%xmm2
pshufd $0xee,%xmm2,%xmm3
pmuludq %xmm2,%xmm3
pshufd $0xaa,%xmm4,%xmm2
pmuludq %xmm4,%xmm2
pmuludq %xmm3,%xmm2
movd   %xmm2,%edx
cmp    %eax,%r8d
je     1248 <midpnt+0x148>
or     $0x1,%eax
not    %eax
add    %esi,%eax
cs nopw 0x0(%rax,%rax,1)
nop
lea    (%rdx,%rdx,2),%edx
add    $0xffffffff,%eax
jne    1240 <midpnt+0x140>
xorps  %xmm3,%xmm3
cvtsi2sd %edx,%xmm3
subsd  %xmm0,%xmm1
movsd  0xdc5(%rip),%xmm2
movsd  %xmm3,0x18(%rsp)
mulsd  %xmm3,%xmm2
movsd  %xmm1,0x20(%rsp)
divsd  %xmm2,%xmm1
movapd %xmm1,%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,0x28(%rsp)
movsd  0xda3(%rip),%xmm2
movsd  %xmm1,0x30(%rsp)
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
cmp    $0x2,%edx
cmovae %edx,%ebp
xorpd  %xmm1,%xmm1
nopl   (%rax)
movsd  %xmm1,0x10(%rsp)
movsd  %xmm2,0x8(%rsp)
movapd %xmm2,%xmm0
call   *%rbx
addsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  0x8(%rsp),%xmm0
addsd  0x28(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x8(%rsp),%xmm0
call   *%rbx
movsd  0x8(%rsp),%xmm2
movapd %xmm0,%xmm1
addsd  0x10(%rsp),%xmm1
addsd  0x30(%rsp),%xmm2
add    $0xffffffff,%ebp
jne    12a0 <midpnt+0x1a0>
movsd  0x20(%rsp),%xmm0
mulsd  %xmm1,%xmm0
divsd  0x18(%rsp),%xmm0
addsd  0x2d1d(%rip),%xmm0
divsd  0xd0d(%rip),%xmm0
movsd  %xmm0,0x2d0d(%rip)
add    $0x38,%rsp
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
@midpnt.s = internal global double 0.000000e+00
define double @midpnt(double (double)* %0, double %1, double %2, i32 %3) {
%5 = alloca double
%6 = alloca double (double)*
%7 = alloca double
%8 = alloca double
%9 = alloca i32
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca i32
%16 = alloca i32
store double (double)* %0, double (double)** %6
store double %1, double* %7
store double %2, double* %8
store i32 %3, i32* %9
%17 = load i32, i32* %9
%18 = icmp eq i32 %17, 1
br i1 %18, label %19, label %30
19: ; preds = %4
%20 = load double, double* %8
%21 = load double, double* %7
%22 = fsub double %20, %21
%23 = load double (double)*, double (double)** %6
%24 = load double, double* %7
%25 = load double, double* %8
%26 = fadd double %24, %25
%27 = fmul double 5.000000e-01, %26
%28 = call double %23(double %27)
%29 = fmul double %22, %28
store double %29, double* @midpnt.s
store double %29, double* %5
br label %93
30: ; preds = %4
store i32 1, i32* %15
store i32 1, i32* %16
br label %31
31: ; preds = %39, %30
%32 = load i32, i32* %16
%33 = load i32, i32* %9
%34 = sub nsw i32 %33, 1
%35 = icmp slt i32 %32, %34
br i1 %35, label %36, label %42
36: ; preds = %31
%37 = load i32, i32* %15
%38 = mul nsw i32 %37, 3
store i32 %38, i32* %15
br label %39
39: ; preds = %36
%40 = load i32, i32* %16
%41 = add nsw i32 %40, 1
store i32 %41, i32* %16
br label %31
42: ; preds = %31
%43 = load i32, i32* %15
%44 = sitofp i32 %43 to double
store double %44, double* %11
%45 = load double, double* %8
%46 = load double, double* %7
%47 = fsub double %45, %46
%48 = load double, double* %11
%49 = fmul double 3.000000e+00, %48
%50 = fdiv double %47, %49
store double %50, double* %13
%51 = load double, double* %13
%52 = load double, double* %13
%53 = fadd double %51, %52
store double %53, double* %14
%54 = load double, double* %7
%55 = load double, double* %13
%56 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %55, double %54)
store double %56, double* %10
store double 0.000000e+00, double* %12
store i32 1, i32* %16
br label %57
57: ; preds = %78, %42
%58 = load i32, i32* %16
%59 = load i32, i32* %15
%60 = icmp sle i32 %58, %59
br i1 %60, label %61, label %81
61: ; preds = %57
%62 = load double (double)*, double (double)** %6
%63 = load double, double* %10
%64 = call double %62(double %63)
%65 = load double, double* %12
%66 = fadd double %65, %64
store double %66, double* %12
%67 = load double, double* %14
%68 = load double, double* %10
%69 = fadd double %68, %67
store double %69, double* %10
%70 = load double (double)*, double (double)** %6
%71 = load double, double* %10
%72 = call double %70(double %71)
%73 = load double, double* %12
%74 = fadd double %73, %72
store double %74, double* %12
%75 = load double, double* %13
%76 = load double, double* %10
%77 = fadd double %76, %75
store double %77, double* %10
br label %78
78: ; preds = %61
%79 = load i32, i32* %16
%80 = add nsw i32 %79, 1
store i32 %80, i32* %16
br label %57
81: ; preds = %57
%82 = load double, double* @midpnt.s
%83 = load double, double* %8
%84 = load double, double* %7
%85 = fsub double %83, %84
%86 = load double, double* %12
%87 = fmul double %85, %86
%88 = load double, double* %11
%89 = fdiv double %87, %88
%90 = fadd double %82, %89
%91 = fdiv double %90, 3.000000e+00
store double %91, double* @midpnt.s
%92 = load double, double* @midpnt.s
store double %92, double* %5
br label %93
93: ; preds = %81, %19
%94 = load double, double* %5
ret double %94
}

LLVM IR (-O2):
@midpnt.s = internal unnamed_addr global double 0.000000e+00
define double @midpnt(double (double)* nocapture %0, double %1, double %2, i32 %3) {
%5 = icmp eq i32 %3, 1
br i1 %5, label %53, label %6
6: ; preds = %4
%7 = add i32 %3, -1
%8 = icmp sgt i32 %3, 2
br i1 %8, label %9, label %65
9: ; preds = %6
%10 = add i32 %3, -2
%11 = icmp ult i32 %10, 8
br i1 %11, label %50, label %12
12: ; preds = %9
%13 = and i32 %10, -8
%14 = or i32 %13, 1
%15 = add i32 %13, -8
%16 = lshr exact i32 %15, 3
%17 = add nuw nsw i32 %16, 1
%18 = and i32 %17, 7
%19 = icmp ult i32 %15, 56
br i1 %19, label %30, label %20
20: ; preds = %12
%21 = and i32 %17, 1073741816
br label %22
22: ; preds = %22, %20
%23 = phi <4 x i32> [ <i32 1, i32 1, i32 1, i32 1>, %20 ], [ %26, %22 ]
%24 = phi <4 x i32> [ <i32 1, i32 1, i32 1, i32 1>, %20 ], [ %27, %22 ]
%25 = phi i32 [ 0, %20 ], [ %28, %22 ]
%26 = mul <4 x i32> %23, <i32 6561, i32 6561, i32 6561, i32 6561>
%27 = mul <4 x i32> %24, <i32 6561, i32 6561, i32 6561, i32 6561>
%28 = add nuw i32 %25, 8
%29 = icmp eq i32 %28, %21
br i1 %29, label %30, label %22
30: ; preds = %22, %12
%31 = phi <4 x i32> [ undef, %12 ], [ %26, %22 ]
%32 = phi <4 x i32> [ undef, %12 ], [ %27, %22 ]
%33 = phi <4 x i32> [ <i32 1, i32 1, i32 1, i32 1>, %12 ], [ %26, %22 ]
%34 = phi <4 x i32> [ <i32 1, i32 1, i32 1, i32 1>, %12 ], [ %27, %22 ]
%35 = icmp eq i32 %18, 0
br i1 %35, label %44, label %36
36: ; preds = %30, %36
%37 = phi <4 x i32> [ %40, %36 ], [ %33, %30 ]
%38 = phi <4 x i32> [ %41, %36 ], [ %34, %30 ]
%39 = phi i32 [ %42, %36 ], [ 0, %30 ]
%40 = mul <4 x i32> %37, <i32 3, i32 3, i32 3, i32 3>
%41 = mul <4 x i32> %38, <i32 3, i32 3, i32 3, i32 3>
%42 = add i32 %39, 1
%43 = icmp eq i32 %42, %18
br i1 %43, label %44, label %36
44: ; preds = %36, %30
%45 = phi <4 x i32> [ %31, %30 ], [ %40, %36 ]
%46 = phi <4 x i32> [ %32, %30 ], [ %41, %36 ]
%47 = mul <4 x i32> %46, %45
%48 = call i32 @llvm.vector.reduce.mul.v4i32(<4 x i32> %47)
%49 = icmp eq i32 %10, %13
br i1 %49, label %65, label %50
50: ; preds = %9, %44
%51 = phi i32 [ 1, %9 ], [ %14, %44 ]
%52 = phi i32 [ 1, %9 ], [ %48, %44 ]
br label %59
53: ; preds = %4
%54 = fsub double %2, %1
%55 = fadd double %1, %2
%56 = fmul double %55, 5.000000e-01
%57 = tail call double %0(double %56)
%58 = fmul double %54, %57
br label %92
59: ; preds = %50, %59
%60 = phi i32 [ %63, %59 ], [ %51, %50 ]
%61 = phi i32 [ %62, %59 ], [ %52, %50 ]
%62 = mul nsw i32 %61, 3
%63 = add nuw nsw i32 %60, 1
%64 = icmp eq i32 %63, %7
br i1 %64, label %65, label %59
65: ; preds = %59, %44, %6
%66 = phi i32 [ 1, %6 ], [ %48, %44 ], [ %62, %59 ]
%67 = sitofp i32 %66 to double
%68 = fsub double %2, %1
%69 = fmul double %67, 3.000000e+00
%70 = fdiv double %68, %69
%71 = fadd double %70, %70
%72 = tail call double @llvm.fmuladd.f64(double %70, double 5.000000e-01, double %1)
%73 = call i32 @llvm.umax.i32(i32 %66, i32 1)
br label %74
74: ; preds = %65, %74
%75 = phi i32 [ %84, %74 ], [ 1, %65 ]
%76 = phi double [ %83, %74 ], [ %72, %65 ]
%77 = phi double [ %82, %74 ], [ 0.000000e+00, %65 ]
%78 = tail call double %0(double %76)
%79 = fadd double %77, %78
%80 = fadd double %71, %76
%81 = tail call double %0(double %80)
%82 = fadd double %79, %81
%83 = fadd double %70, %80
%84 = add nuw i32 %75, 1
%85 = icmp eq i32 %75, %73
br i1 %85, label %86, label %74
86: ; preds = %74
%87 = load double, double* @midpnt.s
%88 = fmul double %68, %82
%89 = fdiv double %88, %67
%90 = fadd double %89, %87
%91 = fdiv double %90, 3.000000e+00
br label %92
92: ; preds = %86, %53
%93 = phi double [ %91, %86 ], [ %58, %53 ]
store double %93, double* @midpnt.s
ret double %93
}