[讨论]FDTD三维程序主解释

关于FDTD三维的程序，请大神给点注解，小弟感激不尽。 <!g]q1  
#include <math.h> EYn?YiVFU  
#include <stdlib.h> w$/lq~zU  
#include <stdio.h> haBmwq(f  
#define ie 26 ; ?f+  
#define je 26 FJDC^@Ne  
#define ke 26 ]}PXN1(  
main() %SOXw8-  
{ o_`6oC"s  
 float gax[ie][je][ke],gay[ie][je][ke],gaz[ie][je][ke]; t$(#$Z,RS  
 float dx[ie][je][ke],dy[ie][je][ke],dz[ie][je][ke]; ij!d-eM/b  
 float ex[ie][je][ke],ey[ie][je][ke],ez[ie][je][ke];  NdRcA  
 float hx[ie][je][ke],hy[ie][je][ke],hz[ie][je][ke]; {N
>ju  
 int l,n,i,j,k,ic,jc,kc,nsteps,npml; @]%eL  
 float ddx,dt,T,epsz,pi,epsilon,sigma,eaf; {PX&#,_  
// float xn,xxn,xnum,xd,curl_e; gy%.+!4>v`  
 float t0,spread,pulse; e,epKtL  
 FILE *fp,*fopen();  u&#>)h  
 ic = ie/2; 7tXy3-~biz  
 jc = je/2; r) x  
 kc = ke/2; eJTU'aX*   
 ddx = 0.01; u+Ix''Fn#%  
 dt = ddx/6e8; <wk  
 epsz = 8.8e-12; ':fp|m)M  
 pi = 3.14159; `v!. ,Yr  
// printf("jc=%d \n",jc); ,IyQmN y  
 for(k=0;k<ke;k++) PkrVQH9^w  
 { WI4<2u;  
  for(j=0;j<je;j++) $/s"It  
  { !Pw*p*z  
   for(i=0;i<ie;i++) x)*Lu">  
   { aSvv(iV  
    ex[j][k] = 0.0; beBG40  
    ey[j][k] = 0.0; kW)3naUf<  
    ez[j][k] = 0.0; =4MTb_  
    dx[j][k] = 0.0; *0xL(  
    dy[j][k] = 0.0; UP^8Yhdo  
    dz[j][k] = 0.0; zI4rAsysL  
    hx[j][k] = 0.0; j{OA%G(I  
    hy[j][k] = 0.0; 4WvW11q8U  
    hz[j][k] = 0.0; KR#,6  
    gax[j][k] = 1.0; ":$4/b6  
    gay[j][k] = 1.0; &'$Bk5D@G  
    gaz[j][k] = 1.0; /Ne#{*z)hO  
   } z )'9[t  
  } h
40;Q<D  
 } ##6\~!P  
 for(k=7;k<21;k++) 3~T ~Bs  
 { xcr2|  
  gaz[ic][jc][k] = 0.0; u
VUU1@  
 } x6`mv8~9Db  
 gaz[ic][jc][kc] = 1.0; [(@K;6o  
 t0 = 20.0; R>O_2`c  
 spread = 6.0; M!'d  
 T = 0; >O?WRCB  
 nsteps = 1; ,fp+nu8,  
 while (nsteps>0) YK6LJv}  
 { i"}z9Ae~.  
  printf("nsteps-->"); (M$0'B
V0  
  scanf("%d",&nsteps); VkFh(Br<{  
  printf("%d \n",nsteps); HpEd$+Mz  
  for(n=1;n<=nsteps;n++) L]H'$~xx*  
  { ;&&<zWq3h  
   T = T+1; Q4m> 3I  
   &N^~=y^`C'  
  //calculate the Dx field _ l|%~  
  for(k=1;k<ke;k++) E xhih^[_  
  { u9"yU:1keb  
   for(j=1;j<je;j++) 20V~?xs~  
   { = g{I`u  
    for(i=1;i<ie;i++) d}4NL:=&  
    { K khuPBd2  
     dx[j][k] = dx[j][k] +0.5*(hz[j][k] - hz[j-1][k] - hy[j][k] + hy[j][k-1]); @4*:qj?  
    } K>DN6{hnV;  
   } 0~fjY^(  
  } qR'FbI  
  //calculate the Dy field ZU`9]7"87B  
  for(k=1;k<ke;k++) GK}
52,NM  
  { #"4ioTL2  
   for(j=1;j<je;j++) ,!l_  
   { aDdGhB  
    for(i=1;i<ie;i++) ri{*\LV*@  
    { \:4*h  
     dy[j][k] = dy[j][k] +0.5*(hx[j][k] - hx[j][k-1] - hz[j][k] + hz[i-1][j][k]); Cjf[]aNJe`  
    } oRCc8&  
   } +r3)\L{U  
  } h.X4x2(.  
   //calculate the Dz field oh
8:1E,I  
  for(k=1;k<ke;k++) gS'{JZu2  
  { $XhMI;h  
   for(j=1;j<je;j++) 8X,6U_>#a  
   { Ma|4nLC}  
    for(i=1;i<ie;i++) Rlewp8?LB  
    { W^[FWFUTY  
     dz[j][k] = dz[j][k] +0.5*(hy[j][k] - hy[i-1][j][k] - hx[j][k] + hx[j-1][k]); K5qCPt`'  
    } 7<xnE]jdq  
   } }n==^2  
  } XLn9NBT4K  
WiytHuUF  
  pulse = exp(-0.5*(pow((t0-T)/spread,2.0))); PT2;%=f
 
  dz[ic][jc][kc] = pulse; L(TM& ps\-  
//  dz[0][0][0] = pulse; P~trxp=k  
  printf("******************************************************\n"); 8:{id>Mm^  
  printf("pulse=%6.3f\n",dz[ic][jc][kc]); \EfX3ghPI  
  printf("******************************************************\n"); `3F#k[IR  
-'$ob~*  
   //calculate the E from D field _iJ~O1qx,w  
  for(k=1;k<ke-1;k++) 8 ??-H0P  
  { 7.h{"xOx{  
   for(j=1;j<je-1;j++) 2%pED xui  
   { '0D$C},^|8  
    for(i=1;i<ie-1;i++) 1(o\GI3:  
    { p/<DR|  
     ex[j][k] = gax[j][k] * dx[j][k]; Q~(Gll;  
     ey[j][k] = gay[j][k] * dy[j][k]; Dx-P]j)4x  
     ez[j][k] = gaz[j][k] * dz[j][k]; M-^I!C
 
    } m;dwt1'Zw  
   } ^Gk`n  
  } R])
Eg&  
XZIapT  
   //calculate the HX field j?2~6W/[  
  for(k=1;k<ke-1;k++) 3s%?)z  
  { >B{NxL3->  
   for(j=1;j<je-1;j++) }Z"iW/?"  
   { }YDi/b7  
    for(i=1;i<ie;i++) &zJI~R  
    { X+;[Gc}(W  
     hx[j][k] = hx[j][k] +0.5*(ey[j][k+1] - ey[j][k] - ez[j+1][k] + ez[j][k]); [7@g*!+d  
    } %
Ln`c.C  
   } L0wT:x*  
  } q}P@}TE  
Wzl/ @CPM  
   //calculate the HY field @)U.Dbm  
  for(k=1;k<ke-1;k++) tY"eoPme  
  { *C*ZmC5  
   for(j=1;j<je;j++) GY]P(NU  
   { ZQfxlzj+X  
    for(i=1;i<ie-1;i++) (GmBv  
    { 6j6CA?|  
     hy[j][k] = hy[j][k] +0.5*(ez[i+1][j][k] - ez[j][k] - ex[j][k+1] + ex[j][k]); kZw"a*6  
    } x?wvS]EBg  
   } wm`<+K  
  } 8S1@,O,  
   //calculate the HZ field ;R*-cm  
  for(k=1;k<ke;k++) (TT=i  
  { D;Qx9^.  
   for(j=1;j<je-1;j++) { ptdOrN  
   { |owr?tC  
    for(i=1;i<ie-1;i++) KO

y{?  
    { Fv~lasW[  
     hz[j][k] = hz[j][k] +0.5*(ex[j+1][k] - ex[j][k] - ey[i+1][j][k] + ey[j][k]); .
==D?#bn  
    } .C^P6S2oJ  
   } aS [[ AL  
  } qi,) l*?f  
  }//end for [{7#IZL  
  printf("EZ \n"); G WIsT\J  
  for(k=1;k<ke;k++) t$!zgUJ  
  { dOeM0_o  
   printf("%2d  \n",k); !zE{`Ha~  
   printf("\n"); arL>{mj  
   .. ~~/
,2^  
RAO+<m  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

这个写得挺清楚的不知道楼主要什么注解？是不是得具体些。

我想要一些关于程序的解释，如 *)r_Y|vg  
for(k=7;k<21;k++) )cN=/i  
 { V13^SVM  
  gaz[ic][jc][k] = 0.0; J{Ij  
 } :qXREF@h  
 gaz[ic][jc][kc] = 1.0; rDaiA
x&  
为什么要从7开始。 还有为什么要引入D等等。 b0f6?s  
不知道你是否懂，能不能给点提示。感激不尽。 f:q2JgX  

从7开始应该是和边界条件有关系吧，D好像是公式中的一个系数还是什么的，具体忘记了，你可以参考一下其它人分享的资料，里面有解释的。