cwp.h File Reference

#include <stdio.h>
#include <stddef.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <ctype.h>
#include <limits.h>
#include <float.h>
#include <fcntl.h>
#include <unistd.h>
#include <sys/types.h>

Go to the source code of this file.

Defines
#define	NULL   ((void *)0)
#define	EXIT_FAILURE   (1)
#define	EXIT_SUCCESS   (0)
#define	SEEK_SET   (0)
#define	SEEK_CUR   (1)
#define	SEEK_END   (2)
#define	PI   (3.141592653589793)
#define	GOLDEN_RATIO   (1.618034)
#define	TRUE   (1)
#define	FALSE   (0)
#define	YES   (1)
#define	NO   (0)
#define	SGN(x)   ((x) < 0 ? -1.0 : 1.0)
#define	ABS(x)   ((x) < 0 ? -(x) : (x))
#define	MAX(x, y)   ((x) > (y) ? (x) : (y))
#define	MIN(x, y)   ((x) < (y) ? (x) : (y))
#define	NINT(x)   ((int)((x)>0.0?(x)+0.5:(x)-0.5))
#define	CLOSETO(x, y)   ((ABS((x) - (y)) <= FLT_EPSILON*ABS(y))?cwp_true:cwp_false)
#define	ISODD(n)   ((n) & 01)
#define	ISIZE   sizeof(int)
#define	FSIZE   sizeof(float)
#define	DSIZE   sizeof(double)
#define	STREQ(s, t)   (strcmp(s,t) == 0)
#define	STRLT(s, t)   (strcmp(s,t) < 0)
#define	STRGT(s, t)   (strcmp(s,t) > 0)
#define	DIM(a)   (sizeof(a)/sizeof(a[0]))
#define	off_t   long
#define	ftello   ftell
#define	fseeko   fseek
Functions
void *	alloc1 (size_t n1, size_t size)
void *	realloc1 (void *v, size_t n1, size_t size)
void **	alloc2 (size_t n1, size_t n2, size_t size)
void ***	alloc3 (size_t n1, size_t n2, size_t n3, size_t size)
void ****	alloc4 (size_t n1, size_t n2, size_t n3, size_t n4, size_t size)
void *****	alloc5 (size_t n1, size_t n2, size_t n3, size_t n4, size_t n5, size_t size)
void ******	alloc6 (size_t n1, size_t n2, size_t n3, size_t n4, size_t n5, size_t n6, size_t size)
void	free1 (void *p)
void	free2 (void **p)
void	free3 (void ***p)
void	free4 (void ****p)
void	free5 (void *****p)
void	free6 (void ******p)
int *	alloc1int (size_t n1)
int *	realloc1int (int *v, size_t n1)
int **	alloc2int (size_t n1, size_t n2)
int ***	alloc3int (size_t n1, size_t n2, size_t n3)
float *	alloc1float (size_t n1)
float *	realloc1float (float *v, size_t n1)
float **	alloc2float (size_t n1, size_t n2)
float ***	alloc3float (size_t n1, size_t n2, size_t n3)
float ****	alloc4float (size_t n1, size_t n2, size_t n3, size_t n4)
void	free4float (float ****p)
float *****	alloc5float (size_t n1, size_t n2, size_t n3, size_t n4, size_t n5)
void	free5float (float *****p)
float ******	alloc6float (size_t n1, size_t n2, size_t n3, size_t n4, size_t n5, size_t n6)
void	free6float (float ******p)
int ****	alloc4int (size_t n1, size_t n2, size_t n3, size_t n4)
void	free4int (int ****p)
int *****	alloc5int (size_t n1, size_t n2, size_t n3, size_t n4, size_t n5)
void	free5int (int *****p)
unsigned short ******	alloc6ushort (size_t n1, size_t n2, size_t n3, size_t n4, size_t n5, size_t n6)
unsigned char *****	alloc5uchar (size_t n1, size_t n2, size_t n3, size_t n4, size_t n5)
void	free5uchar (unsigned char *****p)
unsigned short *****	alloc5ushort (size_t n1, size_t n2, size_t n3, size_t n4, size_t n5)
void	free5ushort (unsigned short *****p)
unsigned char ******	alloc6uchar (size_t n1, size_t n2, size_t n3, size_t n4, size_t n5, size_t n6)
void	free6uchar (unsigned char ******p)
void	free6ushort (unsigned short ******p)
double *	alloc1double (size_t n1)
double *	realloc1double (double *v, size_t n1)
double **	alloc2double (size_t n1, size_t n2)
double ***	alloc3double (size_t n1, size_t n2, size_t n3)
void	free1int (int *p)
void	free2int (int **p)
void	free3int (int ***p)
void	free1float (float *p)
void	free2float (float **p)
void	free3float (float ***p)
void	free1double (double *p)
void	free2double (double **p)
void	free3double (double ***p)
void *	bmalloc (int nbpe, int n1, int n2)
void	bmfree (void *bm)
void	bmread (void *bm, int dir, int k1, int k2, int n, void *v)
void	bmwrite (void *bm, int dir, int k1, int k2, int n, void *v)
float	fsinc (float x)
double	dsinc (double x)
void	mksinc (float d, int lsinc, float sinc[])
void	ints8r (int nxin, float dxin, float fxin, float yin[], float yinl, float yinr, int nxout, float xout[], float yout[])
void	intt8r (int ntable, float table[][8], int nxin, float dxin, float fxin, float yin[], float yinl, float yinr, int nxout, float xout[], float yout[])
void	ress8r (int nxin, float dxin, float fxin, float yin[], float yinl, float yinr, int nxout, float dxout, float fxout, float yout[])
void	shfs8r (float dx, int nxin, float fxin, float yin[], float yinl, float yinr, int nxout, float fxout, float yout[])
void	xindex (int nx, float ax[], float x, int *index)
void	intl2b (int nxin, float dxin, float fxin, int nyin, float dyin, float fyin, unsigned char *zin, int nxout, float dxout, float fxout, int nyout, float dyout, float fyout, unsigned char *zout)
void	intlin (int nin, float xin[], float yin[], float yinl, float yinr, int nout, float xout[], float yout[])
void	intcub (int ideriv, int nin, float xin[], float ydin[][4], int nout, float xout[], float yout[])
void	cakima (int n, float x[], float y[], float yd[][4])
void	cmonot (int n, float x[], float y[], float yd[][4])
void	csplin (int n, float x[], float y[], float yd[][4])
void	yxtoxy (int nx, float dx, float fx, float y[], int ny, float dy, float fy, float xylo, float xyhi, float x[])
void	bfhighpass (int npoles, float f3db, int n, float p[], float q[])
void	bflowpass (int npoles, float f3db, int n, float p[], float q[])
void	bfdesign (float fpass, float apass, float fstop, float astop, int *npoles, float *f3db)
void	mkdiff (int n, float a, float h, int l, int m, float d[])
void	mkhdiff (float h, int l, float d[])
void	holbergd1 (float e, int n, float d[])
void	conv (int lx, int ifx, float *x, int ly, int ify, float *y, int lz, int ifz, float *z)
void	xcor (int lx, int ifx, float *x, int ly, int ify, float *y, int lz, int ifz, float *z)
void	hilbert (int n, float x[], float y[])
void	antialias (float frac, int phase, int n, float p[], float q[])
void *	abelalloc (int n)
void	abelfree (void *at)
void	abel (void *at, float f[], float g[])
void *	hankelalloc (int nfft)
void	hankelfree (void *ht)
void	hankel0 (void *ht, float f[], float h[])
void	hankel1 (void *ht, float f[], float h[])
void	hpsort (int n, float a[])
void	qksort (int n, float a[])
void	qkfind (int m, int n, float a[])
void	qkisort (int n, float a[], int i[])
void	qkifind (int m, int n, float a[], int i[])
float	quest (float p, int n, float x[])
void *	questalloc (float p, int n, float x[])
float	questupdate (void *q, int n, float x[])
void	questfree (void *q)
void	swap_short_2 (short *tni2)
void	swap_u_short_2 (unsigned short *tni2)
void	swap_int_4 (int *tni4)
void	swap_u_int_4 (unsigned int *tni4)
void	swap_long_4 (long *tni4)
void	swap_u_long_4 (unsigned long *tni4)
void	swap_float_4 (float *tnf4)
void	swap_double_8 (double *tndd8)
int	npfa (int nmin)
int	npfao (int nmin, int nmax)
int	npfar (int nmin)
int	npfaro (int nmin, int nmax)
int	npfa_d (int nmin)
int	npfao_d (int nmin, int nmax)
int	npfar_d (int nmin)
int	npfaro_d (int nmin, int nmax)
int	isamax (int n, float *sx, int incx)
float	sasum (int n, float *sx, int incx)
void	saxpy (int n, float sa, float *sx, int incx, float *sy, int incy)
void	scopy (int n, float *sx, int incx, float *sy, int incy)
float	sdot (int n, float *sx, int incx, float *sy, int incy)
float	snrm2 (int n, float *sx, int incx)
void	sscal (int n, float sa, float *sx, int incx)
void	sswap (int n, float *sx, int incx, float *sy, int incy)
int	idamax (int n, double *sx, int incx)
double	dasum (int n, double *sx, int incx)
void	daxpy (int n, double sa, double *sx, int incx, double *sy, int incy)
void	dcopy (int n, double *sx, int incx, double *sy, int incy)
double	ddot (int n, double *sx, int incx, double *sy, int incy)
double	dnrm2 (int n, double *sx, int incx)
void	dscal (int n, double sa, double *sx, int incx)
void	dswap (int n, double *sx, int incx, double *sy, int incy)
void	sgeco (float **a, int n, int *ipvt, float *rcond, float *z)
void	sgefa (float **a, int n, int *ipvt, int *info)
void	sgesl (float **a, int n, int *ipvt, float *b, int job)
void	sqrdc (float **x, int n, int p, float *qraux, int *jpvt, float *work, int job)
void	sqrsl (float **x, int n, int k, float *qraux, float *y, float *qy, float *qty, float *b, float *rsd, float *xb, int job, int *info)
void	sqrst (float **x, int n, int p, float *y, float tol, float *b, float *rsd, int *k, int *jpvt, float *qraux, float *work)
void	dgeco (double **a, int n, int *ipvt, double *rcond, double *z)
void	dgefa (double **a, int n, int *ipvt, int *info)
void	dgesl (double **a, int n, int *ipvt, double *b, int job)
void	stoepd (int n, double r[], double g[], double f[], double a[])
void	stoepf (int n, float r[], float g[], float f[], float a[])
void	vanded (int n, double v[], double b[], double x[])
void	vandef (int n, float v[], float b[], float x[])
void	tridif (int n, float a[], float b[], float c[], float r[], float u[])
void	tridid (int n, double a[], double b[], double c[], double r[], double u[])
void	tripd (float *d, float *e, float *b, int n)
void	tripp (int n, float *d, float *e, float *c, float *b)
int	mnewt (int maxiter, float ftol, float dxtol, int n, float *x, void *aux, void(*fdfdx)(int n, float *x, float *f, float **dfdx, void *aux))
void	recttopolar (int nx, float dx, float fx, int ny, float dy, float fy, float **p, int na, float da, float fa, int nr, float dr, float fr, float **q)
void	polartorect (int na, float da, float fa, int nr, float dr, float fr, float **q, int nx, float dx, float fx, int ny, float dy, float fy, float **p)
void	rfwtva (int n, float z[], float zmin, float zmax, float zbase, int yzmin, int yzmax, int xfirst, int xlast, int wiggle, int nbpr, unsigned char *bits, int endian)
void	rfwtvaint (int n, float z[], float zmin, float zmax, float zbase, int yzmin, int yzmax, int xfirst, int xlast, int wiggle, int nbpr, unsigned char *bits, int endian)
void	scaxis (float x1, float x2, int *nxnum, float *dxnum, float *fxnum)
int	yclip (int nx, float dx, float fx, float y[], float ymin, float ymax, float xc[], float yc[])
float	airya (float x)
float	airyb (float x)
float	airyap (float x)
float	airybp (float x)
float	cpusec (void)
float	cputime (void)
float	wallsec (void)
float	walltime (void)
float	franuni (void)
void	sranuni (int seed)
float	frannor (void)
void	srannor (int seed)
void	LU_decomposition (int nrows, float **matrix, int *index, float *d)
void	backward_substitution (int nrows, float **matrix, int *index, float *b)
void	inverse_matrix (int nrows, float **matrix)
void	inverse_matrix_multiply (int nrows1, float **matrix1, int ncols2, int nrows2, float **matrix2, float **out_matrix)
void	compute_svd (float **a, int m, int n, float w[], float **v)
void	svd_backsubstitute (float **u, float w[], float **v, int m, int n, float b[], float x[])
void	ricker1_wavelet (int nt, float dt, float fpeak, float *wavelet)
void	ricker2_wavelet (int hlw, float dt, float period, float ampl, float distort, float *wavelet)
void	akb_wavelet (int nt, float dt, float fpeak, float *wavelet)
void	spike_wavelet (int nt, int tindex, float *wavelet)
void	unit_wavelet (int nt, float *wavelet)
void	zero_wavelet (int nt, float *wavelet)
void	pp1d (FILE *fp, char *title, int lx, int ifx, float x[])
void	pplot1 (FILE *fp, char *title, int nx, float ax[])
char *	temporary_filename (char *prefix)
void	zasc (char *ainput, char *aoutput, int nchar)
void	zebc (char *ainput, char *aoutput, int nchar)
void	IBMFLT (float *in, float *out, int *nwds, int *idirec)
char *	cwp_strdup (char *str)
void	strchop (char *s, char *t)
int	CWP_Exit ()

---

Define Documentation

#define NULL   ((void *)0)

Definition at line 49 of file cwp.h.

#define EXIT_FAILURE   (1)

Definition at line 52 of file cwp.h.

#define EXIT_SUCCESS   (0)

Definition at line 55 of file cwp.h.

#define SEEK_SET   (0)

Definition at line 58 of file cwp.h.

#define SEEK_CUR   (1)

Definition at line 61 of file cwp.h.

#define SEEK_END   (2)

Definition at line 64 of file cwp.h.

#define PI   (3.141592653589793)

Definition at line 67 of file cwp.h.

#define GOLDEN_RATIO   (1.618034)

Definition at line 70 of file cwp.h.

#define TRUE   (1)

Definition at line 73 of file cwp.h.

#define FALSE   (0)

Definition at line 76 of file cwp.h.

#define YES   (1)

Definition at line 79 of file cwp.h.

#define NO   (0)

Definition at line 82 of file cwp.h.

#define SGN

(

x

)

((x) < 0 ? -1.0 : 1.0)

Definition at line 85 of file cwp.h.

#define ABS

(

x

)

((x) < 0 ? -(x) : (x))

Definition at line 88 of file cwp.h.

#define MAX	(	x,
		y		)	((x) > (y) ? (x) : (y))

Definition at line 91 of file cwp.h.

#define MIN	(	x,
		y		)	((x) < (y) ? (x) : (y))

Definition at line 94 of file cwp.h.

#define NINT

(

x

)

((int)((x)>0.0?(x)+0.5:(x)-0.5))

Definition at line 96 of file cwp.h.

#define CLOSETO	(	x,
		y		)	((ABS((x) - (y)) <= FLT_EPSILON*ABS(y))?cwp_true:cwp_false)

Definition at line 97 of file cwp.h.

#define ISODD

(

n

)

((n) & 01)

Definition at line 98 of file cwp.h.

#define ISIZE   sizeof(int)

Definition at line 99 of file cwp.h.

#define FSIZE   sizeof(float)

Definition at line 100 of file cwp.h.

#define DSIZE   sizeof(double)

Definition at line 101 of file cwp.h.

#define STREQ	(	s,
		t		)	(strcmp(s,t) == 0)

Definition at line 102 of file cwp.h.

#define STRLT	(	s,
		t		)	(strcmp(s,t) < 0)

Definition at line 103 of file cwp.h.

#define STRGT	(	s,
		t		)	(strcmp(s,t) > 0)

Definition at line 104 of file cwp.h.

#define DIM

(

a

)

(sizeof(a)/sizeof(a[0]))

Definition at line 105 of file cwp.h.

#define off_t   long

Definition at line 416 of file cwp.h.

#define ftello   ftell

Definition at line 420 of file cwp.h.

#define fseeko   fseek

Definition at line 424 of file cwp.h.

---

Function Documentation

void* alloc1	(	size_t	n1,
		size_t	size
	)

void* realloc1	(	void *	v,
		size_t	n1,
		size_t	size
	)

void** alloc2	(	size_t	n1,
		size_t	n2,
		size_t	size
	)

void*** alloc3	(	size_t	n1,
		size_t	n2,
		size_t	n3,
		size_t	size
	)

void**** alloc4	(	size_t	n1,
		size_t	n2,
		size_t	n3,
		size_t	n4,
		size_t	size
	)

void***** alloc5	(	size_t	n1,
		size_t	n2,
		size_t	n3,
		size_t	n4,
		size_t	n5,
		size_t	size
	)

void****** alloc6	(	size_t	n1,
		size_t	n2,
		size_t	n3,
		size_t	n4,
		size_t	n5,
		size_t	n6,
		size_t	size
	)

void free1

(

void *

p

)

void free2

(

void **

p

)

void free3

(

void ***

p

)

void free4

(

void ****

p

)

void free5

(

void *****

p

)

void free6

(

void ******

p

)

int* alloc1int

(

size_t

n1

)

int* realloc1int	(	int *	v,
		size_t	n1
	)

int** alloc2int	(	size_t	n1,
		size_t	n2
	)

int*** alloc3int	(	size_t	n1,
		size_t	n2,
		size_t	n3
	)

float* alloc1float

(

size_t

n1

)

float* realloc1float	(	float *	v,
		size_t	n1
	)

float** alloc2float	(	size_t	n1,
		size_t	n2
	)

float*** alloc3float	(	size_t	n1,
		size_t	n2,
		size_t	n3
	)

float**** alloc4float	(	size_t	n1,
		size_t	n2,
		size_t	n3,
		size_t	n4
	)

void free4float

(

float ****

p

)

float***** alloc5float	(	size_t	n1,
		size_t	n2,
		size_t	n3,
		size_t	n4,
		size_t	n5
	)

void free5float

(

float *****

p

)

float****** alloc6float	(	size_t	n1,
		size_t	n2,
		size_t	n3,
		size_t	n4,
		size_t	n5,
		size_t	n6
	)

void free6float

(

float ******

p

)

int**** alloc4int	(	size_t	n1,
		size_t	n2,
		size_t	n3,
		size_t	n4
	)

void free4int

(

int ****

p

)

int***** alloc5int	(	size_t	n1,
		size_t	n2,
		size_t	n3,
		size_t	n4,
		size_t	n5
	)

void free5int

(

int *****

p

)

unsigned short ****** alloc6ushort	(	size_t	n1,
		size_t	n2,
		size_t	n3,
		size_t	n4,
		size_t	n5,
		size_t	n6
	)

unsigned char***** alloc5uchar	(	size_t	n1,
		size_t	n2,
		size_t	n3,
		size_t	n4,
		size_t	n5
	)

void free5uchar

(

unsigned char *****

p

)

unsigned short***** alloc5ushort	(	size_t	n1,
		size_t	n2,
		size_t	n3,
		size_t	n4,
		size_t	n5
	)

void free5ushort

(

unsigned short *****

p

)

unsigned char****** alloc6uchar	(	size_t	n1,
		size_t	n2,
		size_t	n3,
		size_t	n4,
		size_t	n5,
		size_t	n6
	)

void free6uchar

(

unsigned char ******

p

)

void free6ushort

(

unsigned short ******

p

)

double* alloc1double

(

size_t

n1

)

double* realloc1double	(	double *	v,
		size_t	n1
	)

double** alloc2double	(	size_t	n1,
		size_t	n2
	)

double*** alloc3double	(	size_t	n1,
		size_t	n2,
		size_t	n3
	)

void free1int

(

int *

p

)

void free2int

(

int **

p

)

void free3int

(

int ***

p

)

void free1float

(

float *

p

)

void free2float

(

float **

p

)

void free3float

(

float ***

p

)

void free1double

(

double *

p

)

void free2double

(

double **

p

)

void free3double

(

double ***

p

)

void* bmalloc	(	int	nbpe,
		int	n1,
		int	n2
	)

void bmfree

(

void *

bm

)

void bmread	(	void *	bm,
		int	dir,
		int	k1,
		int	k2,
		int	n,
		void *	v
	)

void bmwrite	(	void *	bm,
		int	dir,
		int	k1,
		int	k2,
		int	n,
		void *	v
	)

float fsinc

(

float

x

)

double dsinc

(

double

x

)

void mksinc	(	float	d,
		int	lsinc,
		float	sinc[]
	)

void ints8r	(	int	nxin,
		float	dxin,
		float	fxin,
		float	yin[],
		float	yinl,
		float	yinr,
		int	nxout,
		float	xout[],
		float	yout[]
	)

void intt8r	(	int	ntable,
		float	table[][8],
		int	nxin,
		float	dxin,
		float	fxin,
		float	yin[],
		float	yinl,
		float	yinr,
		int	nxout,
		float	xout[],
		float	yout[]
	)

void ress8r	(	int	nxin,
		float	dxin,
		float	fxin,
		float	yin[],
		float	yinl,
		float	yinr,
		int	nxout,
		float	dxout,
		float	fxout,
		float	yout[]
	)

void shfs8r	(	float	dx,
		int	nxin,
		float	fxin,
		float	yin[],
		float	yinl,
		float	yinr,
		int	nxout,
		float	fxout,
		float	yout[]
	)

void xindex	(	int	nx,
		float	ax[],
		float	x,
		int *	index
	)

void intl2b	(	int	nxin,
		float	dxin,
		float	fxin,
		int	nyin,
		float	dyin,
		float	fyin,
		unsigned char *	zin,
		int	nxout,
		float	dxout,
		float	fxout,
		int	nyout,
		float	dyout,
		float	fyout,
		unsigned char *	zout
	)

void intlin	(	int	nin,
		float	xin[],
		float	yin[],
		float	yinl,
		float	yinr,
		int	nout,
		float	xout[],
		float	yout[]
	)

void intcub	(	int	ideriv,
		int	nin,
		float	xin[],
		float	ydin[][4],
		int	nout,
		float	xout[],
		float	yout[]
	)

void cakima	(	int	n,
		float	x[],
		float	y[],
		float	yd[][4]
	)

void cmonot	(	int	n,
		float	x[],
		float	y[],
		float	yd[][4]
	)

void csplin	(	int	n,
		float	x[],
		float	y[],
		float	yd[][4]
	)

void yxtoxy	(	int	nx,
		float	dx,
		float	fx,
		float	y[],
		int	ny,
		float	dy,
		float	fy,
		float	xylo,
		float	xyhi,
		float	x[]
	)

void bfhighpass	(	int	npoles,
		float	f3db,
		int	n,
		float	p[],
		float	q[]
	)

void bflowpass	(	int	npoles,
		float	f3db,
		int	n,
		float	p[],
		float	q[]
	)

void bfdesign	(	float	fpass,
		float	apass,
		float	fstop,
		float	astop,
		int *	npoles,
		float *	f3db
	)

void mkdiff	(	int	n,
		float	a,
		float	h,
		int	l,
		int	m,
		float	d[]
	)

void mkhdiff	(	float	h,
		int	l,
		float	d[]
	)

void holbergd1	(	float	e,
		int	n,
		float	d[]
	)

void conv	(	int	lx,
		int	ifx,
		float *	x,
		int	ly,
		int	ify,
		float *	y,
		int	lz,
		int	ifz,
		float *	z
	)

void xcor	(	int	lx,
		int	ifx,
		float *	x,
		int	ly,
		int	ify,
		float *	y,
		int	lz,
		int	ifz,
		float *	z
	)

void hilbert	(	int	n,
		float	x[],
		float	y[]
	)

void antialias	(	float	frac,
		int	phase,
		int	n,
		float	p[],
		float	q[]
	)

void* abelalloc

(

int

n

)

void abelfree

(

void *

at

)

void abel	(	void *	at,
		float	f[],
		float	g[]
	)

void* hankelalloc

(

int

nfft

)

void hankelfree

(

void *

ht

)

void hankel0	(	void *	ht,
		float	f[],
		float	h[]
	)

void hankel1	(	void *	ht,
		float	f[],
		float	h[]
	)

void hpsort	(	int	n,
		float	a[]
	)

void qksort	(	int	n,
		float	a[]
	)

void qkfind	(	int	m,
		int	n,
		float	a[]
	)

void qkisort	(	int	n,
		float	a[],
		int	i[]
	)

void qkifind	(	int	m,
		int	n,
		float	a[],
		int	i[]
	)

float quest	(	float	p,
		int	n,
		float	x[]
	)

void* questalloc	(	float	p,
		int	n,
		float	x[]
	)

float questupdate	(	void *	q,
		int	n,
		float	x[]
	)

void questfree

(

void *

q

)

void swap_short_2

(

short *

tni2

)

void swap_u_short_2

(

unsigned short *

tni2

)

void swap_int_4

(

int *

tni4

)

void swap_u_int_4

(

unsigned int *

tni4

)

void swap_long_4

(

long *

tni4

)

void swap_u_long_4

(

unsigned long *

tni4

)

void swap_float_4

(

float *

tnf4

)

void swap_double_8

(

double *

tndd8

)

int npfa

(

int

nmin

)

int npfao	(	int	nmin,
		int	nmax
	)

int npfar

(

int

nmin

)

int npfaro	(	int	nmin,
		int	nmax
	)

int npfa_d

(

int

nmin

)

int npfao_d	(	int	nmin,
		int	nmax
	)

int npfar_d

(

int

nmin

)

int npfaro_d	(	int	nmin,
		int	nmax
	)

int isamax	(	int	n,
		float *	sx,
		int	incx
	)

float sasum	(	int	n,
		float *	sx,
		int	incx
	)

void saxpy	(	int	n,
		float	sa,
		float *	sx,
		int	incx,
		float *	sy,
		int	incy
	)

void scopy	(	int	n,
		float *	sx,
		int	incx,
		float *	sy,
		int	incy
	)

float sdot	(	int	n,
		float *	sx,
		int	incx,
		float *	sy,
		int	incy
	)

float snrm2	(	int	n,
		float *	sx,
		int	incx
	)

void sscal	(	int	n,
		float	sa,
		float *	sx,
		int	incx
	)

void sswap	(	int	n,
		float *	sx,
		int	incx,
		float *	sy,
		int	incy
	)

int idamax	(	int	n,
		double *	sx,
		int	incx
	)

double dasum	(	int	n,
		double *	sx,
		int	incx
	)

void daxpy	(	int	n,
		double	sa,
		double *	sx,
		int	incx,
		double *	sy,
		int	incy
	)

void dcopy	(	int	n,
		double *	sx,
		int	incx,
		double *	sy,
		int	incy
	)

double ddot	(	int	n,
		double *	sx,
		int	incx,
		double *	sy,
		int	incy
	)

double dnrm2	(	int	n,
		double *	sx,
		int	incx
	)

void dscal	(	int	n,
		double	sa,
		double *	sx,
		int	incx
	)

void dswap	(	int	n,
		double *	sx,
		int	incx,
		double *	sy,
		int	incy
	)

void sgeco	(	float **	a,
		int	n,
		int *	ipvt,
		float *	rcond,
		float *	z
	)

void sgefa	(	float **	a,
		int	n,
		int *	ipvt,
		int *	info
	)

void sgesl	(	float **	a,
		int	n,
		int *	ipvt,
		float *	b,
		int	job
	)

void sqrdc	(	float **	x,
		int	n,
		int	p,
		float *	qraux,
		int *	jpvt,
		float *	work,
		int	job
	)

void sqrsl	(	float **	x,
		int	n,
		int	k,
		float *	qraux,
		float *	y,
		float *	qy,
		float *	qty,
		float *	b,
		float *	rsd,
		float *	xb,
		int	job,
		int *	info
	)

void sqrst	(	float **	x,
		int	n,
		int	p,
		float *	y,
		float	tol,
		float *	b,
		float *	rsd,
		int *	k,
		int *	jpvt,
		float *	qraux,
		float *	work
	)

void dgeco	(	double **	a,
		int	n,
		int *	ipvt,
		double *	rcond,
		double *	z
	)

void dgefa	(	double **	a,
		int	n,
		int *	ipvt,
		int *	info
	)

void dgesl	(	double **	a,
		int	n,
		int *	ipvt,
		double *	b,
		int	job
	)

void stoepd	(	int	n,
		double	r[],
		double	g[],
		double	f[],
		double	a[]
	)

void stoepf	(	int	n,
		float	r[],
		float	g[],
		float	f[],
		float	a[]
	)

void vanded	(	int	n,
		double	v[],
		double	b[],
		double	x[]
	)

void vandef	(	int	n,
		float	v[],
		float	b[],
		float	x[]
	)

void tridif	(	int	n,
		float	a[],
		float	b[],
		float	c[],
		float	r[],
		float	u[]
	)

void tridid	(	int	n,
		double	a[],
		double	b[],
		double	c[],
		double	r[],
		double	u[]
	)

void tripd	(	float *	d,
		float *	e,
		float *	b,
		int	n
	)

void tripp	(	int	n,
		float *	d,
		float *	e,
		float *	c,
		float *	b
	)

int mnewt	(	int	maxiter,
		float	ftol,
		float	dxtol,
		int	n,
		float *	x,
		void *	aux,
		void(*)(int n, float *x, float *f, float **dfdx, void *aux)	fdfdx
	)

void recttopolar	(	int	nx,
		float	dx,
		float	fx,
		int	ny,
		float	dy,
		float	fy,
		float **	p,
		int	na,
		float	da,
		float	fa,
		int	nr,
		float	dr,
		float	fr,
		float **	q
	)

void polartorect	(	int	na,
		float	da,
		float	fa,
		int	nr,
		float	dr,
		float	fr,
		float **	q,
		int	nx,
		float	dx,
		float	fx,
		int	ny,
		float	dy,
		float	fy,
		float **	p
	)

void rfwtva	(	int	n,
		float	z[],
		float	zmin,
		float	zmax,
		float	zbase,
		int	yzmin,
		int	yzmax,
		int	xfirst,
		int	xlast,
		int	wiggle,
		int	nbpr,
		unsigned char *	bits,
		int	endian
	)

void rfwtvaint	(	int	n,
		float	z[],
		float	zmin,
		float	zmax,
		float	zbase,
		int	yzmin,
		int	yzmax,
		int	xfirst,
		int	xlast,
		int	wiggle,
		int	nbpr,
		unsigned char *	bits,
		int	endian
	)

void scaxis	(	float	x1,
		float	x2,
		int *	nxnum,
		float *	dxnum,
		float *	fxnum
	)

int yclip	(	int	nx,
		float	dx,
		float	fx,
		float	y[],
		float	ymin,
		float	ymax,
		float	xc[],
		float	yc[]
	)

float airya

(

float

x

)

float airyb

(

float

x

)

float airyap

(

float

x

)

float airybp

(

float

x

)

float cpusec

(

void

)

float cputime

(

void

)

float wallsec

(

void

)

float walltime

(

void

)

float franuni

(

void

)

void sranuni

(

int

seed

)

float frannor

(

void

)

void srannor

(

int

seed

)

void LU_decomposition	(	int	nrows,
		float **	matrix,
		int *	index,
		float *	d
	)

void backward_substitution	(	int	nrows,
		float **	matrix,
		int *	index,
		float *	b
	)

void inverse_matrix	(	int	nrows,
		float **	matrix
	)

void inverse_matrix_multiply	(	int	nrows1,
		float **	matrix1,
		int	ncols2,
		int	nrows2,
		float **	matrix2,
		float **	out_matrix
	)

void compute_svd	(	float **	a,
		int	m,
		int	n,
		float	w[],
		float **	v
	)

void svd_backsubstitute	(	float **	u,
		float	w[],
		float **	v,
		int	m,
		int	n,
		float	b[],
		float	x[]
	)

void ricker1_wavelet	(	int	nt,
		float	dt,
		float	fpeak,
		float *	wavelet
	)

void ricker2_wavelet	(	int	hlw,
		float	dt,
		float	period,
		float	ampl,
		float	distort,
		float *	wavelet
	)

void akb_wavelet	(	int	nt,
		float	dt,
		float	fpeak,
		float *	wavelet
	)

void spike_wavelet	(	int	nt,
		int	tindex,
		float *	wavelet
	)

void unit_wavelet	(	int	nt,
		float *	wavelet
	)

void zero_wavelet	(	int	nt,
		float *	wavelet
	)

void pp1d	(	FILE *	fp,
		char *	title,
		int	lx,
		int	ifx,
		float	x[]
	)

void pplot1	(	FILE *	fp,
		char *	title,
		int	nx,
		float	ax[]
	)

char* temporary_filename

(

char *

prefix

)

void zasc	(	char *	ainput,
		char *	aoutput,
		int	nchar
	)

void zebc	(	char *	ainput,
		char *	aoutput,
		int	nchar
	)

void IBMFLT	(	float *	in,
		float *	out,
		int *	nwds,
		int *	idirec
	)

char* cwp_strdup

(

char *

str

)

void strchop	(	char *	s,
		char *	t
	)

int CWP_Exit

(

)

---