// File: /workspace/github_c_makefile_star2_O0/Alessandro-Barbieri[P]SpaceLib/src/SPACELIB.c
// Function: trasf_miamit
void trasf_miamit(MAT4 A1, MAT4 m, MAT4 A2)
{
	int i,j,h,k;
	real t,t1,t2;
	real tmp[3];

	for (i=0;i<3;i++)             /* build temporary vector */
	{
		for (t=0.,j=0;j<3;j++)
			t += -m[j][i]*m[j][U];
		tmp[i] = t;
	}
	for (i=0;i<3;i++)             /* 3×3 rotational sub-matrix */
		for (j=0;j<3;j++)
		{
			t=0.; t1=0.; t2=0.;
			for (h=0;h<3;h++)
			{
				t1 += A1[U][h] * m[h][j];
				t2 += A1[h][U] * m[h][i];
				for (k=0;k<3;k++)
					t += m[h][i]*A1[h][k]*m[k][j];
			}
			A2[i][j] =(real) t +t1*tmp[i] +t2*tmp[j] +tmp[i]*tmp[j]*A1[U][U];
		}

	for (i=0;i<3;i++)             /* 3*1 & 1*3 parts */
	{
		t1=0.; t2=0;
		for (j=0;j<3;j++)
		{
			t1 += m[j][i]*A1[j][U];
			t2 += m[j][i]*A1[U][j];
		}
		A2[i][U]=(real)t1+tmp[i]*A1[U][U];
		A2[U][i]=(real)t2+tmp[i]*A1[U][U];
	}
	A2[U][U] = A1[U][U];
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movq %rdx, -0x18(%rbp)
movl $0x0, -0x1c(%rbp)
cmpl $0x3, -0x1c(%rbp)
jge 0xd193
xorps %xmm0, %xmm0
movsd %xmm0, -0x30(%rbp)
movl $0x0, -0x20(%rbp)
cmpl $0x3, -0x20(%rbp)
jge 0xd176
movq -0x10(%rbp), %rax
movslq -0x20(%rbp), %rcx
shlq $0x5, %rcx
addq %rcx, %rax
movslq -0x1c(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
movq -0x10(%rbp), %rax
movslq -0x20(%rbp), %rcx
shlq $0x5, %rcx
addq %rcx, %rax
movsd 0x18(%rax), %xmm2
movsd -0x30(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x30(%rbp)
movl -0x20(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x20(%rbp)
jmp 0xd110
movsd -0x30(%rbp), %xmm0
movslq -0x1c(%rbp), %rax
movsd %xmm0, -0x60(%rbp,%rax,8)
movl -0x1c(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x1c(%rbp)
jmp 0xd0f7
movl $0x0, -0x1c(%rbp)
cmpl $0x3, -0x1c(%rbp)
jge 0xd364
movl $0x0, -0x20(%rbp)
cmpl $0x3, -0x20(%rbp)
jge 0xd354
xorps %xmm0, %xmm0
movsd %xmm0, -0x30(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x38(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x40(%rbp)
movl $0x0, -0x24(%rbp)
cmpl $0x3, -0x24(%rbp)
jge 0xd2d6
movq -0x8(%rbp), %rax
movslq -0x24(%rbp), %rcx
movsd 0x60(%rax,%rcx,8), %xmm0
movq -0x10(%rbp), %rax
movslq -0x24(%rbp), %rcx
shlq $0x5, %rcx
addq %rcx, %rax
movslq -0x20(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm2
movsd -0x38(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x38(%rbp)
movq -0x8(%rbp), %rax
movslq -0x24(%rbp), %rcx
shlq $0x5, %rcx
addq %rcx, %rax
movsd 0x18(%rax), %xmm0
movq -0x10(%rbp), %rax
movslq -0x24(%rbp), %rcx
shlq $0x5, %rcx
addq %rcx, %rax
movslq -0x1c(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm2
movsd -0x40(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x40(%rbp)
movl $0x0, -0x28(%rbp)
cmpl $0x3, -0x28(%rbp)
jge 0xd2c6
movq -0x10(%rbp), %rax
movslq -0x24(%rbp), %rcx
shlq $0x5, %rcx
addq %rcx, %rax
movslq -0x1c(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movq -0x8(%rbp), %rax
movslq -0x24(%rbp), %rcx
shlq $0x5, %rcx
addq %rcx, %rax
movslq -0x28(%rbp), %rcx
mulsd (%rax,%rcx,8), %xmm0
movq -0x10(%rbp), %rax
movslq -0x28(%rbp), %rcx
shlq $0x5, %rcx
addq %rcx, %rax
movslq -0x20(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm2
movsd -0x30(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x30(%rbp)
movl -0x28(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x28(%rbp)
jmp 0xd25b
jmp 0xd2c8
movl -0x24(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x24(%rbp)
jmp 0xd1d4
movsd -0x30(%rbp), %xmm1
movsd -0x38(%rbp), %xmm0
movslq -0x1c(%rbp), %rax
movsd -0x60(%rbp,%rax,8), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd -0x40(%rbp), %xmm1
movslq -0x20(%rbp), %rax
movsd -0x60(%rbp,%rax,8), %xmm2
mulsd %xmm2, %xmm1
addsd %xmm0, %xmm1
movslq -0x1c(%rbp), %rax
movsd -0x60(%rbp,%rax,8), %xmm0
movslq -0x20(%rbp), %rax
mulsd -0x60(%rbp,%rax,8), %xmm0
movq -0x8(%rbp), %rax
movsd 0x78(%rax), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq -0x18(%rbp), %rax
movslq -0x1c(%rbp), %rcx
shlq $0x5, %rcx
addq %rcx, %rax
movslq -0x20(%rbp), %rcx
movsd %xmm0, (%rax,%rcx,8)
movl -0x20(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x20(%rbp)
jmp 0xd1ab
jmp 0xd356
movl -0x1c(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x1c(%rbp)
jmp 0xd19a
movl $0x0, -0x1c(%rbp)
cmpl $0x3, -0x1c(%rbp)
jge 0xd48a
xorps %xmm0, %xmm0
movsd %xmm0, -0x38(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x40(%rbp)
movl $0x0, -0x20(%rbp)
cmpl $0x3, -0x20(%rbp)
jge 0xd41a
movq -0x10(%rbp), %rax
movslq -0x20(%rbp), %rcx
shlq $0x5, %rcx
addq %rcx, %rax
movslq -0x1c(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movq -0x8(%rbp), %rax
movslq -0x20(%rbp), %rcx
shlq $0x5, %rcx
addq %rcx, %rax
movsd 0x18(%rax), %xmm2
movsd -0x38(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x38(%rbp)
movq -0x10(%rbp), %rax
movslq -0x20(%rbp), %rcx
shlq $0x5, %rcx
addq %rcx, %rax
movslq -0x1c(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movq -0x8(%rbp), %rax
movslq -0x20(%rbp), %rcx
movsd 0x60(%rax,%rcx,8), %xmm2
movsd -0x40(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x40(%rbp)
movl -0x20(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x20(%rbp)
jmp 0xd38c
movsd -0x38(%rbp), %xmm1
movslq -0x1c(%rbp), %rax
movsd -0x60(%rbp,%rax,8), %xmm0
movq -0x8(%rbp), %rax
movsd 0x78(%rax), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq -0x18(%rbp), %rax
movslq -0x1c(%rbp), %rcx
shlq $0x5, %rcx
addq %rcx, %rax
movsd %xmm0, 0x18(%rax)
movsd -0x40(%rbp), %xmm1
movslq -0x1c(%rbp), %rax
movsd -0x60(%rbp,%rax,8), %xmm0
movq -0x8(%rbp), %rax
movsd 0x78(%rax), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq -0x18(%rbp), %rax
movslq -0x1c(%rbp), %rcx
movsd %xmm0, 0x60(%rax,%rcx,8)
movl -0x1c(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x1c(%rbp)
jmp 0xd36b
movq -0x8(%rbp), %rax
movsd 0x78(%rax), %xmm0
movq -0x18(%rbp), %rax
movsd %xmm0, 0x78(%rax)
popq %rbp
retq
nop