// File: /workspace/github_c_makefile_star2_O0/jiqingxiaoxi[P]RPA/RPA.c
// Function: maxTM
void maxTM(int i,int j,double stackEntropies[],double stackEnthalpies[],double Initdouble[],int Initint[],double enthalpyDPT[],double entropyDPT[],char numSeq1[],char numSeq2[])
{
	double T0,T1,S0,S1,H0,H1;

	S0=entropyDPT[(i-1)*Initint[2]+j-1];
	H0=enthalpyDPT[(i-1)*Initint[2]+j-1];
	T0=(H0+Initdouble[0])/(S0+Initdouble[1]+Initdouble[2]); // at current position 
	if(fabs(enthalpyDPT[(i-2)*Initint[2]+j-2])<999999999&&fabs(Hs(i-1,j-1,1,stackEnthalpies,Initint,numSeq1,numSeq2))<999999999)
	{
		S1=(entropyDPT[(i-2)*Initint[2]+j-2]+Ss(i-1,j-1,1,stackEntropies,Initint,numSeq1,numSeq2));
		H1=(enthalpyDPT[(i-2)*Initint[2]+j-2]+Hs(i-1,j-1,1,stackEnthalpies,Initint,numSeq1,numSeq2));
	}
	else
	{
		S1=-1.0;
		H1=1.0*INFINITY;
	}
	T1=(H1+Initdouble[0])/(S1+Initdouble[1]+Initdouble[2]);

	if(S1<-2500.0)
	{
// to not give dH any value if dS is unreasonable
		S1=-3224.0;
		H1=0.0;
	}
	if(S0<-2500.0)
	{
// to not give dH any value if dS is unreasonable
		S0=-3224.0;
		H0=0.0;
	}
	if((T1>T0)||(S0>0&&H0>0)) // T1 on suurem 
	{
		entropyDPT[(i-1)*Initint[2]+j-1]=S1;
		enthalpyDPT[(i-1)*Initint[2]+j-1]=H1;
	}
	else if(T0>=T1)
	{
		entropyDPT[(i-1)*Initint[2]+j-1]=S0;
		enthalpyDPT[(i-1)*Initint[2]+j-1]=H0;
	}
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x70, %rsp
movq 0x28(%rbp), %rax
movq 0x20(%rbp), %rax
movq 0x18(%rbp), %rax
movq 0x10(%rbp), %rax
movl %edi, -0x4(%rbp)
movl %esi, -0x8(%rbp)
movq %rdx, -0x10(%rbp)
movq %rcx, -0x18(%rbp)
movq %r8, -0x20(%rbp)
movq %r9, -0x28(%rbp)
movq 0x18(%rbp), %rax
movl -0x4(%rbp), %ecx
subl $0x1, %ecx
movq -0x28(%rbp), %rdx
imull 0x8(%rdx), %ecx
addl -0x8(%rbp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, -0x40(%rbp)
movq 0x10(%rbp), %rax
movl -0x4(%rbp), %ecx
subl $0x1, %ecx
movq -0x28(%rbp), %rdx
imull 0x8(%rdx), %ecx
addl -0x8(%rbp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, -0x50(%rbp)
movsd -0x50(%rbp), %xmm0
movq -0x20(%rbp), %rax
addsd (%rax), %xmm0
movsd -0x40(%rbp), %xmm1
movq -0x20(%rbp), %rax
addsd 0x8(%rax), %xmm1
movq -0x20(%rbp), %rax
addsd 0x10(%rax), %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0x30(%rbp)
movq 0x10(%rbp), %rax
movl -0x4(%rbp), %ecx
subl $0x2, %ecx
movq -0x28(%rbp), %rdx
imull 0x8(%rdx), %ecx
addl -0x8(%rbp), %ecx
subl $0x2, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm1
movaps 0xc9e4(%rip), %xmm0     # 0x171f0
pand %xmm0, %xmm1
movsd 0xc920(%rip), %xmm0     # 0x17138
ucomisd %xmm1, %xmm0
jbe 0xa92e
movl -0x4(%rbp), %edi
subl $0x1, %edi
movl -0x8(%rbp), %esi
subl $0x1, %esi
movq -0x18(%rbp), %rcx
movq -0x28(%rbp), %r8
movq 0x20(%rbp), %r9
movq 0x28(%rbp), %rax
movl $0x1, %edx
movq %rax, (%rsp)
callq 0x3c50
movaps %xmm0, %xmm1
movaps 0xc99a(%rip), %xmm0     # 0x171f0
pand %xmm0, %xmm1
movsd 0xc8d6(%rip), %xmm0     # 0x17138
ucomisd %xmm1, %xmm0
jbe 0xa92e
movq 0x18(%rbp), %rax
movl -0x4(%rbp), %ecx
subl $0x2, %ecx
movq -0x28(%rbp), %rdx
imull 0x8(%rdx), %ecx
addl -0x8(%rbp), %ecx
subl $0x2, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, -0x68(%rbp)
movl -0x4(%rbp), %edi
subl $0x1, %edi
movl -0x8(%rbp), %esi
subl $0x1, %esi
movq -0x10(%rbp), %rcx
movq -0x28(%rbp), %r8
movq 0x20(%rbp), %r9
movq 0x28(%rbp), %rax
movl $0x1, %edx
movq %rax, (%rsp)
callq 0x3af0
movaps %xmm0, %xmm1
movsd -0x68(%rbp), %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x48(%rbp)
movq 0x10(%rbp), %rax
movl -0x4(%rbp), %ecx
subl $0x2, %ecx
movq -0x28(%rbp), %rdx
imull 0x8(%rdx), %ecx
addl -0x8(%rbp), %ecx
subl $0x2, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, -0x60(%rbp)
movl -0x4(%rbp), %edi
subl $0x1, %edi
movl -0x8(%rbp), %esi
subl $0x1, %esi
movq -0x18(%rbp), %rcx
movq -0x28(%rbp), %r8
movq 0x20(%rbp), %r9
movq 0x28(%rbp), %rax
movl $0x1, %edx
movq %rax, (%rsp)
callq 0x3c50
movaps %xmm0, %xmm1
movsd -0x60(%rbp), %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x58(%rbp)
jmp 0xa948
movsd 0xc80a(%rip), %xmm0     # 0x17140
movsd %xmm0, -0x48(%rbp)
movsd 0xc7ed(%rip), %xmm0     # 0x17130
movsd %xmm0, -0x58(%rbp)
movsd -0x58(%rbp), %xmm0
movq -0x20(%rbp), %rax
addsd (%rax), %xmm0
movsd -0x48(%rbp), %xmm1
movq -0x20(%rbp), %rax
addsd 0x8(%rax), %xmm1
movq -0x20(%rbp), %rax
addsd 0x10(%rax), %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0x38(%rbp)
movsd 0xc7db(%rip), %xmm0     # 0x17158
ucomisd -0x48(%rbp), %xmm0
jbe 0xa999
movsd 0xc7c4(%rip), %xmm0     # 0x17150
movsd %xmm0, -0x48(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x58(%rbp)
movsd 0xc7b7(%rip), %xmm0     # 0x17158
ucomisd -0x40(%rbp), %xmm0
jbe 0xa9bd
movsd 0xc7a0(%rip), %xmm0     # 0x17150
movsd %xmm0, -0x40(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x50(%rbp)
movsd -0x38(%rbp), %xmm0
ucomisd -0x30(%rbp), %xmm0
ja 0xa9e5
movsd -0x40(%rbp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jbe 0xaa31
movsd -0x50(%rbp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jbe 0xaa31
movsd -0x48(%rbp), %xmm0
movq 0x18(%rbp), %rax
movl -0x4(%rbp), %ecx
subl $0x1, %ecx
movq -0x28(%rbp), %rdx
imull 0x8(%rdx), %ecx
addl -0x8(%rbp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movsd -0x58(%rbp), %xmm0
movq 0x10(%rbp), %rax
movl -0x4(%rbp), %ecx
subl $0x1, %ecx
movq -0x28(%rbp), %rdx
imull 0x8(%rdx), %ecx
addl -0x8(%rbp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
jmp 0xaa89
movsd -0x30(%rbp), %xmm0
ucomisd -0x38(%rbp), %xmm0
jb 0xaa87
movsd -0x40(%rbp), %xmm0
movq 0x18(%rbp), %rax
movl -0x4(%rbp), %ecx
subl $0x1, %ecx
movq -0x28(%rbp), %rdx
imull 0x8(%rdx), %ecx
addl -0x8(%rbp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movsd -0x50(%rbp), %xmm0
movq 0x10(%rbp), %rax
movl -0x4(%rbp), %ecx
subl $0x1, %ecx
movq -0x28(%rbp), %rdx
imull 0x8(%rdx), %ecx
addl -0x8(%rbp), %ecx
subl $0x1, %ecx
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
jmp 0xaa89
addq $0x70, %rsp
popq %rbp
retq
nop