// File: 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
subq $0x120, %rsp            # imm = 0x120
movq %rdi, -0x18(%rbp)
movsd %xmm0, -0x10(%rbp)
movq %rsi, -0x20(%rbp)
movq %rdx, -0x28(%rbp)
movq %rcx, -0x30(%rbp)
leaq -0x70(%rbp), %rdi
leaq 0x2d8f5(%rip), %rsi     # 0xc4a90
movl $0x40, %edx
callq 0x9610
movq -0x28(%rbp), %rax
movsd 0x3038f(%rip), %xmm0    # 0xc7540
ucomisd 0x10(%rax), %xmm0
ja 0x971c6
movq -0x30(%rbp), %rax
xorps %xmm0, %xmm0
ucomisd 0x8(%rax), %xmm0
jb 0x971d3
xorps %xmm0, %xmm0
movsd %xmm0, -0x8(%rbp)
jmp 0x976cc
movq -0x20(%rbp), %rdi
movl $0x8, %esi
callq 0x8c820
movaps %xmm0, %xmm1
xorps %xmm0, %xmm0
ucomisd %xmm1, %xmm0
jb 0x971f5
leaq -0x70(%rbp), %rax
movq %rax, -0x20(%rbp)
movq -0x30(%rbp), %rax
movsd 0x8(%rax), %xmm1
movsd 0x25c52(%rip), %xmm0    # 0xbce58
divsd %xmm0, %xmm1
movsd 0x3034e(%rip), %xmm0    # 0xc7560
addsd %xmm0, %xmm1
movsd 0x3033a(%rip), %xmm0    # 0xc7558
divsd %xmm1, %xmm0
movsd 0x30326(%rip), %xmm1    # 0xc7550
subsd %xmm1, %xmm0
movsd %xmm0, -0x88(%rbp)
movq -0x28(%rbp), %rax
movsd (%rax), %xmm0
movsd 0x25c12(%rip), %xmm1    # 0xbce58
divsd %xmm1, %xmm0
movsd %xmm0, -0xd8(%rbp)
movsd -0x88(%rbp), %xmm0
movsd %xmm0, -0xd0(%rbp)
movq -0x30(%rbp), %rax
movsd (%rax), %xmm0
callq 0x9e00
movsd -0xd8(%rbp), %xmm1
movaps %xmm0, %xmm2
movsd -0xd0(%rbp), %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x90(%rbp)
movsd -0x90(%rbp), %xmm0
movsd 0x302a6(%rip), %xmm1    # 0xc7548
ucomisd %xmm1, %xmm0
jbe 0x972ba
movsd 0x30298(%rip), %xmm0    # 0xc7548
movsd %xmm0, -0x90(%rbp)
jmp 0x972de
movsd 0x302a6(%rip), %xmm0    # 0xc7568
ucomisd -0x90(%rbp), %xmm0
jbe 0x972dc
movsd 0x30294(%rip), %xmm0    # 0xc7568
movsd %xmm0, -0x90(%rbp)
jmp 0x972de
movq -0x28(%rbp), %rax
movsd 0x8(%rax), %xmm0
movsd 0x25b69(%rip), %xmm1    # 0xbce58
divsd %xmm1, %xmm0
movsd %xmm0, -0xe8(%rbp)
movsd -0x88(%rbp), %xmm0
movsd %xmm0, -0xf8(%rbp)
movq -0x30(%rbp), %rax
movsd (%rax), %xmm0
callq 0x91a0
movaps %xmm0, %xmm1
movsd -0xf8(%rbp), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0xf0(%rbp)
movsd 0x25b21(%rip), %xmm0    # 0xbce58
mulsd -0x90(%rbp), %xmm0
callq 0x9e00
movsd -0xf0(%rbp), %xmm1
movaps %xmm0, %xmm2
movsd -0xe8(%rbp), %xmm0
divsd %xmm2, %xmm1
addsd %xmm1, %xmm0
movsd %xmm0, -0x98(%rbp)
movsd -0x98(%rbp), %xmm0
movsd 0x30201(%rip), %xmm1    # 0xc7578
subsd %xmm1, %xmm0
movsd 0x25ad5(%rip), %xmm1    # 0xbce58
mulsd %xmm1, %xmm0
callq 0x9e00
movsd -0x90(%rbp), %xmm1
movsd 0x26a24(%rip), %xmm2    # 0xbddc0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x90(%rbp)
movsd -0x98(%rbp), %xmm0
movsd %xmm0, -0xe0(%rbp)
movq -0x18(%rbp), %rax
movq %rax, -0xc8(%rbp)
movq -0x10(%rbp), %rax
movq %rax, -0xc0(%rbp)
movq -0xc8(%rbp), %rdi
movsd -0xc0(%rbp), %xmm0
leaq -0xb4(%rbp), %rsi
callq 0x8e540
movaps %xmm0, %xmm1
movsd -0xe0(%rbp), %xmm0
movsd 0x26970(%rip), %xmm2    # 0xbdd70
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x78(%rbp)
movsd -0x78(%rbp), %xmm0
movsd 0x25a06(%rip), %xmm1    # 0xbce20
divsd %xmm1, %xmm0
callq 0x9c90
movsd -0x78(%rbp), %xmm1
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
movsd 0x259d9(%rip), %xmm2    # 0xbce20
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x78(%rbp)
movq -0x30(%rbp), %rax
movsd 0x8(%rax), %xmm1
movsd 0x259f3(%rip), %xmm0    # 0xbce58
divsd %xmm0, %xmm1
movsd 0x300ff(%rip), %xmm0    # 0xc7570
subsd %xmm1, %xmm0
movsd 0x259bb(%rip), %xmm1    # 0xbce38
callq 0x9d10
movsd 0x27a36(%rip), %xmm1    # 0xbeec0
mulsd %xmm1, %xmm0
movsd 0x249ca(%rip), %xmm1    # 0xbbe60
addsd %xmm1, %xmm0
movsd %xmm0, -0x80(%rbp)
movq -0x20(%rbp), %rax
movsd (%rax), %xmm1
movsd -0x90(%rbp), %xmm0
movq -0x20(%rbp), %rax
movsd 0x8(%rax), %xmm3
movsd -0x90(%rbp), %xmm2
movq -0x20(%rbp), %rax
movsd 0x10(%rax), %xmm5
movsd -0x90(%rbp), %xmm4
movq -0x20(%rbp), %rax
movsd 0x18(%rax), %xmm6
mulsd %xmm6, %xmm4
addsd %xmm5, %xmm4
mulsd %xmm4, %xmm2
addsd %xmm3, %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0xa0(%rbp)
movq -0x20(%rbp), %rax
movsd 0x20(%rax), %xmm1
movsd -0x90(%rbp), %xmm0
movq -0x20(%rbp), %rax
movsd 0x28(%rax), %xmm3
movsd -0x90(%rbp), %xmm2
movq -0x20(%rbp), %rax
movsd 0x30(%rax), %xmm5
movsd -0x90(%rbp), %xmm4
movq -0x20(%rbp), %rax
movsd 0x38(%rax), %xmm6
mulsd %xmm6, %xmm4
addsd %xmm5, %xmm4
mulsd %xmm4, %xmm2
addsd %xmm3, %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0xa8(%rbp)
xorps %xmm0, %xmm0
ucomisd -0xa0(%rbp), %xmm0
jbe 0x97570
xorps %xmm0, %xmm0
movsd %xmm0, -0x100(%rbp)
jmp 0x97580
movsd -0xa0(%rbp), %xmm0
movsd %xmm0, -0x100(%rbp)
movsd -0x100(%rbp), %xmm0
movsd %xmm0, -0xa0(%rbp)
movsd 0x2ffe8(%rip), %xmm0    # 0xc7580
ucomisd -0xa8(%rbp), %xmm0
jbe 0x975b4
movsd 0x2ffd6(%rip), %xmm0    # 0xc7580
movsd %xmm0, -0x108(%rbp)
jmp 0x975c4
movsd -0xa8(%rbp), %xmm0
movsd %xmm0, -0x108(%rbp)
movsd -0x108(%rbp), %xmm0
movsd %xmm0, -0xa8(%rbp)
movsd -0x78(%rbp), %xmm1
movsd 0x2ffaf(%rip), %xmm0    # 0xc7590
subsd %xmm0, %xmm1
movsd 0x25873(%rip), %xmm0    # 0xbce60
mulsd %xmm1, %xmm0
divsd -0xa8(%rbp), %xmm0
movsd %xmm0, -0xb0(%rbp)
movsd 0x2679f(%rip), %xmm0    # 0xbdda8
mulsd -0x80(%rbp), %xmm0
movsd %xmm0, -0x110(%rbp)
movsd -0xb0(%rbp), %xmm1
movaps 0x25bfb(%rip), %xmm0    # 0xbd220
pand %xmm0, %xmm1
movsd 0x2ff57(%rip), %xmm0    # 0xc7588
ucomisd %xmm1, %xmm0
jbe 0x976a1
movsd -0xa0(%rbp), %xmm0
movsd -0xb0(%rbp), %xmm1
mulsd -0xb0(%rbp), %xmm1
movsd -0xb0(%rbp), %xmm3
mulsd -0xb0(%rbp), %xmm3
movsd 0x2ff39(%rip), %xmm2    # 0xc75a0
divsd %xmm2, %xmm3
movsd 0x258b5(%rip), %xmm2    # 0xbcf28
addsd %xmm3, %xmm2
mulsd %xmm2, %xmm1
movsd 0x247dd(%rip), %xmm2    # 0xbbe60
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movsd 0x2ff05(%rip), %xmm1    # 0xc7598
addsd %xmm1, %xmm0
movsd %xmm0, -0x118(%rbp)
jmp 0x976b3
movsd 0x2feef(%rip), %xmm0    # 0xc7598
movsd %xmm0, -0x118(%rbp)
jmp 0x976b3
movsd -0x110(%rbp), %xmm0
movsd -0x118(%rbp), %xmm1
mulsd %xmm1, %xmm0
movsd %xmm0, -0x8(%rbp)
movsd -0x8(%rbp), %xmm0
addq $0x120, %rsp            # imm = 0x120
popq %rbp
retq
nopw (%rax,%rax)