void set_zero(long L[3][3], const int i, const int j, const long a,
              const long b, const long r, const long s, const long t) {
  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] = s;
  L[i][j] = t;
  L[j][i] = -b / r;
  L[j][j] = a / r;
}


x86 assembly (clang, -O0):
<set_zero>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%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    0x10(%rbp),%rdx
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0xc(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    0x18(%rbp),%rdx
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x10(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
xor    %eax,%eax
sub    -0x20(%rbp),%rax
cqto
idivq  -0x28(%rbp)
mov    %rax,%rdx
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0xc(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
cqto
idivq  -0x28(%rbp)
mov    %rax,%rdx
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x10(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
pop    %rbp
ret

x86 assembly (clang, -O2):
<set_zero>:
push   %rbx
mov    %r8,%rax
mov    %rcx,%r8
mov    0x18(%rsp),%r10
mov    0x10(%rsp),%r11
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,%rcx
lea    (%rcx,%rcx,2),%rsi
lea    (%rdi,%rsi,8),%rbx
mov    %r11,(%rbx,%rcx,8)
movslq %edx,%rsi
mov    %r10,(%rbx,%rsi,8)
mov    %rax,%rdx
or     %r9,%rdx
shr    $0x20,%rdx
je     116c <set_zero+0x6c>
cqto
idiv   %r9
jmp    1171 <set_zero+0x71>
xor    %edx,%edx
div    %r9d
neg    %rax
lea    (%rsi,%rsi,2),%rdx
lea    (%rdi,%rdx,8),%rdi
mov    %rax,(%rdi,%rcx,8)
mov    %r8,%rax
or     %r9,%rax
shr    $0x20,%rax
je     1196 <set_zero+0x96>
mov    %r8,%rax
cqto
idiv   %r9
jmp    119e <set_zero+0x9e>
mov    %r8d,%eax
xor    %edx,%edx
div    %r9d
mov    %rax,(%rdi,%rsi,8)
pop    %rbx
ret

LLVM IR (-O0):
define void @set_zero([3 x i64]* %0, i32 %1, i32 %2, i64 %3, i64 %4, i64 %5, i64 %6, i64 %7) {
%9 = alloca [3 x i64]*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i64
%13 = alloca i64
%14 = alloca i64
%15 = alloca i64
%16 = alloca i64
store [3 x i64]* %0, [3 x i64]** %9
store i32 %1, i32* %10
store i32 %2, i32* %11
store i64 %3, i64* %12
store i64 %4, i64* %13
store i64 %5, i64* %14
store i64 %6, i64* %15
store i64 %7, i64* %16
%17 = load [3 x i64]*, [3 x i64]** %9
%18 = getelementptr inbounds [3 x i64], [3 x i64]* %17, i64 0
%19 = getelementptr inbounds [3 x i64], [3 x i64]* %18, i64 0, i64 0
store i64 1, i64* %19
%20 = load [3 x i64]*, [3 x i64]** %9
%21 = getelementptr inbounds [3 x i64], [3 x i64]* %20, i64 0
%22 = getelementptr inbounds [3 x i64], [3 x i64]* %21, i64 0, i64 1
store i64 0, i64* %22
%23 = load [3 x i64]*, [3 x i64]** %9
%24 = getelementptr inbounds [3 x i64], [3 x i64]* %23, i64 0
%25 = getelementptr inbounds [3 x i64], [3 x i64]* %24, i64 0, i64 2
store i64 0, i64* %25
%26 = load [3 x i64]*, [3 x i64]** %9
%27 = getelementptr inbounds [3 x i64], [3 x i64]* %26, i64 1
%28 = getelementptr inbounds [3 x i64], [3 x i64]* %27, i64 0, i64 0
store i64 0, i64* %28
%29 = load [3 x i64]*, [3 x i64]** %9
%30 = getelementptr inbounds [3 x i64], [3 x i64]* %29, i64 1
%31 = getelementptr inbounds [3 x i64], [3 x i64]* %30, i64 0, i64 1
store i64 1, i64* %31
%32 = load [3 x i64]*, [3 x i64]** %9
%33 = getelementptr inbounds [3 x i64], [3 x i64]* %32, i64 1
%34 = getelementptr inbounds [3 x i64], [3 x i64]* %33, i64 0, i64 2
store i64 0, i64* %34
%35 = load [3 x i64]*, [3 x i64]** %9
%36 = getelementptr inbounds [3 x i64], [3 x i64]* %35, i64 2
%37 = getelementptr inbounds [3 x i64], [3 x i64]* %36, i64 0, i64 0
store i64 0, i64* %37
%38 = load [3 x i64]*, [3 x i64]** %9
%39 = getelementptr inbounds [3 x i64], [3 x i64]* %38, i64 2
%40 = getelementptr inbounds [3 x i64], [3 x i64]* %39, i64 0, i64 1
store i64 0, i64* %40
%41 = load [3 x i64]*, [3 x i64]** %9
%42 = getelementptr inbounds [3 x i64], [3 x i64]* %41, i64 2
%43 = getelementptr inbounds [3 x i64], [3 x i64]* %42, i64 0, i64 2
store i64 1, i64* %43
%44 = load i64, i64* %15
%45 = load [3 x i64]*, [3 x i64]** %9
%46 = load i32, i32* %10
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds [3 x i64], [3 x i64]* %45, i64 %47
%49 = load i32, i32* %10
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds [3 x i64], [3 x i64]* %48, i64 0, i64 %50
store i64 %44, i64* %51
%52 = load i64, i64* %16
%53 = load [3 x i64]*, [3 x i64]** %9
%54 = load i32, i32* %10
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds [3 x i64], [3 x i64]* %53, i64 %55
%57 = load i32, i32* %11
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds [3 x i64], [3 x i64]* %56, i64 0, i64 %58
store i64 %52, i64* %59
%60 = load i64, i64* %13
%61 = sub nsw i64 0, %60
%62 = load i64, i64* %14
%63 = sdiv i64 %61, %62
%64 = load [3 x i64]*, [3 x i64]** %9
%65 = load i32, i32* %11
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds [3 x i64], [3 x i64]* %64, i64 %66
%68 = load i32, i32* %10
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds [3 x i64], [3 x i64]* %67, i64 0, i64 %69
store i64 %63, i64* %70
%71 = load i64, i64* %12
%72 = load i64, i64* %14
%73 = sdiv i64 %71, %72
%74 = load [3 x i64]*, [3 x i64]** %9
%75 = load i32, i32* %11
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds [3 x i64], [3 x i64]* %74, i64 %76
%78 = load i32, i32* %11
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds [3 x i64], [3 x i64]* %77, i64 0, i64 %79
store i64 %73, i64* %80
ret void
}

LLVM IR (-O2):
define void @set_zero([3 x i64]* nocapture writeonly %0, i32 %1, i32 %2, i64 %3, i64 %4, i64 %5, i64 %6, i64 %7) {
%9 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 0, i64 0
store i64 1, i64* %9
%10 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 0, i64 1
%11 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 1, i64 1
%12 = bitcast i64* %10 to i8*
store i64 1, i64* %11
%13 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 1, i64 2
%14 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 2, i64 2
%15 = bitcast i64* %13 to i8*
store i64 1, i64* %14
%16 = sext i32 %1 to i64
%17 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 %16, i64 %16
store i64 %6, i64* %17
%18 = sext i32 %2 to i64
%19 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 %16, i64 %18
store i64 %7, i64* %19
%20 = sdiv i64 %4, %5
%21 = sub nsw i64 0, %20
%22 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 %18, i64 %16
store i64 %21, i64* %22
%23 = sdiv i64 %3, %5
%24 = getelementptr inbounds [3 x i64], [3 x i64]* %0, i64 %18, i64 %18
store i64 %23, i64* %24
ret void
}