void mx_antisym(unsigned int n, double *a) {
  unsigned int i, j;
  double t;
  for (i = 0; i < n; i++) {
    a[i * n + i] = 0.0;
    for (j = i + 1; j < n; j++) {
      t = a[i * n + j] - a[j * n + i];
      a[i * n + j] = t;
      a[j * n + i] = -t;
    }
  }
  return;
}


x86 assembly (clang, -O0):
<mx_antisym>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x4(%rbp),%eax
jae    11e6 <mx_antisym+0xe6>
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x14(%rbp),%ecx
mov    %ecx,%ecx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x4(%rbp),%eax
jae    11d3 <mx_antisym+0xd3>
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x18(%rbp),%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x14(%rbp),%ecx
mov    %ecx,%ecx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x18(%rbp),%ecx
mov    %ecx,%ecx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x14(%rbp),%ecx
mov    %ecx,%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    113f <mx_antisym+0x3f>
jmp    11d8 <mx_antisym+0xd8>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1112 <mx_antisym+0x12>
pop    %rbp
ret

x86 assembly (clang, -O2):
<mx_antisym>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edi,%edi
je     13fd <mx_antisym+0x2fd>
mov    %edi,%r9d
lea    -0x2(%rdi),%r14d
mov    %edi,%r12d
neg    %r12d
lea    (%rdi,%rdi,2),%r13d
lea    0x1(%rdi),%eax
mov    %eax,-0x34(%rsp)
lea    0x0(,%rdi,4),%eax
mov    %eax,-0x44(%rsp)
lea    (%rdi,%rdi,1),%r10d
mov    $0x1,%r15d
movl   $0x3,-0x54(%rsp)
xor    %eax,%eax
mov    %rax,-0x50(%rsp)
movapd 0xeaa(%rip),%xmm0
movapd 0xea2(%rip),%xmm1
mov    %r14d,-0x58(%rsp)
mov    %edi,%eax
mov    %edi,-0x5c(%rsp)
xor    %ecx,%ecx
mov    %rdi,-0x28(%rsp)
mov    %r9,-0x30(%rsp)
mov    %r14d,-0x3c(%rsp)
mov    %r12d,-0x40(%rsp)
jmp    11c2 <mx_antisym+0xc2>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%r15
mov    -0x34(%rsp),%eax
add    %eax,%r13d
add    %eax,-0x54(%rsp)
add    %eax,-0x5c(%rsp)
addl   $0xffffffff,-0x58(%rsp)
mov    -0x50(%rsp),%rax
add    %edi,%eax
mov    %rax,-0x50(%rsp)
mov    -0x8(%rsp),%rcx
cmp    %r9,%rcx
je     13fd <mx_antisym+0x2fd>
mov    %ecx,%ebx
mov    %ecx,%r8d
imul   %edi,%r8d
lea    (%r8,%rbx,1),%eax
movq   $0x0,(%rsi,%rax,8)
add    $0x1,%rcx
mov    %rcx,-0x8(%rsp)
cmp    %r9,%rcx
jae    1190 <mx_antisym+0x90>
mov    %r14d,%ecx
mov    %r15,%r11
sub    %ebx,%ecx
je     1340 <mx_antisym+0x240>
mov    %r15,%r11
cmp    $0x1,%edi
jne    1340 <mx_antisym+0x240>
mov    %ebx,%eax
not    %eax
add    %eax,%eax
mov    %r15,%r11
cmp    %ecx,%eax
jb     1340 <mx_antisym+0x240>
add    $0x1,%rcx
mov    %rcx,%r11
movabs $0x1fffffffe,%rax
and    %rax,%r11
lea    -0x2(%r11),%rax
mov    %rax,%rdx
shr    %rdx
add    $0x1,%rdx
test   %rax,%rax
mov    %r13d,-0x38(%rsp)
mov    %rbx,-0x10(%rsp)
mov    %rcx,-0x18(%rsp)
mov    %rdx,-0x20(%rsp)
je     12dc <mx_antisym+0x1dc>
mov    -0x58(%rsp),%edi
add    $0x1,%rdi
and    $0xfffffffffffffffe,%rdi
add    $0xfffffffffffffffe,%rdi
shr    %rdi
add    $0x1,%rdi
and    $0xfffffffffffffffe,%rdi
mov    -0x5c(%rsp),%eax
mov    -0x54(%rsp),%ecx
mov    %ecx,%ebp
mov    %r13d,%r12d
xor    %r14d,%r14d
mov    -0x44(%rsp),%ebx
xchg   %ax,%ax
lea    -0x2(%rbp),%r13d
movupd (%rsi,%r13,8),%xmm2
mov    %eax,%r9d
movupd (%rsi,%r9,8),%xmm3
subpd  %xmm3,%xmm2
movupd %xmm2,(%rsi,%r13,8)
xorpd  %xmm1,%xmm2
movupd %xmm2,(%rsi,%r9,8)
mov    %ebp,%ecx
movupd (%rsi,%rcx,8),%xmm2
mov    %r12d,%edx
movupd (%rsi,%rdx,8),%xmm3
subpd  %xmm3,%xmm2
movupd %xmm2,(%rsi,%rcx,8)
xorpd  %xmm1,%xmm2
movupd %xmm2,(%rsi,%rdx,8)
add    $0x4,%r14
add    %ebx,%r12d
add    $0x4,%ebp
add    %ebx,%eax
add    $0xfffffffffffffffe,%rdi
jne    1280 <mx_antisym+0x180>
jmp    12df <mx_antisym+0x1df>
xor    %r14d,%r14d
testb  $0x1,-0x20(%rsp)
mov    -0x28(%rsp),%rdi
mov    -0x10(%rsp),%rbx
je     131c <mx_antisym+0x21c>
add    %r15d,%r14d
lea    (%r8,%r14,1),%eax
movupd (%rsi,%rax,8),%xmm2
imul   %edi,%r14d
add    %ebx,%r14d
movupd (%rsi,%r14,8),%xmm3
subpd  %xmm3,%xmm2
movupd %xmm2,(%rsi,%rax,8)
xorpd  %xmm1,%xmm2
movupd %xmm2,(%rsi,%r14,8)
cmp    %r11,-0x18(%rsp)
mov    -0x30(%rsp),%r9
mov    -0x3c(%rsp),%r14d
mov    -0x40(%rsp),%r12d
mov    -0x38(%rsp),%r13d
je     1190 <mx_antisym+0x90>
add    %r15,%r11
xchg   %ax,%ax
mov    %edi,%eax
sub    %r11d,%eax
test   $0x1,%al
jne    1350 <mx_antisym+0x250>
mov    %r11,%r8
jmp    137c <mx_antisym+0x27c>
xchg   %ax,%ax
add    %r11d,%r8d
movsd  (%rsi,%r8,8),%xmm2
mov    %r11d,%eax
imul   %edi,%eax
add    %ebx,%eax
subsd  (%rsi,%rax,8),%xmm2
movsd  %xmm2,(%rsi,%r8,8)
xorpd  %xmm1,%xmm2
movlpd %xmm2,(%rsi,%rax,8)
mov    %r11,%r8
add    $0x1,%r8
not    %r11d
cmp    %r12d,%r11d
je     1190 <mx_antisym+0x90>
mov    -0x50(%rsp),%rax
lea    (%rax,%r8,1),%edx
lea    0x1(%r8),%eax
imul   %edi,%eax
add    %ebx,%eax
mov    %edi,%ebp
sub    %r8d,%ebp
imul   %edi,%r8d
add    %ebx,%r8d
cs nopw 0x0(%rax,%rax,1)
mov    %edx,%ecx
movsd  (%rsi,%rcx,8),%xmm2
mov    %r8d,%ebx
subsd  (%rsi,%rbx,8),%xmm2
movsd  %xmm2,(%rsi,%rcx,8)
xorpd  %xmm0,%xmm2
movlpd %xmm2,(%rsi,%rbx,8)
lea    0x1(%rdx),%ecx
movsd  (%rsi,%rcx,8),%xmm2
mov    %eax,%ebx
subsd  (%rsi,%rbx,8),%xmm2
movsd  %xmm2,(%rsi,%rcx,8)
xorpd  %xmm0,%xmm2
movlpd %xmm2,(%rsi,%rbx,8)
add    $0x2,%edx
add    %r10d,%eax
add    %r10d,%r8d
add    $0xfffffffe,%ebp
jne    13b0 <mx_antisym+0x2b0>
jmp    1190 <mx_antisym+0x90>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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