vectors.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 sq(float ss)
{
	ss=ss*ss;
	return(ss);
}


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

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

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

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

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

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

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

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

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

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

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