并行计算算例.docVIP

并行计算算例 题目： 一维（X）电磁场FDTD模拟，高斯源，一阶Mur吸收边界，阻塞消息通信方式（SSend) 同步通信模式采用c*delta(t)/delta(x)=0.5简化起动进程数必须大于等于2Module parameter_mod* 自定义模型，此处不需要End Module parameter_mod * main program******************************************************************************Program FDTD1D_MPIuse parameter_modimplicit noneinclude mpif.h*定义常量* c:光速，dx:离散间隔， N:x长度，dt:时间间隔（满足稳定条件）， T：时间长度，Nt：时间点数，default：默认参数开关。double precision??? :: c =3.0e8double precision??? :: dx=0.005double precision??? :: N=1.0integer?????? :: Nx0double precision??? :: dtdouble precision??? :: T=1.0e-9integer?????? :: Ntinteger?????? :: default=0 *Ey:用于保存Ey结果(不必要），node_Ey,node_Hz：该进程中的Ey,Hz分量，node_Ey_right,node_Hzleft：进程间电磁场传递增，*node_Eyoldleft：最后一个进程（numberofprocessors-1)中的node_Ey(nx-1)分量，用于计算吸收边界*?node_Eyoldright：第一进程（0，即主进程）中的node_Ey分量，用于计算吸收边界double precision, allocatable :: Ey(:) double precision, allocatable :: node_Ey(:)double precision, allocatable :: node_Hz(:) double precision???? node_Ey_rightdouble precision???? node_Hz_left double precision???? node_Eyoldleftdouble precision???? node_Eyoldright* nx:该进程中的结点数，source_id：源所在进程ID，ksn：源所在进程中的结点号，* n_average:平均每个进程结点数，n_remain：结点余数（分配给从进程）。* displs recvcounts：用于数据收集的控制数组。displs（i）:第i号进程数据偏移量，recvcounts（i）:第i号进程数据量。* num_Ey：该进程Ey总量。integer nxinteger source_id, ksninteger n_average, n_remaininteger, allocatable::displs(:), recvcounts(:)integer???? :: num_Ey* i,nstep：循环变量* nxtemp,kstemp,temp：其它蓝临时变量* flag：用于寻找源的位置integer i, nstepinteger nxtemp, kstemp, templogical flag* w:源脉冲宽度，t0：脉冲峰值出现位置（点号），ks：源的位置（X，点号），pulse：源值。double precision :: wdouble precision :: t0double precision :: ksdouble precision :: pulse* 本进程ID，ierr错误号, numberofprocessors：总进程数，status：消息接收状态，main_id:主进程ID（0），* tagnx,tagEy,tagHz,TagNNt:消息信封标志。integer processorID, ierrinteger numberofprocessorsinteger status(MPI_Status_size)integer,parameter :: main_id=0integer,parameter :: tagnx=1, tagEy=2, tagHz=3, tagNNt=4* 并行初始化，生成通信域（子）(MPI_C