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