// File: jamesliu516[P]amrdg2d/wall4amrdg2013.cpp
// Function: getRadius(pointxy const*, pointxy const*, pointxy const*, pointxy&)
double getRadius(const PXYZ *pt1, const PXYZ *pt2,const PXYZ *pt3, PXYZ &centr)
{
    //Ҫ�����㹲�߻��߽ӽ������������Ŀ��Ҳ��� TODO
    double rd,dls1,dls2,dls3,dls4;
    dls1=(pt1->x-pt2->x)*(pt1->x-pt2->x)+(pt1->y-pt2->y)*(pt1->y-pt2->y);
    dls2=(pt1->x-pt3->x)*(pt1->x-pt3->x)+(pt1->y-pt3->y)*(pt1->y-pt3->y);
    dls3=(pt2->x-pt3->x)*(pt2->x-pt3->x)+(pt2->y-pt3->y)*(pt2->y-pt3->y);
    dls4=2.0*(pt3->x * (pt2->y - pt1->y) + pt2->x * (pt1->y - pt3->y) 
        + pt1->x * (pt3->y-pt2->y));
    rd = sqrt(dls1*dls2*dls3)/(fabs(dls4)+1.0e-18);

    double xc12,yc12;
    xc12=0.5*(pt1->x+pt2->x);
    yc12=0.5*(pt1->y+pt2->y);

    double dis1c, dis;
    dis1c=sqrt((pt1->x-xc12)*(pt1->x-xc12)+(pt1->y-yc12)*(pt1->y-yc12));
    dis=sqrt(rd*rd-dis1c*dis1c);

    //   PXYZ tdir;
    PXYZ nml, cen1;       
    //  tdir.x=(pt1->x-xc12)/dis1c;
    // tdir.y=(pt1->y-yc12)/dis1c;       

    nml.x=(pt1->y-yc12)/dis1c;  
    nml.y=-(pt1->x-xc12)/dis1c;

    cen1.x= xc12+nml.x*dis;
    cen1.y= yc12+nml.y*dis;

    if(fabs(rd*rd-(cen1.x-pt3->x)*(cen1.x-pt3->x)-(cen1.y-pt3->y)*(cen1.y-pt3->y))<ERRS) {
        centr=cen1;
    }
    else {
        centr.x= xc12-nml.x*dis;
        centr.y= yc12-nml.y*dis;        
    }  

    double a1,b1,a2,b2;
    a1=pt2->x-pt1->x;
    b1=pt2->y-pt1->y;

    a2=pt3->x-pt1->x;
    b2=pt3->y-pt1->y;
    if(fabs(a1*b2-b1*a2)<1.0e-12) rd=1.0e15;
    return rd;
}

// Assembly:
subq $0xa8, %rsp
movq %rdi, 0xa0(%rsp)
movq %rsi, 0x98(%rsp)
movq %rdx, 0x90(%rsp)
movq %rcx, 0x88(%rsp)
movq 0xa0(%rsp), %rax
movsd (%rax), %xmm0
movq 0x98(%rsp), %rax
subsd (%rax), %xmm0
movq 0xa0(%rsp), %rax
movsd (%rax), %xmm2
movq 0x98(%rsp), %rax
subsd (%rax), %xmm2
movq 0xa0(%rsp), %rax
movsd 0x8(%rax), %xmm1
movq 0x98(%rsp), %rax
subsd 0x8(%rax), %xmm1
movq 0xa0(%rsp), %rax
movsd 0x8(%rax), %xmm3
movq 0x98(%rsp), %rax
subsd 0x8(%rax), %xmm3
mulsd %xmm3, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x78(%rsp)
movq 0xa0(%rsp), %rax
movsd (%rax), %xmm0
movq 0x90(%rsp), %rax
subsd (%rax), %xmm0
movq 0xa0(%rsp), %rax
movsd (%rax), %xmm2
movq 0x90(%rsp), %rax
subsd (%rax), %xmm2
movq 0xa0(%rsp), %rax
movsd 0x8(%rax), %xmm1
movq 0x90(%rsp), %rax
subsd 0x8(%rax), %xmm1
movq 0xa0(%rsp), %rax
movsd 0x8(%rax), %xmm3
movq 0x90(%rsp), %rax
subsd 0x8(%rax), %xmm3
mulsd %xmm3, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x70(%rsp)
movq 0x98(%rsp), %rax
movsd (%rax), %xmm0
movq 0x90(%rsp), %rax
subsd (%rax), %xmm0
movq 0x98(%rsp), %rax
movsd (%rax), %xmm2
movq 0x90(%rsp), %rax
subsd (%rax), %xmm2
movq 0x98(%rsp), %rax
movsd 0x8(%rax), %xmm1
movq 0x90(%rsp), %rax
subsd 0x8(%rax), %xmm1
movq 0x98(%rsp), %rax
movsd 0x8(%rax), %xmm3
movq 0x90(%rsp), %rax
subsd 0x8(%rax), %xmm3
mulsd %xmm3, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x68(%rsp)
movq 0x90(%rsp), %rax
movsd (%rax), %xmm0
movq 0x98(%rsp), %rax
movsd 0x8(%rax), %xmm2
movq 0xa0(%rsp), %rax
subsd 0x8(%rax), %xmm2
movq 0x98(%rsp), %rax
movsd (%rax), %xmm1
movq 0xa0(%rsp), %rax
movsd 0x8(%rax), %xmm3
movq 0x90(%rsp), %rax
subsd 0x8(%rax), %xmm3
mulsd %xmm3, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0xa0(%rsp), %rax
movsd (%rax), %xmm1
movq 0x90(%rsp), %rax
movsd 0x8(%rax), %xmm2
movq 0x98(%rsp), %rax
subsd 0x8(%rax), %xmm2
mulsd %xmm2, %xmm1
addsd %xmm0, %xmm1
movsd 0x165f9(%rip), %xmm0    # 0x3f0a0
mulsd %xmm1, %xmm0
movsd %xmm0, 0x60(%rsp)
movsd 0x78(%rsp), %xmm0
mulsd 0x70(%rsp), %xmm0
mulsd 0x68(%rsp), %xmm0
callq 0x31c0
movsd 0x60(%rsp), %xmm1
movaps 0x1672b(%rip), %xmm2    # 0x3f200
pand %xmm2, %xmm1
movsd 0x16d1f(%rip), %xmm2    # 0x3f800
addsd %xmm2, %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, 0x80(%rsp)
movq 0xa0(%rsp), %rax
movsd (%rax), %xmm1
movq 0x98(%rsp), %rax
addsd (%rax), %xmm1
movsd 0x1662e(%rip), %xmm0    # 0x3f140
mulsd %xmm1, %xmm0
movsd %xmm0, 0x58(%rsp)
movq 0xa0(%rsp), %rax
movsd 0x8(%rax), %xmm1
movq 0x98(%rsp), %rax
addsd 0x8(%rax), %xmm1
movsd 0x16602(%rip), %xmm0    # 0x3f140
mulsd %xmm1, %xmm0
movsd %xmm0, 0x50(%rsp)
movq 0xa0(%rsp), %rax
movsd (%rax), %xmm0
subsd 0x58(%rsp), %xmm0
movq 0xa0(%rsp), %rax
movsd (%rax), %xmm2
subsd 0x58(%rsp), %xmm2
movq 0xa0(%rsp), %rax
movsd 0x8(%rax), %xmm1
subsd 0x50(%rsp), %xmm1
movq 0xa0(%rsp), %rax
movsd 0x8(%rax), %xmm3
subsd 0x50(%rsp), %xmm3
mulsd %xmm3, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
callq 0x31c0
movsd %xmm0, 0x48(%rsp)
movsd 0x80(%rsp), %xmm0
movsd 0x80(%rsp), %xmm2
movsd 0x48(%rsp), %xmm1
mulsd 0x48(%rsp), %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
callq 0x31c0
movsd %xmm0, 0x40(%rsp)
leaq 0x30(%rsp), %rdi
callq 0x6310
leaq 0x20(%rsp), %rdi
callq 0x6310
movq 0xa0(%rsp), %rax
movsd 0x8(%rax), %xmm0
subsd 0x50(%rsp), %xmm0
divsd 0x48(%rsp), %xmm0
movsd %xmm0, 0x30(%rsp)
movq 0xa0(%rsp), %rax
movsd (%rax), %xmm0
subsd 0x58(%rsp), %xmm0
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
divsd 0x48(%rsp), %xmm0
movsd %xmm0, 0x38(%rsp)
movsd 0x58(%rsp), %xmm1
movsd 0x30(%rsp), %xmm0
movsd 0x40(%rsp), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x20(%rsp)
movsd 0x50(%rsp), %xmm1
movsd 0x38(%rsp), %xmm0
movsd 0x40(%rsp), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x28(%rsp)
movsd 0x80(%rsp), %xmm0
movsd 0x80(%rsp), %xmm2
movsd 0x20(%rsp), %xmm1
movq 0x90(%rsp), %rax
subsd (%rax), %xmm1
movsd 0x20(%rsp), %xmm3
movq 0x90(%rsp), %rax
subsd (%rax), %xmm3
mulsd %xmm3, %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd 0x28(%rsp), %xmm1
movq 0x90(%rsp), %rax
subsd 0x8(%rax), %xmm1
movsd 0x28(%rsp), %xmm2
movq 0x90(%rsp), %rax
subsd 0x8(%rax), %xmm2
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm2, %xmm1
addsd %xmm0, %xmm1
movaps 0x164c2(%rip), %xmm0    # 0x3f200
pand %xmm0, %xmm1
movsd 0x163fe(%rip), %xmm0    # 0x3f148
ucomisd %xmm1, %xmm0
jbe 0x28d6b
movq 0x88(%rsp), %rax
movq 0x20(%rsp), %rcx
movq %rcx, (%rax)
movq 0x28(%rsp), %rcx
movq %rcx, 0x8(%rax)
jmp 0x28de6
movsd 0x58(%rsp), %xmm1
movsd 0x30(%rsp), %xmm0
movsd 0x40(%rsp), %xmm2
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x88(%rsp), %rax
movsd %xmm0, (%rax)
movsd 0x50(%rsp), %xmm1
movsd 0x38(%rsp), %xmm0
movsd 0x40(%rsp), %xmm2
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x88(%rsp), %rax
movsd %xmm0, 0x8(%rax)
movq 0x98(%rsp), %rax
movsd (%rax), %xmm0
movq 0xa0(%rsp), %rax
subsd (%rax), %xmm0
movsd %xmm0, 0x18(%rsp)
movq 0x98(%rsp), %rax
movsd 0x8(%rax), %xmm0
movq 0xa0(%rsp), %rax
subsd 0x8(%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movq 0x90(%rsp), %rax
movsd (%rax), %xmm0
movq 0xa0(%rsp), %rax
subsd (%rax), %xmm0
movsd %xmm0, 0x8(%rsp)
movq 0x90(%rsp), %rax
movsd 0x8(%rax), %xmm0
movq 0xa0(%rsp), %rax
subsd 0x8(%rax), %xmm0
movsd %xmm0, (%rsp)
movsd 0x18(%rsp), %xmm1
movsd (%rsp), %xmm2
movsd 0x10(%rsp), %xmm0
mulsd 0x8(%rsp), %xmm0
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
mulsd %xmm2, %xmm1
addsd %xmm0, %xmm1
movaps 0x16362(%rip), %xmm0    # 0x3f200
pand %xmm0, %xmm1
movsd 0x1694e(%rip), %xmm0    # 0x3f7f8
ucomisd %xmm1, %xmm0
jbe 0x28ec1
movsd 0x16950(%rip), %xmm0    # 0x3f808
movsd %xmm0, 0x80(%rsp)
movsd 0x80(%rsp), %xmm0
addq $0xa8, %rsp
retq
nopw %cs:(%rax,%rax)