a2b.c

seismic software,very useful

/* Copyright (c) Colorado School of Mines, 1990.
/* All rights reserved.                       */

/* a2b - convert ascii floats to binary */

#include "par.h"

/*********************** self documentation **********************/
string sdoc =
" 								\n"
" A2B - convert ascii floats to binary 				\n"
" 								\n"
" a2b <stdin >stdout outpar=/dev/tty 				\n"
" 								\n"
" Required parameters:						\n"
" 	none							\n"
" 								\n"
" Optional parameters:						\n"
" 	n1=2		floats per line in input file		\n"
" 								\n"
" 	outpar=/dev/tty	output parameter file, contains the	\n"
"			number of lines (n=)			\n"
" 								\n"
;
/**************** end self doc ***********************************/

/* Credits:
 *	CWP: Jack, Dave
 */


main(int argc, char **argv)
{
	string outpar;		/* name of file holding output parfile	*/
	FILE *outparfp;		/* ... its file pointer			*/
	int n1;			/* number of floats per line		*/
	int n2 = 0;		/* number of lines in input file 	*/
	char buf[BUFSIZ];	/* buffer for the ascii floats		*/
	float *x;		/* binary floats			*/


	/* Hook up getpar */
	initargs(argc, argv);
	askdoc(1);


	/* Prevent floats from dumping on screen */
	switch(filestat(STDOUT)) {
	case BADFILETYPE:
		warn("stdout is illegal filetype");
		selfdoc();
	break;
	case TTY:
		warn("stdout can't be tty");
		selfdoc();
	break;
	}


	/* Get parameters and do set up */
	if (!getparstring("outpar", &outpar))  outpar = "/dev/tty" ;
	outparfp = efopen(outpar, "w");

	if (!getparint("n1", &n1))  n1 = 2;
	x = ealloc1float(n1);


	/* Loop over data, converting to binary */
	while (gets(buf)) {
		char *p = buf;
		register int i1;

		for (i1 = 0; i1 < n1; ++i1) {
			while (*p == ' ' || *p == '\t')  ++p;
			if (1 != sscanf(p, "%f", x+i1))
				err("line #%d: scan failed on float #%d:\n%s",
					n2+1, i1+1, buf);
			while (*p != ' ' && *p != '\t' && *p != '\0')  ++p;
		}
		efwrite(x, FSIZE, n1, stdout);
		++n2;
	}


	/* Make par file */
	fprintf(outparfp, "n=%d\n", n2);


	return EXIT_SUCCESS;
}