Manuals / IBM / Communications / PDAs & Smartphones

IBM Release 1.93 manual periodstep+=2 i=periodstep%4 if i!=0 periodstep=periodstep+4-i

Models: Release 1.93

1 93

Download 93 pages 50.46 Kb

Page 54

Image 54

TEMPEST uses the following code for choosing the time step (MAXDIMENSION is 3):

dt=(dx/(C/min_index*sqrt((DOUBLE)(MAXDIMENSION)))) ; period_step=ceil((DOUBLE)((2.0*PI)/(omega*dt))) ;

period_step+=2; i=period_step%4;

if (i!=0) period_step=(period_step+4-i) ;

dt=((2.0PI)/(omega(FLOAT)(period_step))) ;

printf("time step after adjustment: %e\n",dt);fflush(stdout);

fprintf(stderr,"period_step: %d\n",period_step);

Page 54

Image 54

IBM Release 1.93 manual periodstep+=2 i=periodstep%4 if i!=0 periodstep=periodstep+4-i, dt=2.0*PI/omega*FLOATperiodstep

Contents

IBM T. J. Watson Research Center Yorktown Heights, NY January 22 Updated January 23 Phil HobbsA Programmable Optimizing Electro-Magnetic Simulator Release ?88,.d88b, d8888b d8888bPage IBM T. J. Watson Research Center Yorktown Heights, NY A Programmable Optimizing Electro-Magnetic Simulator ReleaseAugust 21 Phil Hobbs POEMSChapter Appendix A. V-Antenna Optimization RunChapter Introduction 1.1. Motivation1.1.1. Philosophy 1.1.2. Structure 1.1.3. OptimizationPage FDTD Engine FIDO + CLUSTER SCRIPT or TEMPEST Optimization LoopFront-End Script POEMS.CMD Postprocessor2.2. The Front-End Script poems.cmd 2. HOW POEMS WORKS2.1. Program Organization 2.2.1. Script Operation2.3. The FDTD Engine FIDO/TEMPEST 2.4. The Postprocessor EMPOST2.5. The Visualization System VIS5D 2.6. Cluster Controlempost ordersfile resultsfile 2.6.1. Parallel Processing 3.1.1. poems Command-Line Options 3. USING POEMS3.1. Command Reference 3.1.2. GLOBAL GroupFREQ 200*THz aliases in different SUBDOMAIN statements in the parameters file NB Hosts are identified by their hostname as specified by the$HOSTNAME environment variable, rather than by IP address. This which means that the host’s predefined hostname is not usedMACRO Syntax PRINT string expression 3.1.3. WORLD Group 3.1.4. MATERIAL Group Syntax XRANGE minexpression maxexpressionParameters epsReal epsImag muReal muImag DEFINE matname= BlackGlass type=dielectric n=1.52 k=1e-3DEFINE matname=copper type=conductor conductivity=5.8e7 3.1.5. OBJECT Group xhi1,yhi1,zhi1 and xlo2,ylo2,zlo2, xhi2,yhi2,zhi2. These must beParameters matname orientation xlo xhi ylo yhi zlo zhi widthfunc phasefunc intervalxhi2 ylo2 yhi2 zlo2 zhi2 Parameters mattype xlo xhi ylo yhi zlo zhi thicknessNot implemented in this release Parameters mattype xlo xhi ylo yhi zlo zhiParameters mattype slo shi perp hfunc vfunc Parameters mattype axis lo hi maxradius insidefunc3.1.6. SOURCE Group Parameters x y z width kx ky kz Ex Ey Ez mag phase n k orientationParameters x y z kx ky kz Ex Ey Ez mag phase orientation ParametersNot implemented in this release3.1.8. OUTPUT Group 3.1.7. COMMAND GroupParameters x y z Ex Ey Ez mag phase Parameters reltolerance maxcycles mincycles timestepphase state 3.1.9. POSTPROCESS GroupParameters variable xlo xhi ylo yhi zlo zhi x y z file towards from xInside, yInside, zInside supplied plane, xmin ≤ x ≤ xmax, ymin ≤ y≤ ymax, zmin ≤ z ≤ zmaxParameters name xlo xhi ylo yhi zlo zhi file direction Parameters name xinside yinside zinside xlo xhi ylo yhi zlo zhiParameters file orientation xlo xhi ylo yhi zlo zhi phase indexn The value of coordinates can be polar or rectangular, and arrangementvariable6 variable7 variable8 variable9 xlo xhi ylo yhi zlo zhi file Parametersfile orientation xlo xhi ylo yhi zlo zhi dt frames variableParameters frames variable1 variable2 variable3 variable4 variable5 MOVIEDISSIPATION Parameters name xlo xhi ylo yhi zlo zhi file xinside yinside zinsidelevel level3.1.10. OPTIMIZE Group not be updated at postprocessing time. That means that the Syntax GUESS var1 = expression var2 =, expressionarbitrary boundary. Not implemented in this release Parameters iterations tolerance restarts 3.1.11. SCHEDULE GroupPENALTY expression Stepping is not implemented in this releaseUnfolded 2N-2 Rows 3.2. The Computational Domain3.2.1. Symmetry Actual Domain N rows Image Domain N-2 Rows Next Period3.3. OBJECTS 3.5. SOURCES 3.3.1. Perfectly-Matched Layers3.4. MATERIALS 3.5.1. POINT SOURCES3.5.2. PLANE WAVES Page 3.5.3. BEAM SOURCES 3.5.4. MODEFILE SOURCES3.6. Optimization 3.6.1. Merit Functionsd. RMS Voltage across terminal points c. Phase uniformity across a planee. Power dissipation in a load dissipation f. Reflection from an interface3.6.2. Worked Example Optimizing a V Antenna F i g u r e 2 . 1 2 Optimized V-antennaPage 3.6.3. Worked Example Doped Silica Waveguide Mode 3.6.4. Worked Example Glass Ridge Waveguide to Free Space Coupler3.7. Predefined Constants 3.7.1. Reserved Names i l l u m i n a t i o n p l a n e a f t e r3.8. Predefined Mathematical Functions 3.8.1. Arithmetic Operatorsa. Differences from C/C++ 3.8.4. ACOS 3.8.2. Logical Operators3.8.3. ABS 3.8.5. ACOSH3.8.12. COS 3.8.10. ATAN23.8.11. CEIL 3.8.13. ELINTK3.8.21. MAX 3.8.19. INTEGRAL3.8.20. LN 3.8.22. MIN3.8.28. ROUND 3.8.26. RANDOM3.8.27. ROOT1D 3.8.29. SIGN3.9. Analytical Pupil Functions 3.10. Material Parameter FunctionsA. TEMPEST and General FDTD Information 1. Startup and Steady State2. Time step dt=2.0*PI/omega*FLOATperiodstep printftime step after adjustment %e\n,dtfflushstdoutperiodstep+=2 i=periodstep%4 if i!=0 periodstep=periodstep+4-i fprintfstderr,periodstep %d\n,periodstepAppendix A. V-Antenna Optimization Run The run will produce lots of bitmaps at appropriately chosen locationsA.1. POEMS INPUT DIPOLE2I.PAR set dipoleheight=0.2*um set dipolethickness=0.30*um Fundamental dimensionsset phi=0 /* phase of excitation */ set ky0=step/lambda set kx0=0 /* direction of plane wave set loadlength=2*stepk=kload DEFINE matname= Loadtype=dielectric n=nload END WORLD SUBDOMAIN ALL XRANGE 0 xsize YRANGE 0 ysize ZRANGE 0 zsizezlo2 = Zvertex-DipoleThickness/2 zhi2= Zvertex+DipoleThickness/2 END COMMAND COMPUTE maxcycles = 25 reltolerance = END OUTPUTzhi=zsize phase=0 state=steady file=Dipole2Iexi FAN matname=PCL taper=linear taperpar=0END OPTIMIZE STORE Dipole2I.simplex MERIT diss*1e12POSTPROCESS PHASEEXxlo=0 xhi=xsize ylo=0 yhi=ysize zlo=0 zhi=zsize file=Dipole2IPhaseEx AMPLEXzhi = zsize-Tpml-step xinside=xsize/2 yinside=ysize/2 zinside=zsize/2 zlo=zsize/2+step zhi=zsize/2+step, file=Dipole2IIndexXY.bmp phase=0file=Dipole2IPZXY.bmp phase=0 palette=redblue FLUX name=MiddleFlux , /* crossing the middle z=zsize/2 xLo = Tpmlfile=Dipole2IPXZX.bmp phase=0 palette=redblue file=Dipole2IPZZX.bmp phase=0 palette=redbluefile=Dipole2IPXXY.bmp phase=0 palette=redblue file=Dipole2IAmpExZX.bmp phase=0 palette=flamefile=Dipole2IExIZX.bmp phase=0 palette=redblue SLICE orientation=xy variable=PhaseEx , xlo=0 xhi=xsizeSLICE orientation=zx variable=PhaseEx , xlo=0 xhi=xsize file=Dipole2IDissXY.bmp palette=redbluefile=Dipole2IExQZX.bmp phase=pi/2 palette=redblue END ylo=ysize/2 yhi=ysize/2, zlo=0 zhi=zsize0 100 index 1.000000 A.2. tempest Input File DIPOLE2I.PAR.INrectangle node plot hy steady 0.250000 node 0 99 0 99 0 99 dipole2ihyqWritten by 56 pages of pointsource statements omitted A.3. Postprocessor orders DIPOLE2I.ORDERS empost dipole2i.orders outfile v5dfileDIPOLE2IEYI NULLEXI ARRAY Domain for comparison e.g. modematch 0.00000 SLICE POSTPROC.16.COMPARISONDOMAIN 0.00000 SLICE A.4. Run Results DIPOLE2I.SIMPLEX 0.00000069266248611 90.849532778 0.00000152214633334 Penalty = Contracting by 0.5 in 1D 1-D shrink unsuccessful0.00000064767135 98.83191 0.000001519205 Penalty = 0.000000883423566665 77.7563866665 0.0000017351705 Penalty =Reflecting through centroid. Contracting by 0.5 in 1D Reflecting through centroid. Contracting by 0.5 in 1DReflecting through centroid. Contracting by 0.5 in 1D Reflecting through centroid. Contracting by 0.5 in 1DReflecting through centroid. Contracting by 0.5 in 1D Reflecting through centroid. Contracting by 0.5 in 1DReflecting through centroid. Contracting by 0.5 in 1D Failed to converge in 36 function evaluations1-D shrink unsuccessful Reflecting through centroid. Contracting by 0.5 in 1DReflecting through centroid. Contracting by 0.5 in 1D 0.000000572297880735 102.392565535 0.00000148853488604 Penalty =Reflecting through centroid. Contracting by 0.5 in 1D Reflecting through centroid. Contracting by 0.5 in 1DConverged in 56 function evaluations Appendix B. Configuration B.1. FDTD and TEMPESTB.1.1. tempest patches b. Understanding convergence B.1.3. Advice common to all or most FDTD programsB.1.2. tempest limitations c. Hints and KinksB.3. X Window System Configuration B.3.1. Sample X11 ConfigurationB.2. REXX B.3.2. Running Vis5D B.4. Release NotesB.4.1. Wish list B.4.2. Beta Release Limitations XY plane Z perpendicular +X right, +Y up, +Z out of screen YZ plane X perpendicular +Y right, +Z up, +X out of screenZX plane Y perpendicular +Z right, +X up, +Y out of screen Appendix C. Multiprocessing BASICSTEP 16, 18, 24, 54, 62 Boundary 7, 15, 19, 29, 31 INDEXAiry pattern Antenna 1, 39, 40, 45, 51, 52 C++ 7Periodic boundary conditions 31, 34

period_step+=2; i=period_step%4;

dt=((2.0*PI)/(omega*(FLOAT)(period_step))) ;

printf("time step after adjustment: %e\n",dt);fflush(stdout);

fprintf(stderr,"period_step: %d\n",period_step);

dt=((2.0PI)/(omega(FLOAT)(period_step))) ;