#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 */
double fabs(float);
int /*<<< orphan*/ free_vector(float *, int, int);
void stub1(int, float *, float *);
float *vector(int, int);
void fdjac(int n, float x[], float fvec[], float **df,
           void (*vecfunc)(int, float[], float[])) {
  int i, j;
  float h, temp, *f;
  f = vector(1, n);
  for (j = 1; j <= n; j++) {
    temp = x[j];
    h = 1.0e-4 * fabs(temp);
    if (h == 0.0)
      h = 1.0e-4;
    x[j] = temp + h;
    h = x[j] - temp;
    (*vecfunc)(n, x, f);
    x[j] = temp;
    for (i = 1; i <= n; i++)
      df[i][j] = (f[i] - fvec[i]) / h;
  }
  free_vector(f, 1, n);
}


x86 assembly (clang, -O0):
<fdjac>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x4(%rbp),%esi
mov    $0x1,%edi
call   1050 <vector@plt>
mov    %rax,-0x40(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1285 <fdjac+0x155>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x38(%rbp),%xmm0
call   1040 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xe6a(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11cc <fdjac+0x9c>
jp     11cc <fdjac+0x9c>
movss  0xe41(%rip),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x38(%rbp),%xmm0
addss  -0x34(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x4(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x40(%rbp),%rdx
call   *%rax
movss  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1272 <fdjac+0x142>
mov    -0x40(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
divss  -0x34(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1224 <fdjac+0xf4>
jmp    1277 <fdjac+0x147>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1163 <fdjac+0x33>
mov    -0x40(%rbp),%rdi
mov    -0x4(%rbp),%edx
mov    $0x1,%esi
call   1030 <free_vector@plt>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<fdjac>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r8,0x28(%rsp)
mov    %rcx,%r13
mov    %rdx,%rbp
mov    %rsi,0x20(%rsp)
mov    %edi,%r14d
mov    $0x1,%edi
mov    %r14d,%esi
call   1050 <vector@plt>
mov    %rax,%rbx
mov    %r14d,0x8(%rsp)
test   %r14d,%r14d
jle    12b1 <fdjac+0x181>
mov    0x8(%rsp),%r14d
lea    0x1(%r14),%rdx
mov    %r14,0x18(%rsp)
and    $0xfffffffe,%r14d
neg    %r14
mov    $0x1,%r15d
mov    %rdx,0x10(%rsp)
mov    0x20(%rsp),%r12
jmp    11ad <fdjac+0x7d>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r15
cmp    %rdx,%r15
je     12b1 <fdjac+0x181>
movss  (%r12,%r15,4),%xmm0
movss  %xmm0,0xc(%rsp)
call   1040 <fabs@plt>
mulsd  0xe3a(%rip),%xmm0
cvtsd2ss %xmm0,%xmm0
movaps %xmm0,%xmm2
cmpeqss 0xe36(%rip),%xmm2
movaps %xmm2,%xmm1
andnps %xmm0,%xmm1
movss  0xe24(%rip),%xmm0
andps  %xmm0,%xmm2
orps   %xmm1,%xmm2
addss  0xc(%rsp),%xmm2
movaps %xmm2,0x30(%rsp)
movss  %xmm2,(%r12,%r15,4)
mov    0x8(%rsp),%edi
mov    %r12,%rsi
mov    %rbx,%rdx
call   *0x28(%rsp)
movaps 0x30(%rsp),%xmm1
mov    0x10(%rsp),%rdx
movss  0xc(%rsp),%xmm0
subss  %xmm0,%xmm1
movss  %xmm0,(%r12,%r15,4)
mov    $0x1,%eax
cmp    $0x2,%edx
je     1287 <fdjac+0x157>
mov    $0x1,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%rbx,%rax,4),%xmm0
subss  0x0(%rbp,%rax,4),%xmm0
divss  %xmm1,%xmm0
mov    0x0(%r13,%rax,8),%rcx
movss  %xmm0,(%rcx,%r15,4)
movss  0x4(%rbx,%rax,4),%xmm0
subss  0x4(%rbp,%rax,4),%xmm0
divss  %xmm1,%xmm0
mov    0x8(%r13,%rax,8),%rcx
movss  %xmm0,(%rcx,%r15,4)
lea    (%r14,%rax,1),%rcx
add    $0x2,%rcx
add    $0x2,%rax
cmp    $0x1,%rcx
jne    1240 <fdjac+0x110>
testb  $0x1,0x18(%rsp)
je     11a0 <fdjac+0x70>
movss  (%rbx,%rax,4),%xmm0
subss  0x0(%rbp,%rax,4),%xmm0
divss  %xmm1,%xmm0
mov    0x0(%r13,%rax,8),%rax
movss  %xmm0,(%rax,%r15,4)
jmp    11a0 <fdjac+0x70>
mov    %rbx,%rdi
mov    $0x1,%esi
mov    0x8(%rsp),%edx
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free_vector@plt>

LLVM IR (-O0):
define void @fdjac(i32 %0, float* %1, float* %2, float** %3, void (i32, float*, float*)* %4) {
%6 = alloca i32
%7 = alloca float*
%8 = alloca float*
%9 = alloca float**
%10 = alloca void (i32, float*, float*)*
%11 = alloca i32
%12 = alloca i32
%13 = alloca float
%14 = alloca float
%15 = alloca float*
store i32 %0, i32* %6
store float* %1, float** %7
store float* %2, float** %8
store float** %3, float*** %9
store void (i32, float*, float*)* %4, void (i32, float*, float*)** %10
%16 = load i32, i32* %6
%17 = call float* @vector(i32 1, i32 %16)
store float* %17, float** %15
store i32 1, i32* %12
br label %18
18: ; preds = %90, %5
%19 = load i32, i32* %12
%20 = load i32, i32* %6
%21 = icmp sle i32 %19, %20
br i1 %21, label %22, label %93
22: ; preds = %18
%23 = load float*, float** %7
%24 = load i32, i32* %12
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds float, float* %23, i64 %25
%27 = load float, float* %26
store float %27, float* %14
%28 = load float, float* %14
%29 = call double @fabs(float %28)
%30 = fmul double 1.000000e-04, %29
%31 = fptrunc double %30 to float
store float %31, float* %13
%32 = load float, float* %13
%33 = fpext float %32 to double
%34 = fcmp oeq double %33, 0.000000e+00
br i1 %34, label %35, label %36
35: ; preds = %22
store float 0x3F1A36E2E0000000, float* %13
br label %36
36: ; preds = %35, %22
%37 = load float, float* %14
%38 = load float, float* %13
%39 = fadd float %37, %38
%40 = load float*, float** %7
%41 = load i32, i32* %12
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds float, float* %40, i64 %42
store float %39, float* %43
%44 = load float*, float** %7
%45 = load i32, i32* %12
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds float, float* %44, i64 %46
%48 = load float, float* %47
%49 = load float, float* %14
%50 = fsub float %48, %49
store float %50, float* %13
%51 = load void (i32, float*, float*)*, void (i32, float*, float*)** %10
%52 = load i32, i32* %6
%53 = load float*, float** %7
%54 = load float*, float** %15
%55 = load float, float* %14
%56 = load float*, float** %7
%57 = load i32, i32* %12
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds float, float* %56, i64 %58
store float %55, float* %59
store i32 1, i32* %11
br label %60
60: ; preds = %86, %36
%61 = load i32, i32* %11
%62 = load i32, i32* %6
%63 = icmp sle i32 %61, %62
br i1 %63, label %64, label %89
64: ; preds = %60
%65 = load float*, float** %15
%66 = load i32, i32* %11
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds float, float* %65, i64 %67
%69 = load float, float* %68
%70 = load float*, float** %8
%71 = load i32, i32* %11
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds float, float* %70, i64 %72
%74 = load float, float* %73
%75 = fsub float %69, %74
%76 = load float, float* %13
%77 = fdiv float %75, %76
%78 = load float**, float*** %9
%79 = load i32, i32* %11
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds float*, float** %78, i64 %80
%82 = load float*, float** %81
%83 = load i32, i32* %12
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float, float* %82, i64 %84
store float %77, float* %85
br label %86
86: ; preds = %64
%87 = load i32, i32* %11
%88 = add nsw i32 %87, 1
store i32 %88, i32* %11
br label %60
89: ; preds = %60
br label %90
90: ; preds = %89
%91 = load i32, i32* %12
%92 = add nsw i32 %91, 1
store i32 %92, i32* %12
br label %18
93: ; preds = %18
%94 = load float*, float** %15
%95 = load i32, i32* %6
%96 = call i32 @free_vector(float* %94, i32 1, i32 %95)
ret void
}

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