1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
|
<!DOCTYPE HTML SYSTEM "html.dtd">
<HTML>
<HEAD>
<TITLE>Section VII - file Storage Requirements </TITLE>
</HEAD>
<BODY>
<H2>Section VII - file Storage Requirements</H2>
<P>Depending on the requirements of a run, NEC-2 may use the following files.
</P>
<DL>
<DD>11,12,13,14,15,16 - scratch files for matrix manipulations. </DD>
<DD>20 - <A HREF="ngf.html">NGF</A> file </DD>
<DD>21 - <A HREF="somnec.html">Sommerfeld/Norton</A> data. </DD>
</DL>
<P>The scratch files are used only when the matrix will not fit into core
stroage. Fro a case that does not use the <A HREF="ngf.html">NGF</A> (but may
write a <A HREF="ngf.html">NGF</A> file), there are five options for matrix
storage. IF </P>
<DL>
<DD>N=the number of equations (number of segments plus twice the number of
patches), </DD>
<DD>N<SUB>x</SUB>=the number of equations for a symmetric section, and </DD>
<DD>I<SUB>R</SUB>=number of complex numbers for the matrix in core storage
(4000), </DD>
</DL>
<P>then the case, indicated by the value of ICASE in the code, are</P>
<OL>
<LI>matrix in core, no symmetry (N<SUP>2</SUP> <=I<SUB>R</SUB>). </LI>
<LI>matrix in core,symmetry (NN<SUB>x</SUB> <=I<SUB>R</SUB>). </LI>
<LI>matrix in core, no symmetry (N<SUP>2</SUP> >=I<SUB>R</SUB>). </LI>
<LI>matrix in core,symmetry blocks fit in core (NN<SUB>x</SUB> >
I<SUB>R</SUB>,N<SUB>x</SUB><SUP>2</SUP><=I<SUB>R</SUB>). </LI>
<LI>matrix out of core,symmetry blocks do not fit in core
N<SUB>x</SUB><SUP>2</SUP>>I<SUB>R</SUB>). </LI>
</OL>
<P>File storage is used in cases 3,4 and 5. Only the four files (11,12,13 and
14) are used when the <A HREF="ngf.html">NGF</A> is not in use. The size of
each file is approximately @NN<SUB>x</SUB> words. If the computer system
requires that the user specify the file size, a safety margin should be
included in the request. A more accurate estimate of the file size is </P>
<DL>
<DD>L=2NN<SUB>c</SUB>[N<SUB>x</SUB>/N<SUB>c</SUB>], </DD>
</DL>
<P>Where </P>
<DL>
<DD>N<SUB>c</SUB>=[I<SUB>R</SUB>)/(2N)] </DD>
</DL>
<P>and [ ] indicates truncation. N<SUB>c</SUB>, which is the number of matrix
columns in I/O block, must be at least 1.</P>
<P>When the <A HREF="ngf.html">NGF</A> is used, all six scratch files may be
required. For the <A HREF="ngf.html">NGF</A> the matrix is partitioned into
four sections as<BR>
</P>
<PRE> _ _
|A B|
|C D|
- -
</PRE>
<P>A is the matrix for the <A HREF="ngf.html">NGF</A> structure and is factored
before the <A HREF="ngf.html">NGF</A>file is written. The stroage case for A is
indicated in the <A HREF="ngf.html">NGF</A>lable by the value of ICASE (see
exampel 10 in section IV). When the <A HREF="ngf.html">NGF</A> is used, matrix
A is read from file 20 and, if ICASE is 3,4, or 5, is stored on file 13. The
size of file 13 when the <A HREF="ngf.html">NGF</A> is used is approximately:
</P>
<PRE> <U>ICASE</U> <U>Length of file 13</U>
3 4N</PRE>
<P><SUP>2</SUP> 4 2NN<SUB>x</SUB> 5 4NN<SUB>x</SUB> There are four options for
storage of the matrices B, C, and D. these are associated with the integer
ICASX as follows:<BR>
<U>ICASE</U> </P>
<OL>
<LI>A<SUB>R</SUB>, B, C, and D fit in core together where<BR>
A<SUB>R</SUB>=A for ICASE=1 or 2,<BR>
=one I/O block of A for ICASE=3 or 5,<BR>
=one submatrix for ICASE=4 </LI>
<LI>B, C, an dD fit in core but not with A<SUB>R</SUB>. This is possible only
for ICASE=3, 4, or 5 when A does not need dedicated space in core.
A<SUB>R</SUB> an dB must also fit in core together. </LI>
<LI>B, C, and D do not fit in core, but D fits in core alone.<BR>
A and D must fit together if ICASE=1 or 2. </LI>
<LI>D does not fit in core. </LI>
</OL>
<P>The sizees of matrices B, C, and D depend on the number of new unknowns
N<SUB>n</SUB> where </P>
<DL>
<DT>N<SUB>n</SUB>=N<SUB>n</SUB>+N<SUB>t</SUB>+2N<SUB>p</SUB>+10N<SUB>q</SUB>
</DT>
<DD>N<SUB>s</SUB>=number of new segments added to <A HREF="ngf.html">NGF</A>,
</DD>
<DD>N<SUB>t</SUB>=number of <A HREF="ngf.html">NGF</A> segments connected to
new segments or patches, </DD>
<DD>N<SUB>p</SUB>=number of new patches, </DD>
<DD>N<SUB>q</SUB>=number of <A HREF="ngf.html">NGF</A> patches connected to new
segments. </DD>
</DL>
<P>The sizes of matrices B and C are 2NN<SUB>n</SUB> and the size of D is
2N<SUB>n</SUB><SUP>2</SUP> words. The file lengths are approximately
2N<SUB>n</SUB><SUP>2</SUP> words for files 11 and 12 and 2NN<SUB>n</SUB> for
files 14, 15 and 16. when ICASX is 1 these files are not used, and when ICASX
is 2 file 16 is not used. </P>
<P>The length of the <A HREF="ngf.html">NGF</A> file (20) is approximately
4N(N<SUB>x</SUB>+3). The length of the Summerfeld/Norton data file (21) is
about 2200 words. </P>
<HR>
</BODY>
</HTML>
|
