// File: /workspace/github_c_makefile_star2_O2/halcy[P]quasi_crystals_bot/Vector/MatrixDouble.c
// Function: dmat4x4inverse
dmat4x4_t dmat4x4inverse(dmat4x4_t m)
{
	dmat4x4_t res;

	double a0=m.m[0]*m.m[5]-m.m[1]*m.m[4];
	double a1=m.m[0]*m.m[6]-m.m[2]*m.m[4];
	double a2=m.m[0]*m.m[7]-m.m[3]*m.m[4];
	double a3=m.m[1]*m.m[6]-m.m[2]*m.m[5];
	double a4=m.m[1]*m.m[7]-m.m[3]*m.m[5];
	double a5=m.m[2]*m.m[7]-m.m[3]*m.m[6];
	double b0=m.m[8]*m.m[13]-m.m[9]*m.m[12];
	double b1=m.m[8]*m.m[14]-m.m[10]*m.m[12];
	double b2=m.m[8]*m.m[15]-m.m[11]*m.m[12];
	double b3=m.m[9]*m.m[14]-m.m[10]*m.m[13];
	double b4=m.m[9]*m.m[15]-m.m[11]*m.m[13];
	double b5=m.m[10]*m.m[15]-m.m[11]*m.m[14];
	double det=a0*b5-a1*b4+a2*b3+a3*b2-a4*b1+a5*b0;
	// singular if det==0

	res.m[0]=(m.m[5]*b5-m.m[6]*b4+m.m[7]*b3)/det;
	res.m[4]=-(m.m[4]*b5-m.m[6]*b2+m.m[7]*b1)/det;
	res.m[8]=(m.m[4]*b4-m.m[5]*b2+m.m[7]*b0)/det;
	res.m[12]=-(m.m[4]*b3-m.m[5]*b1+m.m[6]*b0)/det;

	res.m[1]=-(m.m[1]*b5-m.m[2]*b4+m.m[3]*b3)/det;
	res.m[5]=(m.m[0]*b5-m.m[2]*b2+m.m[3]*b1)/det;
	res.m[9]=-(m.m[0]*b4-m.m[1]*b2+m.m[3]*b0)/det;
	res.m[13]=(m.m[0]*b3-m.m[1]*b1+m.m[2]*b0)/det;

	res.m[2]=(m.m[13]*a5-m.m[14]*a4+m.m[15]*a3)/det;
	res.m[6]=-(m.m[12]*a5-m.m[14]*a2+m.m[15]*a1)/det;
	res.m[10]=(m.m[12]*a4-m.m[13]*a2+m.m[15]*a0)/det;
	res.m[14]=-(m.m[12]*a3-m.m[13]*a1+m.m[14]*a0)/det;

	res.m[3]=-(m.m[9]*a5-m.m[10]*a4+m.m[11]*a3)/det;
	res.m[7]=(m.m[8]*a5-m.m[10]*a2+m.m[11]*a1)/det;
	res.m[11]=-(m.m[8]*a4-m.m[9]*a2+m.m[11]*a0)/det;
	res.m[15]=(m.m[8]*a3-m.m[9]*a1+m.m[10]*a0)/det;

	return res;
}

// Assembly:
leaq -0x128(%rsp), %rsp
movsd 0x168(%rsp), %xmm9
movsd 0x130(%rsp), %xmm8
movsd 0x140(%rsp), %xmm5
movsd 0x180(%rsp), %xmm4
movsd 0x1a8(%rsp), %xmm6
movsd 0x170(%rsp), %xmm1
movsd 0x190(%rsp), %xmm2
movsd 0x1a0(%rsp), %xmm7
movsd 0x178(%rsp), %xmm3
movupd 0x148(%rsp), %xmm15
movq %rdi, %rax
movapd %xmm9, %xmm0
movapd %xmm5, -0x60(%rsp)
unpcklpd %xmm9, %xmm5            # xmm5 = xmm5[0],xmm9[0]
movapd %xmm6, %xmm11
movapd %xmm1, 0x80(%rsp)
movapd %xmm8, %xmm14
movapd %xmm8, -0x70(%rsp)
movapd %xmm2, 0xf0(%rsp)
movapd %xmm9, (%rsp)
movupd 0x198(%rsp), %xmm9
movapd %xmm6, -0x20(%rsp)
movapd %xmm7, 0x30(%rsp)
movapd %xmm3, 0x110(%rsp)
movapd %xmm4, %xmm10
movapd %xmm4, 0x10(%rsp)
unpcklpd %xmm8, %xmm0            # xmm0 = xmm0[0],xmm8[0]
unpcklpd %xmm1, %xmm11           # xmm11 = xmm11[0],xmm1[0]
movapd %xmm2, %xmm1
movapd %xmm2, %xmm8
movupd 0x188(%rsp), %xmm2
mulpd %xmm0, %xmm5
movapd %xmm4, %xmm0
unpcklpd %xmm4, %xmm1            # xmm1 = xmm1[0],xmm4[0]
movapd %xmm15, -0x50(%rsp)
unpcklpd %xmm6, %xmm0            # xmm0 = xmm0[0],xmm6[0]
movsd 0x188(%rsp), %xmm6
movapd %xmm9, %xmm13
mulpd %xmm0, %xmm11
movapd %xmm7, %xmm0
movapd %xmm6, 0x50(%rsp)
movapd %xmm13, 0xe0(%rsp)
unpcklpd %xmm6, %xmm0            # xmm0 = xmm0[0],xmm6[0]
mulsd %xmm13, %xmm6
mulpd %xmm0, %xmm2
movapd %xmm3, %xmm0
unpcklpd %xmm9, %xmm0            # xmm0 = xmm0[0],xmm9[0]
movupd 0x138(%rsp), %xmm9
mulpd %xmm1, %xmm0
movupd 0x158(%rsp), %xmm1
movapd %xmm9, -0x40(%rsp)
movapd %xmm5, %xmm12
movapd %xmm7, %xmm5
movapd %xmm3, %xmm7
movapd %xmm14, %xmm3
movapd %xmm1, %xmm4
movapd %xmm1, -0x30(%rsp)
unpcklpd %xmm9, %xmm3            # xmm3 = xmm3[0],xmm9[0]
shufpd $0x1, %xmm15, %xmm4     # xmm4 = xmm4[1],xmm15[0]
mulpd %xmm1, %xmm3
subpd %xmm2, %xmm11
movapd %xmm11, %xmm2
movapd %xmm3, %xmm14
movapd %xmm15, %xmm3
shufpd $0x1, %xmm1, %xmm3      # xmm3 = xmm3[1],xmm1[0]
movapd %xmm15, %xmm1
movapd 0x170(%rsp), %xmm15
mulpd %xmm13, %xmm15
mulpd %xmm1, %xmm4
mulpd %xmm9, %xmm3
movsd 0x160(%rsp), %xmm1
movapd %xmm1, %xmm13
movapd %xmm1, 0x20(%rsp)
subpd %xmm0, %xmm15
movapd %xmm12, %xmm0
subpd %xmm3, %xmm14
movapd %xmm10, %xmm3
subpd %xmm4, %xmm0
mulsd %xmm8, %xmm3
movapd %xmm15, %xmm12
mulpd %xmm14, %xmm2
movapd %xmm15, %xmm8
movapd %xmm14, 0x40(%rsp)
movapd (%rsp), %xmm14
mulpd %xmm0, %xmm12
movapd %xmm0, 0x100(%rsp)
movapd %xmm5, %xmm0
movapd -0x20(%rsp), %xmm5
unpckhpd %xmm15, %xmm8           # xmm8 = xmm8[1],xmm15[1]
mulsd 0x80(%rsp), %xmm0
mulsd %xmm7, %xmm5
movapd -0x50(%rsp), %xmm7
movapd %xmm8, -0x10(%rsp)
subsd %xmm3, %xmm0
subsd %xmm6, %xmm5
movsd 0x150(%rsp), %xmm6
movapd %xmm8, %xmm3
movsd %xmm6, -0x78(%rsp)
mulsd %xmm7, %xmm3
haddpd %xmm2, %xmm12
movapd %xmm11, %xmm2
mulsd %xmm9, %xmm2
movapd %xmm0, %xmm4
movddup %xmm1, %xmm0            # xmm0 = xmm1[0,0]
movapd %xmm0, %xmm1
movapd %xmm0, 0x60(%rsp)
movapd -0x60(%rsp), %xmm0
movddup %xmm3, %xmm10           # xmm10 = xmm3[0,0]
movapd %xmm4, 0x70(%rsp)
mulsd %xmm6, %xmm0
movapd %xmm14, %xmm6
mulsd -0x70(%rsp), %xmm1
mulsd %xmm9, %xmm6
movapd -0x30(%rsp), %xmm9
mulsd %xmm9, %xmm7
movddup %xmm2, %xmm2            # xmm2 = xmm2[0,0]
addpd %xmm2, %xmm10
movapd %xmm5, %xmm2
haddpd %xmm12, %xmm12
mulpd %xmm13, %xmm2
movapd %xmm5, %xmm13
movapd %xmm1, %xmm3
movapd %xmm5, %xmm1
subsd %xmm0, %xmm3
movapd %xmm11, %xmm0
subsd %xmm7, %xmm6
movapd %xmm4, %xmm7
mulsd %xmm3, %xmm1
movapd %xmm3, 0xc0(%rsp)
movapd %xmm4, %xmm3
movapd %xmm8, %xmm4
movsd -0x78(%rsp), %xmm8
mulsd %xmm6, %xmm7
movapd %xmm6, 0xb0(%rsp)
movapd %xmm11, %xmm6
mulsd %xmm14, %xmm4
mulsd %xmm8, %xmm0
mulsd %xmm8, %xmm13
movapd -0x40(%rsp), %xmm8
mulsd %xmm9, %xmm6
movapd %xmm14, %xmm9
movapd %xmm3, %xmm14
movddup %xmm5, %xmm3            # xmm3 = xmm5[0,0]
blendpd $0x2, %xmm10, %xmm2     # xmm2 = xmm2[0],xmm10[1]
movapd 0x20(%rsp), %xmm10
mulpd %xmm11, %xmm8
mulpd -0x40(%rsp), %xmm3
mulsd %xmm9, %xmm14
mulsd %xmm15, %xmm9
addsd %xmm1, %xmm7
movsd 0x2216(%rip), %xmm1     # 0xc050
addsd %xmm6, %xmm4
movapd -0x70(%rsp), %xmm6
addpd %xmm0, %xmm14
movapd %xmm11, %xmm0
addpd %xmm13, %xmm9
movapd -0x50(%rsp), %xmm13
mulsd %xmm6, %xmm5
mulsd %xmm6, %xmm0
subsd %xmm7, %xmm12
movapd 0x70(%rsp), %xmm7
movsd %xmm4, %xmm3            # xmm3 = xmm4[0],xmm3[1]
movapd %xmm7, %xmm4
mulsd %xmm13, %xmm4
subpd %xmm2, %xmm3
movddup %xmm0, %xmm0            # xmm0 = xmm0[0,0]
movsd %xmm14, %xmm8           # xmm8 = xmm14[0],xmm8[1]
movapd 0xe0(%rsp), %xmm14
divsd %xmm12, %xmm1
movapd %xmm13, %xmm12
mulsd %xmm15, %xmm12
movddup %xmm4, %xmm2            # xmm2 = xmm4[0,0]
movddup 0x138(%rsp), %xmm4      # xmm4 = mem[0,0]
addpd %xmm0, %xmm2
movapd %xmm11, %xmm0
movapd %xmm4, -0x50(%rsp)
movapd %xmm3, (%rsp)
movapd -0x30(%rsp), %xmm3
unpckhpd %xmm11, %xmm0           # xmm0 = xmm0[1],xmm11[1]
mulpd %xmm4, %xmm11
mulpd %xmm0, %xmm10
movddup %xmm12, %xmm4           # xmm4 = xmm12[0,0]
mulpd %xmm3, %xmm0
movapd 0x10(%rsp), %xmm12
mulpd %xmm7, %xmm3
movsd %xmm9, %xmm11           # xmm11 = xmm9[0],xmm11[1]
movapd %xmm3, -0x30(%rsp)
movapd %xmm12, %xmm3
movapd %xmm1, 0xd0(%rsp)
movddup %xmm5, %xmm1            # xmm1 = xmm5[0,0]
movapd %xmm10, %xmm5
movapd %xmm14, %xmm10
addpd %xmm1, %xmm4
movapd 0x60(%rsp), %xmm1
blendpd $0x2, %xmm2, %xmm5      # xmm5 = xmm5[0],xmm2[1]
movapd -0x10(%rsp), %xmm2
mulsd %xmm15, %xmm1
mulsd -0x60(%rsp), %xmm15
subpd %xmm8, %xmm5
movapd %xmm1, 0x60(%rsp)
movsd -0x78(%rsp), %xmm1
blendpd $0x2, %xmm4, %xmm0      # xmm0 = xmm0[0],xmm4[1]
movapd 0xb0(%rsp), %xmm4
mulsd %xmm2, %xmm1
mulsd %xmm6, %xmm2
subpd %xmm0, %xmm11
movapd %xmm5, 0x20(%rsp)
movapd -0x20(%rsp), %xmm5
mulsd %xmm4, %xmm3
movapd %xmm5, %xmm6
movapd %xmm5, %xmm9
movsd %xmm1, -0x78(%rsp)
movddup %xmm15, %xmm1           # xmm1 = xmm15[0,0]
movddup %xmm2, %xmm0            # xmm0 = xmm2[0,0]
movapd 0x110(%rsp), %xmm15
addpd %xmm0, %xmm1
movapd 0x30(%rsp), %xmm0
movapd %xmm11, -0x40(%rsp)
movapd %xmm0, %xmm11
movapd %xmm1, -0x10(%rsp)
movapd 0x40(%rsp), %xmm1
mulsd %xmm1, %xmm5
movapd %xmm1, %xmm13
unpckhpd %xmm13, %xmm13          # xmm13 = xmm13[1,1]
mulsd %xmm13, %xmm6
movapd %xmm5, -0x20(%rsp)
movapd %xmm0, %xmm5
mulsd %xmm1, %xmm0
mulsd %xmm4, %xmm5
movapd %xmm6, 0x10(%rsp)
movapd 0xf0(%rsp), %xmm6
movapd %xmm6, %xmm8
movapd %xmm0, 0x30(%rsp)
movapd 0x50(%rsp), %xmm0
movapd %xmm0, %xmm2
movapd %xmm0, -0x60(%rsp)
movapd %xmm0, -0x70(%rsp)
movapd %xmm1, %xmm0
movapd %xmm5, 0x90(%rsp)
movapd 0xc0(%rsp), %xmm5
mulsd %xmm12, %xmm0
mulsd %xmm1, %xmm2
mulsd %xmm5, %xmm9
movapd 0x100(%rsp), %xmm1
movapd %xmm1, %xmm7
mulsd %xmm1, %xmm10
mulsd %xmm1, %xmm8
unpckhpd %xmm1, %xmm7            # xmm7 = xmm7[1],xmm1[1]
mulsd %xmm7, %xmm11
mulsd %xmm7, %xmm12
subsd 0x90(%rsp), %xmm10
movapd %xmm0, 0x40(%rsp)
movapd -0x70(%rsp), %xmm0
movapd %xmm2, 0x50(%rsp)
movapd %xmm14, %xmm2
movapd %xmm9, 0xa0(%rsp)
movapd %xmm6, %xmm9
mulsd %xmm7, %xmm14
mulsd %xmm15, %xmm7
mulsd %xmm13, %xmm6
mulsd %xmm5, %xmm0
mulsd %xmm5, %xmm2
mulsd %xmm4, %xmm9
mulsd %xmm15, %xmm5
mulsd %xmm1, %xmm15
subsd %xmm8, %xmm11
movapd 0x10(%rsp), %xmm8
addsd %xmm10, %xmm8
movapd 0x40(%rsp), %xmm10
subsd 0xa0(%rsp), %xmm11
movapd %xmm0, -0x70(%rsp)
movapd -0x60(%rsp), %xmm0
subsd %xmm14, %xmm9
subsd %xmm6, %xmm2
movapd -0x20(%rsp), %xmm6
subsd %xmm15, %xmm3
mulsd %xmm13, %xmm0
subsd 0x30(%rsp), %xmm2
addsd %xmm9, %xmm6
movapd -0x70(%rsp), %xmm9
movapd %xmm0, -0x60(%rsp)
movapd 0x80(%rsp), %xmm0
mulsd %xmm0, %xmm1
mulsd %xmm0, %xmm13
mulsd %xmm0, %xmm4
movddup 0x70(%rsp), %xmm0       # xmm0 = mem[0,0]
subsd -0x60(%rsp), %xmm3
mulpd -0x50(%rsp), %xmm0
subsd %xmm4, %xmm7
movapd 0x60(%rsp), %xmm4
subsd %xmm5, %xmm13
movapd -0x30(%rsp), %xmm5
subsd %xmm12, %xmm1
unpcklpd %xmm3, %xmm8            # xmm8 = xmm8[0],xmm3[0]
movapd (%rsp), %xmm3
addsd -0x78(%rsp), %xmm4
blendpd $0x2, -0x10(%rsp), %xmm5 # xmm5 = xmm5[0],mem[1]
subsd 0x50(%rsp), %xmm7
movapd %xmm5, %xmm12
movapd 0x20(%rsp), %xmm5
addsd %xmm1, %xmm9
addsd %xmm13, %xmm10
movapd %xmm6, %xmm1
movapd -0x40(%rsp), %xmm6
movsd %xmm4, %xmm0            # xmm0 = xmm4[0],xmm0[1]
subpd %xmm0, %xmm12
movddup 0xd0(%rsp), %xmm0       # xmm0 = mem[0,0]
unpcklpd %xmm7, %xmm1            # xmm1 = xmm1[0],xmm7[0]
mulpd %xmm0, %xmm3
mulpd %xmm0, %xmm5
mulpd %xmm0, %xmm6
mulpd %xmm0, %xmm8
mulpd %xmm0, %xmm1
unpcklpd %xmm9, %xmm11           # xmm11 = xmm11[0],xmm9[0]
unpcklpd %xmm10, %xmm2           # xmm2 = xmm2[0],xmm10[0]
mulpd %xmm0, %xmm11
mulpd %xmm0, %xmm2
mulpd %xmm0, %xmm12
movupd %xmm3, (%rdi)
movupd %xmm5, 0x20(%rdi)
movupd %xmm6, 0x40(%rdi)
movupd %xmm12, 0x60(%rdi)
movupd %xmm8, 0x10(%rdi)
movupd %xmm11, 0x30(%rdi)
movupd %xmm1, 0x50(%rdi)
movupd %xmm2, 0x70(%rdi)
leaq 0x128(%rsp), %rsp
retq