Los Alamos National Laboratory Poisson Superfish Program Pandira written by Ron F. Holsinger in collaboration with Klaus Halbach The original Poisson Superfish codes were developed by Ron F. Holsinger in collaboration with Klaus Halbach. These programs are provided as a service to the accelerator community by the Los Alamos Accelerator Code Group (LAACG). (c) Copyright 1985-2005, by the Regents of the University of California. This software was produced under U. S. Government contract W-7405-ENG-36 by Los Alamos National Laboratory, which is operated by the University of California for the U. S. Department of Energy. Neither the Government nor the University makes any warranty, express or implied, or assumes any liability or responsibility for its use, or represents that use of this software would not infringe privately owned rights. Unpublished - rights reserved under Copyright Laws of the United States. Program Pandira 7.17 released 1-13-2006 Starting from file N:\FE\POISSON\MESH_DIEL_NO_RG.T35. Program file: C:\LANL\PANDIRA.EXE SF.INI file: C:\LANL\SF.INI 7-17-2012 14:55:56 Memory used for MESH_DIEL_NO_RG.T35 arrays: 39.006 M Memory used for other Pandira arrays: 55.068 M Problem description: based on FILENAME Problem file: N:\FE\POISSON\MESH_DIEL_NO_RG.AM 8-06-2012 23:40:18 Problem file length: 857 bytes Originating program: Automesh Problem type: Electrostatic Coordinates and lengths have dimensions of centimeters. Problem constants and variables. Letter A in the code column indicates a value supplied in the Automesh input file. Variable Code Value Description ANGLE 0.0 Extent of arc for interpolating potential ANGLZ 0.0 Initial point on arc for interpolation CLIGHT 2.997924580E+10 Exact speed of light in cm/sec CONV A 1.00000000 Length conversion (number of units per cm) DIAGDLL 0 If 1, DLL writes diagnostics to DiagDLL.txt DX1 3.012048193E-04 First X mesh interval (at XMIN) DXMIN 3.000000000E-04 Minimum X mesh interval (found by Automesh) DYMIN 2.999850007E-04 Minimum Y mesh interval (found by Automesh) EPS0 8.854187818E-12 Permittivity of free space EPSILA 5.000000000E-07 Convergence PARAMETER: air, interface points EPSILI 5.000000000E-07 Convergence parameter for iron points EPSO 1.000000000E-05 Convergence parameter in mesh optimization FIXEPS 9.00000000 Permittivity when MODE=-1 FIXGAM 4.000000000E-03 Reluctivity when MODE=-1 FMU0 1.256637061E-06 Permeability of free space ICAL 0 Type of formula for currents near coils ICCC 1 1 in Poisson/Pandira problems (all real arrays) ICYLIN A 0 0 for X,Y problems, 1 for Z,R problems ICYSEN 0 If zero, do not print boundary integrals IENERGY 0 If 1, calculate stored energy IHDL 100000 Cycles between H*dl integrals IMAX 86 KMAX+2 IPERM 0 0 = real currents, 1 = no real currents IPIVOT 1 Pivoting in matrix inversion routines IPRFQ 0 Print frequency during Poisson iterations IRMAX 25 Used in optimization of RHOXY ISKIP 1 Cycles between calculation of reluctance ITOT 573620 (KMAX+2)*(LMAX+2) IVERG 10 Number of cycles between convergence tests KBZERO 1 Logical K coordinate where B = BDES KMAX 84 Number of horizontal logical mesh points KMIN 1 Lower K bound for computing field KPROB A 0 Problem type indicator (Poisson) KTOP 84 Upper K bound for computing field KTYPE 2 Symmetry indicator for harmonic analysis LAST35 1 Code for last program to update T35 file LBZERO 1 Logical L coordinate where B = BDES LCYCLE 87 Iteration number LMAX 6668 Number of vertical logical mesh points LMIN 1 Lower L bound for computing field LTOP 1 Upper L bound for computing field MAP 1 Conformal mapping parameter MAXCY 20 Maximum number of cycles (-1: use default) MAXPPR 13502 Maximum points per region MODE A -1 Material property indicator NAIR 195888 Number of air points NAMAX 0 Number of points for recalculating couplings NBND 0 Number of Dirichlet boundary points NBSLF A 1 Left-side boundary condition NBSLO 1 Lower boundary condition NBSRT A 1 Right-side boundary condition NBSUP 0 Upper boundary condition NEGAT 0 Zero-area triangle indicator NFE 364056 Number of iron points NGMAX 359805 Number of points for recalculating gammas NGSAM 359556 Number of points for reluctivity when NM6=NM1 NINTER 84 Number of interface points NMATR 1 Number of material records in T35 file NOTE 1 Determines relaxation order NPBOUND 22671 Total number of boundary points in the mesh NPINP 560112 Total points in problem NPONTS 546448 Number of unknown relaxation points NPTC 0 Number of arc points for interpolation NREG 5 Number of regions NSEG 18 Number of boundary segments NSPL 84 Number of special-potential points NTERM 0 Number of coefficients in harmonic analysis NWMAX 361972 Number of points for recalculating couplings OMEGAM 1.000000000E-03 Used in optimization of RHOXY OMEGAP 1.000000000E-03 Optimization parameter for RHOAIR in Poisson PI 3.14159265 The number pi to machine precision RATIO 0.0 BZERO/XJFACT for air solution RHOAIR 1.90000000 Over-relaxation factor: air, interface points RHOFE 1.00000000 Over-relaxation factor for iron points RHOGAM 8.000000000E-02 Under-relaxation parameter for reluctivity RHOPT1 1.90000000 If = RHOAIR, causes optimization of RHOAIR RHOXY 1.60000000 Over-relaxation factor for mesh optimization RINT 0.0 Radius of arc for interpolating potential RNORM 0.0 Aperture radius in harmonic analysis RZERO 1.00000000 Scaling factor for conformal transformation SNEGG 0.0 Total negative current at generation SPOSG 0.0 Total positive current at generation STACK 1.00000000 Stacking factor STOTG 0.0 Total current in problem at generation TNEGC 0.0 Negative current after conformal transformation TPOSC 0.0 Positive current after conformal transformation TRIAVG 4.517846397E-08 Average area of all triangles TRIMAX 5.957084330E-08 Area of the largest positive-area triangle TRIMIN 2.711171151E-08 Area of the smallest positive-area triangle VOLUME 0.0 Cavity volume (cylindrical symmetry only) XJFACT A 0.0 Factor multiplying currents XJFEND 1.00000000 Ending XJFACT for a range of multipliers XJSTEPS 0 Number of steps from XJFACT to XJFEND XJTOL 5.000000000E-07 Tolerance on XJFACT XMAXF 0.0 Upper X bound for computing field XMAXG 2.500000000E-02 Upper X bound of the problem geometry XMINF 0.0 Lower X bound for computing field XMING 0.0 Lower X bound of the problem geometry XORG 0.0 X coordinate center of harmonic analysis circle XYAREA 5.000000000E-02 Total cross sectional area YMAXF 0.0 Upper Y bound for computing field YMAXG 2.00000000 Upper Y bound of the problem geometry YMINF 0.0 Lower Y bound for computing field YMING 0.0 Lower Y bound of the problem geometry YORG 0.0 Y coordinate center of harmonic analysis circle Material properties: Material 1 corresponds to air regions. Material 3 uses the fixed permittivity Epsilon = 3.7 from MT namelist ID 3. Material 4 uses the fixed permittivity Epsilon = 9 from MT namelist ID 4. Kmax = 84, Lmax = 6668, Total number of mesh points = 573620 ---------------------------------------------------------------- Writing temporary data arrays to memory. Memory used: 371.893 M ---------------------------------------------------------------- Solution matrix minimum = 7.4000E+03 maximum = 7.4000E+03 Convergence rate in iron = 1.000 Residual = 1.0000E+00 Elapsed time = 22.19 sec End of cycle 1 ---------------------------------------------------------------- Solution converged at the end of cycle 1. Elapsed time = 22.19 sec XJFACT = 0.0 Field interpolation at cycle 1. K L X Y Ex Ey |E| V Chi^2/d.f. Fn (cm) (cm) (V/cm) (V/cm) (V/cm) (V) 1 1 0.00000 0.00000 0.000000E+00 0.000000E+00 0.000000E+00 7.400000E+03 1.60E-31 20 2 1 3.012048E-04 0.00000 2.872513E-09 0.000000E+00 2.872513E-09 7.400000E+03 1.49E-31 20 3 1 6.024096E-04 0.00000 5.026898E-09 0.000000E+00 5.026898E-09 7.400000E+03 1.58E-31 20 4 1 9.036145E-04 0.00000 2.793968E-09 0.000000E+00 2.793968E-09 7.400000E+03 4.94E-31 19 5 1 1.204819E-03 0.00000 1.862645E-09 0.000000E+00 1.862645E-09 7.400000E+03 4.29E-31 19 6 1 1.506024E-03 0.00000 -5.587935E-09 0.000000E+00 5.587935E-09 7.400000E+03 4.08E-31 19 7 1 1.807229E-03 0.00000 -1.234002E-08 0.000000E+00 1.234002E-08 7.400000E+03 3.31E-31 19 8 1 2.108434E-03 0.00000 -9.778887E-09 0.000000E+00 9.778887E-09 7.400000E+03 4.19E-31 19 9 1 2.409639E-03 0.00000 -1.001172E-08 0.000000E+00 1.001172E-08 7.400000E+03 3.90E-31 19 10 1 2.710843E-03 0.00000 -6.984919E-09 0.000000E+00 6.984919E-09 7.400000E+03 4.84E-31 19 11 1 3.012048E-03 0.00000 -5.355105E-09 0.000000E+00 5.355105E-09 7.400000E+03 5.03E-31 19 12 1 3.313253E-03 0.00000 -3.259629E-09 0.000000E+00 3.259629E-09 7.400000E+03 5.38E-31 19 13 1 3.614458E-03 0.00000 -1.629815E-09 0.000000E+00 1.629815E-09 7.400000E+03 5.58E-31 19 14 1 3.915663E-03 0.00000 -1.629815E-09 0.000000E+00 1.629815E-09 7.400000E+03 5.67E-31 19 15 1 4.216867E-03 0.00000 -2.095476E-09 0.000000E+00 2.095476E-09 7.400000E+03 6.07E-31 19 16 1 4.518072E-03 0.00000 -6.984919E-10 0.000000E+00 6.984919E-10 7.400000E+03 8.74E-31 19 17 1 4.819277E-03 0.00000 9.313226E-10 0.000000E+00 9.313226E-10 7.400000E+03 8.01E-31 19 18 1 5.120482E-03 0.00000 5.820766E-09 0.000000E+00 5.820766E-09 7.400000E+03 5.65E-31 19 19 1 5.421687E-03 0.00000 1.001172E-08 0.000000E+00 1.001172E-08 7.400000E+03 3.50E-31 19 20 1 5.722892E-03 0.00000 1.257285E-08 0.000000E+00 1.257285E-08 7.400000E+03 2.34E-31 19 21 1 6.024096E-03 0.00000 2.002344E-08 0.000000E+00 2.002344E-08 7.400000E+03 3.05E-31 19 22 1 6.325301E-03 0.00000 1.979060E-08 0.000000E+00 1.979060E-08 7.400000E+03 2.56E-31 19 23 1 6.626506E-03 0.00000 2.165325E-08 0.000000E+00 2.165325E-08 7.400000E+03 2.32E-31 19 24 1 6.927711E-03 0.00000 1.839362E-08 0.000000E+00 1.839362E-08 7.400000E+03 4.21E-31 19 25 1 7.228916E-03 0.00000 1.583248E-08 0.000000E+00 1.583248E-08 7.400000E+03 2.67E-31 19 26 1 7.530120E-03 0.00000 1.466833E-08 0.000000E+00 1.466833E-08 7.400000E+03 1.79E-31 19 27 1 7.831325E-03 0.00000 1.792796E-08 0.000000E+00 1.792796E-08 7.400000E+03 2.97E-31 19 28 1 8.132530E-03 0.00000 2.072193E-08 0.000000E+00 2.072193E-08 7.400000E+03 2.99E-31 19 29 1 8.433735E-03 0.00000 2.072193E-08 0.000000E+00 2.072193E-08 7.400000E+03 2.86E-31 19 30 1 8.734940E-03 0.00000 1.979060E-08 0.000000E+00 1.979060E-08 7.400000E+03 3.19E-31 19 31 1 9.036145E-03 0.00000 2.165325E-08 0.000000E+00 2.165325E-08 7.400000E+03 3.27E-31 19 32 1 9.337349E-03 0.00000 2.072193E-08 0.000000E+00 2.072193E-08 7.400000E+03 2.51E-31 19 33 1 9.638554E-03 0.00000 2.188608E-08 0.000000E+00 2.188608E-08 7.400000E+03 1.35E-31 19 34 1 9.939759E-03 0.00000 1.979060E-08 0.000000E+00 1.979060E-08 7.400000E+03 1.90E-31 19 35 1 1.024096E-02 0.00000 2.095476E-08 0.000000E+00 2.095476E-08 7.400000E+03 6.13E-32 19 36 1 1.054217E-02 0.00000 1.839362E-08 0.000000E+00 1.839362E-08 7.400000E+03 8.71E-32 19 37 1 1.084337E-02 0.00000 2.072193E-08 0.000000E+00 2.072193E-08 7.400000E+03 2.29E-31 19 38 1 1.114458E-02 0.00000 2.328306E-08 0.000000E+00 2.328306E-08 7.400000E+03 2.71E-31 19 39 1 1.144578E-02 0.00000 2.025627E-08 0.000000E+00 2.025627E-08 7.400000E+03 2.42E-31 19 40 1 1.174699E-02 0.00000 1.513399E-08 0.000000E+00 1.513399E-08 7.400000E+03 2.92E-31 19 41 1 1.204819E-02 0.00000 1.140870E-08 0.000000E+00 1.140870E-08 7.400000E+03 2.52E-31 19 42 1 1.234940E-02 0.00000 1.071021E-08 0.000000E+00 1.071021E-08 7.400000E+03 1.80E-31 19 43 1 1.265060E-02 0.00000 1.303852E-08 0.000000E+00 1.303852E-08 7.400000E+03 2.67E-31 19 44 1 1.295181E-02 0.00000 1.606531E-08 0.000000E+00 1.606531E-08 7.400000E+03 4.51E-31 19 45 1 1.325301E-02 0.00000 2.072193E-08 0.000000E+00 2.072193E-08 7.400000E+03 4.87E-31 19 46 1 1.355422E-02 0.00000 2.584420E-08 0.000000E+00 2.584420E-08 7.400000E+03 5.24E-31 19 47 1 1.385542E-02 0.00000 2.677552E-08 0.000000E+00 2.677552E-08 7.400000E+03 4.30E-31 19 48 1 1.415663E-02 0.00000 2.025627E-08 0.000000E+00 2.025627E-08 7.400000E+03 3.25E-31 19 49 1 1.445783E-02 0.00000 1.280569E-08 0.000000E+00 1.280569E-08 7.400000E+03 1.49E-31 19 50 1 1.475904E-02 0.00000 7.916242E-09 0.000000E+00 7.916242E-09 7.400000E+03 1.02E-31 19 51 1 1.506024E-02 0.00000 4.423782E-09 0.000000E+00 4.423782E-09 7.400000E+03 1.03E-31 19 52 1 1.536145E-02 0.00000 1.629815E-09 0.000000E+00 1.629815E-09 7.400000E+03 2.16E-31 19 53 1 1.566265E-02 0.00000 9.313226E-10 0.000000E+00 9.313226E-10 7.400000E+03 8.80E-32 19 54 1 1.596386E-02 0.00000 0.000000E+00 0.000000E+00 0.000000E+00 7.400000E+03 1.03E-31 19 55 1 1.626506E-02 0.00000 -2.095476E-09 0.000000E+00 2.095476E-09 7.400000E+03 1.53E-31 19 56 1 1.656627E-02 0.00000 2.095476E-09 0.000000E+00 2.095476E-09 7.400000E+03 1.16E-31 19 57 1 1.686747E-02 0.00000 7.217750E-09 0.000000E+00 7.217750E-09 7.400000E+03 1.39E-31 19 58 1 1.716867E-02 0.00000 1.117587E-08 0.000000E+00 1.117587E-08 7.400000E+03 1.52E-31 19 59 1 1.746988E-02 0.00000 1.490116E-08 0.000000E+00 1.490116E-08 7.400000E+03 1.44E-31 19 60 1 1.777108E-02 0.00000 2.095476E-08 0.000000E+00 2.095476E-08 7.400000E+03 2.01E-31 19 61 1 1.807229E-02 0.00000 2.817251E-08 0.000000E+00 2.817251E-08 7.400000E+03 2.94E-31 19 62 1 1.837349E-02 0.00000 3.376044E-08 0.000000E+00 3.376044E-08 7.400000E+03 2.67E-31 19 63 1 1.867470E-02 0.00000 3.748573E-08 0.000000E+00 3.748573E-08 7.400000E+03 1.75E-31 19 64 1 1.897590E-02 0.00000 3.725290E-08 0.000000E+00 3.725290E-08 7.400000E+03 1.95E-31 19 65 1 1.927711E-02 0.00000 3.282912E-08 0.000000E+00 3.282912E-08 7.400000E+03 2.08E-31 19 66 1 1.957831E-02 0.00000 2.980232E-08 0.000000E+00 2.980232E-08 7.400000E+03 1.06E-31 19 67 1 1.987952E-02 0.00000 2.770685E-08 0.000000E+00 2.770685E-08 7.400000E+03 5.15E-32 19 68 1 2.018072E-02 0.00000 2.724119E-08 0.000000E+00 2.724119E-08 7.400000E+03 2.31E-32 19 69 1 2.048193E-02 0.00000 2.840534E-08 0.000000E+00 2.840534E-08 7.400000E+03 6.22E-32 19 70 1 2.078313E-02 0.00000 2.724119E-08 0.000000E+00 2.724119E-08 7.400000E+03 9.86E-32 19 71 1 2.108434E-02 0.00000 2.561137E-08 0.000000E+00 2.561137E-08 7.400000E+03 1.59E-31 19 72 1 2.138554E-02 0.00000 2.235174E-08 0.000000E+00 2.235174E-08 7.400000E+03 1.99E-31 19 73 1 2.168675E-02 0.00000 1.885928E-08 0.000000E+00 1.885928E-08 7.400000E+03 1.62E-31 19 74 1 2.198795E-02 0.00000 1.769513E-08 0.000000E+00 1.769513E-08 7.400000E+03 1.32E-31 19 75 1 2.228916E-02 0.00000 1.699664E-08 0.000000E+00 1.699664E-08 7.400000E+03 8.44E-32 19 76 1 2.259036E-02 0.00000 1.722947E-08 0.000000E+00 1.722947E-08 7.400000E+03 1.18E-31 19 77 1 2.289157E-02 0.00000 1.862645E-08 0.000000E+00 1.862645E-08 7.400000E+03 2.00E-31 19 78 1 2.319277E-02 0.00000 1.722947E-08 0.000000E+00 1.722947E-08 7.400000E+03 1.58E-31 19 79 1 2.349398E-02 0.00000 1.490116E-08 0.000000E+00 1.490116E-08 7.400000E+03 1.35E-31 19 80 1 2.379518E-02 0.00000 1.303852E-08 0.000000E+00 1.303852E-08 7.400000E+03 1.12E-31 19 81 1 2.409639E-02 0.00000 6.894031E-09 0.000000E+00 6.894031E-09 7.400000E+03 1.51E-31 20 82 1 2.439759E-02 0.00000 3.231577E-09 0.000000E+00 3.231577E-09 7.400000E+03 4.77E-31 20 83 1 2.469880E-02 0.00000 1.149005E-09 0.000000E+00 1.149005E-09 7.400000E+03 3.20E-31 20 84 1 2.500000E-02 0.00000 0.000000E+00 0.000000E+00 0.000000E+00 7.400000E+03 3.99E-24 20 Problem variables computed by this code. Variable Code Value Description ICYCLE 1 Present iteration number SNEGA 0.0 Total negative current, this solution SPOSA 0.0 total positive current, this solution STOTA 0.0 Total current in problem, this solution