// File: /workspace/github_cpp_makefile_star2_O3/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 0x20726(%rip), %xmm1    # 0x74cc8
ucomisd 0x10(%rdx), %xmm1
ja 0x548a1
movq %rcx, %rbx
ucomisd 0x8(%rcx), %xmm0
jae 0x548a1
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 0x545c9
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
movsd %xmm2, -0x50(%rbp)
jb 0x545f6
sqrtsd %xmm0, %xmm0
jmp 0x545ff
callq 0x84e0
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
leaq 0x1e126(%rip), %r13     # 0x72730
cmovbq %r12, %r13
movsd (%rbx), %xmm0
movsd 0x8(%rbx), %xmm1
movsd 0x162f1(%rip), %xmm2    # 0x6a910
divsd %xmm2, %xmm1
addsd 0x206a5(%rip), %xmm1    # 0x74cd0
movsd 0x206a5(%rip), %xmm3    # 0x74cd8
divsd %xmm1, %xmm3
addsd 0x206a1(%rip), %xmm3    # 0x74ce0
movsd %xmm3, -0x40(%rbp)
movsd (%r15), %xmm1
divsd %xmm2, %xmm1
movsd %xmm1, -0x30(%rbp)
callq 0x8bf0
mulsd -0x40(%rbp), %xmm0
addsd -0x30(%rbp), %xmm0
movsd 0x2067f(%rip), %xmm2    # 0x74ce8
ucomisd %xmm2, %xmm0
ja 0x54685
movsd 0x20679(%rip), %xmm1    # 0x74cf0
ucomisd %xmm0, %xmm1
movapd %xmm0, %xmm2
jbe 0x54685
movapd %xmm1, %xmm2
movsd %xmm2, -0x30(%rbp)
movsd 0x8(%r15), %xmm1
movsd 0x16278(%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 0x1624b(%rip), %xmm0    # 0x6a910
callq 0x8bf0
movsd -0x40(%rbp), %xmm1
divsd %xmm0, %xmm1
addsd -0x48(%rbp), %xmm1
movsd %xmm1, -0x40(%rbp)
movsd 0x20613(%rip), %xmm0    # 0x74cf8
addsd %xmm1, %xmm0
mulsd 0x1621f(%rip), %xmm0    # 0x6a910
callq 0x8bf0
mulsd 0x17692(%rip), %xmm0    # 0x6bd90
addsd -0x30(%rbp), %xmm0
movsd %xmm0, -0x30(%rbp)
leaq 0x1df91(%rip), %rdi     # 0x726a0
callq 0x4f027
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 0x176b5(%rip), %xmm1    # 0x6be00
mulsd 0x1618d(%rip), %xmm1    # 0x6a8e0
addsd %xmm0, %xmm1
addsd -0x50(%rbp), %xmm1
movsd -0x40(%rbp), %xmm0
mulsd 0x1768f(%rip), %xmm0    # 0x6bdf8
addsd %xmm1, %xmm0
movsd %xmm0, -0x40(%rbp)
divsd 0x16156(%rip), %xmm0    # 0x6a8d0
callq 0x8ac0
mulsd 0x17679(%rip), %xmm0    # 0x6be00
addsd -0x40(%rbp), %xmm0
movapd %xmm0, -0x40(%rbp)
movsd 0x8(%rbx), %xmm1
divsd 0x16172(%rip), %xmm1    # 0x6a910
movsd 0x2055a(%rip), %xmm0    # 0x74d00
subsd %xmm1, %xmm0
movsd 0x16136(%rip), %xmm1    # 0x6a8e8
callq 0x8b20
movapd -0x40(%rbp), %xmm4
movsd -0x30(%rbp), %xmm3
mulsd 0x185af(%rip), %xmm0    # 0x6cd78
addsd 0x1514f(%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 0x2050b(%rip), %xmm2    # 0x74d08
maxsd %xmm1, %xmm2
addsd 0x20507(%rip), %xmm4    # 0x74d10
mulsd 0x16107(%rip), %xmm4    # 0x6a918
divsd %xmm2, %xmm4
mulsd 0x17583(%rip), %xmm0    # 0x6bda0
movapd 0x1510b(%rip), %xmm1    # 0x69930
andpd %xmm4, %xmm1
movsd 0x204ef(%rip), %xmm2    # 0x74d20
ucomisd %xmm1, %xmm2
jbe 0x54895
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 0x204bd(%rip), %xmm1    # 0x74d28
addsd 0x1617d(%rip), %xmm1    # 0x6a9f0
mulsd %xmm4, %xmm1
xorpd %xmm3, %xmm3
addsd 0x1509d(%rip), %xmm1    # 0x69920
maxsd %xmm2, %xmm3
mulsd %xmm3, %xmm1
addsd 0x20485(%rip), %xmm1    # 0x74d18
jmp 0x5489d
movsd 0x2047b(%rip), %xmm1    # 0x74d18
mulsd %xmm1, %xmm0
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq