void jachs51(double *p, double *jac, int m, int n, void *data) {
  register int i, j;
  double t1, t2, t3, t4;
  t1 = p[0] - p[1];
  t2 = p[1] + p[2] - 2.0;
  t3 = p[3] - 1.0;
  t4 = p[4] - 1.0;
  for (i = j = 0; i < n; ++i) {
    jac[j++] = 2.0 * t1;
    jac[j++] = 2.0 * (t2 - t1);
    jac[j++] = 2.0 * t2;
    jac[j++] = 2.0 * t3;
    jac[j++] = 2.0 * t4;
  }
}


x86 assembly (clang, -O0):
<jachs51>:
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)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
subsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x8(%rbp),%rax
addsd  0x10(%rax),%xmm0
movsd  0xec2(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x18(%rax),%xmm0
movsd  0xea0(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x20(%rax),%xmm0
movsd  0xe86(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1268 <jachs51+0x168>
movsd  0xe63(%rip),%xmm0
mulsd  -0x30(%rbp),%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)
movsd  -0x38(%rbp),%xmm1
subsd  -0x30(%rbp),%xmm1
movsd  0xe35(%rip),%xmm0
mulsd  %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)
movsd  0xe12(%rip),%xmm0
mulsd  -0x38(%rbp),%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)
movsd  0xdee(%rip),%xmm0
mulsd  -0x40(%rbp),%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)
movsd  0xdca(%rip),%xmm0
mulsd  -0x48(%rbp),%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    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1191 <jachs51+0x91>
pop    %rbp
ret

x86 assembly (clang, -O2):
<jachs51>:
test   %ecx,%ecx
jle    11b2 <jachs51+0xb2>
movsd  (%rdi),%xmm1
movsd  0x8(%rdi),%xmm0
subsd  %xmm0,%xmm1
addsd  0x10(%rdi),%xmm0
addsd  0xede(%rip),%xmm0
movsd  0xede(%rip),%xmm3
movsd  0x20(%rdi),%xmm2
addsd  %xmm3,%xmm2
addsd  0x18(%rdi),%xmm3
movapd %xmm0,%xmm4
subsd  %xmm1,%xmm4
unpcklpd %xmm4,%xmm1
addpd  %xmm1,%xmm1
unpcklpd %xmm3,%xmm0
addpd  %xmm0,%xmm0
addsd  %xmm2,%xmm2
cmp    $0x1,%ecx
jne    115d <jachs51+0x5d>
xor    %eax,%eax
jmp    119c <jachs51+0x9c>
mov    %ecx,%edx
and    $0xfffffffe,%edx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd %xmm1,(%rsi,%rax,8)
movupd %xmm0,0x10(%rsi,%rax,8)
movsd  %xmm2,0x20(%rsi,%rax,8)
movupd %xmm1,0x28(%rsi,%rax,8)
movupd %xmm0,0x38(%rsi,%rax,8)
movsd  %xmm2,0x48(%rsi,%rax,8)
add    $0xa,%rax
add    $0xfffffffe,%edx
jne    1170 <jachs51+0x70>
test   $0x1,%cl
je     11b2 <jachs51+0xb2>
movupd %xmm1,(%rsi,%rax,8)
movupd %xmm0,0x10(%rsi,%rax,8)
movsd  %xmm2,0x20(%rsi,%rax,8)
ret

LLVM IR (-O0):
define void @jachs51(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
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double
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
%17 = load double*, double** %6
%18 = getelementptr inbounds double, double* %17, i64 0
%19 = load double, double* %18
%20 = load double*, double** %6
%21 = getelementptr inbounds double, double* %20, i64 1
%22 = load double, double* %21
%23 = fsub double %19, %22
store double %23, double* %13
%24 = load double*, double** %6
%25 = getelementptr inbounds double, double* %24, i64 1
%26 = load double, double* %25
%27 = load double*, double** %6
%28 = getelementptr inbounds double, double* %27, i64 2
%29 = load double, double* %28
%30 = fadd double %26, %29
%31 = fsub double %30, 2.000000e+00
store double %31, double* %14
%32 = load double*, double** %6
%33 = getelementptr inbounds double, double* %32, i64 3
%34 = load double, double* %33
%35 = fsub double %34, 1.000000e+00
store double %35, double* %15
%36 = load double*, double** %6
%37 = getelementptr inbounds double, double* %36, i64 4
%38 = load double, double* %37
%39 = fsub double %38, 1.000000e+00
store double %39, double* %16
store i32 0, i32* %12
store i32 0, i32* %11
br label %40
40: ; preds = %82, %5
%41 = load i32, i32* %11
%42 = load i32, i32* %9
%43 = icmp slt i32 %41, %42
br i1 %43, label %44, label %85
44: ; preds = %40
%45 = load double, double* %13
%46 = fmul double 2.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
%52 = load double, double* %14
%53 = load double, double* %13
%54 = fsub double %52, %53
%55 = fmul double 2.000000e+00, %54
%56 = load double*, double** %7
%57 = load i32, i32* %12
%58 = add nsw i32 %57, 1
store i32 %58, i32* %12
%59 = sext i32 %57 to i64
%60 = getelementptr inbounds double, double* %56, i64 %59
store double %55, double* %60
%61 = load double, double* %14
%62 = fmul double 2.000000e+00, %61
%63 = load double*, double** %7
%64 = load i32, i32* %12
%65 = add nsw i32 %64, 1
store i32 %65, i32* %12
%66 = sext i32 %64 to i64
%67 = getelementptr inbounds double, double* %63, i64 %66
store double %62, double* %67
%68 = load double, double* %15
%69 = fmul double 2.000000e+00, %68
%70 = load double*, double** %7
%71 = load i32, i32* %12
%72 = add nsw i32 %71, 1
store i32 %72, i32* %12
%73 = sext i32 %71 to i64
%74 = getelementptr inbounds double, double* %70, i64 %73
store double %69, double* %74
%75 = load double, double* %16
%76 = fmul double 2.000000e+00, %75
%77 = load double*, double** %7
%78 = load i32, i32* %12
%79 = add nsw i32 %78, 1
store i32 %79, i32* %12
%80 = sext i32 %78 to i64
%81 = getelementptr inbounds double, double* %77, i64 %80
store double %76, double* %81
br label %82
82: ; preds = %44
%83 = load i32, i32* %11
%84 = add nsw i32 %83, 1
store i32 %84, i32* %11
br label %40
85: ; preds = %40
ret void
}

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