// File: /workspace/github_c_makefile_star2_O3/halcy[P]quasi_crystals_bot/Vector/MatrixDouble.c
// Function: dmat4x4affinemul
dmat4x4_t dmat4x4affinemul(dmat4x4_t a,dmat4x4_t b)
{
	return dmat4x4(a.m[0]*b.m[0]+a.m[4]*b.m[1]+a.m[8]*b.m[2],
	              a.m[0]*b.m[4]+a.m[4]*b.m[5]+a.m[8]*b.m[6],
	              a.m[0]*b.m[8]+a.m[4]*b.m[9]+a.m[8]*b.m[10],
	              a.m[0]*b.m[12]+a.m[4]*b.m[13]+a.m[8]*b.m[14]+a.m[12],

	              a.m[1]*b.m[0]+a.m[5]*b.m[1]+a.m[9]*b.m[2],
	              a.m[1]*b.m[4]+a.m[5]*b.m[5]+a.m[9]*b.m[6],
	              a.m[1]*b.m[8]+a.m[5]*b.m[9]+a.m[9]*b.m[10],
	              a.m[1]*b.m[12]+a.m[5]*b.m[13]+a.m[9]*b.m[14]+a.m[13],

	              a.m[2]*b.m[0]+a.m[6]*b.m[1]+a.m[10]*b.m[2],
	              a.m[2]*b.m[4]+a.m[6]*b.m[5]+a.m[10]*b.m[6],
	              a.m[2]*b.m[8]+a.m[6]*b.m[9]+a.m[10]*b.m[10],
	              a.m[2]*b.m[12]+a.m[6]*b.m[13]+a.m[10]*b.m[14]+a.m[14],

	              0,0,0,1);
}

// Assembly:
movsd 0x88(%rsp), %xmm6
movsd 0x90(%rsp), %xmm2
movapd 0x8(%rsp), %xmm12
movapd 0x28(%rsp), %xmm8
movsd 0x98(%rsp), %xmm15
movapd 0x48(%rsp), %xmm7
movsd 0xb0(%rsp), %xmm11
movsd 0xa8(%rsp), %xmm3
movsd 0xe8(%rsp), %xmm0
xorl %ecx, %ecx
movq %rdi, %rax
movddup %xmm6, %xmm1            # xmm1 = xmm6[0,0]
movddup %xmm2, %xmm4            # xmm4 = xmm2[0,0]
movddup %xmm15, %xmm5           # xmm5 = xmm15[0,0]
movddup %xmm11, %xmm9           # xmm9 = xmm11[0,0]
movapd %xmm7, %xmm10
movddup %xmm0, %xmm14           # xmm14 = xmm0[0,0]
movapd %xmm0, -0x58(%rsp)
mulpd %xmm12, %xmm1
mulpd %xmm8, %xmm4
mulpd %xmm7, %xmm5
mulpd %xmm8, %xmm9
mulpd %xmm12, %xmm14
movapd %xmm10, -0x78(%rsp)
addpd %xmm1, %xmm4
movsd 0xc8(%rsp), %xmm1
movapd %xmm5, %xmm7
movddup %xmm3, %xmm5            # xmm5 = xmm3[0,0]
movapd %xmm1, -0x48(%rsp)
movddup %xmm1, %xmm1            # xmm1 = xmm1[0,0]
mulpd %xmm12, %xmm5
mulpd %xmm12, %xmm1
addpd %xmm4, %xmm7
addpd %xmm5, %xmm9
movsd 0x18(%rsp), %xmm5
movapd %xmm1, -0x38(%rsp)
movsd 0xd0(%rsp), %xmm1
movapd %xmm5, %xmm4
movddup %xmm1, %xmm13           # xmm13 = xmm1[0,0]
movapd %xmm1, %xmm0
movsd 0x38(%rsp), %xmm1
mulsd %xmm6, %xmm4
movapd %xmm7, -0x18(%rsp)
movsd 0xb8(%rsp), %xmm7
movapd %xmm1, %xmm12
mulpd %xmm8, %xmm13
mulsd %xmm2, %xmm12
movsd 0xf0(%rsp), %xmm2
movddup %xmm2, %xmm6            # xmm6 = xmm2[0,0]
movapd %xmm2, -0x68(%rsp)
movsd 0x58(%rsp), %xmm2
mulpd %xmm8, %xmm6
movapd %xmm9, -0x28(%rsp)
movddup %xmm7, %xmm9            # xmm9 = xmm7[0,0]
mulpd %xmm10, %xmm9
addsd %xmm4, %xmm12
movapd %xmm2, %xmm4
addpd -0x38(%rsp), %xmm13
mulsd %xmm15, %xmm4
movapd %xmm5, %xmm15
mulsd %xmm3, %xmm15
movapd %xmm1, %xmm3
mulsd %xmm11, %xmm3
movsd 0xd8(%rsp), %xmm11
movddup %xmm11, %xmm8           # xmm8 = xmm11[0,0]
addpd -0x28(%rsp), %xmm9
mulpd %xmm10, %xmm8
movsd 0xf8(%rsp), %xmm10
addsd %xmm12, %xmm4
movapd %xmm2, %xmm12
addpd %xmm14, %xmm6
mulsd %xmm7, %xmm12
movapd %xmm5, %xmm7
mulsd -0x58(%rsp), %xmm5
addsd %xmm15, %xmm3
movapd %xmm1, %xmm15
mulsd -0x68(%rsp), %xmm1
mulsd -0x48(%rsp), %xmm7
mulsd %xmm0, %xmm15
movddup %xmm10, %xmm0           # xmm0 = xmm10[0,0]
mulpd -0x78(%rsp), %xmm0
addsd %xmm5, %xmm1
addsd %xmm3, %xmm12
addsd %xmm7, %xmm15
movapd %xmm2, %xmm7
mulsd %xmm10, %xmm2
addpd %xmm13, %xmm8
mulsd %xmm11, %xmm7
addsd %xmm1, %xmm2
movaps -0x18(%rsp), %xmm1
addpd %xmm6, %xmm0
addsd %xmm15, %xmm7
movups %xmm1, (%rdi)
movsd %xmm4, 0x10(%rdi)
movq %rcx, 0x18(%rdi)
movupd %xmm9, 0x20(%rdi)
movsd %xmm12, 0x30(%rdi)
movq %rcx, 0x38(%rdi)
movupd %xmm8, 0x40(%rdi)
addsd 0x78(%rsp), %xmm2
addpd 0x68(%rsp), %xmm0
movsd %xmm7, 0x50(%rdi)
movq %rcx, 0x58(%rdi)
movabsq $0x3ff0000000000000, %rcx # imm = 0x3FF0000000000000
movupd %xmm0, 0x60(%rdi)
movsd %xmm2, 0x70(%rdi)
movq %rcx, 0x78(%rdi)
retq