node47.html
来自「htmdoc for html coding」· HTML 代码 · 共 456 行
HTML
456 行
<!DOCTYPE HTML PUBLIC "-//IETF//DTD HTML 3.2 Final//FR"><!-- Converted with LaTeX2HTML 95.1 (Fri Jan 20 1995) --><!-- by Nikos Drakos (nikos@cbl.leeds.ac.uk), CBLU, University of Leeds --><!-- Modified Simulog 03/97 --><HTML><HEAD><TITLE>1.11 DS MUA</TITLE><LINK REL=STYLESHEET TYPE="text/css" HREF="./Modulef.css" TITLE="Modulef CSS"><meta name="description" value="1.11 DS MUA"><meta name="keywords" value="Guide2"><meta name="resource-type" value="document"><meta name="distribution" value="global"></HEAD><BODY BGCOLOR="#FFFFFF"><P> <IMG SRC="../icons/smallmod.gif" WIDTH=211 HEIGHT=50 ALIGN=BOTTOM ALT="Modulef"><A NAME=tex2html1001 HREF="node46.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif" ALT="previous"></A><A NAME=tex2html1007 HREF="node35.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif" ALT="up"></A><A NAME=tex2html1009 HREF="node48.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif" ALT="next"></A><A NAME=tex2html1011 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif" ALT="contents"></A><A NAME=tex2html1012 HREF="node58.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/index_motif.gif" ALT="index"></A><A HREF="../Guide2-18/node47.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM ALT="[BIG]"></A><A HREF="../Guide2-14/node47.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM ALT="[Normal]"></A><A HREF="../Guide2-10/node47.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM ALT="[small]"></A><BR><B> Next: </B> <A NAME=tex2html1010 HREF="node48.html">1.12 DS NDL1</A><B>Up: </B> <A NAME=tex2html1008 HREF="node35.html">1 Description of DS by type</A><B> Prev: </B> <A NAME=tex2html1002 HREF="node46.html">1.10 DS MPLR</A><B><A HREF="node58.html" >Index</A></B><B><A HREF="node2.html" >Contents</A></B><HR SIZE=3 WIDTH="75%"><H1><A NAME=SECTION051120000000000000000>1.11 DS MUA</A></H1><A NAME=6022> </A><P><P><P><H2><A NAME=SECTION051121000000000000000>1.11.1 Contents</A></H2><P><P><P>This DS stores, for each line, the coefficients lying between the first column with a <em> a priori</em> non-zerocoefficient and the diagonal of a sparse finite element matrix. This type of storage is called <em> profile</em>or <em> skyline</em> storage.<P>DS <b> MUA</b> consists of 6 arrays of predefined order.<P><DL COMPACT><DT>Array MUA0:<DD> General information. <BR><P>This integer array contains 32 variables, consisting of a general descriptionof the job (title, date, name), of DS <b> MUA</b> (type, level, ...), andindicates the presence or absence of array <b> MUA1</b>.<P><DL COMPACT><DT> 1:20 <tt> TITRE</tt><DD><P>the job title in 20 words of 4 characters,<P> </DL><P><DL COMPACT><DT> 21:22 <tt> DATE</tt><DD><P>the data of creation in 2 words of 4 characters,<P> </DL><P><DL COMPACT><DT> 23:28 <tt> NOMCRE</tt><DD><P>the creator's name in 6 words of 4 characters,<P> </DL><P><DL COMPACT><DT> 29 <tt> 'MUA'</tt><DD><P>the DS type,<P> </DL><P><DL COMPACT><DT> 30 <tt> NIVEAU</tt><DD><P>the DS level,<P> </DL><P><DL COMPACT><DT> 31 <tt> ETAT</tt><DD><P>a reserved parameter,<P> </DL><P><DL COMPACT><DT> 32 <tt> NTACM</tt><DD><P>the number of supplementary arrays associated with the DS(they are described in array <b> MUA1</b>).<P> </DL><P><DT>Array MUA1:<DD> Description of any supplementary arrays. <BR><P>This array is analogous to array <b> B1</b> of DS <b> B</b> (see this DS).<P><DT>Array MUA2:<DD> General description of the matrix. <BR><P>This integer array contains 12 values.<P><DL COMPACT><DT> 1 <tt> NTYP</tt><DD><P>the type of matrix coefficients,<P> </DL><P><DL COMPACT><DT> 2 <tt> NBLOC</tt><DD><P>the number of blocks or pages in the matrix,<P> </DL><P><DL COMPACT><DT> 3 <tt> NTCOL</tt><DD><P>the number of columns of the largest block,<P> </DL><P><DL COMPACT><DT> 4 <tt> NMOPB</tt><DD><P>the number of words required in m.m. to store the largest block,<P> </DL><P><DL COMPACT><DT> 5 <tt> NMATA</tt><DD><P>the number of matrices on file,<P> </DL><P><DL COMPACT><DT> 6 <tt> NCODSA</tt><DD><P>the type of matrix storage:<UL><LI> 1: symmetric matrix, only the lower triangle is stored, line by line,<LI> 0: diagonal matrix, or<LI> -1: non-symmetric matrix, ordered line after line, with the diagonal coefficient placed at the end. <BR> </UL> </DL><P><DL COMPACT><DT> 7 <tt> LBDP</tt><DD><P>the largest difference + 1 between the 2 node numbers in the same element ifNCODSA is non-zero, and the half band-width in terms of nodes if NCODSA=0,<P> </DL><P><DL COMPACT><DT> 8 <tt> LBDPDL</tt><DD><P>the largest difference + 1 between the degree of freedom numbers of 2 nodes belonging to the same element ifNCODSA is non-zero, and the half band-width in terms of degrees of freedom if NCODSA=0,<P> </DL><P><DL COMPACT><DT> 9 <tt> ND</tt><DD><P>the number of degrees of freedom per node if it is constant, 0 otherwise,<P> </DL><P><DL COMPACT><DT> 10 <tt> NOE</tt><DD><P>the number of nodes,<P> </DL><P><DL COMPACT><DT> 11 <tt> NTDL</tt><DD><P>the matrix order,<P> </DL><P><DL COMPACT><DT> 12 <tt> MUA5AR</tt><DD><P>the assembly is already performed if 1,<BR> the assembly is not yet performed if 0.<P> </DL><P><DT>Array MUA3:<DD> Pointer to the last degree of freedom of each block. <BR><P>This integer array , of length NBLOC+1, contains:<P><UL><LI> MUA3(1) = 0,<LI> MUA3(i+1) = the number of the last degree of freedom of block i in the matrix.</UL><P><DT>Array MUA4:<DD> Pointer to the diagonal coefficient of each line. <BR><P>This integer array, of length NTDL+1 or 2 depending on NCODSA, contains:<P><UL><LI> NCODSA = 0<UL><LI> MUA4(1) = 0<LI> MUA4(2) = NTDL,</UL><LI> NCODSA <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img96.gif"> 0<UL><LI> MUA4(1) = 0<LI> MUA4(i+1) = the address of diagonal coefficient i.</UL></UL><P><DT>Array MUA5:<DD> The matrix coefficients. <BR><P>This NTYP type array contains the matrix coefficients:<P><UL><LI> MUA5(k) = value of the k-th matrix coefficient (if it is not paginated),<LI> or MUA5(k-MUA3(i)) the k-th coefficient of block i (if it is paginated).</UL><P> </DL><P><P><H2><A NAME=SECTION051122000000000000000>1.11.2 Corresponding tools</A></H2><P><P><P><H3><A NAME=SECTION051122100000000000000> Reading and writing a DS MUA</A></H3><P>This DS is a category 2 DS. A DS residing in main memory is written partially (only its 5generic arrays and any associated arrays described in array <b> MUA1</b>) on sequential or direct access file via module <b> SDECRI</b>.<P>Similarly, a DS residing on file is read and transferred partially (only its 5generic arrays and any associated arrays described in array <b> MUA1</b>)to main memory via module <b> SDLECT</b>; only array <b> MUA5</b> is addressed:<P><UL><LI> When reading, array <b> MUA5</b> is only addressed and the read is performed by the program via a WRITE type command:<P><blockquote> READ(NFMUA) LE,(M(IAMUA5-1+I),I=1,LE) </blockquote><P> or if it is a direct access file:<P><blockquote> READ(NFMUA,NEMUA) LE,(M(IAMUA5-1+I),I=1,LE)</blockquote><P><LI> When writing, array <b> MUA5</b> is transferred via a WRITE type command.</UL><P> <P><P><H3><A NAME=SECTION051122200000000000000> Printing a DS MUA and plotting its profile</A></H3><P><P><P>The contents (total or partial) of a DS <b> MUA</b> can be printed by module<b> IMMUA</b>. Preprocessor <b> IMAGXX</b> is used when a conversational call of <b> IMMUA</b> is desired.<P>The matrix profile contained in a file can be plotted via preprocessor <b> TRPRXX</b> (see [<A HREF="node57.html#mod_96">96</A>]).<P><H3><A NAME=SECTION051122300000000000000> Creation and modification module for DS MUA</A></H3><P>DS <b> MUA</b> is created using a DS <b> TAE</b> by assembly of the element matrices, for example via module <b> ASSMUA</b> and <b> ASMAPS</b>.<P><P><P><H2><A NAME=SECTION051123000000000000000>1.11.3 Storage structure</A></H2><P><P><P>Consider the following 2 matrix examples:<P><DL COMPACT><DT>Example 1:<DD> Symmetric matrix.<P>Consider the 12-th order matrix given below (the values correspond to the the rows in the storage and not to the matrix coefficients) constituting 3 pages:<P><DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img130.gif"></DIV> <P><P><P><P>Array <b> MUA5</b>, of length MUA4(13), is subdivided into 3 pages, line 1 to 5, 6 to 10 and11 to 12. Each page contains a maximum of 22 words. Thus, we have:<P><UL><LI> NCODSA = 1 <LI> array <b> MUA3</b>:<DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img131.gif"></DIV> <P> <LI> array <b> MUA4</b>:<DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img132.gif"></DIV> <P> </UL><P><P><DT>Example 2:<DD> Non-symmetric matrix.<P>Consider the 12-th order matrix given below (the values correspond to the the rows in the storage and not to the matrix coefficients) constituting 3 pages:<P><DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img133.gif"></DIV> <P><P><P><P>Here, each page contains a maximum of 32 words. Thus, we have:<P><UL><LI> NCODSA = -1 <LI> array <b> MUA3</b>:<DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img134.gif"></DIV> <P> <LI> array <b> MUA4</b>:<DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img135.gif"></DIV> <P> </UL> </DL><P> <HR SIZE=3 WIDTH="75%"><IMG SRC="../icons/smallmod.gif" WIDTH=211 HEIGHT=50 ALIGN=BOTTOM ALT="Modulef"><A NAME=tex2html1001 HREF="node46.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif" ALT="previous"></A><A NAME=tex2html1007 HREF="node35.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif" ALT="up"></A><A NAME=tex2html1009 HREF="node48.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif" ALT="next"></A><A NAME=tex2html1011 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif" ALT="contents"></A><A NAME=tex2html1012 HREF="node58.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/index_motif.gif" ALT="index"></A><A HREF="../Guide2-18/node47.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM ALT="[BIG]"></A><A HREF="../Guide2-14/node47.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM ALT="[Normal]"></A><A HREF="../Guide2-10/node47.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM ALT="[small]"></A><BR><B> Next: </B> <A NAME=tex2html1010 HREF="node48.html">1.12 DS NDL1</A><B>Up: </B> <A NAME=tex2html1008 HREF="node35.html">1 Description of DS by type</A><B> Prev: </B> <A NAME=tex2html1002 HREF="node46.html">1.10 DS MPLR</A><B><A HREF="node58.html" >Index</A></B><B><A HREF="node2.html" >Contents</A></B><BR> <HR><P><ADDRESS></ADDRESS></BODY></HTML>⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?