// File: jamesliu516[P]amrdg2d/wall4amrdg2013.cpp
// Function: getGridIndex(bodygridclass*, pointxy const*, bodygridclass*&)
void getGridIndex(OctCell *parent00, const PXYZ *p_sm, ptrOctCell &pcell)
{
   OctCell *pp[8]={NULL};
  find_sol_refp(pp,bodygrid+parent00->NOparent);
/*
    OctCell *pp[24]={NULL}, *parent; 
    parent=bodygrid+parent00->NOparent;
    pp[0]=EastNeighbor(parent);
    pp[1]=NorthNeighbor(pp[0]);
    pp[2]=NorthNeighbor(parent);
    pp[3]=WestNeighbor(pp[2]);
    pp[4]=WestNeighbor(parent);
    pp[5]=SouthNeighbor(pp[4]);
    pp[6]=SouthNeighbor(parent);
    pp[7]=EastNeighbor(pp[6]);

    pp[8]=NorthNeighbor(pp[1]);

    pp[9]=WestNeighbor(pp[8]);
    pp[10]=WestNeighbor(pp[9]);
    pp[11]=WestNeighbor(pp[10]);
    pp[12]=WestNeighbor(pp[3]);
    pp[13]=WestNeighbor(pp[4]);
    pp[14]=WestNeighbor(pp[5]);
    pp[15]=SouthNeighbor(pp[14]);
    pp[16]=EastNeighbor(pp[15]);
    pp[17]=EastNeighbor(pp[16]);
    pp[18]=EastNeighbor(pp[17]);
    pp[19]=EastNeighbor(pp[18]);
    pp[20]=NorthNeighbor(pp[19]);
    pp[21]=NorthNeighbor(pp[20]);
    pp[22]=NorthNeighbor(pp[21]);
    pp[23]=NorthNeighbor(pp[22]);
    */
    /*
    //================
    pp[24]= NorthNeighbor(pp[23]);
    pp[25]=WestNeighbor(pp[24]);
    pp[26]=WestNeighbor(pp[25]);
    pp[27]=WestNeighbor(pp[26]);  
    pp[28]=WestNeighbor(pp[27]);  
    pp[29]=WestNeighbor(pp[28]);    

    pp[30]=WestNeighbor(pp[11]);
    pp[31]=WestNeighbor(pp[12]);
    pp[32]=WestNeighbor(pp[13]);  
    pp[33]=WestNeighbor(pp[14]);  
    pp[34]=WestNeighbor(pp[15]);    

    pp[35]=SouthNeighbor(pp[34]);

    pp[36]=SouthNeighbor(pp[15]);
    pp[37]=SouthNeighbor(pp[16]);
    pp[38]=SouthNeighbor(pp[17]);
    pp[39]=SouthNeighbor(pp[18]);
    pp[40]=SouthNeighbor(pp[19]);
    pp[41]=EastNeighbor(pp[40]);
    pp[42]=NorthNeighbor(pp[41]);
    pp[43]=NorthNeighbor(pp[42]);  
    pp[44]=NorthNeighbor(pp[43]);  
    pp[45]=NorthNeighbor(pp[44]);  
    pp[46]=NorthNeighbor(pp[45]); 
    pp[47]=NorthNeighbor(pp[46]);  */

    OctCell *smInCell=NULL;//�ԳƵ����ڵ�cell

    smInCell=findPointInCell(p_sm, bodygrid+parent00->NOparent);
    if(smInCell==NULL){ 
        for(int ii=0; ii<8; ++ii) {
            if(pp[ii]!=NULL) {
                smInCell=findPointInCell(p_sm, pp[ii]);
                if(smInCell!=NULL) break;
            }
        } 
    }
#ifdef DEBUG  
    if(smInCell==NULL){
        cout<<"010wrong find symitric point cell\n";
        cout<<parent00->xc1<<","<<parent00->yc1<<endl;
        cout<<p_sm->x<<","<<p_sm->y<<endl;
    }

#endif   
    pcell= smInCell;
//    if(smInCell==NULL)pcell=parent; //2012/3/36
}

// Assembly:
subq $0x68, %rsp
movq %rdi, 0x60(%rsp)
movq %rsi, 0x58(%rsp)
movq %rdx, 0x50(%rsp)
leaq 0x10(%rsp), %rdi
xorl %esi, %esi
movl $0x40, %edx
callq 0x3150
leaq 0x10(%rsp), %rdi
leaq 0x21570(%rip), %rax     # 0x4b370
movq (%rax), %rsi
movq 0x60(%rsp), %rax
movl 0x200(%rax), %eax
imulq $0x360, %rax, %rax      # imm = 0x360
addq %rax, %rsi
callq 0x3b130
movq $0x0, 0x8(%rsp)
movq 0x58(%rsp), %rdi
leaq 0x2153e(%rip), %rax     # 0x4b370
movq (%rax), %rsi
movq 0x60(%rsp), %rax
movl 0x200(%rax), %eax
imulq $0x360, %rax, %rax      # imm = 0x360
addq %rax, %rsi
callq 0x29c10
movq %rax, 0x8(%rsp)
cmpq $0x0, 0x8(%rsp)
jne 0x29eae
movl $0x0, 0x4(%rsp)
cmpl $0x8, 0x4(%rsp)
jge 0x29eac
movslq 0x4(%rsp), %rax
cmpq $0x0, 0x10(%rsp,%rax,8)
je 0x29e9d
movq 0x58(%rsp), %rdi
movslq 0x4(%rsp), %rax
movq 0x10(%rsp,%rax,8), %rsi
callq 0x29c10
movq %rax, 0x8(%rsp)
cmpq $0x0, 0x8(%rsp)
je 0x29e9b
jmp 0x29eac
jmp 0x29e9d
jmp 0x29e9f
movl 0x4(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x4(%rsp)
jmp 0x29e64
jmp 0x29eae
cmpq $0x0, 0x8(%rsp)
jne 0x29f5e
movq 0x2110f(%rip), %rdi     # 0x4afd0
leaq 0x15950(%rip), %rsi     # 0x3f818
callq 0x3240
movq 0x60(%rsp), %rax
movsd 0x278(%rax), %xmm0
movq 0x210ef(%rip), %rdi     # 0x4afd0
callq 0x3480
movq %rax, %rdi
leaq 0x1519f(%rip), %rsi     # 0x3f08f
callq 0x3240
movq %rax, %rdi
movq 0x60(%rsp), %rax
movsd 0x280(%rax), %xmm0
callq 0x3480
movq %rax, %rdi
movq 0x210a4(%rip), %rsi     # 0x4afb8
callq 0x3290
movq 0x58(%rsp), %rax
movsd (%rax), %xmm0
movq 0x210a7(%rip), %rdi     # 0x4afd0
callq 0x3480
movq %rax, %rdi
leaq 0x15157(%rip), %rsi     # 0x3f08f
callq 0x3240
movq %rax, %rdi
movq 0x58(%rsp), %rax
movsd 0x8(%rax), %xmm0
callq 0x3480
movq %rax, %rdi
movq 0x2105f(%rip), %rsi     # 0x4afb8
callq 0x3290
movq 0x8(%rsp), %rcx
movq 0x50(%rsp), %rax
movq %rcx, (%rax)
addq $0x68, %rsp
retq