// File: /workspace/github_cpp_makefile_star2_O1/DoubleString[P]RtConverter/src/Rtklib/rtkcmn.cpp
// Function: bamboo::ionmodel(bamboo::gtime_t, double const*, double const*, double const*)
extern double ionmodel(gtime_t t, const double *ion, const double *pos,
                       const double *azel)
{
    const double ion_default[]={ /* 2004/1/1 */
        0.1118E-07,-0.7451E-08,-0.5961E-07, 0.1192E-06,
        0.1167E+06,-0.2294E+06,-0.1311E+06, 0.1049E+07
    };
    double tt,f,psi,phi,lam,amp,per,x;
    int week;
    
    if (pos[2]<-1E3||azel[1]<=0) return 0.0;
    if (norm(ion,8)<=0.0) ion=ion_default;
    
    /* earth centered angle (semi-circle) */
    psi=0.0137/(azel[1]/pi+0.11)-0.022;
    
    /* subionospheric latitude/longitude (semi-circle) */
    phi=pos[0]/pi+psi*cos(azel[0]);
    if      (phi> 0.416) phi= 0.416;
    else if (phi<-0.416) phi=-0.416;
    lam=pos[1]/pi+psi*sin(azel[0])/cos(phi*pi);
    
    /* geomagnetic latitude (semi-circle) */
    phi+=0.064*cos((lam-1.617)*pi);
    
    /* local time (s) */
    tt=43200.0*lam+time2gpst(t,&week);
    tt-=floor(tt/86400.0)*86400.0; /* 0<=tt<86400 */
    
    /* slant factor */
    f=1.0+16.0*pow(0.53-azel[1]/pi,3.0);
    
    /* ionospheric delay */
    amp=ion[0]+phi*(ion[1]+phi*(ion[2]+phi*ion[3]));
    per=ion[4]+phi*(ion[5]+phi*(ion[6]+phi*ion[7]));
    amp=amp<    0.0?    0.0:amp;
    per=per<72000.0?72000.0:per;
    x=2.0*pi*(tt-50400.0)/per;
    
    return CLIGHT*f*(fabs(x)<1.57?5E-9+amp*(1.0+x*x*(-0.5+x*x/24.0)):5E-9);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movapd %xmm0, %xmm2
xorpd %xmm0, %xmm0
movsd 0x1fca2(%rip), %xmm1    # 0x74ce8
ucomisd 0x10(%rdx), %xmm1
ja 0x55345
movq %rcx, %rbx
ucomisd 0x8(%rcx), %xmm0
jae 0x55345
movq %rdx, %r15
movq %rsi, %r12
movq %rdi, %r14
movl $0x9, %eax
movsd -0x10(%r12,%rax,8), %xmm1
mulsd %xmm1, %xmm1
addsd %xmm1, %xmm0
decq %rax
cmpq $0x1, %rax
ja 0x5506d
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
movsd %xmm2, -0x50(%rbp)
jb 0x5509a
sqrtsd %xmm0, %xmm0
jmp 0x550a3
callq 0x84e0
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
leaq 0x1d6a2(%rip), %r13     # 0x72750
cmovbq %r12, %r13
movsd (%rbx), %xmm0
movsd 0x8(%rbx), %xmm1
movsd 0x1584d(%rip), %xmm2    # 0x6a910
divsd %xmm2, %xmm1
addsd 0x1fc21(%rip), %xmm1    # 0x74cf0
movsd 0x1fc21(%rip), %xmm3    # 0x74cf8
divsd %xmm1, %xmm3
addsd 0x1fc1d(%rip), %xmm3    # 0x74d00
movsd %xmm3, -0x40(%rbp)
movsd (%r15), %xmm1
divsd %xmm2, %xmm1
movsd %xmm1, -0x30(%rbp)
callq 0x8c00
mulsd -0x40(%rbp), %xmm0
addsd -0x30(%rbp), %xmm0
movsd 0x1fbfb(%rip), %xmm2    # 0x74d08
ucomisd %xmm2, %xmm0
ja 0x55129
movsd 0x1fbf5(%rip), %xmm1    # 0x74d10
ucomisd %xmm0, %xmm1
movapd %xmm0, %xmm2
jbe 0x55129
movapd %xmm1, %xmm2
movsd %xmm2, -0x30(%rbp)
movsd 0x8(%r15), %xmm1
movsd 0x157d4(%rip), %xmm0    # 0x6a910
divsd %xmm0, %xmm1
movsd %xmm1, -0x48(%rbp)
movsd (%rbx), %xmm0
callq 0x8150
movsd -0x40(%rbp), %xmm1
mulsd %xmm0, %xmm1
movsd %xmm1, -0x40(%rbp)
movsd -0x30(%rbp), %xmm0
mulsd 0x157a7(%rip), %xmm0    # 0x6a910
callq 0x8c00
movsd -0x40(%rbp), %xmm1
divsd %xmm0, %xmm1
addsd -0x48(%rbp), %xmm1
movsd %xmm1, -0x40(%rbp)
movsd 0x1fb8f(%rip), %xmm0    # 0x74d18
addsd %xmm1, %xmm0
mulsd 0x1577b(%rip), %xmm0    # 0x6a910
callq 0x8c00
mulsd 0x16bee(%rip), %xmm0    # 0x6bd90
addsd -0x30(%rbp), %xmm0
movsd %xmm0, -0x30(%rbp)
leaq 0x1d50d(%rip), %rdi     # 0x726c0
callq 0x4f9a4
subq %rax, %r14
movabsq $0x377aef2669de1559, %rcx # imm = 0x377AEF2669DE1559
movq %r14, %rax
imulq %rcx
movq %rdx, %rax
shrq $0x3f, %rax
shrq $0x11, %rdx
addl %eax, %edx
xorps %xmm0, %xmm0
cvtsi2sd %r14, %xmm0
xorps %xmm1, %xmm1
cvtsi2sd %edx, %xmm1
mulsd 0x16c11(%rip), %xmm1    # 0x6be00
mulsd 0x156e9(%rip), %xmm1    # 0x6a8e0
addsd %xmm0, %xmm1
addsd -0x50(%rbp), %xmm1
movsd -0x40(%rbp), %xmm0
mulsd 0x16beb(%rip), %xmm0    # 0x6bdf8
addsd %xmm1, %xmm0
movsd %xmm0, -0x40(%rbp)
divsd 0x156b2(%rip), %xmm0    # 0x6a8d0
callq 0x8ad0
mulsd 0x16bd5(%rip), %xmm0    # 0x6be00
addsd -0x40(%rbp), %xmm0
movapd %xmm0, -0x40(%rbp)
movsd 0x8(%rbx), %xmm1
divsd 0x156ce(%rip), %xmm1    # 0x6a910
movsd 0x1fad6(%rip), %xmm0    # 0x74d20
subsd %xmm1, %xmm0
movsd 0x15692(%rip), %xmm1    # 0x6a8e8
callq 0x8b30
movapd -0x40(%rbp), %xmm4
movsd -0x30(%rbp), %xmm3
mulsd 0x17b0b(%rip), %xmm0    # 0x6cd78
addsd 0x146ab(%rip), %xmm0    # 0x69920
movsd 0x38(%r13), %xmm1
mulsd %xmm3, %xmm1
addsd 0x30(%r13), %xmm1
mulsd %xmm3, %xmm1
addsd 0x28(%r13), %xmm1
mulsd %xmm3, %xmm1
addsd 0x20(%r13), %xmm1
movsd 0x1fa87(%rip), %xmm2    # 0x74d28
maxsd %xmm1, %xmm2
addsd 0x1fa83(%rip), %xmm4    # 0x74d30
mulsd 0x15663(%rip), %xmm4    # 0x6a918
divsd %xmm2, %xmm4
mulsd 0x16adf(%rip), %xmm0    # 0x6bda0
movapd 0x14667(%rip), %xmm1    # 0x69930
andpd %xmm4, %xmm1
movsd 0x1fa6b(%rip), %xmm2    # 0x74d40
ucomisd %xmm1, %xmm2
jbe 0x55339
movsd 0x18(%r13), %xmm2
mulsd %xmm3, %xmm2
addsd 0x10(%r13), %xmm2
mulsd %xmm3, %xmm2
addsd 0x8(%r13), %xmm2
mulsd %xmm3, %xmm2
addsd (%r13), %xmm2
mulsd %xmm4, %xmm4
movapd %xmm4, %xmm1
divsd 0x1fa39(%rip), %xmm1    # 0x74d48
addsd 0x156d9(%rip), %xmm1    # 0x6a9f0
mulsd %xmm4, %xmm1
xorpd %xmm3, %xmm3
addsd 0x145f9(%rip), %xmm1    # 0x69920
maxsd %xmm2, %xmm3
mulsd %xmm3, %xmm1
addsd 0x1fa01(%rip), %xmm1    # 0x74d38
jmp 0x55341
movsd 0x1f9f7(%rip), %xmm1    # 0x74d38
mulsd %xmm1, %xmm0
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq