void jacpowell(double *p, double *jac, int m, int n, void *data) {
  register int i, j;
  for (i = j = 0; i < n; i += 2) {
    jac[j++] = 1.0;
    jac[j++] = 0.0;
    jac[j++] = 1.0 / ((p[0] + 0.1) * (p[0] + 0.1));
    jac[j++] = 4.0 * p[1];
  }
}


x86 assembly (clang, -O0):
<jacpowell>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11e6 <jacpowell+0xe6>
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x28(%rbp)
movslq %ecx,%rcx
movsd  0xebe(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x28(%rbp)
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movsd  0xe9b(%rip),%xmm1
addsd  (%rax),%xmm1
mov    -0x8(%rbp),%rax
movsd  0xe8b(%rip),%xmm0
addsd  (%rax),%xmm0
mulsd  %xmm0,%xmm1
movsd  0xe73(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x28(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movsd  0xe44(%rip),%xmm0
mulsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x28(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x24(%rbp)
jmp    1124 <jacpowell+0x24>
pop    %rbp
ret

x86 assembly (clang, -O2):
<jacpowell>:
test   %ecx,%ecx
jle    11ad <jacpowell+0xad>
lea    -0x1(%rcx),%r10d
cmp    $0x2,%r10d
jb     1135 <jacpowell+0x35>
shr    %r10d
lea    0x10(%rdi),%rax
cmp    %rsi,%rax
jbe    11ae <jacpowell+0xae>
mov    %r10,%rax
shl    $0x5,%rax
add    %rsi,%rax
add    $0x20,%rax
cmp    %rdi,%rax
jbe    11ae <jacpowell+0xae>
xor    %r8d,%r8d
xor    %eax,%eax
lea    (%rsi,%r8,8),%rdx
add    $0x18,%rdx
movabs $0x3ff0000000000000,%rsi
movsd  0xeac(%rip),%xmm0
movsd  0xeac(%rip),%xmm1
movsd  0xeac(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %rsi,-0x18(%rdx)
movq   $0x0,-0x10(%rdx)
movsd  (%rdi),%xmm3
addsd  %xmm0,%xmm3
mulsd  %xmm3,%xmm3
movapd %xmm1,%xmm4
divsd  %xmm3,%xmm4
movsd  %xmm4,-0x8(%rdx)
movsd  0x8(%rdi),%xmm3
mulsd  %xmm2,%xmm3
movsd  %xmm3,(%rdx)
add    $0x2,%eax
add    $0x20,%rdx
cmp    %ecx,%eax
jl     1170 <jacpowell+0x70>
ret
add    $0x1,%r10d
mov    %r10d,%r9d
and    $0xfffffffe,%r9d
lea    0x0(,%r9,4),%r8
lea    (%r9,%r9,1),%eax
movsd  0xe33(%rip),%xmm0
movsd  0xe33(%rip),%xmm1
movsd  0xe33(%rip),%xmm2
movsd  0xe23(%rip),%xmm3
mov    %r9,%r11
mov    %rsi,%rdx
nopl   0x0(%rax,%rax,1)
movsd  (%rdi),%xmm4
addsd  %xmm0,%xmm4
mulsd  %xmm4,%xmm4
movapd %xmm1,%xmm5
divsd  %xmm4,%xmm5
movsd  0x8(%rdi),%xmm4
mulsd  %xmm2,%xmm4
unpcklpd %xmm4,%xmm5
movupd %xmm3,0x20(%rdx)
movupd %xmm3,(%rdx)
movupd %xmm5,0x30(%rdx)
movupd %xmm5,0x10(%rdx)
add    $0x40,%rdx
add    $0xfffffffffffffffe,%r11
jne    11f0 <jacpowell+0xf0>
cmp    %r10,%r9
jne    113a <jacpowell+0x3a>
jmp    11ad <jacpowell+0xad>

LLVM IR (-O0):
define void @jacpowell(double* %0, double* %1, i32 %2, i32 %3, i8* %4) {
%6 = alloca double*
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i8*
%11 = alloca i32
%12 = alloca i32
store double* %0, double** %6
store double* %1, double** %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store i8* %4, i8** %10
store i32 0, i32* %12
store i32 0, i32* %11
br label %13
13: ; preds = %52, %5
%14 = load i32, i32* %11
%15 = load i32, i32* %9
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %55
17: ; preds = %13
%18 = load double*, double** %7
%19 = load i32, i32* %12
%20 = add nsw i32 %19, 1
store i32 %20, i32* %12
%21 = sext i32 %19 to i64
%22 = getelementptr inbounds double, double* %18, i64 %21
store double 1.000000e+00, double* %22
%23 = load double*, double** %7
%24 = load i32, i32* %12
%25 = add nsw i32 %24, 1
store i32 %25, i32* %12
%26 = sext i32 %24 to i64
%27 = getelementptr inbounds double, double* %23, i64 %26
store double 0.000000e+00, double* %27
%28 = load double*, double** %6
%29 = getelementptr inbounds double, double* %28, i64 0
%30 = load double, double* %29
%31 = fadd double %30, 1.000000e-01
%32 = load double*, double** %6
%33 = getelementptr inbounds double, double* %32, i64 0
%34 = load double, double* %33
%35 = fadd double %34, 1.000000e-01
%36 = fmul double %31, %35
%37 = fdiv double 1.000000e+00, %36
%38 = load double*, double** %7
%39 = load i32, i32* %12
%40 = add nsw i32 %39, 1
store i32 %40, i32* %12
%41 = sext i32 %39 to i64
%42 = getelementptr inbounds double, double* %38, i64 %41
store double %37, double* %42
%43 = load double*, double** %6
%44 = getelementptr inbounds double, double* %43, i64 1
%45 = load double, double* %44
%46 = fmul double 4.000000e+00, %45
%47 = load double*, double** %7
%48 = load i32, i32* %12
%49 = add nsw i32 %48, 1
store i32 %49, i32* %12
%50 = sext i32 %48 to i64
%51 = getelementptr inbounds double, double* %47, i64 %50
store double %46, double* %51
br label %52
52: ; preds = %17
%53 = load i32, i32* %11
%54 = add nsw i32 %53, 2
store i32 %54, i32* %11
br label %13
55: ; preds = %13
ret void
}

LLVM IR (-O2):
define void @jacpowell(double* nocapture %0, double* nocapture writeonly %1, i32 %2, i32 %3, i8* nocapture readnone %4) {
%6 = icmp sgt i32 %3, 0
br i1 %6, label %7, label %75
7: ; preds = %5
%8 = getelementptr inbounds double, double* %0, i64 1
%9 = add i32 %3, -1
%10 = lshr i32 %9, 1
%11 = add nuw i32 %10, 1
%12 = zext i32 %11 to i64
%13 = icmp ult i32 %9, 2
br i1 %13, label %53, label %14
14: ; preds = %7
%15 = add i32 %3, -1
%16 = lshr i32 %15, 1
%17 = zext i32 %16 to i64
%18 = shl nuw nsw i64 %17, 2
%19 = add nuw nsw i64 %18, 4
%20 = getelementptr double, double* %1, i64 %19
%21 = getelementptr double, double* %0, i64 2
%22 = icmp ugt double* %21, %1
%23 = icmp ugt double* %20, %0
%24 = and i1 %22, %23
br i1 %24, label %53, label %25
25: ; preds = %14
%26 = and i64 %12, 4294967294
%27 = shl nuw nsw i64 %26, 2
%28 = trunc i64 %26 to i32
%29 = shl i32 %28, 1
%30 = getelementptr inbounds double, double* %1, i64 -3
br label %31
31: ; preds = %31, %25
%32 = phi i64 [ 0, %25 ], [ %49, %31 ]
%33 = shl i64 %32, 2
%34 = load double, double* %0
%35 = insertelement <2 x double> poison, double %34, i64 0
%36 = fadd <2 x double> %35, <double 1.000000e-01, double poison>
%37 = fmul <2 x double> %36, %36
%38 = fdiv <2 x double> <double 1.000000e+00, double poison>, %37
%39 = shufflevector <2 x double> %38, <2 x double> poison, <2 x i32> zeroinitializer
%40 = or i64 %33, 3
%41 = load double, double* %8
%42 = insertelement <2 x double> poison, double %41, i64 0
%43 = fmul <2 x double> %42, <double 4.000000e+00, double poison>
%44 = shufflevector <2 x double> %43, <2 x double> poison, <2 x i32> zeroinitializer
%45 = getelementptr inbounds double, double* %30, i64 %40
%46 = bitcast double* %45 to <8 x double>*
%47 = shufflevector <2 x double> %39, <2 x double> %44, <4 x i32> <i32 0, i32 1, i32 2, i32 3>
%48 = shufflevector <4 x double> <double 1.000000e+00, double 1.000000e+00, double 0.000000e+00, double 0.000000e+00>, <4 x double> %47, <8 x i32> <i32 0, i32 2, i32 4, i32 6, i32 1, i32 3, i32 5, i32 7>
store <8 x double> %48, <8 x double>* %46
%49 = add nuw i64 %32, 2
%50 = icmp eq i64 %49, %26
br i1 %50, label %51, label %31
51: ; preds = %31
%52 = icmp eq i64 %26, %12
br i1 %52, label %75, label %53
53: ; preds = %14, %7, %51
%54 = phi i64 [ 0, %14 ], [ 0, %7 ], [ %27, %51 ]
%55 = phi i32 [ 0, %14 ], [ 0, %7 ], [ %29, %51 ]
br label %56
56: ; preds = %53, %56
%57 = phi i64 [ %71, %56 ], [ %54, %53 ]
%58 = phi i32 [ %73, %56 ], [ %55, %53 ]
%59 = or i64 %57, 1
%60 = getelementptr inbounds double, double* %1, i64 %57
store double 1.000000e+00, double* %60
%61 = or i64 %57, 2
%62 = getelementptr inbounds double, double* %1, i64 %59
store double 0.000000e+00, double* %62
%63 = load double, double* %0
%64 = fadd double %63, 1.000000e-01
%65 = fmul double %64, %64
%66 = fdiv double 1.000000e+00, %65
%67 = or i64 %57, 3
%68 = getelementptr inbounds double, double* %1, i64 %61
store double %66, double* %68
%69 = load double, double* %8
%70 = fmul double %69, 4.000000e+00
%71 = add nuw nsw i64 %57, 4
%72 = getelementptr inbounds double, double* %1, i64 %67
store double %70, double* %72
%73 = add nuw nsw i32 %58, 2
%74 = icmp slt i32 %73, %3
br i1 %74, label %56, label %75
75: ; preds = %56, %51, %5
ret void
}