void
test_wsac3() {
float x[1024], y[1024];
int err, max, leven;
x[0] = 0;
y[0] = 0;
leven = TRUE;
max = 1;
/* File with one data point */
newhdr();
setnhv("npts", &max, &err, -1);
setlhv("leven", &leven, &err, -1);
setfhv("b", &(x[0]), &err, -1);
setfhv("e", &(x[0]), &err, -1);
setihv("iftype", "itime", &err, -1, -1);
/* Filename Length Special < 0 */
wsac3(FILE_WRITE, &(x[0]), &(y[0]), &err, -1);
ok(err == SAC_OK, "wsac3 writing sac file with 1 point");
file_check(FILE_WRITE, SAC_HEADER_SIZE + 1 * sizeof(float), "1 data point");
unlink(FILE_WRITE);
/* Filename Correct Length */
wsac3(FILE_WRITE, &(x[0]), &(y[0]), &err, strlen(FILE_WRITE));
ok(err == SAC_OK, "wsac3 writing sac file with 1 point");
file_check(FILE_WRITE, SAC_HEADER_SIZE + 1 * sizeof(float), "1 data point");
unlink(FILE_WRITE);
/* Filename Length too long */
wsac3(FILE_WRITE, &(x[0]), &(y[0]), &err, strlen(FILE_WRITE) * 3);
ok(err == SAC_OK, "wsac3 writing sac file with 1 point");
file_check(FILE_WRITE, SAC_HEADER_SIZE + 1 * sizeof(float), "1 data point");
unlink(FILE_WRITE);
/* Filename Length too short*/
wsac3(FILE_WRITE, &(x[0]), &(y[0]), &err, strlen(FILE_WRITE) - 4);
ok(err == SAC_OK, "wsac3 writing sac file with 1 point");
file_check("test_write", SAC_HEADER_SIZE + 1 * sizeof(float), "1 data point");
unlink("test_write");
/* Filename Length = 0 */
wsac3(FILE_WRITE, &(x[0]), &(y[0]), &err, 0);
ok(err == 101, "wsac3 writing sac file with 1 point value %d expected %d",
err, 101);
/* Increase the number of data points to 1024 */
{
int i;
for(i = 0; i < 1024; i++) {
x[i] = y[i] = i;
}
}
/* But do not tell the sacio library about the increased size*/
wsac3(FILE_WRITE, &(x[0]), &(y[0]), &err, -1);
ok(err == SAC_OK, "wsac3 writing sac file with 1 / 1024 data points "
"value %d expected %d", err, SAC_OK);
file_check(FILE_WRITE, SAC_HEADER_SIZE + (1 * sizeof(float)), "1 data point");
unlink(FILE_WRITE);
/* Inform sacio that the file is now longer */
max = 1024;
setnhv("npts", &max, &err, -1);
wsac3(FILE_WRITE, &(x[0]), &(y[0]), &err, -1);
ok(err == SAC_OK, "wsac3 writing sac file with 1024 data points "
"value %d expected %d", err, SAC_OK);
file_check(FILE_WRITE, SAC_HEADER_SIZE + (1024 * sizeof(float)),
"1024 data points");
unlink(FILE_WRITE);
/* Write a file with < 0 points */
max = -1;
setnhv("npts", &max, &err, -1);
wsac3(FILE_WRITE, &(x[0]), &(y[0]), &err, -1);
ok(err == 115, "wsac3 writing sac file with -1 data points "
"value %d expected %d", err, 115);
ok(!file_exists(FILE_WRITE), "-1 data points");
unlink(FILE_WRITE);
/* Write the other types of sac files */
max = 2;
setnhv("npts", &max, &err, -1);
/* Spectral Real / Imaginary */
setihv("iftype", "irlim", &err, -1, -1);
wsac3(FILE_WRITE, &(x[0]), &(y[0]), &err, -1);
ok(err == SAC_OK, "wsac3 writing sac file (real/imag) with 2 data points "
"value %d expected %d", err, SAC_OK);
file_check_not_ok(FILE_WRITE, SAC_HEADER_SIZE + max * 2 * sizeof(float),
"real/imag");
unlink(FILE_WRITE);
/* Spectral Amplitude / Phase */
setihv("iftype", "iamph", &err, -1, -1);
wsac3(FILE_WRITE, &(x[0]), &(y[0]), &err, -1);
ok(err == SAC_OK, "wsac3 writing sac file (amp/phase)with 2 data points "
"value %d expected %d", err, SAC_OK);
file_check_not_ok(FILE_WRITE, SAC_HEADER_SIZE + max * 2 * sizeof(float),
"amp/phase");
unlink(FILE_WRITE);
/* IXY spaced file */
setihv("iftype", "ixy", &err, -1, -1);
wsac3(FILE_WRITE, &(x[0]), &(y[0]), &err, -1);
ok(err == SAC_OK, "wsac3 writing sac file (ixy) with 2 data points "
"value %d expected %d", err, SAC_OK);
file_check_not_ok(FILE_WRITE, SAC_HEADER_SIZE + (max * 2 * sizeof(float)),
"ixy");
unlink(FILE_WRITE);
/* unevenly spaced file */
leven = 0;
setlhv("leven", &leven, &err, -1);
setihv("iftype", "itime", &err, -1, -1);
wsac3(FILE_WRITE, &(x[0]), &(y[0]), &err, -1);
ok(err == SAC_OK, "wsac3 writing sac file (uneven/time) with 2 data points "
"value %d expected %d", err, SAC_OK);
file_check(FILE_WRITE, SAC_HEADER_SIZE + max * 2 * sizeof(float),
"uneven/time");
unlink(FILE_WRITE);
}
void asc2sac(char *ascfn) {
sac_int_t npts, max_npts, nerr, nconv, itmp;
float ftmp;
char *sacfn;
float *time, *data, a, b;
FILE *f;
sacfn = (char *)malloc(strlen(ascfn) + strlen(".sac") + 1);
if (!sacfn) {
fprintf(stderr, "Out of memory\n");
exit(1);
}
strcpy(sacfn, ascfn);
strcat(sacfn, ".sac");
/* reading ascii file */
f = fopen(ascfn, "r");
if(f == NULL) {
fprintf(stderr, "Cannot open file '%s' to read\n", ascfn);
exit(1);
}
time = data = 0;
max_npts = 0;
for (npts = 0; (nconv = fscanf(f, "%f %f\n", &a, &b)) == 2; ++npts) {
if (nconv != 2) {
fprintf(stderr, "error while reading file '%s'\n", ascfn);
exit(1);
}
if (npts >= max_npts) {
max_npts = max_npts ? 2 * max_npts : 1024;
time = (float*) realloc(time, max_npts * sizeof(float));
data = (float*) realloc(data, max_npts * sizeof(float));
if(time == NULL || data == NULL) {
fprintf(stderr, "Out of memory\n");
exit(1);
}
}
time[npts] = a;
data[npts] = b;
}
if (nconv != EOF || ferror(f)) {
fprintf(stderr, "error while reading file '%s' (on or near line %ld)\n", ascfn, (long)npts);
exit(1);
}
fclose(f);
/* finished reading ascii file */
/* write SAC data usng SAC IO library */
nerr = 0;
newhdr();
setnhv("npts", &npts, &nerr, strlen("npts"));
itmp = 1;
setlhv("leven", &itmp, &nerr, strlen("leven"));
ftmp = time[1] - time[0];
setfhv("delta", &ftmp, &nerr, strlen("delta"));
setfhv("b", &(time[0]), &nerr, strlen("b"));
setfhv("e", &(time[npts-1]), &nerr, strlen("e"));
setihv("iftype", "itime", &nerr, strlen("iftype"), strlen("itime"));
setihv("idep", "idisp", &nerr, strlen("idep"), strlen("idisp"));
if(nerr) {
fprintf(stderr, "error when setting header for '%s'\n", sacfn);
exit(1);
}
wsac0(sacfn, time, data, &nerr, strlen(sacfn));
if(nerr) {
fprintf(stderr, "error when writing '%s'\n", sacfn);
exit(1);
}
free(time);
free(data);
free(sacfn);
return;
}