buffers.h.txt

FDTD的并行算法

#define c 2.99795637e8
#define eps0 8.854e-12
#define mu0 1.2566306e-6

#define dx 1e-3
#define dy 1e-3
#define ny 100
#define nmax 500

#define w_pulse 15e-12

float ** exbufl1a(int nx)
{
	int i;
	float **a;
	a=(float**)calloc(nx+2,sizeof(float));
	a[0]=(float*)calloc((nx+2)*(2),sizeof(float));
	for(i=1;i<=nx+2;i++)
	{	
		a[i]=a[0]+i*(2);
	}
	return(a);
}

float ** exbufl2a(int nx)
{
	int i;
	float **a;
	a=(float**)calloc(nx+2,sizeof(float));
	a[0]=(float*)calloc((nx+2)*(2),sizeof(float));
	for(i=1;i<=nx+2;i++)
	{	
		a[i]=a[0]+i*(2);
	}
	return(a);
}

float ** exbufr1a(int nx)
{
	int i;
	float **a;
	a=(float**)calloc(nx+2,sizeof(float));
	a[0]=(float*)calloc((nx+2)*(2),sizeof(float));
	for(i=1;i<=nx+2;i++)
	{	
		a[i]=a[0]+i*(2);
	}
	return(a);
}

float ** exbufr2a(int nx)
{
	int i;
	float **a;
	a=(float**)calloc(nx+2,sizeof(float));
	a[0]=(float*)calloc((nx+2)*(2),sizeof(float));
	for(i=1;i<=nx+2;i++)
	{	
		a[i]=a[0]+i*(2);
	}
	return(a);
}

float ** eybufl1a(int nx)
{
	int i;
	float **a;
	a=(float**)calloc(2,sizeof(float));
	a[0]=(float*)calloc((2)*(ny+2),sizeof(float));
	for(i=1;i<=2;i++)
	{	
		a[i]=a[0]+i*(ny+2);
	}
	return(a);
}

float ** eybufl2a(int nx)
{
	int i;
	float **a;
	a=(float**)calloc(2,sizeof(float));
	a[0]=(float*)calloc((2)*(ny+2),sizeof(float));
	for(i=1;i<=2;i++)
	{	
		a[i]=a[0]+i*(ny+2);
	}
	return(a);
}

float ** eybufr1a(int nx)
{
	int i;
	float **a;
	a=(float**)calloc(2,sizeof(float));
	a[0]=(float*)calloc((2)*(ny+2),sizeof(float));
	for(i=1;i<=2;i++)
	{	
		a[i]=a[0]+i*(ny+2);
	}
	return(a);
}

float ** eybufr2a(int nx)
{
	int i;
	float **a;
	a=(float**)calloc(2,sizeof(float));
	a[0]=(float*)calloc((2)*(ny+2),sizeof(float));
	for(i=1;i<=2;i++)
	{	
		a[i]=a[0]+i*(ny+2);
	}
	return(a);
}