// File: /workspace/github_cpp_makefile_star2_O2/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 $0x38, %rsp
movapd %xmm0, %xmm2
xorpd %xmm0, %xmm0
movsd 0x1f057(%rip), %xmm1    # 0x6acc8
ucomisd 0x10(%rdx), %xmm1
ja 0x4bf0e
movq %rcx, %rbx
ucomisd 0x8(%rcx), %xmm0
jae 0x4bf0e
movq %rdx, %r15
movq %rsi, %r12
movq %rdi, %r14
movsd %xmm2, -0x50(%rbp)
pushq $0x8
popq %rsi
movq %r12, %rdi
callq 0x45e31
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
leaq 0x1ca7e(%rip), %r13     # 0x68730
cmovbq %r12, %r13
movsd (%rbx), %xmm0
movsd 0x8(%rbx), %xmm1
movsd 0x14c39(%rip), %xmm2    # 0x60900
divsd %xmm2, %xmm1
addsd 0x1effd(%rip), %xmm1    # 0x6acd0
movsd 0x1effd(%rip), %xmm3    # 0x6acd8
divsd %xmm1, %xmm3
addsd 0x1eff9(%rip), %xmm3    # 0x6ace0
movapd %xmm3, -0x40(%rbp)
movsd (%r15), %xmm1
divsd %xmm2, %xmm1
movsd %xmm1, -0x30(%rbp)
callq 0x8bf0
mulsd -0x40(%rbp), %xmm0
addsd -0x30(%rbp), %xmm0
movsd 0x1efd7(%rip), %xmm2    # 0x6ace8
ucomisd %xmm2, %xmm0
ja 0x4bd2d
movsd 0x1efd1(%rip), %xmm1    # 0x6acf0
ucomisd %xmm0, %xmm1
movapd %xmm0, %xmm2
jbe 0x4bd2d
movapd %xmm1, %xmm2
movsd %xmm2, -0x30(%rbp)
movsd 0x8(%r15), %xmm1
movsd 0x14bc0(%rip), %xmm0    # 0x60900
divsd %xmm0, %xmm1
movsd %xmm1, -0x48(%rbp)
movsd (%rbx), %xmm0
callq 0x8160
movapd -0x40(%rbp), %xmm1
mulsd %xmm0, %xmm1
movapd %xmm1, -0x40(%rbp)
movsd -0x30(%rbp), %xmm0
mulsd 0x14b93(%rip), %xmm0    # 0x60900
callq 0x8bf0
movapd -0x40(%rbp), %xmm1
divsd %xmm0, %xmm1
addsd -0x48(%rbp), %xmm1
movapd %xmm1, -0x40(%rbp)
movsd 0x1ef6b(%rip), %xmm0    # 0x6acf8
addsd %xmm1, %xmm0
mulsd 0x14b67(%rip), %xmm0    # 0x60900
callq 0x8bf0
mulsd 0x15fea(%rip), %xmm0    # 0x61d90
addsd -0x30(%rbp), %xmm0
movsd %xmm0, -0x30(%rbp)
leaq -0x54(%rbp), %rsi
movq %r14, %rdi
movsd -0x50(%rbp), %xmm0
callq 0x470f1
movapd -0x40(%rbp), %xmm1
mulsd 0x1601a(%rip), %xmm1    # 0x61de8
addsd %xmm0, %xmm1
movapd %xmm1, -0x40(%rbp)
movapd %xmm1, %xmm0
divsd 0x14add(%rip), %xmm0    # 0x608c0
callq 0x8ac0
mulsd 0x14ad0(%rip), %xmm0    # 0x608c0
movapd -0x40(%rbp), %xmm1
subsd %xmm0, %xmm1
movapd %xmm1, -0x40(%rbp)
movsd 0x8(%rbx), %xmm1
divsd 0x14af5(%rip), %xmm1    # 0x60900
movsd 0x1eeed(%rip), %xmm0    # 0x6ad00
subsd %xmm1, %xmm0
movsd 0x14ab9(%rip), %xmm1    # 0x608d8
callq 0x8b20
movsd -0x30(%rbp), %xmm4
movapd -0x40(%rbp), %xmm3
mulsd 0x16f42(%rip), %xmm0    # 0x62d78
addsd 0x13ae2(%rip), %xmm0    # 0x5f920
movsd 0x38(%r13), %xmm1
mulsd %xmm4, %xmm1
addsd 0x30(%r13), %xmm1
mulsd %xmm4, %xmm1
addsd 0x28(%r13), %xmm1
mulsd %xmm4, %xmm1
addsd 0x20(%r13), %xmm1
movsd 0x1ee9e(%rip), %xmm2    # 0x6ad08
maxsd %xmm1, %xmm2
addsd 0x1ee9a(%rip), %xmm3    # 0x6ad10
mulsd 0x14a8a(%rip), %xmm3    # 0x60908
divsd %xmm2, %xmm3
mulsd 0x15f16(%rip), %xmm0    # 0x61da0
movapd 0x13a9e(%rip), %xmm1    # 0x5f930
andpd %xmm3, %xmm1
movsd 0x1ee82(%rip), %xmm2    # 0x6ad20
ucomisd %xmm1, %xmm2
jbe 0x4bf02
movsd 0x18(%r13), %xmm2
mulsd %xmm4, %xmm2
addsd 0x10(%r13), %xmm2
mulsd %xmm4, %xmm2
addsd 0x8(%r13), %xmm2
mulsd %xmm4, %xmm2
addsd (%r13), %xmm2
mulsd %xmm3, %xmm3
movapd %xmm3, %xmm1
divsd 0x1ee50(%rip), %xmm1    # 0x6ad28
addsd 0x14b00(%rip), %xmm1    # 0x609e0
mulsd %xmm3, %xmm1
xorpd %xmm3, %xmm3
addsd 0x13a30(%rip), %xmm1    # 0x5f920
maxsd %xmm2, %xmm3
mulsd %xmm3, %xmm1
addsd 0x1ee18(%rip), %xmm1    # 0x6ad18
jmp 0x4bf0a
movsd 0x1ee0e(%rip), %xmm1    # 0x6ad18
mulsd %xmm1, %xmm0
addq $0x38, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq