void hs51(double *p, double *x, int m, int n, void *data) {
  register int i;
  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 = 0; i < n; ++i)
    x[i] = t1 * t1 + t2 * t2 + t3 * t3 + t4 * t4;
}


x86 assembly (clang, -O0):
<hs51>:
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,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11f1 <hs51+0xf1>
movsd  -0x30(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x48(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    118a <hs51+0x8a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<hs51>:
test   %ecx,%ecx
jle    126e <hs51+0x16e>
movsd  0x8(%rdi),%xmm1
movsd  0xeeb(%rip),%xmm2
movsd  0x20(%rdi),%xmm0
addsd  %xmm2,%xmm0
addsd  0x18(%rdi),%xmm2
movsd  0x10(%rdi),%xmm3
addsd  %xmm1,%xmm3
addsd  0xed4(%rip),%xmm3
movsd  (%rdi),%xmm4
subsd  %xmm1,%xmm4
mulsd  %xmm3,%xmm3
mulsd  %xmm4,%xmm4
addsd  %xmm3,%xmm4
mulsd  %xmm2,%xmm2
addsd  %xmm4,%xmm2
mulsd  %xmm0,%xmm0
addsd  %xmm2,%xmm0
mov    %ecx,%eax
cmp    $0x4,%ecx
jae    1166 <hs51+0x66>
xor    %ecx,%ecx
jmp    1260 <hs51+0x160>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%rcx),%rdx
mov    %rdx,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
mov    %edi,%r8d
and    $0x7,%r8d
cmp    $0x1c,%rdx
jae    1196 <hs51+0x96>
xor    %edx,%edx
jmp    1225 <hs51+0x125>
and    $0xfffffffffffffff8,%rdi
xor    %edx,%edx
nopl   0x0(%rax)
movupd %xmm1,(%rsi,%rdx,8)
movupd %xmm1,0x10(%rsi,%rdx,8)
movupd %xmm1,0x20(%rsi,%rdx,8)
movupd %xmm1,0x30(%rsi,%rdx,8)
movupd %xmm1,0x40(%rsi,%rdx,8)
movupd %xmm1,0x50(%rsi,%rdx,8)
movupd %xmm1,0x60(%rsi,%rdx,8)
movupd %xmm1,0x70(%rsi,%rdx,8)
movupd %xmm1,0x80(%rsi,%rdx,8)
movupd %xmm1,0x90(%rsi,%rdx,8)
movupd %xmm1,0xa0(%rsi,%rdx,8)
movupd %xmm1,0xb0(%rsi,%rdx,8)
movupd %xmm1,0xc0(%rsi,%rdx,8)
movupd %xmm1,0xd0(%rsi,%rdx,8)
movupd %xmm1,0xe0(%rsi,%rdx,8)
movupd %xmm1,0xf0(%rsi,%rdx,8)
add    $0x20,%rdx
add    $0xfffffffffffffff8,%rdi
jne    11a0 <hs51+0xa0>
test   %r8,%r8
je     1254 <hs51+0x154>
lea    (%rsi,%rdx,8),%rdx
add    $0x10,%rdx
shl    $0x5,%r8
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
movupd %xmm1,-0x10(%rdx,%rdi,1)
movupd %xmm1,(%rdx,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%r8
jne    1240 <hs51+0x140>
cmp    %rax,%rcx
je     126e <hs51+0x16e>
nopl   0x0(%rax)
movsd  %xmm0,(%rsi,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1260 <hs51+0x160>
ret

LLVM IR (-O0):
define void @hs51(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 double
%13 = alloca double
%14 = alloca double
%15 = 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
%16 = load double*, double** %6
%17 = getelementptr inbounds double, double* %16, i64 0
%18 = load double, double* %17
%19 = load double*, double** %6
%20 = getelementptr inbounds double, double* %19, i64 1
%21 = load double, double* %20
%22 = fsub double %18, %21
store double %22, double* %12
%23 = load double*, double** %6
%24 = getelementptr inbounds double, double* %23, i64 1
%25 = load double, double* %24
%26 = load double*, double** %6
%27 = getelementptr inbounds double, double* %26, i64 2
%28 = load double, double* %27
%29 = fadd double %25, %28
%30 = fsub double %29, 2.000000e+00
store double %30, double* %13
%31 = load double*, double** %6
%32 = getelementptr inbounds double, double* %31, i64 3
%33 = load double, double* %32
%34 = fsub double %33, 1.000000e+00
store double %34, double* %14
%35 = load double*, double** %6
%36 = getelementptr inbounds double, double* %35, i64 4
%37 = load double, double* %36
%38 = fsub double %37, 1.000000e+00
store double %38, double* %15
store i32 0, i32* %11
br label %39
39: ; preds = %60, %5
%40 = load i32, i32* %11
%41 = load i32, i32* %9
%42 = icmp slt i32 %40, %41
br i1 %42, label %43, label %63
43: ; preds = %39
%44 = load double, double* %12
%45 = load double, double* %12
%46 = load double, double* %13
%47 = load double, double* %13
%48 = fmul double %46, %47
%49 = call double @llvm.fmuladd.f64(double %44, double %45, double %48)
%50 = load double, double* %14
%51 = load double, double* %14
%52 = call double @llvm.fmuladd.f64(double %50, double %51, double %49)
%53 = load double, double* %15
%54 = load double, double* %15
%55 = call double @llvm.fmuladd.f64(double %53, double %54, double %52)
%56 = load double*, double** %7
%57 = load i32, i32* %11
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %56, i64 %58
store double %55, double* %59
br label %60
60: ; preds = %43
%61 = load i32, i32* %11
%62 = add nsw i32 %61, 1
store i32 %62, i32* %11
br label %39
63: ; preds = %39
ret void
}

LLVM IR (-O2):
define void @hs51(double* nocapture %0, double* nocapture writeonly %1, i32 %2, i32 %3, i8* nocapture readnone %4) {
%6 = getelementptr inbounds double, double* %0, i64 1
%7 = load double, double* %6
%8 = icmp sgt i32 %3, 0
br i1 %8, label %9, label %108
9: ; preds = %5
%10 = getelementptr inbounds double, double* %0, i64 4
%11 = load double, double* %10
%12 = fadd double %11, -1.000000e+00
%13 = getelementptr inbounds double, double* %0, i64 3
%14 = load double, double* %13
%15 = fadd double %14, -1.000000e+00
%16 = getelementptr inbounds double, double* %0, i64 2
%17 = load double, double* %16
%18 = fadd double %7, %17
%19 = fadd double %18, -2.000000e+00
%20 = load double, double* %0
%21 = fsub double %20, %7
%22 = fmul double %19, %19
%23 = tail call double @llvm.fmuladd.f64(double %21, double %21, double %22)
%24 = tail call double @llvm.fmuladd.f64(double %15, double %15, double %23)
%25 = tail call double @llvm.fmuladd.f64(double %12, double %12, double %24)
%26 = zext i32 %3 to i64
%27 = icmp ult i32 %3, 4
br i1 %27, label %101, label %28
28: ; preds = %9
%29 = and i64 %26, 4294967292
%30 = insertelement <2 x double> poison, double %25, i64 0
%31 = shufflevector <2 x double> %30, <2 x double> poison, <2 x i32> zeroinitializer
%32 = insertelement <2 x double> poison, double %25, i64 0
%33 = shufflevector <2 x double> %32, <2 x double> poison, <2 x i32> zeroinitializer
%34 = add nsw i64 %29, -4
%35 = lshr exact i64 %34, 2
%36 = add nuw nsw i64 %35, 1
%37 = and i64 %36, 7
%38 = icmp ult i64 %34, 28
br i1 %38, label %86, label %39
39: ; preds = %28
%40 = and i64 %36, 9223372036854775800
br label %41
41: ; preds = %41, %39
%42 = phi i64 [ 0, %39 ], [ %83, %41 ]
%43 = phi i64 [ 0, %39 ], [ %84, %41 ]
%44 = getelementptr inbounds double, double* %1, i64 %42
%45 = bitcast double* %44 to <2 x double>*
store <2 x double> %31, <2 x double>* %45
%46 = getelementptr inbounds double, double* %44, i64 2
%47 = bitcast double* %46 to <2 x double>*
store <2 x double> %33, <2 x double>* %47
%48 = or i64 %42, 4
%49 = getelementptr inbounds double, double* %1, i64 %48
%50 = bitcast double* %49 to <2 x double>*
store <2 x double> %31, <2 x double>* %50
%51 = getelementptr inbounds double, double* %49, i64 2
%52 = bitcast double* %51 to <2 x double>*
store <2 x double> %33, <2 x double>* %52
%53 = or i64 %42, 8
%54 = getelementptr inbounds double, double* %1, i64 %53
%55 = bitcast double* %54 to <2 x double>*
store <2 x double> %31, <2 x double>* %55
%56 = getelementptr inbounds double, double* %54, i64 2
%57 = bitcast double* %56 to <2 x double>*
store <2 x double> %33, <2 x double>* %57
%58 = or i64 %42, 12
%59 = getelementptr inbounds double, double* %1, i64 %58
%60 = bitcast double* %59 to <2 x double>*
store <2 x double> %31, <2 x double>* %60
%61 = getelementptr inbounds double, double* %59, i64 2
%62 = bitcast double* %61 to <2 x double>*
store <2 x double> %33, <2 x double>* %62
%63 = or i64 %42, 16
%64 = getelementptr inbounds double, double* %1, i64 %63
%65 = bitcast double* %64 to <2 x double>*
store <2 x double> %31, <2 x double>* %65
%66 = getelementptr inbounds double, double* %64, i64 2
%67 = bitcast double* %66 to <2 x double>*
store <2 x double> %33, <2 x double>* %67
%68 = or i64 %42, 20
%69 = getelementptr inbounds double, double* %1, i64 %68
%70 = bitcast double* %69 to <2 x double>*
store <2 x double> %31, <2 x double>* %70
%71 = getelementptr inbounds double, double* %69, i64 2
%72 = bitcast double* %71 to <2 x double>*
store <2 x double> %33, <2 x double>* %72
%73 = or i64 %42, 24
%74 = getelementptr inbounds double, double* %1, i64 %73
%75 = bitcast double* %74 to <2 x double>*
store <2 x double> %31, <2 x double>* %75
%76 = getelementptr inbounds double, double* %74, i64 2
%77 = bitcast double* %76 to <2 x double>*
store <2 x double> %33, <2 x double>* %77
%78 = or i64 %42, 28
%79 = getelementptr inbounds double, double* %1, i64 %78
%80 = bitcast double* %79 to <2 x double>*
store <2 x double> %31, <2 x double>* %80
%81 = getelementptr inbounds double, double* %79, i64 2
%82 = bitcast double* %81 to <2 x double>*
store <2 x double> %33, <2 x double>* %82
%83 = add nuw i64 %42, 32
%84 = add i64 %43, 8
%85 = icmp eq i64 %84, %40
br i1 %85, label %86, label %41
86: ; preds = %41, %28
%87 = phi i64 [ 0, %28 ], [ %83, %41 ]
%88 = icmp eq i64 %37, 0
br i1 %88, label %99, label %89
89: ; preds = %86, %89
%90 = phi i64 [ %96, %89 ], [ %87, %86 ]
%91 = phi i64 [ %97, %89 ], [ 0, %86 ]
%92 = getelementptr inbounds double, double* %1, i64 %90
%93 = bitcast double* %92 to <2 x double>*
store <2 x double> %31, <2 x double>* %93
%94 = getelementptr inbounds double, double* %92, i64 2
%95 = bitcast double* %94 to <2 x double>*
store <2 x double> %33, <2 x double>* %95
%96 = add nuw i64 %90, 4
%97 = add i64 %91, 1
%98 = icmp eq i64 %97, %37
br i1 %98, label %99, label %89
99: ; preds = %89, %86
%100 = icmp eq i64 %29, %26
br i1 %100, label %108, label %101
101: ; preds = %9, %99
%102 = phi i64 [ 0, %9 ], [ %29, %99 ]
br label %103
103: ; preds = %101, %103
%104 = phi i64 [ %106, %103 ], [ %102, %101 ]
%105 = getelementptr inbounds double, double* %1, i64 %104
store double %25, double* %105
%106 = add nuw nsw i64 %104, 1
%107 = icmp eq i64 %106, %26
br i1 %107, label %108, label %103
108: ; preds = %103, %99, %5
ret void
}