void sfdif_cent_jac_approx(void (*func)(float *p, float *hx, int m, int n,
                                        void *adata),
                           float *p, float *hxm, float *hxp, float delta,
                           float *jac, int m, int n, void *adata) {
  register int i, j;
  float tmp;
  register float d;
  for (j = 0; j < m; ++j) {
    d = 1E-04F * p[j];
    d = (((d) >= 0.0) ? (d) : -(d));
    if (d < delta)
      d = delta;
    tmp = p[j];
    p[j] -= d;
    (*func)(p, hxm, m, n, adata);
    p[j] = tmp + d;
    (*func)(p, hxp, m, n, adata);
    p[j] = tmp;
    d = 0.5F / d;
    for (i = 0; i < n; ++i) {
      jac[i * m + j] = (hxp[i] - hxm[i]) * d;
    }
  }
}


x86 assembly (clang, -O0):
<sfdif_cent_jac_approx>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movss  %xmm0,-0x24(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9d,-0x34(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    12bd <sfdif_cent_jac_approx+0x1bd>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  0xeb1(%rip),%xmm0
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x44(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     117e <sfdif_cent_jac_approx+0x7e>
movss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
jmp    1195 <sfdif_cent_jac_approx+0x95>
movss  -0x44(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x44(%rbp),%xmm1
movss  -0x24(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    11bc <sfdif_cent_jac_approx+0xbc>
movss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x44(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x34(%rbp),%edx
mov    0x10(%rbp),%ecx
mov    0x18(%rbp),%r8
call   *%rax
movss  -0x40(%rbp),%xmm0
addss  -0x44(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x34(%rbp),%edx
mov    0x10(%rbp),%ecx
mov    0x18(%rbp),%r8
call   *%rax
movss  -0x40(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  0xdba(%rip),%xmm0
divss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    12aa <sfdif_cent_jac_approx+0x1aa>
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
mulss  -0x44(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x34(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    125b <sfdif_cent_jac_approx+0x15b>
jmp    12af <sfdif_cent_jac_approx+0x1af>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1133 <sfdif_cent_jac_approx+0x33>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sfdif_cent_jac_approx>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r8,0x10(%rsp)
movss  %xmm0,0xc(%rsp)
mov    %rsi,0x38(%rsp)
mov    %rdi,0x30(%rsp)
mov    %r9d,(%rsp)
test   %r9d,%r9d
jle    12e9 <sfdif_cent_jac_approx+0x1e9>
mov    %rcx,%rbx
mov    %rdx,%rbp
mov    (%rsp),%edx
mov    0x80(%rsp),%eax
mov    %eax,%esi
and    $0xfffffffe,%esi
lea    0x0(,%rdx,8),%rdi
xor    %r12d,%r12d
mov    0x10(%rsp),%r8
mov    %rdx,0x28(%rsp)
mov    %rsi,0x20(%rsp)
mov    %rdi,0x18(%rsp)
mov    0x38(%rsp),%r15
mov    0x30(%rsp),%r14
jmp    1181 <sfdif_cent_jac_approx+0x81>
add    $0x1,%r12
add    $0x4,%r8
cmp    %rdx,%r12
je     12e9 <sfdif_cent_jac_approx+0x1e9>
mov    %r8,0x40(%rsp)
movss  (%r15,%r12,4),%xmm1
movaps %xmm1,%xmm0
movaps %xmm1,%xmm3
movss  %xmm1,0x8(%rsp)
mulss  0xe60(%rip),%xmm0
movaps %xmm0,%xmm1
xorps  0xe66(%rip),%xmm1
xorps  %xmm2,%xmm2
cmpless %xmm0,%xmm2
andps  %xmm2,%xmm0
andnps %xmm1,%xmm2
orps   %xmm0,%xmm2
movss  0xc(%rsp),%xmm1
maxss  %xmm2,%xmm1
movss  %xmm1,0x4(%rsp)
movaps %xmm3,%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%r15,%r12,4)
mov    %r15,%rdi
mov    %rbp,%rsi
mov    (%rsp),%r13d
mov    %r13d,%edx
mov    0x80(%rsp),%ecx
mov    0x88(%rsp),%r8
call   *%r14
movss  0x8(%rsp),%xmm0
addss  0x4(%rsp),%xmm0
movss  %xmm0,(%r15,%r12,4)
mov    %r15,%rdi
mov    %rbx,%rsi
mov    %r13d,%edx
mov    0x80(%rsp),%r13d
mov    %r13d,%ecx
mov    0x88(%rsp),%r8
call   *%r14
mov    0x40(%rsp),%r8
mov    0x18(%rsp),%rdi
mov    0x20(%rsp),%rsi
mov    0x28(%rsp),%rdx
movss  0x8(%rsp),%xmm0
movss  %xmm0,(%r15,%r12,4)
test   %r13d,%r13d
jle    1170 <sfdif_cent_jac_approx+0x70>
movss  0xdab(%rip),%xmm0
divss  0x4(%rsp),%xmm0
cmpl   $0x1,0x80(%rsp)
01
jne    1270 <sfdif_cent_jac_approx+0x170>
xor    %eax,%eax
jmp    12b4 <sfdif_cent_jac_approx+0x1b4>
nopl   (%rax)
mov    %r8,%rcx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%rbx,%rax,4),%xmm1
subss  0x0(%rbp,%rax,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%rcx)
movss  0x4(%rbx,%rax,4),%xmm1
subss  0x4(%rbp,%rax,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%rcx,%rdx,4)
add    $0x2,%rax
add    %rdi,%rcx
cmp    %rax,%rsi
jne    1280 <sfdif_cent_jac_approx+0x180>
mov    0x80(%rsp),%ecx
test   $0x1,%cl
je     1170 <sfdif_cent_jac_approx+0x70>
movss  (%rbx,%rax,4),%xmm1
subss  0x0(%rbp,%rax,4),%xmm1
mulss  %xmm1,%xmm0
imul   %rdx,%rax
add    %r12,%rax
mov    0x10(%rsp),%rcx
movss  %xmm0,(%rcx,%rax,4)
jmp    1170 <sfdif_cent_jac_approx+0x70>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @sfdif_cent_jac_approx(void (float*, float*, i32, i32, i8*)* %0, float* %1, float* %2, float* %3, float %4, float* %5, i32 %6, i32 %7, i8* %8) {
%10 = alloca void (float*, float*, i32, i32, i8*)*
%11 = alloca float*
%12 = alloca float*
%13 = alloca float*
%14 = alloca float
%15 = alloca float*
%16 = alloca i32
%17 = alloca i32
%18 = alloca i8*
%19 = alloca i32
%20 = alloca i32
%21 = alloca float
%22 = alloca float
store void (float*, float*, i32, i32, i8*)* %0, void (float*, float*, i32, i32, i8*)** %10
store float* %1, float** %11
store float* %2, float** %12
store float* %3, float** %13
store float %4, float* %14
store float* %5, float** %15
store i32 %6, i32* %16
store i32 %7, i32* %17
store i8* %8, i8** %18
store i32 0, i32* %20
br label %23
23: ; preds = %118, %9
%24 = load i32, i32* %20
%25 = load i32, i32* %16
%26 = icmp slt i32 %24, %25
br i1 %26, label %27, label %121
27: ; preds = %23
%28 = load float*, float** %11
%29 = load i32, i32* %20
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds float, float* %28, i64 %30
%32 = load float, float* %31
%33 = fmul float 0x3F1A36E2E0000000, %32
store float %33, float* %22
%34 = load float, float* %22
%35 = fpext float %34 to double
%36 = fcmp oge double %35, 0.000000e+00
br i1 %36, label %37, label %39
37: ; preds = %27
%38 = load float, float* %22
br label %42
39: ; preds = %27
%40 = load float, float* %22
%41 = fneg float %40
br label %42
42: ; preds = %39, %37
%43 = phi float [ %38, %37 ], [ %41, %39 ]
store float %43, float* %22
%44 = load float, float* %22
%45 = load float, float* %14
%46 = fcmp olt float %44, %45
br i1 %46, label %47, label %49
47: ; preds = %42
%48 = load float, float* %14
store float %48, float* %22
br label %49
49: ; preds = %47, %42
%50 = load float*, float** %11
%51 = load i32, i32* %20
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds float, float* %50, i64 %52
%54 = load float, float* %53
store float %54, float* %21
%55 = load float, float* %22
%56 = load float*, float** %11
%57 = load i32, i32* %20
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds float, float* %56, i64 %58
%60 = load float, float* %59
%61 = fsub float %60, %55
store float %61, float* %59
%62 = load void (float*, float*, i32, i32, i8*)*, void (float*, float*, i32, i32, i8*)** %10
%63 = load float*, float** %11
%64 = load float*, float** %12
%65 = load i32, i32* %16
%66 = load i32, i32* %17
%67 = load i8*, i8** %18
%68 = load float, float* %21
%69 = load float, float* %22
%70 = fadd float %68, %69
%71 = load float*, float** %11
%72 = load i32, i32* %20
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds float, float* %71, i64 %73
store float %70, float* %74
%75 = load void (float*, float*, i32, i32, i8*)*, void (float*, float*, i32, i32, i8*)** %10
%76 = load float*, float** %11
%77 = load float*, float** %13
%78 = load i32, i32* %16
%79 = load i32, i32* %17
%80 = load i8*, i8** %18
%81 = load float, float* %21
%82 = load float*, float** %11
%83 = load i32, i32* %20
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float, float* %82, i64 %84
store float %81, float* %85
%86 = load float, float* %22
%87 = fdiv float 5.000000e-01, %86
store float %87, float* %22
store i32 0, i32* %19
br label %88
88: ; preds = %114, %49
%89 = load i32, i32* %19
%90 = load i32, i32* %17
%91 = icmp slt i32 %89, %90
br i1 %91, label %92, label %117
92: ; preds = %88
%93 = load float*, float** %13
%94 = load i32, i32* %19
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds float, float* %93, i64 %95
%97 = load float, float* %96
%98 = load float*, float** %12
%99 = load i32, i32* %19
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds float, float* %98, i64 %100
%102 = load float, float* %101
%103 = fsub float %97, %102
%104 = load float, float* %22
%105 = fmul float %103, %104
%106 = load float*, float** %15
%107 = load i32, i32* %19
%108 = load i32, i32* %16
%109 = mul nsw i32 %107, %108
%110 = load i32, i32* %20
%111 = add nsw i32 %109, %110
%112 = sext i32 %111 to i64
%113 = getelementptr inbounds float, float* %106, i64 %112
store float %105, float* %113
br label %114
114: ; preds = %92
%115 = load i32, i32* %19
%116 = add nsw i32 %115, 1
store i32 %116, i32* %19
br label %88
117: ; preds = %88
br label %118
118: ; preds = %117
%119 = load i32, i32* %20
%120 = add nsw i32 %119, 1
store i32 %120, i32* %20
br label %23
121: ; preds = %23
ret void
}

LLVM IR (-O2):
define void @sfdif_cent_jac_approx(void (float*, float*, i32, i32, i8*)* nocapture %0, float* %1, float* %2, float* %3, float %4, float* nocapture writeonly %5, i32 %6, i32 %7, i8* %8) {
%10 = icmp sgt i32 %6, 0
br i1 %10, label %11, label %74
11: ; preds = %9
%12 = icmp sgt i32 %7, 0
%13 = zext i32 %6 to i64
%14 = zext i32 %6 to i64
%15 = zext i32 %7 to i64
%16 = and i64 %15, 1
%17 = icmp eq i32 %7, 1
%18 = and i64 %15, 4294967294
%19 = icmp eq i64 %16, 0
br label %20
20: ; preds = %11, %71
%21 = phi i64 [ 0, %11 ], [ %72, %71 ]
%22 = getelementptr inbounds float, float* %1, i64 %21
%23 = load float, float* %22
%24 = fmul float %23, 0x3F1A36E2E0000000
%25 = fcmp oge float %24, 0.000000e+00
%26 = fneg float %24
%27 = select i1 %25, float %24, float %26
%28 = fcmp olt float %27, %4
%29 = select i1 %28, float %4, float %27
%30 = fsub float %23, %29
store float %30, float* %22
tail call void %0(float* %1, float* %2, i32 %6, i32 %7, i8* %8)
%31 = fadd float %23, %29
store float %31, float* %22
tail call void %0(float* %1, float* %3, i32 %6, i32 %7, i8* %8)
store float %23, float* %22
%32 = fdiv float 5.000000e-01, %29
br i1 %12, label %33, label %71
33: ; preds = %20
br i1 %17, label %59, label %34
34: ; preds = %33, %34
%35 = phi i64 [ %56, %34 ], [ 0, %33 ]
%36 = phi i64 [ %57, %34 ], [ 0, %33 ]
%37 = getelementptr inbounds float, float* %3, i64 %35
%38 = load float, float* %37
%39 = getelementptr inbounds float, float* %2, i64 %35
%40 = load float, float* %39
%41 = fsub float %38, %40
%42 = fmul float %32, %41
%43 = mul nsw i64 %35, %13
%44 = add nuw nsw i64 %43, %21
%45 = getelementptr inbounds float, float* %5, i64 %44
store float %42, float* %45
%46 = or i64 %35, 1
%47 = getelementptr inbounds float, float* %3, i64 %46
%48 = load float, float* %47
%49 = getelementptr inbounds float, float* %2, i64 %46
%50 = load float, float* %49
%51 = fsub float %48, %50
%52 = fmul float %32, %51
%53 = mul nsw i64 %46, %13
%54 = add nuw nsw i64 %53, %21
%55 = getelementptr inbounds float, float* %5, i64 %54
store float %52, float* %55
%56 = add nuw nsw i64 %35, 2
%57 = add i64 %36, 2
%58 = icmp eq i64 %57, %18
br i1 %58, label %59, label %34
59: ; preds = %34, %33
%60 = phi i64 [ 0, %33 ], [ %56, %34 ]
br i1 %19, label %71, label %61
61: ; preds = %59
%62 = getelementptr inbounds float, float* %3, i64 %60
%63 = load float, float* %62
%64 = getelementptr inbounds float, float* %2, i64 %60
%65 = load float, float* %64
%66 = fsub float %63, %65
%67 = fmul float %32, %66
%68 = mul nsw i64 %60, %13
%69 = add nuw nsw i64 %68, %21
%70 = getelementptr inbounds float, float* %5, i64 %69
store float %67, float* %70
br label %71
71: ; preds = %61, %59, %20
%72 = add nuw nsw i64 %21, 1
%73 = icmp eq i64 %72, %14
br i1 %73, label %74, label %20
74: ; preds = %71, %9
ret void
}