// File: /workspace/github_cpp_makefile_star2_O1/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 0xd0(%rsp), %r14
xorl %r15d, %r15d
movl $0xc0, %edx
movq %r14, %rdi
xorl %esi, %esi
callq 0x3120
xorps %xmm0, %xmm0
movaps %xmm0, 0xc0(%rsp)
movaps %xmm0, 0xb0(%rsp)
movaps %xmm0, 0xa0(%rsp)
movaps %xmm0, 0x90(%rsp)
movaps %xmm0, 0x80(%rsp)
movaps %xmm0, 0x70(%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, 0xf0(%rsp,%r15)
movups %xmm1, 0xe0(%rsp,%r15)
movups %xmm0, 0xd0(%rsp,%r15)
addq $0x30, %r15
cmpq $0xc0, %r15
jne 0x194f2
movaps 0x559c(%rip), %xmm0     # 0x1ead0
movaps %xmm0, 0x70(%rsp)
movabsq $0x3ff0000000000000, %rax # imm = 0x3FF0000000000000
movq %rax, 0x80(%rsp)
movsd 0x4bf5(%rip), %xmm0     # 0x1e148
movapd %xmm4, %xmm1
mulsd %xmm0, %xmm1
mulsd %xmm5, %xmm0
xorl %eax, %eax
leaq 0x5588(%rip), %rcx      # 0x1eaf0
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 0x19568
xorl %eax, %eax
movsd 0x70(%rsp,%rax,8), %xmm6
movapd %xmm1, %xmm7
mulsd %xmm6, %xmm7
addsd %xmm2, %xmm7
movsd %xmm7, 0xb0(%rsp,%rax,8)
mulsd %xmm0, %xmm6
addsd %xmm3, %xmm6
movsd %xmm6, 0x90(%rsp,%rax,8)
incq %rax
cmpq $0x3, %rax
jne 0x1959f
leaq 0xbc3d(%rip), %r12      # 0x25218
movsd (%r12), %xmm0
ucomisd 0xd0(%rsp), %xmm0
movsd %xmm4, 0x28(%rsp)
movsd %xmm5, 0x20(%rsp)
movsd %xmm2, 0x18(%rsp)
movsd %xmm3, 0x10(%rsp)
jbe 0x1962b
xorpd %xmm0, %xmm0
movupd %xmm0, 0xe8(%rsp)
movupd %xmm0, 0xd8(%rsp)
movq $0x0, 0xf8(%rsp)
jmp 0x1975c
movsd 0x548d(%rip), %xmm0     # 0x1eac0
xorl %r13d, %r13d
leaq 0xd0(%rsp), %r15
movsd 0x190(%rsp,%r13,8), %xmm1
movsd %xmm1, (%rsp)
xorl %ebp, %ebp
movsd %xmm0, 0x8(%rsp)
movsd 0xb0(%rsp,%rbp,8), %xmm0
movsd (%rsp), %xmm1
movsd 0x28(%rsp), %xmm4
movsd 0x20(%rsp), %xmm5
movq %r15, %rdi
movl $0x6, %esi
callq 0xdcc9
movsd 0x10(%rsp), %xmm3
movsd 0x18(%rsp), %xmm2
minsd 0x8(%rsp), %xmm0
incq %rbp
cmpq $0x3, %rbp
jne 0x1964f
incq %r13
cmpq $0x3, %r13
jne 0x1963e
xorl %r13d, %r13d
leaq 0xd0(%rsp), %r15
movsd 0x1b0(%rsp,%r13,8), %xmm1
movsd %xmm1, (%rsp)
xorl %ebp, %ebp
movsd %xmm0, 0x8(%rsp)
movsd 0x90(%rsp,%rbp,8), %xmm1
movsd (%rsp), %xmm0
movsd 0x28(%rsp), %xmm4
movsd 0x20(%rsp), %xmm5
movq %r15, %rdi
movl $0x6, %esi
callq 0xdcc9
movsd 0x10(%rsp), %xmm3
movsd 0x18(%rsp), %xmm2
minsd 0x8(%rsp), %xmm0
incq %rbp
cmpq $0x3, %rbp
jne 0x196bc
incq %r13
cmpq $0x3, %r13
jne 0x196ab
movsd 0xd0(%rsp), %xmm4
movapd %xmm4, %xmm1
subsd (%r12), %xmm1
subsd %xmm0, %xmm4
divsd %xmm4, %xmm1
andpd 0x4a90(%rip), %xmm1     # 0x1e1c0
minsd 0x4a68(%rip), %xmm1     # 0x1e1a0
movl $0x1, %eax
movsd 0xd0(%rsp,%rax,8), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, 0xd0(%rsp,%rax,8)
incq %rax
cmpq $0x6, %rax
jne 0x1973d
movsd 0xd0(%rsp), %xmm4
movsd 0x100(%rsp), %xmm0
divsd %xmm4, %xmm0
movsd 0x130(%rsp), %xmm1
divsd %xmm4, %xmm1
mulsd %xmm1, %xmm1
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
mulsd 0x49ad(%rip), %xmm4     # 0x1e140
mulsd %xmm0, %xmm4
addsd 0x160(%rsp), %xmm4
mulsd 0x4a00(%rip), %xmm4     # 0x1e1a8
movsd (%r12), %xmm0
movsd %xmm4, 0x38(%rsp)
ucomisd %xmm4, %xmm0
jbe 0x197f2
movl $0x8, %eax
xorpd %xmm0, %xmm0
movupd %xmm0, 0xe0(%rsp,%rax)
movupd %xmm0, 0xd0(%rsp,%rax)
movq $0x0, 0xf0(%rsp,%rax)
addq $0x30, %rax
cmpq $0xc8, %rax
jne 0x197c3
jmp 0x19a67
movq %rbx, 0x68(%rsp)
movsd 0x52c1(%rip), %xmm0     # 0x1eac0
xorl %r13d, %r13d
movsd 0x1b0(%rsp,%r13,8), %xmm6
xorl %ebp, %ebp
movapd %xmm0, %xmm1
movsd %xmm6, (%rsp)
movsd %xmm1, 0x30(%rsp)
movsd 0x90(%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 0xdcc9
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 0x19831
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 0x489a(%rip), %xmm1     # 0x1e140
mulsd %xmm0, %xmm1
addsd 0x58(%rsp), %xmm1
mulsd 0x48f0(%rip), %xmm1     # 0x1e1a8
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 0x48b7(%rip), %xmm4     # 0x1e1a0
andpd %xmm4, %xmm0
orpd %xmm1, %xmm0
minsd 0x30(%rsp), %xmm0
incq %rbp
movapd %xmm0, %xmm1
cmpq $0x3, %rbp
jne 0x19817
incq %r13
cmpq $0x3, %r13
jne 0x19802
xorl %r15d, %r15d
leaq 0xd0(%rsp), %r13
movsd 0x190(%rsp,%r15,8), %xmm1
xorl %ebp, %ebp
movapd %xmm0, %xmm4
movsd %xmm1, (%rsp)
movsd %xmm4, 0x30(%rsp)
movsd 0xb0(%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 0xdcc9
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 0x1994f
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 0x4780(%rip), %xmm6     # 0x1e140
mulsd %xmm0, %xmm6
addsd 0x58(%rsp), %xmm6
mulsd 0x47d6(%rip), %xmm6     # 0x1e1a8
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 0x479d(%rip), %xmm4     # 0x1e1a0
andpd %xmm4, %xmm0
orpd %xmm5, %xmm0
minsd 0x30(%rsp), %xmm0
incq %rbp
movapd %xmm0, %xmm4
cmpq $0x3, %rbp
jne 0x19935
incq %r15
cmpq $0x3, %r15
jne 0x19920
leaq 0xd0(%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 0x19a43
incq %rcx
addq $0x30, %rax
cmpq $0x4, %rcx
jne 0x19a3e
addq $0x38, %rbx
xorl %eax, %eax
movq 0xf8(%rsp,%rax), %rcx
movq %rcx, 0x20(%rbx,%rax)
movupd 0xd8(%rsp,%rax), %xmm0
movupd 0xe8(%rsp,%rax), %xmm1
movupd %xmm1, 0x10(%rbx,%rax)
movupd %xmm0, (%rbx,%rax)
addq $0x30, %rax
cmpq $0xc0, %rax
jne 0x19a6d
addq $0x1c8, %rsp            # imm = 0x1C8
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq