// File: /workspace/github_cpp_makefile_star2_O3/jamesliu516[P]amrdg2d/tvbm_limiter4all.cpp
// Function: pos_limiter(bodygridclass*)
void pos_limiter(OctCell *pc)
{
    double dx,dy,xc,yc;
    double dof[4][6]={0.0};
    double cgsx[3];
    double cgsy[3];
    
    double cgslbx[4]={0.0};
    double cgslby[4]={0.0};

    double gslbwei[4]={0.0};
    double gslbpt[4]={0.0};
    double qmin=0.0, rtp,tht1,tht2, pre0;
    double u, v, gma1;
    double pre1;
    double ush[4],ujb[4];

    dx=pc->dx;
    dy=pc->dy;
    xc=pc->xc1;
    yc=pc->yc1;

    for(int i=0;i<4;i++) 
        for(int j=0; j< nDOF; j++) dof[i][j]=pc->dof[i][j];


    if(NLBpt==3){
        for(int i=0; i<3;++i){
            gslbwei[i]=weiLB3[i];
            gslbpt[i]=gsLB3pt[i];
        }
    }
    else{ 
        for(int i=0; i<4;++i){
            gslbwei[i]=weiLB4[i];
            gslbpt[i]=gsLB4pt[i];
        }
    }

    for(int i=0; i<3;++i){
        cgsx[i]=xc+0.5*dx*gspt[i];
        cgsy[i]=yc+0.5*dy*gspt[i];
    }

    for(int i=0; i<NLBpt;++i){
        cgslbx[i]=xc+0.5*dx*gslbpt[i];
        cgslby[i]=yc+0.5*dy*gslbpt[i];
    }

//limit density
    if(dof[0][0]<epsPos){
    //    dof[0][0]=epsPos;
        for(int mm=1;mm<nDOF;++mm)
            dof[0][mm]=0.0;
    }
    else{
        qmin=1.0e20;
        for(int i=0;i<3;++i){
            for(int j=0;j<NLBpt;++j){
                rtp=n_approx_sol(cgslbx[j],cgsy[i],
                    xc,yc,dx,dy, dof[0],nDOF);
                if(rtp<qmin) qmin=rtp;
            }
        }

        for(int i=0;i<3;++i){
            for(int j=0;j<NLBpt;++j){
                rtp=n_approx_sol(cgsx[i],cgslby[j],
                    xc,yc,dx,dy, dof[0],nDOF);
                if(rtp<qmin) qmin=rtp;
            }
        }

        tht1=min2(fabs((dof[0][0]-epsPos)/(dof[0][0]-qmin)), 1.0);

        for(int i=1;i<nDOF;++i)
            dof[0][i]*=tht1;
    }


    gma1=GAMMA-1.0;

    u=dof[1][0]/dof[0][0];
    v=dof[2][0]/dof[0][0];

    pre0=gma1*(dof[3][0]-0.5*dof[0][0]*(u*u+v*v));

    if(pre0<epsPos){
        for(int i=0;i<4;++i){
            for(int j=1;j<nDOF;j++){
                dof[i][j]=0.0;
            }
        }
    }
    else{
        tht2=1.0e20;
        for(int i=0;i<3;++i){
            for(int j=0;j<NLBpt;++j){
                for(int a=0;a<4;++a){
                    ush[a]=n_approx_sol(cgsx[i],cgslby[j],
                        xc,yc,dx,dy, dof[a],nDOF);
                }
                shbl2jbbl(ush,ujb);
                pre1=ujb[3];
                if(pre1>epsPos){
                    rtp=1.0;
                }
                else{
                    rtp=(pre0-epsPos)/(pre0-pre1);
                }

                if(rtp<tht2) tht2=rtp;
            }
        }


        for(int i=0;i<3;++i){
            for(int j=0;j<NLBpt;++j){
                for(int a=0;a<4;++a){
                    ush[a]=n_approx_sol(cgslbx[j],cgsy[i],
                        xc,yc,dx,dy, dof[a],nDOF);
                }
                shbl2jbbl(ush,ujb);
                pre1=ujb[3];
                if(pre1>epsPos){
                    rtp=1.0;
                }
                else{
                    rtp=(pre0-epsPos)/(pre0-pre1);
                }

                if(rtp<tht2) tht2=rtp;
            }
        }

        for(int j=0;j<4;++j)
            for(int i=1;i<nDOF;++i)
                dof[j][i]*=tht2;
    }

    for(int i=0;i<4;i++) 
        for(int j=1; j< nDOF; j++) pc->dof[i][j]= dof[i][j];
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x1c8, %rsp            # imm = 0x1C8
movq %rdi, %rbx
leaq 0x70(%rsp), %r14
xorl %r15d, %r15d
movl $0xc0, %edx
movq %r14, %rdi
xorl %esi, %esi
callq 0x3120
xorps %xmm0, %xmm0
movaps %xmm0, 0x160(%rsp)
movaps %xmm0, 0x150(%rsp)
movaps %xmm0, 0x140(%rsp)
movaps %xmm0, 0x130(%rsp)
movaps %xmm0, 0x180(%rsp)
movsd 0x288(%rbx), %xmm4
movsd 0x290(%rbx), %xmm5
movsd 0x278(%rbx), %xmm2
movsd 0x280(%rbx), %xmm3
leaq 0x30(%rbx), %rax
movups (%rax,%r15), %xmm0
movups 0x10(%rax,%r15), %xmm1
movupd 0x20(%rax,%r15), %xmm6
movupd %xmm6, 0x90(%rsp,%r15)
movups %xmm1, 0x80(%rsp,%r15)
movups %xmm0, 0x70(%rsp,%r15)
addq $0x30, %r15
cmpq $0xc0, %r15
jne 0x18fee
movaps 0x4aa3(%rip), %xmm0     # 0x1dad0
movaps %xmm0, 0x170(%rsp)
movabsq $0x3ff0000000000000, %rax # imm = 0x3FF0000000000000
movq %rax, 0x180(%rsp)
movsd 0x40f9(%rip), %xmm0     # 0x1d148
movapd %xmm4, %xmm1
mulsd %xmm0, %xmm1
mulsd %xmm5, %xmm0
xorl %eax, %eax
leaq 0x4a8c(%rip), %rcx      # 0x1daf0
movsd (%rax,%rcx), %xmm6
movapd %xmm1, %xmm7
mulsd %xmm6, %xmm7
addsd %xmm2, %xmm7
movsd %xmm7, 0x1b0(%rsp,%rax)
mulsd %xmm0, %xmm6
addsd %xmm3, %xmm6
movsd %xmm6, 0x190(%rsp,%rax)
addq $0x8, %rax
cmpq $0x18, %rax
jne 0x19064
xorl %eax, %eax
movsd 0x170(%rsp,%rax,8), %xmm6
movapd %xmm1, %xmm7
mulsd %xmm6, %xmm7
addsd %xmm2, %xmm7
movsd %xmm7, 0x150(%rsp,%rax,8)
mulsd %xmm0, %xmm6
addsd %xmm3, %xmm6
movsd %xmm6, 0x130(%rsp,%rax,8)
incq %rax
cmpq $0x3, %rax
jne 0x1909b
movsd 0x70(%rsp), %xmm0
leaq 0xb138(%rip), %r12      # 0x24218
movsd (%r12), %xmm1
ucomisd %xmm0, %xmm1
movsd %xmm4, 0x28(%rsp)
movsd %xmm5, 0x20(%rsp)
movsd %xmm2, 0x18(%rsp)
movsd %xmm3, 0x10(%rsp)
jbe 0x1912c
movapd %xmm0, %xmm6
xorpd %xmm0, %xmm0
movupd %xmm0, 0x88(%rsp)
movupd %xmm0, 0x78(%rsp)
movq $0x0, 0x98(%rsp)
jmp 0x19256
movsd 0x498c(%rip), %xmm0     # 0x1dac0
xorl %r13d, %r13d
leaq 0x70(%rsp), %r15
movsd 0x190(%rsp,%r13,8), %xmm1
movsd %xmm1, (%rsp)
xorl %ebp, %ebp
movsd %xmm0, 0x8(%rsp)
movsd 0x150(%rsp,%rbp,8), %xmm0
movsd (%rsp), %xmm1
movsd 0x28(%rsp), %xmm4
movsd 0x20(%rsp), %xmm5
movq %r15, %rdi
movl $0x6, %esi
callq 0xd981
movsd 0x10(%rsp), %xmm3
movsd 0x18(%rsp), %xmm2
minsd 0x8(%rsp), %xmm0
incq %rbp
cmpq $0x3, %rbp
jne 0x1914d
incq %r13
cmpq $0x3, %r13
jne 0x1913c
xorl %r13d, %r13d
leaq 0x70(%rsp), %r15
movsd 0x1b0(%rsp,%r13,8), %xmm1
movsd %xmm1, (%rsp)
xorl %ebp, %ebp
movsd %xmm0, 0x8(%rsp)
movsd 0x130(%rsp,%rbp,8), %xmm1
movsd (%rsp), %xmm0
movsd 0x28(%rsp), %xmm4
movsd 0x20(%rsp), %xmm5
movq %r15, %rdi
movl $0x6, %esi
callq 0xd981
movsd 0x10(%rsp), %xmm3
movsd 0x18(%rsp), %xmm2
minsd 0x8(%rsp), %xmm0
incq %rbp
cmpq $0x3, %rbp
jne 0x191b7
incq %r13
cmpq $0x3, %r13
jne 0x191a6
movsd 0x70(%rsp), %xmm5
movsd (%r12), %xmm1
movapd %xmm5, %xmm4
subsd %xmm1, %xmm4
movapd %xmm5, %xmm6
subsd %xmm0, %xmm5
divsd %xmm5, %xmm4
andpd 0x3f90(%rip), %xmm4     # 0x1d1c0
minsd 0x3f68(%rip), %xmm4     # 0x1d1a0
movl $0x1, %eax
movsd 0x70(%rsp,%rax,8), %xmm0
mulsd %xmm4, %xmm0
movsd %xmm0, 0x70(%rsp,%rax,8)
incq %rax
cmpq $0x6, %rax
jne 0x1923d
movsd 0xa0(%rsp), %xmm0
movsd 0xd0(%rsp), %xmm4
divsd %xmm6, %xmm0
divsd %xmm6, %xmm4
mulsd %xmm4, %xmm4
mulsd %xmm0, %xmm0
mulsd 0x3ec0(%rip), %xmm6     # 0x1d140
addsd %xmm4, %xmm0
mulsd %xmm0, %xmm6
addsd 0x100(%rsp), %xmm6
mulsd 0x3f0f(%rip), %xmm6     # 0x1d1a8
movsd %xmm6, 0x38(%rsp)
ucomisd %xmm6, %xmm1
jbe 0x19315
xorpd %xmm0, %xmm0
movupd %xmm0, 0x88(%rsp)
movupd %xmm0, 0x78(%rsp)
xorl %eax, %eax
movq %rax, 0x98(%rsp)
movupd %xmm0, 0xa8(%rsp)
movupd %xmm0, 0xb8(%rsp)
movq %rax, 0xc8(%rsp)
movupd %xmm0, 0xd8(%rsp)
movupd %xmm0, 0xe8(%rsp)
movq %rax, 0xf8(%rsp)
movupd %xmm0, 0x108(%rsp)
movupd %xmm0, 0x118(%rsp)
movq %rax, 0x128(%rsp)
jmp 0x19584
movq %rbx, 0x68(%rsp)
movsd 0x479e(%rip), %xmm0     # 0x1dac0
xorl %r13d, %r13d
movsd 0x1b0(%rsp,%r13,8), %xmm6
xorl %ebp, %ebp
movapd %xmm0, %xmm1
movsd %xmm6, (%rsp)
movsd %xmm1, 0x30(%rsp)
movsd 0x130(%rsp,%rbp,8), %xmm0
movsd %xmm0, 0x8(%rsp)
movq %r14, %r15
xorl %ebx, %ebx
movapd %xmm6, %xmm0
movsd 0x8(%rsp), %xmm1
movsd 0x28(%rsp), %xmm4
movsd 0x20(%rsp), %xmm5
movq %r15, %rdi
movl $0x6, %esi
callq 0xd981
movsd (%rsp), %xmm6
movsd 0x10(%rsp), %xmm3
movsd 0x18(%rsp), %xmm2
movsd %xmm0, 0x40(%rsp,%rbx,8)
incq %rbx
addq $0x30, %r15
cmpq $0x4, %rbx
jne 0x19354
movsd 0x40(%rsp), %xmm1
movsd 0x48(%rsp), %xmm0
divsd %xmm1, %xmm0
movsd 0x50(%rsp), %xmm4
divsd %xmm1, %xmm4
mulsd %xmm4, %xmm4
mulsd %xmm0, %xmm0
addsd %xmm4, %xmm0
mulsd 0x3d77(%rip), %xmm1     # 0x1d140
mulsd %xmm0, %xmm1
addsd 0x58(%rsp), %xmm1
mulsd 0x3dcd(%rip), %xmm1     # 0x1d1a8
movsd (%r12), %xmm0
movsd 0x38(%rsp), %xmm5
movapd %xmm5, %xmm4
subsd %xmm0, %xmm4
subsd %xmm1, %xmm5
divsd %xmm5, %xmm4
cmpltsd %xmm1, %xmm0
movapd %xmm0, %xmm1
andnpd %xmm4, %xmm1
movsd 0x3d94(%rip), %xmm4     # 0x1d1a0
andpd %xmm4, %xmm0
orpd %xmm1, %xmm0
minsd 0x30(%rsp), %xmm0
incq %rbp
movapd %xmm0, %xmm1
cmpq $0x3, %rbp
jne 0x1933a
incq %r13
cmpq $0x3, %r13
jne 0x19325
xorl %r15d, %r15d
leaq 0x70(%rsp), %r13
movsd 0x190(%rsp,%r15,8), %xmm1
xorl %ebp, %ebp
movapd %xmm0, %xmm4
movsd %xmm1, (%rsp)
movsd %xmm4, 0x30(%rsp)
movsd 0x150(%rsp,%rbp,8), %xmm0
movsd %xmm0, 0x8(%rsp)
movq %r13, %r14
xorl %ebx, %ebx
movsd 0x8(%rsp), %xmm0
movsd 0x28(%rsp), %xmm4
movsd 0x20(%rsp), %xmm5
movq %r14, %rdi
movl $0x6, %esi
callq 0xd981
movsd (%rsp), %xmm1
movsd 0x10(%rsp), %xmm3
movsd 0x18(%rsp), %xmm2
movsd %xmm0, 0x40(%rsp,%rbx,8)
incq %rbx
addq $0x30, %r14
cmpq $0x4, %rbx
jne 0x1946f
movsd 0x40(%rsp), %xmm6
movsd 0x48(%rsp), %xmm0
divsd %xmm6, %xmm0
movsd 0x50(%rsp), %xmm4
divsd %xmm6, %xmm4
mulsd %xmm4, %xmm4
mulsd %xmm0, %xmm0
addsd %xmm4, %xmm0
mulsd 0x3c60(%rip), %xmm6     # 0x1d140
mulsd %xmm0, %xmm6
addsd 0x58(%rsp), %xmm6
mulsd 0x3cb6(%rip), %xmm6     # 0x1d1a8
movsd (%r12), %xmm0
movsd 0x38(%rsp), %xmm5
movapd %xmm5, %xmm4
subsd %xmm0, %xmm4
subsd %xmm6, %xmm5
divsd %xmm5, %xmm4
cmpltsd %xmm6, %xmm0
movapd %xmm0, %xmm5
andnpd %xmm4, %xmm5
movsd 0x3c7d(%rip), %xmm4     # 0x1d1a0
andpd %xmm4, %xmm0
orpd %xmm5, %xmm0
minsd 0x30(%rsp), %xmm0
incq %rbp
movapd %xmm0, %xmm4
cmpq $0x3, %rbp
jne 0x19455
incq %r15
cmpq $0x3, %r15
jne 0x19440
leaq 0x70(%rsp), %rax
xorl %ecx, %ecx
movq 0x68(%rsp), %rbx
movl $0x1, %edx
movsd (%rax,%rdx,8), %xmm1
mulsd %xmm0, %xmm1
movsd %xmm1, (%rax,%rdx,8)
incq %rdx
cmpq $0x6, %rdx
jne 0x19560
incq %rcx
addq $0x30, %rax
cmpq $0x4, %rcx
jne 0x1955b
addq $0x38, %rbx
xorl %eax, %eax
movq 0x98(%rsp,%rax), %rcx
movq %rcx, 0x20(%rbx,%rax)
movupd 0x78(%rsp,%rax), %xmm0
movupd 0x88(%rsp,%rax), %xmm1
movupd %xmm1, 0x10(%rbx,%rax)
movupd %xmm0, (%rbx,%rax)
addq $0x30, %rax
cmpq $0xc0, %rax
jne 0x1958a
addq $0x1c8, %rsp            # imm = 0x1C8
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq