void ComputeMc(double **block, int size, double *x, double *y, int s) {
  double row = 0.0;
  double col = 0.0;
  double a = 0.0;
  double b = 0.0;
  double mass = 0.0;
  register int i, j;
  for (i = 0; i < size; i++) {
    a = 0.0;
    b = 0.0;
    for (j = 0; j < size; j++) {
      a += block[i][j];
      b += block[j][i];
      mass += block[j][i];
    }
    row += a * (i + 1);
    col += b * (i + 1);
  }
  if (mass != 0.0) {
    *y = (size + 1) / 2.0 - row / mass;
    *x = col / mass - (size + 1) / 2.0;
  } else {
    *x = 0.0;
    *y = 0.0;
  }
}


x86 assembly (clang, -O0):
<ComputeMc>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1230 <ComputeMc+0x130>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11e0 <ComputeMc+0xe0>
mov    -0x8(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1169 <ComputeMc+0x69>
movsd  -0x40(%rbp),%xmm0
mov    -0x54(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x54(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1146 <ComputeMc+0x46>
movsd  -0x50(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    124d <ComputeMc+0x14d>
jp     124d <ComputeMc+0x14d>
jmp    12aa <ComputeMc+0x1aa>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  0xda1(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  -0x30(%rbp),%xmm1
divsd  -0x50(%rbp),%xmm1
subsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x38(%rbp),%xmm0
divsd  -0x50(%rbp),%xmm0
mov    -0xc(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
movsd  0xd6b(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    12c0 <ComputeMc+0x1c0>
mov    -0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<ComputeMc>:
push   %r14
push   %rbx
test   %esi,%esi
jle    1226 <ComputeMc+0x126>
mov    %esi,%r8d
mov    %r8d,%r9d
and    $0xfffffffe,%r9d
xorpd  %xmm0,%xmm0
xorpd  %xmm1,%xmm1
xor    %r10d,%r10d
jmp    1151 <ComputeMc+0x51>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r10
xorps  %xmm3,%xmm3
cvtsi2sd %r10d,%xmm3
unpcklpd %xmm3,%xmm3
mulpd  %xmm2,%xmm3
addpd  %xmm3,%xmm0
cmp    %r8,%r10
je     11e6 <ComputeMc+0xe6>
mov    (%rdi,%r10,8),%r11
cmp    $0x1,%esi
jne    1170 <ComputeMc+0x70>
xor    %eax,%eax
xorpd  %xmm2,%xmm2
jmp    11bf <ComputeMc+0xbf>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xorpd  %xmm2,%xmm2
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
mov    (%rdi,%rax,8),%rbx
mov    0x8(%rdi,%rax,8),%r14
movsd  (%rbx,%r10,8),%xmm3
addsd  %xmm3,%xmm1
movhpd (%r11,%rax,8),%xmm3
addpd  %xmm2,%xmm3
movsd  (%r14,%r10,8),%xmm4
movapd %xmm4,%xmm2
movhpd 0x8(%r11,%rax,8),%xmm2
addsd  %xmm4,%xmm1
addpd  %xmm3,%xmm2
add    $0x2,%rax
cmp    %rax,%r9
jne    1180 <ComputeMc+0x80>
test   $0x1,%r8b
je     1130 <ComputeMc+0x30>
mov    (%rdi,%rax,8),%rbx
movsd  (%rbx,%r10,8),%xmm3
addsd  %xmm3,%xmm1
movhpd (%r11,%rax,8),%xmm3
addpd  %xmm3,%xmm2
jmp    1130 <ComputeMc+0x30>
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm1
jne    11f2 <ComputeMc+0xf2>
jnp    1226 <ComputeMc+0x126>
add    $0x1,%esi
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
mulsd  0xdfc(%rip),%xmm2
movapd %xmm0,%xmm3
unpckhpd %xmm0,%xmm3
divsd  %xmm1,%xmm3
divsd  %xmm1,%xmm0
subsd  %xmm2,%xmm0
subsd  %xmm3,%xmm2
movsd  %xmm2,(%rcx)
movsd  %xmm0,(%rdx)
jmp    1234 <ComputeMc+0x134>
movq   $0x0,(%rdx)
movq   $0x0,(%rcx)
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @ComputeMc(double** %0, i32 %1, double* %2, double* %3, i32 %4) {
%6 = alloca double**
%7 = alloca i32
%8 = alloca double*
%9 = alloca double*
%10 = alloca i32
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca i32
%17 = alloca i32
store double** %0, double*** %6
store i32 %1, i32* %7
store double* %2, double** %8
store double* %3, double** %9
store i32 %4, i32* %10
store double 0.000000e+00, double* %11
store double 0.000000e+00, double* %12
store double 0.000000e+00, double* %13
store double 0.000000e+00, double* %14
store double 0.000000e+00, double* %15
store i32 0, i32* %16
br label %18
18: ; preds = %77, %5
%19 = load i32, i32* %16
%20 = load i32, i32* %7
%21 = icmp slt i32 %19, %20
br i1 %21, label %22, label %80
22: ; preds = %18
store double 0.000000e+00, double* %13
store double 0.000000e+00, double* %14
store i32 0, i32* %17
br label %23
23: ; preds = %61, %22
%24 = load i32, i32* %17
%25 = load i32, i32* %7
%26 = icmp slt i32 %24, %25
br i1 %26, label %27, label %64
27: ; preds = %23
%28 = load double**, double*** %6
%29 = load i32, i32* %16
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double*, double** %28, i64 %30
%32 = load double*, double** %31
%33 = load i32, i32* %17
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds double, double* %32, i64 %34
%36 = load double, double* %35
%37 = load double, double* %13
%38 = fadd double %37, %36
store double %38, double* %13
%39 = load double**, double*** %6
%40 = load i32, i32* %17
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds double*, double** %39, i64 %41
%43 = load double*, double** %42
%44 = load i32, i32* %16
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds double, double* %43, i64 %45
%47 = load double, double* %46
%48 = load double, double* %14
%49 = fadd double %48, %47
store double %49, double* %14
%50 = load double**, double*** %6
%51 = load i32, i32* %17
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds double*, double** %50, i64 %52
%54 = load double*, double** %53
%55 = load i32, i32* %16
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds double, double* %54, i64 %56
%58 = load double, double* %57
%59 = load double, double* %15
%60 = fadd double %59, %58
store double %60, double* %15
br label %61
61: ; preds = %27
%62 = load i32, i32* %17
%63 = add nsw i32 %62, 1
store i32 %63, i32* %17
br label %23
64: ; preds = %23
%65 = load double, double* %13
%66 = load i32, i32* %16
%67 = add nsw i32 %66, 1
%68 = sitofp i32 %67 to double
%69 = load double, double* %11
%70 = call double @llvm.fmuladd.f64(double %65, double %68, double %69)
store double %70, double* %11
%71 = load double, double* %14
%72 = load i32, i32* %16
%73 = add nsw i32 %72, 1
%74 = sitofp i32 %73 to double
%75 = load double, double* %12
%76 = call double @llvm.fmuladd.f64(double %71, double %74, double %75)
store double %76, double* %12
br label %77
77: ; preds = %64
%78 = load i32, i32* %16
%79 = add nsw i32 %78, 1
store i32 %79, i32* %16
br label %18
80: ; preds = %18
%81 = load double, double* %15
%82 = fcmp une double %81, 0.000000e+00
br i1 %82, label %83, label %102
83: ; preds = %80
%84 = load i32, i32* %7
%85 = add nsw i32 %84, 1
%86 = sitofp i32 %85 to double
%87 = fdiv double %86, 2.000000e+00
%88 = load double, double* %11
%89 = load double, double* %15
%90 = fdiv double %88, %89
%91 = fsub double %87, %90
%92 = load double*, double** %9
store double %91, double* %92
%93 = load double, double* %12
%94 = load double, double* %15
%95 = fdiv double %93, %94
%96 = load i32, i32* %7
%97 = add nsw i32 %96, 1
%98 = sitofp i32 %97 to double
%99 = fdiv double %98, 2.000000e+00
%100 = fsub double %95, %99
%101 = load double*, double** %8
store double %100, double* %101
br label %105
102: ; preds = %80
%103 = load double*, double** %8
store double 0.000000e+00, double* %103
%104 = load double*, double** %9
store double 0.000000e+00, double* %104
br label %105
105: ; preds = %102, %83
ret void
}

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