void h1(float *c, float *a, double b[], int size) {
  int t = 0;
#pragma omp simd safelen(16) linear(t) aligned(c : 32) aligned(a, b)
  for (int i = 0; i < size; ++i) {
    c[i] = a[i] * a[i] + b[i] * b[t];
    ++t;
  }
}


x86 assembly (clang, -O0):
<h1>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1195 <h1+0x95>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1121 <h1+0x21>
pop    %rbp
ret

x86 assembly (clang, -O2):
<h1>:
test   %ecx,%ecx
jle    1235 <h1+0x135>
mov    %ecx,%r8d
cmp    $0x4,%ecx
jb     112a <h1+0x2a>
lea    (%rsi,%r8,4),%rax
cmp    %rdi,%rax
jbe    11ca <h1+0xca>
lea    (%rdi,%r8,4),%rax
cmp    %rsi,%rax
jbe    11ca <h1+0xca>
xor    %ecx,%ecx
mov    %rcx,%rax
not    %rax
test   $0x1,%r8b
je     1162 <h1+0x62>
movss  (%rsi,%rcx,4),%xmm0
mulss  %xmm0,%xmm0
cvtss2sd %xmm0,%xmm0
movsd  (%rdx,%rcx,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm0,(%rdi,%rcx,4)
or     $0x1,%rcx
add    %r8,%rax
je     1235 <h1+0x135>
nopl   0x0(%rax,%rax,1)
movss  (%rsi,%rcx,4),%xmm0
mulss  %xmm0,%xmm0
cvtss2sd %xmm0,%xmm0
movsd  (%rdx,%rcx,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm0,(%rdi,%rcx,4)
movss  0x4(%rsi,%rcx,4),%xmm0
mulss  %xmm0,%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0x8(%rdx,%rcx,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm0,0x4(%rdi,%rcx,4)
add    $0x2,%rcx
cmp    %rcx,%r8
jne    1170 <h1+0x70>
jmp    1235 <h1+0x135>
mov    %r8d,%ecx
and    $0xfffffffc,%ecx
lea    0x0(,%r8,4),%r9
and    $0xfffffffffffffff0,%r9
xor    %eax,%eax
xchg   %ax,%ax
movsd  (%rsi,%rax,1),%xmm0
movsd  0x8(%rsi,%rax,1),%xmm1
mulps  %xmm0,%xmm0
mulps  %xmm1,%xmm1
cvtps2pd %xmm0,%xmm0
cvtps2pd %xmm1,%xmm1
movupd (%rdx,%rax,2),%xmm2
movupd 0x10(%rdx,%rax,2),%xmm3
mulpd  %xmm2,%xmm2
addpd  %xmm0,%xmm2
mulpd  %xmm3,%xmm3
addpd  %xmm1,%xmm3
cvtpd2ps %xmm2,%xmm0
cvtpd2ps %xmm3,%xmm1
unpcklpd %xmm1,%xmm0
movupd %xmm0,(%rdi,%rax,1)
add    $0x10,%rax
cmp    %rax,%r9
jne    11e0 <h1+0xe0>
cmp    %r8,%rcx
jne    112c <h1+0x2c>
ret

LLVM IR (-O0):
define void @h1(float* %0, float* %1, double* %2, i32 %3) {
%5 = alloca float*
%6 = alloca float*
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store float* %0, float** %5
store float* %1, float** %6
store double* %2, double** %7
store i32 %3, i32* %8
store i32 0, i32* %9
store i32 0, i32* %10
br label %11
11: ; preds = %46, %4
%12 = load i32, i32* %10
%13 = load i32, i32* %8
%14 = icmp slt i32 %12, %13
br i1 %14, label %15, label %49
15: ; preds = %11
%16 = load float*, float** %6
%17 = load i32, i32* %10
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds float, float* %16, i64 %18
%20 = load float, float* %19
%21 = load float*, float** %6
%22 = load i32, i32* %10
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds float, float* %21, i64 %23
%25 = load float, float* %24
%26 = fmul float %20, %25
%27 = fpext float %26 to double
%28 = load double*, double** %7
%29 = load i32, i32* %10
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %28, i64 %30
%32 = load double, double* %31
%33 = load double*, double** %7
%34 = load i32, i32* %9
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double, double* %33, i64 %35
%37 = load double, double* %36
%38 = call double @llvm.fmuladd.f64(double %32, double %37, double %27)
%39 = fptrunc double %38 to float
%40 = load float*, float** %5
%41 = load i32, i32* %10
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds float, float* %40, i64 %42
store float %39, float* %43
%44 = load i32, i32* %9
%45 = add nsw i32 %44, 1
store i32 %45, i32* %9
br label %46
46: ; preds = %15
%47 = load i32, i32* %10
%48 = add nsw i32 %47, 1
store i32 %48, i32* %10
br label %11
49: ; preds = %11
ret void
}

LLVM IR (-O2):
define void @h1(float* nocapture writeonly %0, float* nocapture %1, double* nocapture %2, i32 %3) {
%5 = icmp sgt i32 %3, 0
br i1 %5, label %6, label %67
6: ; preds = %4
%7 = zext i32 %3 to i64
%8 = icmp ult i32 %3, 4
br i1 %8, label %47, label %9
9: ; preds = %6
%10 = getelementptr float, float* %0, i64 %7
%11 = getelementptr float, float* %1, i64 %7
%12 = icmp ugt float* %11, %0
%13 = icmp ugt float* %10, %1
%14 = and i1 %12, %13
br i1 %14, label %47, label %15
15: ; preds = %9
%16 = and i64 %7, 4294967292
br label %17
17: ; preds = %17, %15
%18 = phi i64 [ 0, %15 ], [ %43, %17 ]
%19 = getelementptr inbounds float, float* %1, i64 %18
%20 = bitcast float* %19 to <2 x float>*
%21 = load <2 x float>, <2 x float>* %20
%22 = getelementptr inbounds float, float* %19, i64 2
%23 = bitcast float* %22 to <2 x float>*
%24 = load <2 x float>, <2 x float>* %23
%25 = fmul <2 x float> %21, %21
%26 = fmul <2 x float> %24, %24
%27 = fpext <2 x float> %25 to <2 x double>
%28 = fpext <2 x float> %26 to <2 x double>
%29 = getelementptr inbounds double, double* %2, i64 %18
%30 = bitcast double* %29 to <2 x double>*
%31 = load <2 x double>, <2 x double>* %30
%32 = getelementptr inbounds double, double* %29, i64 2
%33 = bitcast double* %32 to <2 x double>*
%34 = load <2 x double>, <2 x double>* %33
%35 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %31, <2 x double> %31, <2 x double> %27)
%36 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %34, <2 x double> %34, <2 x double> %28)
%37 = fptrunc <2 x double> %35 to <2 x float>
%38 = fptrunc <2 x double> %36 to <2 x float>
%39 = getelementptr inbounds float, float* %0, i64 %18
%40 = bitcast float* %39 to <2 x float>*
store <2 x float> %37, <2 x float>* %40
%41 = getelementptr inbounds float, float* %39, i64 2
%42 = bitcast float* %41 to <2 x float>*
store <2 x float> %38, <2 x float>* %42
%43 = add nuw i64 %18, 4
%44 = icmp eq i64 %43, %16
br i1 %44, label %45, label %17
45: ; preds = %17
%46 = icmp eq i64 %16, %7
br i1 %46, label %67, label %47
47: ; preds = %9, %6, %45
%48 = phi i64 [ 0, %9 ], [ 0, %6 ], [ %16, %45 ]
%49 = xor i64 %48, -1
%50 = and i64 %7, 1
%51 = icmp eq i64 %50, 0
br i1 %51, label %63, label %52
52: ; preds = %47
%53 = getelementptr inbounds float, float* %1, i64 %48
%54 = load float, float* %53
%55 = fmul float %54, %54
%56 = fpext float %55 to double
%57 = getelementptr inbounds double, double* %2, i64 %48
%58 = load double, double* %57
%59 = tail call double @llvm.fmuladd.f64(double %58, double %58, double %56)
%60 = fptrunc double %59 to float
%61 = getelementptr inbounds float, float* %0, i64 %48
store float %60, float* %61
%62 = or i64 %48, 1
br label %63
63: ; preds = %52, %47
%64 = phi i64 [ %48, %47 ], [ %62, %52 ]
%65 = sub nsw i64 0, %7
%66 = icmp eq i64 %49, %65
br i1 %66, label %67, label %68
67: ; preds = %63, %68, %45, %4
ret void
68: ; preds = %63, %68
%69 = phi i64 [ %89, %68 ], [ %64, %63 ]
%70 = getelementptr inbounds float, float* %1, i64 %69
%71 = load float, float* %70
%72 = fmul float %71, %71
%73 = fpext float %72 to double
%74 = getelementptr inbounds double, double* %2, i64 %69
%75 = load double, double* %74
%76 = tail call double @llvm.fmuladd.f64(double %75, double %75, double %73)
%77 = fptrunc double %76 to float
%78 = getelementptr inbounds float, float* %0, i64 %69
store float %77, float* %78
%79 = add nuw nsw i64 %69, 1
%80 = getelementptr inbounds float, float* %1, i64 %79
%81 = load float, float* %80
%82 = fmul float %81, %81
%83 = fpext float %82 to double
%84 = getelementptr inbounds double, double* %2, i64 %79
%85 = load double, double* %84
%86 = tail call double @llvm.fmuladd.f64(double %85, double %85, double %83)
%87 = fptrunc double %86 to float
%88 = getelementptr inbounds float, float* %0, i64 %79
store float %87, float* %88
%89 = add nuw nsw i64 %69, 2
%90 = icmp eq i64 %89, %7
br i1 %90, label %67, label %68
}