void disturb_rows(long L[3][3], const int i, const int j) {
  L[0][0] = 1;
  L[0][1] = 0;
  L[0][2] = 0;
  L[1][0] = 0;
  L[1][1] = 1;
  L[1][2] = 0;
  L[2][0] = 0;
  L[2][1] = 0;
  L[2][2] = 1;
  L[i][i] = 1;
  L[i][j] = 1;
  L[j][i] = 0;
  L[j][j] = 1;
}


x86 assembly (clang, -O0):
<disturb_rows>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    -0x8(%rbp),%rax
movq   $0x1,(%rax)
mov    -0x8(%rbp),%rax
movq   $0x0,0x8(%rax)
mov    -0x8(%rbp),%rax
movq   $0x0,0x10(%rax)
mov    -0x8(%rbp),%rax
movq   $0x0,0x18(%rax)
mov    -0x8(%rbp),%rax
movq   $0x1,0x20(%rax)
mov    -0x8(%rbp),%rax
movq   $0x0,0x28(%rax)
mov    -0x8(%rbp),%rax
movq   $0x0,0x30(%rax)
mov    -0x8(%rbp),%rax
movq   $0x0,0x38(%rax)
mov    -0x8(%rbp),%rax
movq   $0x1,0x40(%rax)
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0xc(%rbp),%rcx
movq   $0x1,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x10(%rbp),%rcx
movq   $0x1,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0xc(%rbp),%rcx
movq   $0x0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x10(%rbp),%rcx
movq   $0x1,(%rax,%rcx,8)
pop    %rbp
ret

x86 assembly (clang, -O2):
<disturb_rows>:
movq   $0x1,(%rdi)
xorps  %xmm0,%xmm0
movups %xmm0,0x8(%rdi)
movq   $0x0,0x18(%rdi)
movq   $0x1,0x20(%rdi)
movups %xmm0,0x28(%rdi)
movq   $0x0,0x38(%rdi)
movq   $0x1,0x40(%rdi)
movslq %esi,%rax
lea    (%rax,%rax,2),%rcx
lea    (%rdi,%rcx,8),%rcx
movq   $0x1,(%rcx,%rax,8)
movslq %edx,%rdx
movq   $0x1,(%rcx,%rdx,8)
lea    (%rdx,%rdx,2),%rcx
lea    (%rdi,%rcx,8),%rcx
movq   $0x0,(%rcx,%rax,8)
movq   $0x1,(%rcx,%rdx,8)
ret

LLVM IR (-O0):
define void @disturb_rows([3 x i64]* %0, i32 %1, i32 %2) {
%4 = alloca [3 x i64]*
%5 = alloca i32
%6 = alloca i32
store [3 x i64]* %0, [3 x i64]** %4
store i32 %1, i32* %5
store i32 %2, i32* %6
%7 = load [3 x i64]*, [3 x i64]** %4
%8 = getelementptr inbounds [3 x i64], [3 x i64]* %7, i64 0
%9 = getelementptr inbounds [3 x i64], [3 x i64]* %8, i64 0, i64 0
store i64 1, i64* %9
%10 = load [3 x i64]*, [3 x i64]** %4
%11 = getelementptr inbounds [3 x i64], [3 x i64]* %10, i64 0
%12 = getelementptr inbounds [3 x i64], [3 x i64]* %11, i64 0, i64 1
store i64 0, i64* %12
%13 = load [3 x i64]*, [3 x i64]** %4
%14 = getelementptr inbounds [3 x i64], [3 x i64]* %13, i64 0
%15 = getelementptr inbounds [3 x i64], [3 x i64]* %14, i64 0, i64 2
store i64 0, i64* %15
%16 = load [3 x i64]*, [3 x i64]** %4
%17 = getelementptr inbounds [3 x i64], [3 x i64]* %16, i64 1
%18 = getelementptr inbounds [3 x i64], [3 x i64]* %17, i64 0, i64 0
store i64 0, i64* %18
%19 = load [3 x i64]*, [3 x i64]** %4
%20 = getelementptr inbounds [3 x i64], [3 x i64]* %19, i64 1
%21 = getelementptr inbounds [3 x i64], [3 x i64]* %20, i64 0, i64 1
store i64 1, i64* %21
%22 = load [3 x i64]*, [3 x i64]** %4
%23 = getelementptr inbounds [3 x i64], [3 x i64]* %22, i64 1
%24 = getelementptr inbounds [3 x i64], [3 x i64]* %23, i64 0, i64 2
store i64 0, i64* %24
%25 = load [3 x i64]*, [3 x i64]** %4
%26 = getelementptr inbounds [3 x i64], [3 x i64]* %25, i64 2
%27 = getelementptr inbounds [3 x i64], [3 x i64]* %26, i64 0, i64 0
store i64 0, i64* %27
%28 = load [3 x i64]*, [3 x i64]** %4
%29 = getelementptr inbounds [3 x i64], [3 x i64]* %28, i64 2
%30 = getelementptr inbounds [3 x i64], [3 x i64]* %29, i64 0, i64 1
store i64 0, i64* %30
%31 = load [3 x i64]*, [3 x i64]** %4
%32 = getelementptr inbounds [3 x i64], [3 x i64]* %31, i64 2
%33 = getelementptr inbounds [3 x i64], [3 x i64]* %32, i64 0, i64 2
store i64 1, i64* %33
%34 = load [3 x i64]*, [3 x i64]** %4
%35 = load i32, i32* %5
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds [3 x i64], [3 x i64]* %34, i64 %36
%38 = load i32, i32* %5
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds [3 x i64], [3 x i64]* %37, i64 0, i64 %39
store i64 1, i64* %40
%41 = load [3 x i64]*, [3 x i64]** %4
%42 = load i32, i32* %5
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds [3 x i64], [3 x i64]* %41, i64 %43
%45 = load i32, i32* %6
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds [3 x i64], [3 x i64]* %44, i64 0, i64 %46
store i64 1, i64* %47
%48 = load [3 x i64]*, [3 x i64]** %4
%49 = load i32, i32* %6
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds [3 x i64], [3 x i64]* %48, i64 %50
%52 = load i32, i32* %5
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds [3 x i64], [3 x i64]* %51, i64 0, i64 %53
store i64 0, i64* %54
%55 = load [3 x i64]*, [3 x i64]** %4
%56 = load i32, i32* %6
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds [3 x i64], [3 x i64]* %55, i64 %57
%59 = load i32, i32* %6
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds [3 x i64], [3 x i64]* %58, i64 0, i64 %60
store i64 1, i64* %61
ret void
}

LLVM IR (-O2):
define void @disturb_rows([3 x i64]* nocapture writeonly %0, i32 %1, i32 %2) {
%4 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 0, i64 0
store i64 1, i64* %4
%5 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 0, i64 1
%6 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 1, i64 1
%7 = bitcast i64* %5 to i8*
store i64 1, i64* %6
%8 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 1, i64 2
%9 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 2, i64 2
%10 = bitcast i64* %8 to i8*
store i64 1, i64* %9
%11 = sext i32 %1 to i64
%12 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 %11, i64 %11
store i64 1, i64* %12
%13 = sext i32 %2 to i64
%14 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 %11, i64 %13
store i64 1, i64* %14
%15 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 %13, i64 %11
store i64 0, i64* %15
%16 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 %13, i64 %13
store i64 1, i64* %16
ret void
}