地震波克西霍夫叠后偏移程序源代码

//#define MPICH_SKIP_MPICXX
//#include “mpi.h“
#include 
#include 
#include 
#include 
#include 
//#include 
using namespace std;

//#define INFINITY 1.e+20
#define FMM_OUT ‘\0‘
#define FMM_IN ‘i‘
#define FMM_FRONT ‘f‘
#define CVTIME（p1p2p3） （slow0*hypot（hypot（（p1）*dd1-s1（p2）*dd2-s2）（p3）*dd3-s3））

static float **x **xn **x1;

void fastmarch （int orderfloat s1 float s2 float s3int b1 int b2 int b3 
		        int nz int ny int nxfloat dz float dy float dx float slow0
				float *slow float *ttime）;

void update （int order int i1 int i2 int i3 float* tjunsigned char* mj float s）;

void cvupdate （double cvtimefloat* tjunsigned char* mj）;

void HL（unsigned char *pint n）;

float ibm2