void _dmsum(float *def, int nxdef, int nydef, int nzdef, float *coefs,
            float *dmshape) {
  int i, k;
  int n = (nxdef * nydef);
  float co;
  for (i = 0; i < n; i++) {
    dmshape[i] = 0.0f;
  }
  for (k = 0; k < nzdef; k++) {
    co = coefs[k];
    for (i = 0; i < n; i++) {
      dmshape[i] += co * def[i + k * n];
    }
  }
}


x86 assembly (clang, -O0):
<_dmsum>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    -0xc(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1154 <_dmsum+0x54>
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    112a <_dmsum+0x2a>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11e4 <_dmsum+0xe4>
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x38(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    11d1 <_dmsum+0xd1>
movss  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1180 <_dmsum+0x80>
jmp    11d6 <_dmsum+0xd6>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    115b <_dmsum+0x5b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<_dmsum>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r9,%rbx
mov    %ecx,%r12d
mov    %esi,%r13d
imul   %edx,%r13d
mov    %r13d,%ebp
test   %r13d,%r13d
mov    %r8,0x8(%rsp)
mov    %rdi,(%rsp)
jle    1157 <_dmsum+0x47>
lea    0x0(,%rbp,4),%rdx
mov    %rbx,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
mov    (%rsp),%rdi
mov    0x8(%rsp),%r8
test   %r12d,%r12d
jle    1386 <_dmsum+0x276>
movslq %r13d,%rcx
mov    %r12d,%r11d
lea    (%rbx,%rbp,4),%rax
mov    %rax,0x40(%rsp)
mov    %ebp,%eax
and    $0xfffffff8,%eax
mov    %rax,0x30(%rsp)
add    $0xfffffffffffffff8,%rax
mov    %rax,0x28(%rsp)
mov    %rax,%rdx
shr    $0x3,%rdx
add    $0x1,%rdx
mov    %rdx,%rax
mov    %rdx,0x20(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,0x10(%rsp)
mov    %rbp,%r12
neg    %r12
lea    0x30(%rdi),%r15
lea    0x0(,%rcx,4),%rsi
lea    0x4(%rdi),%r9
xor    %r10d,%r10d
mov    %r11,0x38(%rsp)
mov    %r12,0x18(%rsp)
jmp    11e3 <_dmsum+0xd3>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r10
add    %rsi,%r15
add    %rsi,%r9
cmp    %r11,%r10
je     1386 <_dmsum+0x276>
test   %r13d,%r13d
jle    11d0 <_dmsum+0xc0>
mov    %r10,%r14
imul   %rcx,%r14
movss  (%r8,%r10,4),%xmm0
cmp    $0x8,%r13d
jb     1217 <_dmsum+0x107>
lea    (%r14,%rbp,1),%rax
lea    (%rdi,%rax,4),%rax
cmp    %rbx,%rax
jbe    1289 <_dmsum+0x179>
lea    (%rdi,%r14,4),%rax
cmp    0x40(%rsp),%rax
jae    1289 <_dmsum+0x179>
xor    %edx,%edx
mov    %rdx,%rax
test   $0x1,%bpl
je     1240 <_dmsum+0x130>
add    %rdx,%r14
movss  (%rdi,%r14,4),%xmm1
mulss  %xmm0,%xmm1
addss  (%rbx,%rdx,4),%xmm1
movss  %xmm1,(%rbx,%rdx,4)
mov    %rdx,%rax
or     $0x1,%rax
not    %rdx
cmp    %r12,%rdx
je     11d0 <_dmsum+0xc0>
nopl   0x0(%rax,%rax,1)
movss  -0x4(%r9,%rax,4),%xmm1
mulss  %xmm0,%xmm1
addss  (%rbx,%rax,4),%xmm1
movss  %xmm1,(%rbx,%rax,4)
movss  (%r9,%rax,4),%xmm1
mulss  %xmm0,%xmm1
addss  0x4(%rbx,%rax,4),%xmm1
movss  %xmm1,0x4(%rbx,%rax,4)
add    $0x2,%rax
cmp    %rax,%rbp
jne    1250 <_dmsum+0x140>
jmp    11d0 <_dmsum+0xc0>
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
cmpq   $0x0,0x28(%rsp)
je     137a <_dmsum+0x26a>
mov    0x10(%rsp),%r12
xor    %r11d,%r11d
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movups -0x30(%r15,%r11,4),%xmm2
movups -0x20(%r15,%r11,4),%xmm3
mulps  %xmm1,%xmm2
movups (%rbx,%r11,4),%xmm4
addps  %xmm2,%xmm4
movups 0x10(%rbx,%r11,4),%xmm2
movups 0x20(%rbx,%r11,4),%xmm5
movups 0x30(%rbx,%r11,4),%xmm6
mulps  %xmm1,%xmm3
addps  %xmm2,%xmm3
movups %xmm4,(%rbx,%r11,4)
movups %xmm3,0x10(%rbx,%r11,4)
movups -0x10(%r15,%r11,4),%xmm2
movups (%r15,%r11,4),%xmm3
mulps  %xmm1,%xmm2
addps  %xmm5,%xmm2
mulps  %xmm1,%xmm3
addps  %xmm6,%xmm3
movups %xmm2,0x20(%rbx,%r11,4)
movups %xmm3,0x30(%rbx,%r11,4)
add    $0x10,%r11
add    $0xfffffffffffffffe,%r12
jne    12b0 <_dmsum+0x1a0>
testb  $0x1,0x20(%rsp)
je     1351 <_dmsum+0x241>
lea    (%r11,%r14,1),%rax
mov    (%rsp),%rdx
movups (%rdx,%rax,4),%xmm2
movups 0x10(%rdx,%rax,4),%xmm3
movups (%rbx,%r11,4),%xmm4
movups 0x10(%rbx,%r11,4),%xmm5
mulps  %xmm1,%xmm2
addps  %xmm4,%xmm2
mulps  %xmm1,%xmm3
addps  %xmm5,%xmm3
movups %xmm2,(%rbx,%r11,4)
movups %xmm3,0x10(%rbx,%r11,4)
mov    0x30(%rsp),%rax
mov    %rax,%rdx
cmp    %rbp,%rax
mov    0x8(%rsp),%r8
mov    (%rsp),%rdi
mov    0x38(%rsp),%r11
mov    0x18(%rsp),%r12
je     11d0 <_dmsum+0xc0>
jmp    1219 <_dmsum+0x109>
xor    %r11d,%r11d
testb  $0x1,0x20(%rsp)
jne    131e <_dmsum+0x20e>
jmp    1351 <_dmsum+0x241>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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