mvmpy(r, c, A, V, Y) int r, c;
double *A, *V, *Y;
{
  register double s;
  double *pA, *pV, *pY;
  int i, j;
  pA = A;
  pY = Y;
  for (i = 0; i < r; i++) {
    pV = V;
    s = 0.0;
    for (j = 0; j < c; j++) {
      s += *pA++ * *pV++;
    }
    *pY++ = s;
  }
}


x86 assembly (clang, -O0):
<mvmpy>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x28(%rbp),%rax
mov    %rax,-0x48(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11c8 <mvmpy+0xc8>
mov    -0x20(%rbp),%rax
mov    %rax,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11a2 <mvmpy+0xa2>
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1150 <mvmpy+0x50>
movsd  -0x30(%rbp),%xmm0
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
movsd  %xmm0,(%rax)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    112d <mvmpy+0x2d>
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<mvmpy>:
push   %r14
push   %rbx
test   %edi,%edi
jle    11ef <mvmpy+0xef>
lea    -0x1(%rsi),%r9d
mov    %esi,%r10d
and    $0xfffffffc,%r10d
mov    %esi,%r14d
and    $0x3,%r14d
shl    $0x3,%r14
xor    %r11d,%r11d
jmp    1146 <mvmpy+0x46>
cs nopw 0x0(%rax,%rax,1)
movsd  %xmm0,(%r8)
add    $0x8,%r8
add    $0x1,%r11d
cmp    %edi,%r11d
je     11ef <mvmpy+0xef>
xorpd  %xmm0,%xmm0
test   %esi,%esi
jle    1130 <mvmpy+0x30>
cmp    $0x3,%r9d
jae    1160 <mvmpy+0x60>
mov    %rcx,%rbx
test   $0x3,%sil
jne    11ce <mvmpy+0xce>
jmp    1130 <mvmpy+0x30>
nop
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movsd  (%rdx,%rbx,8),%xmm1
mulsd  (%rcx,%rbx,8),%xmm1
movsd  0x8(%rdx,%rbx,8),%xmm2
addsd  %xmm0,%xmm1
mulsd  0x8(%rcx,%rbx,8),%xmm2
addsd  %xmm1,%xmm2
movsd  0x10(%rdx,%rbx,8),%xmm1
mulsd  0x10(%rcx,%rbx,8),%xmm1
addsd  %xmm2,%xmm1
movsd  0x18(%rdx,%rbx,8),%xmm0
mulsd  0x18(%rcx,%rbx,8),%xmm0
addsd  %xmm1,%xmm0
add    $0xffffffffffffffe0,%rax
add    $0x4,%rbx
cmp    %ebx,%r10d
jne    1170 <mvmpy+0x70>
mov    %rcx,%rbx
sub    %rax,%rbx
sub    %rax,%rdx
test   $0x3,%sil
je     1130 <mvmpy+0x30>
xor    %eax,%eax
movsd  (%rdx,%rax,1),%xmm1
mulsd  (%rbx,%rax,1),%xmm1
addsd  %xmm1,%xmm0
add    $0x8,%rax
cmp    %eax,%r14d
jne    11d0 <mvmpy+0xd0>
add    %rax,%rdx
jmp    1130 <mvmpy+0x30>
pop    %rbx
pop    %r14
ret

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

LLVM IR (-O2):
define i32 @mvmpy(i32 %0, i32 %1, double* nocapture %2, double* nocapture %3, double* nocapture writeonly %4) {
%6 = icmp sgt i32 %0, 0
br i1 %6, label %7, label %70
7: ; preds = %5
%8 = icmp sgt i32 %1, 0
%9 = add i32 %1, -1
%10 = and i32 %1, 3
%11 = icmp ult i32 %9, 3
%12 = and i32 %1, -4
%13 = icmp eq i32 %10, 0
br label %14
14: ; preds = %7, %64
%15 = phi i32 [ 0, %7 ], [ %68, %64 ]
%16 = phi double* [ %4, %7 ], [ %67, %64 ]
%17 = phi double* [ %2, %7 ], [ %66, %64 ]
br i1 %8, label %18, label %64
18: ; preds = %14
br i1 %11, label %46, label %19
19: ; preds = %18, %19
%20 = phi double* [ %41, %19 ], [ %3, %18 ]
%21 = phi double* [ %39, %19 ], [ %17, %18 ]
%22 = phi double [ %43, %19 ], [ 0.000000e+00, %18 ]
%23 = phi i32 [ %44, %19 ], [ 0, %18 ]
%24 = getelementptr inbounds double, double* %21, i64 1
%25 = load double, double* %21
%26 = getelementptr inbounds double, double* %20, i64 1
%27 = load double, double* %20
%28 = tail call double @llvm.fmuladd.f64(double %25, double %27, double %22)
%29 = getelementptr inbounds double, double* %21, i64 2
%30 = load double, double* %24
%31 = getelementptr inbounds double, double* %20, i64 2
%32 = load double, double* %26
%33 = tail call double @llvm.fmuladd.f64(double %30, double %32, double %28)
%34 = getelementptr inbounds double, double* %21, i64 3
%35 = load double, double* %29
%36 = getelementptr inbounds double, double* %20, i64 3
%37 = load double, double* %31
%38 = tail call double @llvm.fmuladd.f64(double %35, double %37, double %33)
%39 = getelementptr inbounds double, double* %21, i64 4
%40 = load double, double* %34
%41 = getelementptr inbounds double, double* %20, i64 4
%42 = load double, double* %36
%43 = tail call double @llvm.fmuladd.f64(double %40, double %42, double %38)
%44 = add i32 %23, 4
%45 = icmp eq i32 %44, %12
br i1 %45, label %46, label %19
46: ; preds = %19, %18
%47 = phi double* [ undef, %18 ], [ %39, %19 ]
%48 = phi double [ undef, %18 ], [ %43, %19 ]
%49 = phi double* [ %3, %18 ], [ %41, %19 ]
%50 = phi double* [ %17, %18 ], [ %39, %19 ]
%51 = phi double [ 0.000000e+00, %18 ], [ %43, %19 ]
br i1 %13, label %64, label %52
52: ; preds = %46, %52
%53 = phi double* [ %59, %52 ], [ %49, %46 ]
%54 = phi double* [ %57, %52 ], [ %50, %46 ]
%55 = phi double [ %61, %52 ], [ %51, %46 ]
%56 = phi i32 [ %62, %52 ], [ 0, %46 ]
%57 = getelementptr inbounds double, double* %54, i64 1
%58 = load double, double* %54
%59 = getelementptr inbounds double, double* %53, i64 1
%60 = load double, double* %53
%61 = tail call double @llvm.fmuladd.f64(double %58, double %60, double %55)
%62 = add i32 %56, 1
%63 = icmp eq i32 %62, %10
br i1 %63, label %64, label %52
64: ; preds = %46, %52, %14
%65 = phi double [ 0.000000e+00, %14 ], [ %48, %46 ], [ %61, %52 ]
%66 = phi double* [ %17, %14 ], [ %47, %46 ], [ %57, %52 ]
%67 = getelementptr inbounds double, double* %16, i64 1
store double %65, double* %16
%68 = add nuw nsw i32 %15, 1
%69 = icmp eq i32 %68, %0
br i1 %69, label %70, label %14
70: ; preds = %64, %5
ret i32 undef
}