static bool parse_listen_faces(char str[], ccarray_t * faces)
{
static const char delims[] = " \t\n,;";
char * tok;
tok = strtok(str, delims);
while ( tok && ccarray_size(faces) < ccarray_capacity(faces) ) {
uint32_t address = 0;
uint16_t port = 0;
if ( getifaddr(tok, &address, &port) == -1 ) {
fprintf(stderr, "FATAL: Can't get address for '%s': %s\n", tok, strerror(errno));
fprintf(stderr, "Check if device name is valid and device is up\n");
return false;
}
if ( !port ) {
fprintf(stderr, "No port specified in '%s'\n", tok);
return false;
}
ccarray_push_back(faces, &(struct sockaddr_in ) {
.sin_family = AF_INET,
.sin_addr.s_addr = htonl(address),
.sin_port = htons(port),
.sin_zero = { 0 }
});
tok = strtok(NULL, delims);
}
/** main() */
int main(int argc, char *argv[])
{
const char * fname[2] =
{ NULL, NULL };
FILE * fp[2] =
{ NULL, NULL };
compression_t compression[2] =
{ compression_unknown, compression_unknown };
int rc[2] =
{ -1, -1 };
int dc[2] =
{ -1, -1 };
int ru[2] =
{ radians, radians };
int du[2] =
{ radians, radians };
ccarray_t * list[2] =
{ NULL, NULL };
size_t capacity[2] =
{ 15000000, 15000000 };
char head[2][MAX_HEADER_LENGTH];
int beverbose = 0;
int invert_match = 0;
double r = -1;
int i;
size_t size1, size2, pos1, pos2;
const obj_t * obj1, * obj2;
//double rcd;
/* parse command line */
for ( i = 1; i < argc; ++i )
{
if ( strcmp(argv[i], "--help") == 0 || strcmp(argv[i], "-help") == 0 ) {
show_usage(stdout, argc, argv);
return 0;
}
if ( strncmp(argv[i], "rc1=", 4) == 0 )
{
if ( sscanf(argv[i] + 4, "%d", &rc[0]) != 1 || rc[0] < 1 ) {
fprintf(stderr,"Invalid value of %s\n", argv[i]);
return -1;
}
}
else if ( strncmp(argv[i], "rc2=", 4) == 0 )
{
if ( sscanf(argv[i] + 4, "%d", &rc[1]) != 1 || rc[1] < 1 ) {
fprintf(stderr,"Invalid value of %s\n", argv[i]);
return -1;
}
}
else if ( strncmp(argv[i], "dc1=", 4) == 0 )
{
if ( sscanf(argv[i] + 4, "%d", &dc[0]) != 1 || dc[0] < 1 ) {
fprintf(stderr,"Invalid value of %s\n", argv[i]);
return -1;
}
}
else if ( strncmp(argv[i], "dc2=", 4) == 0 )
{
if ( sscanf(argv[i] + 4, "%d", &dc[1]) != 1 || dc[1] < 1 ) {
fprintf(stderr,"Invalid value of %s\n", argv[i]);
return -1;
}
}
else if ( strncmp(argv[i], "ru1=", 4) == 0 )
{
if ( (ru[0] = parse_unit(argv[i] + 4)) == -1 ) {
fprintf(stderr,"Invalid value of %s\n", argv[i]);
return -1;
}
}
else if ( strncmp(argv[i], "ru2=", 4) == 0 )
{
if ( (ru[1] = parse_unit(argv[i] + 4)) == -1 ) {
fprintf(stderr,"Invalid value of %s\n", argv[i]);
return -1;
}
}
else if ( strncmp(argv[i], "du1=", 4) == 0 )
{
if ( (du[0] = parse_unit(argv[i] + 4)) == -1 ) {
fprintf(stderr,"Invalid value of %s\n", argv[i]);
return -1;
}
}
else if ( strncmp(argv[i], "du2=", 4) == 0 )
{
if ( (du[1] = parse_unit(argv[i] + 4)) == -1 ) {
fprintf(stderr,"Invalid value of %s\n", argv[i]);
return -1;
}
}
else if ( strncmp(argv[i], "r=", 2) == 0 )
{
if ( sscanf(argv[i] + 2, "%lf", &r) != 1 || r <= 0 ) {
fprintf(stderr,"Invalid value of %s\n", argv[i]);
return -1;
}
}
else if ( strcmp(argv[i], "-i") == 0 ) {
invert_match = 1;
}
else if ( strcmp(argv[i], "-v") == 0 ) {
beverbose = 1;
}
else if ( !fname[0] ) {
fname[0] = argv[i];
}
else if ( !fname[1] ) {
fname[1] = argv[i];
}
else {
fprintf(stderr, "Invalid argument '%s'. Note that only 2 input files are allowed\n", argv[i]);
return -1;
}
}
/* check command line inputs */
if ( !fname[0] || !fname[1] ) {
fprintf(stderr,"Two input file names expected\n");
show_usage(stderr, argc, argv);
return -1;
}
if ( rc[0] < 1 ) {
fprintf(stderr,"rc1 argument is mandatory\n");
show_usage(stderr, argc, argv);
return -1;
}
if ( rc[1] < 1 ) {
fprintf(stderr,"rc2 argument is mandatory\n");
show_usage(stderr, argc, argv);
return -1;
}
if ( dc[0] < 1 ) {
fprintf(stderr,"dc1 argument is mandatory\n");
show_usage(stderr, argc, argv);
return -1;
}
if ( dc[1] < 1 ) {
fprintf(stderr,"dc2 argument is mandatory\n");
show_usage(stderr, argc, argv);
return -1;
}
if ( (r *= (PI / (180 * 3600))) <= 0 ) {
fprintf(stderr,"r argument is mandatory\n");
show_usage(stderr, argc, argv);
return -1;
}
/* check if input files are readable */
for ( i = 0; i < 2; ++i )
{
if ( access(fname[i], R_OK) != 0 ) {
fprintf(stderr, "Can't read %s: %s\n", fname[i], strerror(errno));
return -1;
}
}
/* allocate memory storage */
for ( i = 0; i < 2; ++i )
{
if ( !(list[i] = ccarray_create(capacity[i], sizeof(obj_t))) ) {
fprintf(stderr, "ccarray_create(capacity=%zu) fails: %s\n", capacity[i], strerror(errno));
return -1;
}
}
/* load input files */
for ( i = 0; i < 2; ++i )
{
if ( beverbose ) {
fprintf(stderr,"loading %s....\n", fname[i]);
}
if ( !(fp[i] = open_file(fname[i], &compression[i])) ) {
fprintf(stderr, "Can't read '%s': %s\n", fname[i], strerror(errno));
return -1;
}
if ( load_objects(fp[i], rc[i], dc[i], ru[i], du[i], head[i], list[i]) != 0 ) {
fprintf(stderr, "Can't load %s\n", fname[i]);
return -1;
}
close_file( fp[i], compression[i] );
if ( beverbose ) {
fprintf(stderr,"%s: %zu objects\n", fname[i], ccarray_size(list[i]));
}
if ( ccarray_size(list[i]) == capacity[i] && !feof(fp[i]) ) {
fprintf(stderr, "load_objects() fails for %s: Too many objects. Try increase capacity1 (currently %zu)\n",
fname[i], capacity[i]);
return -1;
}
}
/* sort lists */
for ( i = 0; i < 2; ++i )
{
if ( beverbose ) {
fprintf(stderr, "sort %s....\n", fname[i]);
}
//ccarray_sort(list[i], 0, ccarray_size(list[i]), cmpradec);
ccarray_sort(list[i], 0, ccarray_size(list[i]), cmpdecra);
}
/* search pairs */
if ( beverbose ) {
fprintf(stderr,"search pairs...\n");
}
printf("%s\t%s\tdra\tddec\n", head[0], head[1]);
size1 = ccarray_size(list[0]);
size2 = ccarray_size(list[1]);
double dra_closest=0;
double ddec_closest=0;
size_t imin=0;
size_t imax=size2;
for ( pos1 = 0; pos1 < size1; ++pos1 )
{
char line_closest [MAX_INPUT_LINE_LENGTH];
line_closest[MAX_INPUT_LINE_LENGTH-1]=0;
double decmin, decmax, dra, ddec;
int objcount = 0;
obj1 = ccarray_peek(list[0], pos1);
decmin = obj1->dec - r;
decmax = obj1->dec + r;
get_index(list[1],&imin,&imax,size2,decmin,decmax);
double rmin=r;
double dr;
for ( pos2=imin; pos2 < imax; ++pos2 )
{
obj2 = ccarray_peek(list[1], pos2);
dra = (obj2->ra - obj1->ra) * cos((obj1->dec + obj2->dec) / 2);
ddec = (obj2->dec - obj1->dec);
dr=hypot(dra, ddec);
if ( dr< rmin )
{
rmin=dr;
dra_closest=dra;
ddec_closest=ddec;
strcpy(line_closest,obj2->line);
if ( line_closest[MAX_INPUT_LINE_LENGTH - 1] != 0 ) {
fprintf(stderr,"too long input line in this file\n");
return -1;
}
++objcount;
}
}
if ( !invert_match && objcount>0) {
printf("%s\t%s\t%+9.3f\t%+9.3f\n", obj1->line, line_closest, dra_closest * 180 * 3600 / PI, ddec_closest * 180 * 3600 / PI);
}
if ( invert_match && objcount == 0 ) {
printf("%s\n", obj1->line);
}
}
return 0;
}
/* load lines from input file */
static int load_objects(FILE * input, int rc, int dc, int ru, int du, char header[MAX_HEADER_LENGTH],
ccarray_t * objects)
{
char line[MAX_INPUT_LINE_LENGTH];
int ic;
char * pc;
size_t size;
size_t capacity;
obj_t * obj;
header[MAX_HEADER_LENGTH-1] = 0;
if ( !fgets(header, MAX_HEADER_LENGTH, input) ) {
fprintf(stderr,"fgets(header) fails: %d (%s)\n", errno, strerror(errno));
return -1;
}
if ( header[MAX_HEADER_LENGTH - 1] != 0 ) {
fprintf(stderr,"too long header line in this file\n");
return -1;
}
/* remove trailing new line */
header[strlen(header) - 1] = 0;
size = ccarray_size(objects);
capacity = ccarray_capacity(objects);
while ( size < capacity && !feof(input) )
{
obj = ccarray_peek(objects, size);
line[MAX_INPUT_LINE_LENGTH-1] = 0;
if ( !fgets(line, MAX_INPUT_LINE_LENGTH, input) ) {
break;
}
if ( line[MAX_INPUT_LINE_LENGTH - 1] != 0 ) {
fprintf(stderr,"too long input line in this file\n");
return -1;
}
/* remove trailing new line */
line[strlen(line) - 1] = 0;
for ( pc = line, ic = 1; ic < rc; ++ic ) {
if ( !(pc = strchr(pc + 1, '\t')) ) {
break;
}
}
if ( ic != rc || sscanf(pc, " %lf", &obj->ra) != 1 ) {
continue;
}
for ( pc = line, ic = 1; ic < dc; ++ic ) {
if ( !(pc = strchr(pc + 1, '\t')) ) {
break;
}
}
if ( ic != dc || sscanf(pc, " %lf", &obj->dec) != 1 ) {
continue;
}
if ( ru == degrees ) {
obj->ra *= PI / 180;
}
if ( du == degrees ) {
obj->dec *= PI / 180;
}
obj->line = strdup(line);
ccarray_set_size(objects, ++size );
}
return 0;
}
