#include <config.h>
#ifdef SIM
#include "ntpd.h"
#include "ntp_config.h"
int determine_event_ordering(const Event *e1, const Event *e2);
int determine_recv_buf_ordering(const struct recvbuf *b1,
const struct recvbuf *b2);
void create_server_associations(void);
void init_sim_io(void);
sim_info simulation;
local_clock_info simclock;
queue *event_queue;
queue *recv_queue;
static double sys_residual = 0;
void (*event_ptr[]) (Event *) = {
sim_event_beep, sim_update_clocks, sim_event_timer, sim_event_recv_packet
};
int
determine_event_ordering(
const Event *e1,
const Event *e2
)
{
return (e1->time - e2->time);
}
int
determine_recv_buf_ordering(
const struct recvbuf *b1,
const struct recvbuf *b2
)
{
double recv_time1;
double recv_time2;
LFPTOD(&b1->recv_time, recv_time1);
LFPTOD(&b2->recv_time, recv_time2);
return (int)(recv_time1 - recv_time2);
}
void create_server_associations(void)
{
int i;
for (i = 0; i < simulation.num_of_servers; ++i) {
printf("%s\n", stoa(simulation.servers[i].addr));
if (peer_config(simulation.servers[i].addr,
NULL,
loopback_interface,
MODE_CLIENT,
NTP_VERSION,
NTP_MINDPOLL,
NTP_MAXDPOLL,
0,
0,
0,
NULL ) == 0) {
fprintf(stderr,
"ERROR!! Could not create association for: %s\n",
stoa(simulation.servers[i].addr));
}
}
}
int
ntpsim(
int argc,
char * argv[]
)
{
Event * curr_event;
struct timeval seed;
simclock.local_time = 0;
simclock.adj = 0;
simclock.slew = 500e-6;
simulation.num_of_servers = 0;
simulation.beep_delay = BEEP_DLY;
simulation.sim_time = 0;
simulation.end_time = SIM_TIME;
initializing = TRUE;
msyslog_term = TRUE;
init_sim_io();
init_auth();
init_util();
init_restrict();
init_mon();
init_timer();
init_lib();
init_request();
init_control();
init_peer();
init_proto();
init_loopfilter();
mon_start(MON_OFF);
getconfig(argc, argv);
loop_config(LOOP_DRIFTINIT, 0);
initializing = FALSE;
gettimeofday(&seed, NULL);
ntp_srandom(seed.tv_usec);
event_queue = create_priority_queue((q_order_func)
determine_event_ordering);
recv_queue = create_priority_queue((q_order_func)
determine_recv_buf_ordering);
enqueue(event_queue, event(0, BEEP));
enqueue(event_queue, event(simulation.sim_time + 1.0, TIMER));
while (simulation.sim_time <= simulation.end_time &&
(!empty(event_queue))) {
curr_event = dequeue(event_queue);
sim_update_clocks(curr_event);
(*event_ptr[curr_event->function])(curr_event);
free_node(curr_event);
}
printf("sys_received: %lu\n", sys_received);
printf("sys_badlength: %lu\n", sys_badlength);
printf("sys_declined: %lu\n", sys_declined);
printf("sys_restricted: %lu\n", sys_restricted);
printf("sys_newversion: %lu\n", sys_newversion);
printf("sys_oldversion: %lu\n", sys_oldversion);
printf("sys_limitrejected: %lu\n", sys_limitrejected);
printf("sys_badauth: %lu\n", sys_badauth);
return (0);
}
void
init_sim_io(void)
{
loopback_interface = emalloc_zero(sizeof(*loopback_interface));
ep_list = loopback_interface;
strlcpy(loopback_interface->name, "IPv4loop",
sizeof(loopback_interface->name));
loopback_interface->flags = INT_UP | INT_LOOPBACK;
loopback_interface->fd = -1;
loopback_interface->bfd = -1;
loopback_interface->ifnum = 1;
loopback_interface->family = AF_INET;
AF(&loopback_interface->sin) = AF_INET;
SET_ADDR4(&loopback_interface->sin, LOOPBACKADR);
SET_PORT(&loopback_interface->sin, NTP_PORT);
AF(&loopback_interface->mask) = AF_INET;
SET_ADDR4(&loopback_interface->mask, LOOPNETMASK);
}
Event *event(double t, funcTkn f)
{
Event *e;
if ((e = get_node(sizeof(*e))) == NULL)
abortsim("get_node failed in event");
e->time = t;
e->function = f;
return (e);
}
void sim_event_timer(Event *e)
{
struct recvbuf *rbuf;
timer();
while (!empty(recv_queue)) {
rbuf = (struct recvbuf *)dequeue(recv_queue);
(*rbuf->receiver)(rbuf);
free_node(rbuf);
}
enqueue(event_queue,
event(simulation.sim_time + (1 << EVENT_TIMEOUT), TIMER));
}
int simulate_server(
sockaddr_u *serv_addr,
endpt * inter,
struct pkt *rpkt
)
{
struct pkt xpkt;
struct recvbuf rbuf;
Event *e;
server_info *server;
script_info *curr_script;
int i;
double d1, d2, d3;
double t1, t2, t3, t4;
l_fp lfp_host;
ZERO(xpkt);
ZERO(rbuf);
server = NULL;
for (i = 0; i < simulation.num_of_servers; ++i) {
if (memcmp(simulation.servers[i].addr, serv_addr,
sizeof(*serv_addr)) == 0) {
server = &simulation.servers[i];
break;
}
}
fprintf(stderr, "Received packet from %s on %s\n",
stoa(serv_addr), latoa(inter));
if (server == NULL)
abortsim("Server with specified address not found!!!");
curr_script = server->curr_script;
xpkt.li_vn_mode = PKT_LI_VN_MODE(LEAP_NOWARNING, NTP_VERSION,
MODE_SERVER);
xpkt.stratum = STRATUM_TO_PKT(((u_char)1));
memcpy(&xpkt.refid, "GPS", 4);
xpkt.ppoll = rpkt->ppoll;
xpkt.precision = rpkt->precision;
xpkt.rootdelay = 0;
xpkt.rootdisp = 0;
d1 = poisson(curr_script->prop_delay, curr_script->jitter);
d2 = poisson(curr_script->proc_delay, 0);
d3 = poisson(curr_script->prop_delay, curr_script->jitter);
NTOHL_FP(&rpkt->xmt, &lfp_host);
LFPTOD(&lfp_host, t1);
t2 = server->server_time + d1;
t3 = server->server_time + d1 + d2;
t4 = t1 + d1 + d2 + d3;
xpkt.org = rpkt->xmt;
DTOLFP(t2, &lfp_host);
HTONL_FP(&lfp_host, &xpkt.rec);
DTOLFP(t3, &lfp_host);
HTONL_FP(&lfp_host, &xpkt.xmt);
xpkt.reftime = xpkt.xmt;
rbuf.used = 1;
rbuf.receiver = &receive;
rbuf.recv_length = LEN_PKT_NOMAC;
rbuf.recv_pkt = xpkt;
rbuf.dstadr = inter;
rbuf.fd = inter->fd;
memcpy(&rbuf.srcadr, serv_addr, sizeof(rbuf.srcadr));
memcpy(&rbuf.recv_srcadr, serv_addr, sizeof(rbuf.recv_srcadr));
e = event(t4, PACKET);
e->rcv_buf = rbuf;
enqueue(event_queue, e);
if (curr_script->duration > simulation.sim_time &&
NULL == HEAD_PFIFO(server->script)) {
printf("Hello\n");
UNLINK_FIFO(curr_script, *server->script, link);
}
return (0);
}
void sim_update_clocks(Event *e)
{
double time_gap;
double adj;
int i;
time_gap = e->time - simulation.sim_time;
if (time_gap < 0)
printf("WARNING: e->time %.6g comes before sim_time %.6g (gap %+.6g)\n",
e->time, simulation.sim_time, time_gap);
if (e->time + time_gap < simclock.local_time)
printf("WARNING: e->time + gap %.6g comes before local_time %.6g\n",
e->time + time_gap, simclock.local_time);
simclock.local_time = e->time + time_gap;
simulation.sim_time = e->time;
for (i = 0; i < simulation.num_of_servers; ++i) {
simulation.servers[i].curr_script->freq_offset +=
gauss(0, time_gap * simulation.servers[i].curr_script->wander);
simulation.servers[i].server_time += time_gap *
(1 + simulation.servers[i].curr_script->freq_offset);
}
adj = time_gap * simclock.slew;
if (adj < fabs(simclock.adj)) {
if (simclock.adj < 0) {
simclock.adj += adj;
simclock.local_time -= adj;
} else {
simclock.adj -= adj;
simclock.local_time += adj;
}
} else {
simclock.local_time += simclock.adj;
simclock.adj = 0;
}
}
void sim_event_recv_packet(Event *e)
{
struct recvbuf *rbuf;
if ((rbuf = get_node(sizeof(*rbuf))) == NULL)
abortsim("get_node failed in sim_event_recv_packet");
memcpy(rbuf, &e->rcv_buf, sizeof(*rbuf));
DTOLFP(simclock.local_time, &rbuf->recv_time);
enqueue(recv_queue, rbuf);
}
void sim_event_beep(Event *e)
{
#if 0
static int first_time = 1;
char *dash = "-----------------";
#endif
fprintf(stderr, "BEEP!!!\n");
enqueue(event_queue, event(e->time + simulation.beep_delay, BEEP));
#if 0
if(simulation.beep_delay > 0) {
if (first_time) {
printf("\t%4c T %4c\t%4c T+ERR %3c\t%5cT+ERR+NTP\n",
' ', ' ', ' ', ' ',' ');
printf("\t%s\t%s\t%s\n", dash, dash, dash);
first_time = 0;
printf("\t%16.6f\t%16.6f\t%16.6f\n",
n->time, n->clk_time, n->ntp_time);
return;
}
printf("\t%16.6f\t%16.6f\t%16.6f\n",
simclock.local_time,
n->time, n->clk_time, n->ntp_time);
#endif
}
void abortsim(char *errmsg)
{
perror(errmsg);
exit(1);
}
void
get_systime(
l_fp *now )
{
if (simclock.local_time == simclock.last_read_time)
simclock.local_time += 200e-9;
simclock.last_read_time = simclock.local_time;
DTOLFP(simclock.local_time, now);
}
int
adj_systime(
double now
)
{
struct timeval adjtv;
double dtemp;
long ticks;
int isneg = 0;
dtemp = now + sys_residual;
if (dtemp < 0) {
isneg = 1;
dtemp = -dtemp;
}
adjtv.tv_sec = (long)dtemp;
dtemp -= adjtv.tv_sec;
ticks = (long)(dtemp / sys_tick + .5);
adjtv.tv_usec = (long)(ticks * sys_tick * 1e6);
dtemp -= adjtv.tv_usec / 1e6;
sys_residual = dtemp;
if (isneg) {
adjtv.tv_sec = -adjtv.tv_sec;
adjtv.tv_usec = -adjtv.tv_usec;
sys_residual = -sys_residual;
}
simclock.adj = now;
return (1);
}
int
step_systime(
double now
)
{
#ifdef DEBUG
if (debug)
printf("step_systime: time %.6f adj %.6f\n",
simclock.local_time, now);
#endif
simclock.local_time += now;
return (1);
}
double
gauss(
double m,
double s
)
{
double q1, q2;
if (s == 0)
return (m);
while ((q1 = drand48()) == 0)
;
q2 = drand48();
return (m + s * sqrt(-2. * log(q1)) * cos(2. * PI * q2));
}
double
poisson(
double m,
double s
)
{
double q1;
if (s == 0)
return (m);
while ((q1 = drand48()) == 0)
;
return (m - s * log(q1 * s));
}
#endif