double dd(int k, int n, int m, const double *x) {
  double d = 1;
  double q = x[k];
  for (int l = 0; l < m; l++) {
    for (int ll = l; ll < n; ll += m) {
      double xx = 2 * (q - x[ll]);
      if (xx != 0)
        d *= xx;
    }
  }
  return 1.0 / d;
}


x86 assembly (clang, -O0):
<dd>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %rcx,-0x18(%rbp)
movsd  0xee7(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11ca <dd+0xca>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11b7 <dd+0xb7>
movsd  -0x28(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
movsd  0xe99(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1195 <dd+0x95>
jp     1195 <dd+0x95>
jmp    11a4 <dd+0xa4>
movsd  -0x38(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    11a9 <dd+0xa9>
mov    -0xc(%rbp),%eax
add    -0x30(%rbp),%eax
mov    %eax,-0x30(%rbp)
jmp    1149 <dd+0x49>
jmp    11bc <dd+0xbc>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1137 <dd+0x37>
movsd  0xe2e(%rip),%xmm0
divsd  -0x20(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<dd>:
movsd  0xef8(%rip),%xmm0
movapd %xmm0,%xmm1
test   %edx,%edx
jle    1190 <dd+0x90>
movslq %edi,%rax
movsd  (%rcx,%rax,8),%xmm2
mov    %edx,%eax
movslq %esi,%rdx
movsd  0xed7(%rip),%xmm1
xor    %r8d,%r8d
xorpd  %xmm3,%xmm3
movsd  0xec8(%rip),%xmm4
jmp    1149 <dd+0x49>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r8
cmp    %rax,%r8
je     1190 <dd+0x90>
mov    %r8,%rdi
cmp    %esi,%r8d
jge    1140 <dd+0x40>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movapd %xmm2,%xmm5
subsd  (%rcx,%rdi,8),%xmm5
addsd  %xmm5,%xmm5
movapd %xmm5,%xmm6
cmpneqsd %xmm3,%xmm6
andpd  %xmm6,%xmm5
andnpd %xmm4,%xmm6
orpd   %xmm5,%xmm6
mulsd  %xmm6,%xmm1
add    %rax,%rdi
cmp    %rdx,%rdi
jl     1160 <dd+0x60>
jmp    1140 <dd+0x40>
divsd  %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @dd(i32 %0, i32 %1, i32 %2, double* %3) {
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca double*
%9 = alloca double
%10 = alloca double
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
store i32 %0, i32* %5
store i32 %1, i32* %6
store i32 %2, i32* %7
store double* %3, double** %8
store double 1.000000e+00, double* %9
%14 = load double*, double** %8
%15 = load i32, i32* %5
%16 = sext i32 %15 to i64
%17 = getelementptr inbounds double, double* %14, i64 %16
%18 = load double, double* %17
store double %18, double* %10
store i32 0, i32* %11
br label %19
19: ; preds = %50, %4
%20 = load i32, i32* %11
%21 = load i32, i32* %7
%22 = icmp slt i32 %20, %21
br i1 %22, label %23, label %53
23: ; preds = %19
%24 = load i32, i32* %11
store i32 %24, i32* %12
br label %25
25: ; preds = %45, %23
%26 = load i32, i32* %12
%27 = load i32, i32* %6
%28 = icmp slt i32 %26, %27
br i1 %28, label %29, label %49
29: ; preds = %25
%30 = load double, double* %10
%31 = load double*, double** %8
%32 = load i32, i32* %12
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds double, double* %31, i64 %33
%35 = load double, double* %34
%36 = fsub double %30, %35
%37 = fmul double 2.000000e+00, %36
store double %37, double* %13
%38 = load double, double* %13
%39 = fcmp une double %38, 0.000000e+00
br i1 %39, label %40, label %44
40: ; preds = %29
%41 = load double, double* %13
%42 = load double, double* %9
%43 = fmul double %42, %41
store double %43, double* %9
br label %44
44: ; preds = %40, %29
br label %45
45: ; preds = %44
%46 = load i32, i32* %7
%47 = load i32, i32* %12
%48 = add nsw i32 %47, %46
store i32 %48, i32* %12
br label %25
49: ; preds = %25
br label %50
50: ; preds = %49
%51 = load i32, i32* %11
%52 = add nsw i32 %51, 1
store i32 %52, i32* %11
br label %19
53: ; preds = %19
%54 = load double, double* %9
%55 = fdiv double 1.000000e+00, %54
ret double %55
}

LLVM IR (-O2):
define double @dd(i32 %0, i32 %1, i32 %2, double* nocapture %3) {
%5 = sext i32 %0 to i64
%6 = getelementptr inbounds double, double* %3, i64 %5
%7 = load double, double* %6
%8 = icmp sgt i32 %2, 0
br i1 %8, label %9, label %18
9: ; preds = %4
%10 = zext i32 %2 to i64
%11 = sext i32 %1 to i64
%12 = zext i32 %2 to i64
br label %13
13: ; preds = %9, %21
%14 = phi i64 [ 0, %9 ], [ %23, %21 ]
%15 = phi double [ 1.000000e+00, %9 ], [ %22, %21 ]
%16 = trunc i64 %14 to i32
%17 = icmp slt i32 %16, %1
br i1 %17, label %25, label %21
18: ; preds = %21, %4
%19 = phi double [ 1.000000e+00, %4 ], [ %22, %21 ]
%20 = fdiv double 1.000000e+00, %19
ret double %20
21: ; preds = %25, %13
%22 = phi double [ %15, %13 ], [ %34, %25 ]
%23 = add nuw nsw i64 %14, 1
%24 = icmp eq i64 %23, %12
br i1 %24, label %18, label %13
25: ; preds = %13, %25
%26 = phi i64 [ %35, %25 ], [ %14, %13 ]
%27 = phi double [ %34, %25 ], [ %15, %13 ]
%28 = getelementptr inbounds double, double* %3, i64 %26
%29 = load double, double* %28
%30 = fsub double %7, %29
%31 = fmul double %30, 2.000000e+00
%32 = fcmp une double %31, 0.000000e+00
%33 = select i1 %32, double %31, double 1.000000e+00
%34 = fmul double %27, %33
%35 = add i64 %26, %10
%36 = icmp slt i64 %35, %11
br i1 %36, label %25, label %21
}