void mult3drm(float *array1, float *array2, int Ntot) {
  int i;
  float tmp, *pt1, *pt2;
  pt1 = array1;
  pt2 = array2;
  for (i = 0; i < Ntot; i++) {
    tmp = *pt1 * *pt2;
    *pt1 = tmp;
    pt1++;
    pt2++;
  }
}


x86 assembly (clang, -O0):
<mult3drm>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x30(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    117a <mult3drm+0x7a>
mov    -0x28(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x30(%rbp),%rax
mulss  (%rax),%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x1c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x28(%rbp)
mov    -0x30(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x30(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1126 <mult3drm+0x26>
pop    %rbp
ret

x86 assembly (clang, -O2):
<mult3drm>:
test   %edx,%edx
jle    128d <mult3drm+0x18d>
lea    -0x1(%rdx),%r10d
xor    %r8d,%r8d
cmp    $0x7,%r10d
jb     11f7 <mult3drm+0xf7>
lea    (%rsi,%r10,4),%rax
add    $0x4,%rax
cmp    %rdi,%rax
jbe    1137 <mult3drm+0x37>
lea    (%rdi,%r10,4),%rax
add    $0x4,%rax
cmp    %rsi,%rax
ja     11f7 <mult3drm+0xf7>
add    $0x1,%r10
mov    %r10,%r8
and    $0xfffffffffffffff8,%r8
lea    -0x8(%r8),%rax
mov    %rax,%r9
shr    $0x3,%r9
add    $0x1,%r9
test   %rax,%rax
je     128e <mult3drm+0x18e>
mov    %r9,%rcx
and    $0xfffffffffffffffe,%rcx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movups (%rsi,%rax,4),%xmm0
movups 0x10(%rsi,%rax,4),%xmm1
movups (%rdi,%rax,4),%xmm2
mulps  %xmm0,%xmm2
movups 0x10(%rdi,%rax,4),%xmm0
mulps  %xmm1,%xmm0
movups 0x20(%rdi,%rax,4),%xmm1
movups 0x30(%rdi,%rax,4),%xmm3
movups %xmm2,(%rdi,%rax,4)
movups %xmm0,0x10(%rdi,%rax,4)
movups 0x20(%rsi,%rax,4),%xmm0
mulps  %xmm1,%xmm0
movups 0x30(%rsi,%rax,4),%xmm1
mulps  %xmm3,%xmm1
movups %xmm0,0x20(%rdi,%rax,4)
movups %xmm1,0x30(%rdi,%rax,4)
add    $0x10,%rax
add    $0xfffffffffffffffe,%rcx
jne    1170 <mult3drm+0x70>
test   $0x1,%r9b
je     11e6 <mult3drm+0xe6>
movups (%rdi,%rax,4),%xmm0
movups 0x10(%rdi,%rax,4),%xmm1
movups (%rsi,%rax,4),%xmm2
mulps  %xmm0,%xmm2
movups 0x10(%rsi,%rax,4),%xmm0
mulps  %xmm1,%xmm0
movups %xmm2,(%rdi,%rax,4)
movups %xmm0,0x10(%rdi,%rax,4)
cmp    %r8,%r10
je     128d <mult3drm+0x18d>
lea    (%rsi,%r8,4),%rsi
lea    (%rdi,%r8,4),%rdi
mov    %edx,%ecx
sub    %r8d,%ecx
mov    %r8d,%r9d
not    %r9d
add    %edx,%r9d
and    $0x3,%ecx
je     122e <mult3drm+0x12e>
neg    %ecx
xor    %eax,%eax
xchg   %ax,%ax
movss  (%rdi),%xmm0
mulss  (%rsi),%xmm0
movss  %xmm0,(%rdi)
add    $0x4,%rdi
add    $0x4,%rsi
add    $0xffffffff,%eax
cmp    %eax,%ecx
jne    1210 <mult3drm+0x110>
sub    %eax,%r8d
cmp    $0x3,%r9d
jb     128d <mult3drm+0x18d>
sub    %r8d,%edx
xor    %eax,%eax
nopl   0x0(%rax)
movss  (%rdi,%rax,4),%xmm0
mulss  (%rsi,%rax,4),%xmm0
movss  0x4(%rdi,%rax,4),%xmm1
movss  %xmm0,(%rdi,%rax,4)
mulss  0x4(%rsi,%rax,4),%xmm1
movss  %xmm1,0x4(%rdi,%rax,4)
movss  0x8(%rdi,%rax,4),%xmm0
mulss  0x8(%rsi,%rax,4),%xmm0
movss  %xmm0,0x8(%rdi,%rax,4)
movss  0xc(%rdi,%rax,4),%xmm0
mulss  0xc(%rsi,%rax,4),%xmm0
movss  %xmm0,0xc(%rdi,%rax,4)
add    $0x4,%rax
cmp    %eax,%edx
jne    1240 <mult3drm+0x140>
ret
xor    %eax,%eax
test   $0x1,%r9b
jne    11c5 <mult3drm+0xc5>
jmp    11e6 <mult3drm+0xe6>

LLVM IR (-O0):
define void @mult3drm(float* %0, float* %1, i32 %2) {
%4 = alloca float*
%5 = alloca float*
%6 = alloca i32
%7 = alloca i32
%8 = alloca float
%9 = alloca float*
%10 = alloca float*
store float* %0, float** %4
store float* %1, float** %5
store i32 %2, i32* %6
%11 = load float*, float** %4
store float* %11, float** %9
%12 = load float*, float** %5
store float* %12, float** %10
store i32 0, i32* %7
br label %13
13: ; preds = %29, %3
%14 = load i32, i32* %7
%15 = load i32, i32* %6
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %32
17: ; preds = %13
%18 = load float*, float** %9
%19 = load float, float* %18
%20 = load float*, float** %10
%21 = load float, float* %20
%22 = fmul float %19, %21
store float %22, float* %8
%23 = load float, float* %8
%24 = load float*, float** %9
store float %23, float* %24
%25 = load float*, float** %9
%26 = getelementptr inbounds float, float* %25, i32 1
store float* %26, float** %9
%27 = load float*, float** %10
%28 = getelementptr inbounds float, float* %27, i32 1
store float* %28, float** %10
br label %29
29: ; preds = %17
%30 = load i32, i32* %7
%31 = add nsw i32 %30, 1
store i32 %31, i32* %7
br label %13
32: ; preds = %13
ret void
}

LLVM IR (-O2):
define void @mult3drm(float* nocapture %0, float* nocapture %1, i32 %2) {
%4 = icmp sgt i32 %2, 0
br i1 %4, label %5, label %145
5: ; preds = %3
%6 = add i32 %2, -1
%7 = zext i32 %6 to i64
%8 = add nuw nsw i64 %7, 1
%9 = icmp ult i32 %6, 7
br i1 %9, label %92, label %10
10: ; preds = %5
%11 = add i32 %2, -1
%12 = zext i32 %11 to i64
%13 = add nuw nsw i64 %12, 1
%14 = getelementptr float, float* %0, i64 %13
%15 = getelementptr float, float* %1, i64 %13
%16 = icmp ugt float* %15, %0
%17 = icmp ugt float* %14, %1
%18 = and i1 %16, %17
br i1 %18, label %92, label %19
19: ; preds = %10
%20 = and i64 %8, 8589934584
%21 = getelementptr float, float* %1, i64 %20
%22 = getelementptr float, float* %0, i64 %20
%23 = trunc i64 %20 to i32
%24 = add nsw i64 %20, -8
%25 = lshr exact i64 %24, 3
%26 = add nuw nsw i64 %25, 1
%27 = and i64 %26, 1
%28 = icmp eq i64 %24, 0
br i1 %28, label %70, label %29
29: ; preds = %19
%30 = and i64 %26, 4611686018427387902
br label %31
31: ; preds = %31, %29
%32 = phi i64 [ 0, %29 ], [ %67, %31 ]
%33 = phi i64 [ 0, %29 ], [ %68, %31 ]
%34 = getelementptr float, float* %1, i64 %32
%35 = getelementptr float, float* %0, i64 %32
%36 = bitcast float* %35 to <4 x float>*
%37 = load <4 x float>, <4 x float>* %36
%38 = getelementptr float, float* %35, i64 4
%39 = bitcast float* %38 to <4 x float>*
%40 = load <4 x float>, <4 x float>* %39
%41 = bitcast float* %34 to <4 x float>*
%42 = load <4 x float>, <4 x float>* %41
%43 = getelementptr float, float* %34, i64 4
%44 = bitcast float* %43 to <4 x float>*
%45 = load <4 x float>, <4 x float>* %44
%46 = fmul <4 x float> %37, %42
%47 = fmul <4 x float> %40, %45
%48 = bitcast float* %35 to <4 x float>*
store <4 x float> %46, <4 x float>* %48
%49 = bitcast float* %38 to <4 x float>*
store <4 x float> %47, <4 x float>* %49
%50 = or i64 %32, 8
%51 = getelementptr float, float* %1, i64 %50
%52 = getelementptr float, float* %0, i64 %50
%53 = bitcast float* %52 to <4 x float>*
%54 = load <4 x float>, <4 x float>* %53
%55 = getelementptr float, float* %52, i64 4
%56 = bitcast float* %55 to <4 x float>*
%57 = load <4 x float>, <4 x float>* %56
%58 = bitcast float* %51 to <4 x float>*
%59 = load <4 x float>, <4 x float>* %58
%60 = getelementptr float, float* %51, i64 4
%61 = bitcast float* %60 to <4 x float>*
%62 = load <4 x float>, <4 x float>* %61
%63 = fmul <4 x float> %54, %59
%64 = fmul <4 x float> %57, %62
%65 = bitcast float* %52 to <4 x float>*
store <4 x float> %63, <4 x float>* %65
%66 = bitcast float* %55 to <4 x float>*
store <4 x float> %64, <4 x float>* %66
%67 = add nuw i64 %32, 16
%68 = add i64 %33, 2
%69 = icmp eq i64 %68, %30
br i1 %69, label %70, label %31
70: ; preds = %31, %19
%71 = phi i64 [ 0, %19 ], [ %67, %31 ]
%72 = icmp eq i64 %27, 0
br i1 %72, label %90, label %73
73: ; preds = %70
%74 = getelementptr float, float* %1, i64 %71
%75 = getelementptr float, float* %0, i64 %71
%76 = bitcast float* %75 to <4 x float>*
%77 = load <4 x float>, <4 x float>* %76
%78 = getelementptr float, float* %75, i64 4
%79 = bitcast float* %78 to <4 x float>*
%80 = load <4 x float>, <4 x float>* %79
%81 = bitcast float* %74 to <4 x float>*
%82 = load <4 x float>, <4 x float>* %81
%83 = getelementptr float, float* %74, i64 4
%84 = bitcast float* %83 to <4 x float>*
%85 = load <4 x float>, <4 x float>* %84
%86 = fmul <4 x float> %77, %82
%87 = fmul <4 x float> %80, %85
%88 = bitcast float* %75 to <4 x float>*
store <4 x float> %86, <4 x float>* %88
%89 = bitcast float* %78 to <4 x float>*
store <4 x float> %87, <4 x float>* %89
br label %90
90: ; preds = %70, %73
%91 = icmp eq i64 %8, %20
br i1 %91, label %145, label %92
92: ; preds = %10, %5, %90
%93 = phi float* [ %1, %10 ], [ %1, %5 ], [ %21, %90 ]
%94 = phi float* [ %0, %10 ], [ %0, %5 ], [ %22, %90 ]
%95 = phi i32 [ 0, %10 ], [ 0, %5 ], [ %23, %90 ]
%96 = sub i32 %2, %95
%97 = xor i32 %95, -1
%98 = add i32 %97, %2
%99 = and i32 %96, 3
%100 = icmp eq i32 %99, 0
br i1 %100, label %114, label %101
101: ; preds = %92, %101
%102 = phi float* [ %110, %101 ], [ %93, %92 ]
%103 = phi float* [ %109, %101 ], [ %94, %92 ]
%104 = phi i32 [ %111, %101 ], [ %95, %92 ]
%105 = phi i32 [ %112, %101 ], [ 0, %92 ]
%106 = load float, float* %103
%107 = load float, float* %102
%108 = fmul float %106, %107
store float %108, float* %103
%109 = getelementptr inbounds float, float* %103, i64 1
%110 = getelementptr inbounds float, float* %102, i64 1
%111 = add nuw nsw i32 %104, 1
%112 = add i32 %105, 1
%113 = icmp eq i32 %112, %99
br i1 %113, label %114, label %101
114: ; preds = %101, %92
%115 = phi float* [ %93, %92 ], [ %110, %101 ]
%116 = phi float* [ %94, %92 ], [ %109, %101 ]
%117 = phi i32 [ %95, %92 ], [ %111, %101 ]
%118 = icmp ult i32 %98, 3
br i1 %118, label %145, label %119
119: ; preds = %114, %119
%120 = phi float* [ %142, %119 ], [ %115, %114 ]
%121 = phi float* [ %141, %119 ], [ %116, %114 ]
%122 = phi i32 [ %143, %119 ], [ %117, %114 ]
%123 = load float, float* %121
%124 = load float, float* %120
%125 = fmul float %123, %124
store float %125, float* %121
%126 = getelementptr inbounds float, float* %121, i64 1
%127 = getelementptr inbounds float, float* %120, i64 1
%128 = load float, float* %126
%129 = load float, float* %127
%130 = fmul float %128, %129
store float %130, float* %126
%131 = getelementptr inbounds float, float* %121, i64 2
%132 = getelementptr inbounds float, float* %120, i64 2
%133 = load float, float* %131
%134 = load float, float* %132
%135 = fmul float %133, %134
store float %135, float* %131
%136 = getelementptr inbounds float, float* %121, i64 3
%137 = getelementptr inbounds float, float* %120, i64 3
%138 = load float, float* %136
%139 = load float, float* %137
%140 = fmul float %138, %139
store float %140, float* %136
%141 = getelementptr inbounds float, float* %121, i64 4
%142 = getelementptr inbounds float, float* %120, i64 4
%143 = add nuw nsw i32 %122, 4
%144 = icmp eq i32 %143, %2
br i1 %144, label %145, label %119
145: ; preds = %114, %119, %90, %3
ret void
}