// File: Luthaf[P]zeopp/geometry.cc
// Function: invertMatrix(double (*) [3], double (*) [3])
void invertMatrix(double matrix [3][3], double newMatrix[3][3]){
    double determinant = calcDeterminant(matrix);
    if(determinant == 0){
        cerr << "Determinant of provided matrix is 0. Matrix is not invertible."
        << "\n" << "Exiting ..." << "\n";
        exit(1);
    }
    
    double invDet = 1/determinant;
    newMatrix[0][0] = invDet*   (matrix[2][2]*matrix[1][1]-matrix[2][1]*matrix[1][2]);
    newMatrix[0][1] = invDet*-1*(matrix[2][2]*matrix[0][1]-matrix[2][1]*matrix[0][2]);
    newMatrix[0][2] = invDet*   (matrix[1][2]*matrix[0][1]-matrix[1][1]*matrix[0][2]);
    newMatrix[1][0] = invDet*-1*(matrix[2][2]*matrix[1][0]-matrix[2][0]*matrix[1][2]);
    newMatrix[1][1] = invDet*   (matrix[2][2]*matrix[0][0]-matrix[2][0]*matrix[0][2]);
    newMatrix[1][2] = invDet*-1*(matrix[1][2]*matrix[0][0]-matrix[1][0]*matrix[0][2]);
    newMatrix[2][0] = invDet*   (matrix[2][1]*matrix[1][0]-matrix[2][0]*matrix[1][1]);
    newMatrix[2][1] = invDet*-1*(matrix[2][1]*matrix[0][0]-matrix[2][0]*matrix[0][1]);
    newMatrix[2][2] = invDet*   (matrix[1][1]*matrix[0][0]-matrix[1][0]*matrix[0][1]);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x20, %rsp
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movq -0x8(%rbp), %rdi
callq 0xe9390
movsd %xmm0, -0x18(%rbp)
movsd -0x18(%rbp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jne 0xe95b8
jp 0xe95b8
movq 0x189a7b(%rip), %rdi    # 0x272ff0
leaq 0x112774(%rip), %rsi    # 0x1fbcf0
callq 0x5540
movq %rax, %rdi
leaq 0x11009a(%rip), %rsi    # 0x1f9625
callq 0x5540
movq %rax, %rdi
leaq 0x109669(%rip), %rsi    # 0x1f2c03
callq 0x5540
movq %rax, %rdi
leaq 0x11007c(%rip), %rsi    # 0x1f9625
callq 0x5540
movl $0x1, %edi
callq 0x67c0
movsd 0x108a90(%rip), %xmm0   # 0x1f2050
divsd -0x18(%rbp), %xmm0
movsd %xmm0, -0x20(%rbp)
movsd -0x20(%rbp), %xmm0
movq -0x8(%rbp), %rax
movsd 0x40(%rax), %xmm1
movq -0x8(%rbp), %rax
movsd 0x20(%rax), %xmm3
movq -0x8(%rbp), %rax
movsd 0x38(%rax), %xmm2
movq -0x8(%rbp), %rax
mulsd 0x28(%rax), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movq -0x10(%rbp), %rax
movsd %xmm0, (%rax)
movsd 0x109152(%rip), %xmm0   # 0x1f2778
mulsd -0x20(%rbp), %xmm0
movq -0x8(%rbp), %rax
movsd 0x40(%rax), %xmm1
movq -0x8(%rbp), %rax
movsd 0x8(%rax), %xmm3
movq -0x8(%rbp), %rax
movsd 0x38(%rax), %xmm2
movq -0x8(%rbp), %rax
mulsd 0x10(%rax), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movq -0x10(%rbp), %rax
movsd %xmm0, 0x8(%rax)
movsd -0x20(%rbp), %xmm0
movq -0x8(%rbp), %rax
movsd 0x28(%rax), %xmm1
movq -0x8(%rbp), %rax
movsd 0x8(%rax), %xmm3
movq -0x8(%rbp), %rax
movsd 0x20(%rax), %xmm2
movq -0x8(%rbp), %rax
mulsd 0x10(%rax), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movq -0x10(%rbp), %rax
movsd %xmm0, 0x10(%rax)
movsd 0x1090a0(%rip), %xmm0   # 0x1f2778
mulsd -0x20(%rbp), %xmm0
movq -0x8(%rbp), %rax
movsd 0x40(%rax), %xmm1
movq -0x8(%rbp), %rax
movsd 0x18(%rax), %xmm3
movq -0x8(%rbp), %rax
movsd 0x30(%rax), %xmm2
movq -0x8(%rbp), %rax
mulsd 0x28(%rax), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movq -0x10(%rbp), %rax
movsd %xmm0, 0x18(%rax)
movsd -0x20(%rbp), %xmm0
movq -0x8(%rbp), %rax
movsd 0x40(%rax), %xmm1
movq -0x8(%rbp), %rax
movsd (%rax), %xmm3
movq -0x8(%rbp), %rax
movsd 0x30(%rax), %xmm2
movq -0x8(%rbp), %rax
mulsd 0x10(%rax), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movq -0x10(%rbp), %rax
movsd %xmm0, 0x20(%rax)
movsd 0x108fef(%rip), %xmm0   # 0x1f2778
mulsd -0x20(%rbp), %xmm0
movq -0x8(%rbp), %rax
movsd 0x28(%rax), %xmm1
movq -0x8(%rbp), %rax
movsd (%rax), %xmm3
movq -0x8(%rbp), %rax
movsd 0x18(%rax), %xmm2
movq -0x8(%rbp), %rax
mulsd 0x10(%rax), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movq -0x10(%rbp), %rax
movsd %xmm0, 0x28(%rax)
movsd -0x20(%rbp), %xmm0
movq -0x8(%rbp), %rax
movsd 0x38(%rax), %xmm1
movq -0x8(%rbp), %rax
movsd 0x18(%rax), %xmm3
movq -0x8(%rbp), %rax
movsd 0x30(%rax), %xmm2
movq -0x8(%rbp), %rax
mulsd 0x20(%rax), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movq -0x10(%rbp), %rax
movsd %xmm0, 0x30(%rax)
movsd 0x108f3e(%rip), %xmm0   # 0x1f2778
mulsd -0x20(%rbp), %xmm0
movq -0x8(%rbp), %rax
movsd 0x38(%rax), %xmm1
movq -0x8(%rbp), %rax
movsd (%rax), %xmm3
movq -0x8(%rbp), %rax
movsd 0x30(%rax), %xmm2
movq -0x8(%rbp), %rax
mulsd 0x8(%rax), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movq -0x10(%rbp), %rax
movsd %xmm0, 0x38(%rax)
movsd -0x20(%rbp), %xmm0
movq -0x8(%rbp), %rax
movsd 0x20(%rax), %xmm1
movq -0x8(%rbp), %rax
movsd (%rax), %xmm3
movq -0x8(%rbp), %rax
movsd 0x18(%rax), %xmm2
movq -0x8(%rbp), %rax
mulsd 0x8(%rax), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movq -0x10(%rbp), %rax
movsd %xmm0, 0x40(%rax)
addq $0x20, %rsp
popq %rbp
retq
nopl (%rax,%rax)