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

关于FDTD三维的程序，请大神给点注解，小弟感激不尽。 *-E'$  
#include <math.h> \9}DAM_  
#include <stdlib.h> o,/wE  
#include <stdio.h> e`K{  
#define ie 26 o76{;Bl\O  
#define je 26 :xY9eq=  
#define ke 26 ghTue*A  
main() E+$D$a  
{ e^N}(Kpy  
 float gax[ie][je][ke],gay[ie][je][ke],gaz[ie][je][ke]; iAt&927  
 float dx[ie][je][ke],dy[ie][je][ke],dz[ie][je][ke]; CbOCL~ "  
 float ex[ie][je][ke],ey[ie][je][ke],ez[ie][je][ke]; ss`Sl$  
 float hx[ie][je][ke],hy[ie][je][ke],hz[ie][je][ke]; )kYDN_W  
 int l,n,i,j,k,ic,jc,kc,nsteps,npml; fGWXUJ
 
 float ddx,dt,T,epsz,pi,epsilon,sigma,eaf; a&7uRR26  
// float xn,xxn,xnum,xd,curl_e; nW oh(a  
 float t0,spread,pulse; {/qQ=$t  
 FILE *fp,*fopen(); StiWa<"c  
 ic = ie/2; eAjsMED  
 jc = je/2; ~"8r=8|  
 kc = ke/2; :BB=E'293  
 ddx = 0.01; mrig5{  
 dt = ddx/6e8; jVqpokWH  
 epsz = 8.8e-12; Ml'lZ)  
 pi = 3.14159;
dZ4c!3'F  
// printf("jc=%d \n",jc); msQ?V&+<  
 for(k=0;k<ke;k++) iV!V!0- @  
 { YdN]Tqc  
  for(j=0;j<je;j++) Zg|l:^E  
  { {vQ:4O!:  
   for(i=0;i<ie;i++) }L=Qp
=4  
   { 7;pQ'FmZJ  
    ex[j][k] = 0.0; QZuKM'D+  
    ey[j][k] = 0.0; `A
-  
    ez[j][k] = 0.0; 1*.*\4xo  
    dx[j][k] = 0.0; PZI6{KOis  
    dy[j][k] = 0.0; ` }B,w-,io  
    dz[j][k] = 0.0; (k_9<Yb3  
    hx[j][k] = 0.0; TIK'A<  
    hy[j][k] = 0.0; hS&.-5v  
    hz[j][k] = 0.0; AS'+p%(  
    gax[j][k] = 1.0; P(
a}OlG  
    gay[j][k] = 1.0; `$q
0fTz  
    gaz[j][k] = 1.0; 3$.#\*s_4  
   } RiFUa $  
  } T1'\!6_5  
 } p1W6s0L  
 for(k=7;k<21;k++) Y~?Z'uR  
 { $EzWUt  
  gaz[ic][jc][k] = 0.0; )JJF}m=  
 } 4um^7Ns)7  
 gaz[ic][jc][kc] = 1.0; %/qwqo`Q  
 t0 = 20.0; ~wO-Hgd  
 spread = 6.0; EV7lgKM^  
 T = 0; Ys|SacWC  
 nsteps = 1; kI2+&  
 while (nsteps>0) YOUB%N9+  
 { C~C`K%7  
  printf("nsteps-->"); +O.qYX  
  scanf("%d",&nsteps); |kTq &^$  
  printf("%d \n",nsteps); u4bVp+  
  for(n=1;n<=nsteps;n++) HJP~ lg  
  { T\bpeky~  
   T = T+1; =^\?{oV  
   JpxQS~VX  
  //calculate the Dx field t0Jqr)9}6  
  for(k=1;k<ke;k++) >8b%*f8R  
  { mI]gDL1  
   for(j=1;j<je;j++) 'sE["eC  
   { ;>S|?M4GZ  
    for(i=1;i<ie;i++) *||Q_tlz  
    { 9ExI,  
     dx[j][k] = dx[j][k] +0.5*(hz[j][k] - hz[j-1][k] - hy[j][k] + hy[j][k-1]); &I%E8E  
    } \jmZt*c  
   }
R1eWPtWs  
  }  1'F!C  
  //calculate the Dy field IAmZ_2  
  for(k=1;k<ke;k++) WW2hwB(  
  { )lz~Rt;1i  
   for(j=1;j<je;j++) 7[!dm_  
   { lJdYR'/Wd  
    for(i=1;i<ie;i++) U3`?Z`i(  
    { `uZMln @  
     dy[j][k] = dy[j][k] +0.5*(hx[j][k] - hx[j][k-1] - hz[j][k] + hz[i-1][j][k]); $15H_X*!  
    } R[)bGl6#  
   } ?%Ww3cU+J  
  } UEhFId  
   //calculate the Dz field c{KJNH%7  
  for(k=1;k<ke;k++) Y=g]\%-PB  
  { 6 jm@`pYbE  
   for(j=1;j<je;j++) !l Egta[Ql  
   { e`q*'u1?  
    for(i=1;i<ie;i++) gNc;P[  
    { "z;R"sv\  
     dz[j][k] = dz[j][k] +0.5*(hy[j][k] - hy[i-1][j][k] - hx[j][k] + hx[j-1][k]); gVI`&W__,  
    } t\TxK7i  
   } _N)&<'lB<  
  } `w`F-ke]I  
x+;y0`oL  
  pulse = exp(-0.5*(pow((t0-T)/spread,2.0))); +l.LwA
 
  dz[ic][jc][kc] = pulse; WglpWp)  
//  dz[0][0][0] = pulse; 08D:2 z1z  
  printf("******************************************************\n"); -ZyY95E<  
  printf("pulse=%6.3f\n",dz[ic][jc][kc]);
m l@%H  
  printf("******************************************************\n"); 8FZC0j.^DH  
MLg{Y?@  
   //calculate the E from D field VJFFH\!`  
  for(k=1;k<ke-1;k++) z Tz_"NI  
  { }/,Rp/+7]  
   for(j=1;j<je-1;j++) R!lug;u#  
   { jzGK(%sw"  
    for(i=1;i<ie-1;i++) 6}!1a?X  
    { nMfR<%r  
     ex[j][k] = gax[j][k] * dx[j][k]; }6<5mq)%  
     ey[j][k] = gay[j][k] * dy[j][k]; 4*Gv0#dga  
     ez[j][k] = gaz[j][k] * dz[j][k]; 41s\^'^&  
    } T93st<F=R  
   } &[_@f
#  
  } 4d}=g]P  
2|d^#8)ZC  
   //calculate the HX field &q9=0So4\  
  for(k=1;k<ke-1;k++) S'kgpF"bm  
  { }f14# y;  
   for(j=1;j<je-1;j++) 9V[}#(f$  
   { q\|RI;W  
    for(i=1;i<ie;i++) F",TP,X  
    {
0a^bAEP  
     hx[j][k] = hx[j][k] +0.5*(ey[j][k+1] - ey[j][k] - ez[j+1][k] + ez[j][k]); 4>LaA7)v  
    } N)AlQ'Lwx  
   } VZ=:`)  
  } LHXR7Fjc  
G2?#MO  
   //calculate the HY field VJT /9O)Z|  
  for(k=1;k<ke-1;k++) nE*S3  
  { Xt$P!~Lu  
   for(j=1;j<je;j++) hITYBPqRO  
   { R! On  
    for(i=1;i<ie-1;i++) o 9/,@Ri\5  
    { \AH5zdK  
     hy[j][k] = hy[j][k] +0.5*(ez[i+1][j][k] - ez[j][k] - ex[j][k+1] + ex[j][k]); ]8DTk!  
    } FfrC/"N  
   } _DT,iF*6  
  } c$^v~lQS  
   //calculate the HZ field o3GkTn O  
  for(k=1;k<ke;k++) BRskxyL&,  
  { 6M_:D  
   for(j=1;j<je-1;j++) bgKC^Q/F  
   {  -gS9I^  
    for(i=1;i<ie-1;i++) K/ &`  
    { II|;_j  
     hz[j][k] = hz[j][k] +0.5*(ex[j+1][k] - ex[j][k] - ey[i+1][j][k] + ey[j][k]); UcOP 0_/  
    } "~._G5i.  
   } k ,+,,W  
  } .%e>>U>F  
  }//end for 6 U[VoUU   
  printf("EZ \n"); Mj6,VD9L  
  for(k=1;k<ke;k++) !4=_l6kg~+  
  { zrG  
   printf("%2d  \n",k); Z`
MpH  
   printf("\n"); g?Nk-cg
  
   .. #asi%&3pP  
( $2M
"n  

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

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

我想要一些关于程序的解释，如 AA^K/y  
for(k=7;k<21;k++) M| Gl&   
 { )-[$m%  
  gaz[ic][jc][k] = 0.0; 9yTdbpY  
 } AYY(<b  
 gaz[ic][jc][kc] = 1.0; QObVJg,GD  
为什么要从7开始。 还有为什么要引入D等等。 Xn=yC Pi  
不知道你是否懂，能不能给点提示。感激不尽。 vv,<#4d
 

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