📄 node14.html
字号:
<!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>2.3 Manipulation of a D.S. B resulting from a transient problem</TITLE><LINK REL=STYLESHEET TYPE="text/css" HREF="./Modulef.css" TITLE="Modulef CSS"><meta name="description" value="2.3 Manipulation of a D.S. B resulting from a transient problem"><meta name="keywords" value="Guide6"><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=tex2html650 HREF="node13.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif" ALT="previous"></A><A NAME=tex2html656 HREF="node11.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif" ALT="up"></A><A NAME=tex2html658 HREF="node15.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif" ALT="next"></A><A NAME=tex2html660 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif" ALT="contents"></A><A HREF="../Guide6-18/node14.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM ALT="[BIG]"></A><A HREF="../Guide6-14/node14.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM ALT="[Normal]"></A><A HREF="../Guide6-10/node14.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM ALT="[small]"></A><BR><B> Next: </B> <A NAME=tex2html659 HREF="node15.html">2.4 Changing the type of a </A><B>Up: </B> <A NAME=tex2html657 HREF="node11.html">2 Manipulation of solutions</A><B> Prev: </B> <A NAME=tex2html651 HREF="node13.html">2.2 Symmetrization and "gluing together" of </A><B><A HREF="node2.html" >Contents</A></B><HR SIZE=3 WIDTH="75%"><H1><A NAME=SECTION03330000000000000000>2.3 Manipulation of a D.S. B resulting from a transient problem</A></H1><P><P><P>The solution of a parabolic problem by a variable time-step method leads to computation of the solutionsat times <em> a priori</em> unknown. The aim of module INTERQ is therefore to interpolatethese solutions to obtain a <b> B</b> structure containing the solutions at those times selected.<P><PRE> SUBROUTINE INTERQ(M,XM,NQ,NFB,NIB,PAS,NPAS,T0,TSAV,NFBS,NIBS)C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++C SUBROUTINE INTERQ . SINGLE PRECISIONC AIM : INTERPOLATE TO THE ORDER NQ THE SOLUTIONS RESULTING FROMC EVO*GE(RK)C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++</PRE><P>where<P><UL><LI> M and XM represent the super array,<LI> NQ is the order of interpolation, it also encodes the manner in which the solutions is saved:<UL><LI> 1 <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img11.gif"> NQ <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img11.gif"> 10: interpolation of the order NQ and saves the solution at a specific time TSAV,<LI> 11 <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img11.gif"> NQ <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img11.gif"> 20: interpolation of the order NQ - 10 and saves the solution at all time steps,</UL><LI> NFB, NIB are the file numbers and the level of structure <b> B</b> resulting from the transient module,<LI> PAS is the time step,<LI> NPAS is the number of time steps (initial solution included),<LI> T0 is the initial time,<LI> TSAV is not used if NQ <b>></b> 10, otherwise it is the time at which we want to save the solution. The value oftime, obtained after T0, PAS and NPAS, is defined by:<P><DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img12.gif"></DIV><P><LI> NFBS, NIBS are the file number and level of structure <b> B</b> on exit.</UL><P>The solution of a parabolic problem by an arbitrary method can, in certain cases, lead tocomputing the solutions in two overlapping time intervals <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img13.gif"> and <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img14.gif">.The aim of module UNIONB is therefore to form a solution in the interval <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img15.gif"> by removing the intersection.<P>Several options are available. To describe them, let us consider the following example:<P><UL><LI> the first structure <b> B</b> contains the solutions at times<DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img16.gif"></DIV> while an associated array contains the values of <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img17.gif"> <LI> the second structure contains the solutions at times<DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img18.gif"></DIV> while an associated array contains the values of <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img19.gif"></UL> therefore<P><UL><LI> the concatenated structure contains the solutions at times:<DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img20.gif"></DIV> if we choose the truncation option in increasing order. Its associated arraycontains the values of all these times.<LI> the concatenated structure contains the solutions at times:<DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img21.gif"></DIV> if we choose the truncation option in decreasing order. Its associated arraycontains the values of all these times.<LI> the concatenated structure contains the solutions placed end-to-end, if this option ischosen</UL><P>Preprocessor <b> UNIOXX</b> calls module UNIONB which performs the concatenation.<P><PRE> SUBROUTINE UNIONB(M,XM,NFB1,NIB1,NFB2,NIB2,JOPTIO,NFBS,NIBS)C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++C SUBROUTINE UNIONB . SINGLE PRECISIONC AIM : GLUE TOGETHER 2 I.D.S. B INTO A SINGLE O.D.S. BC +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++</PRE><P>where<P><UL><LI> M and XM represent the super array,<LI> NFB1, NIB1 are the file number and level of the first <b> B</b> structure,<LI> NFB2, NIB2 are the file number and level of the second <b> B</b> structure,<LI> JOPTIO is the truncation option:<DL COMPACT><DT>0<DD>: place end-to-end<DT>1<DD>: increasing order according to the last associated array,<DT>2<DD>: decreasing order according to the last associated array,<P> </DL><LI> NFBS, NIBS are the file number and level of the output structure <b> B</b>.</UL><P><P><P><HR SIZE=3 WIDTH="75%"><IMG SRC="../icons/smallmod.gif" WIDTH=211 HEIGHT=50 ALIGN=BOTTOM ALT="Modulef"><A NAME=tex2html650 HREF="node13.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif" ALT="previous"></A><A NAME=tex2html656 HREF="node11.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif" ALT="up"></A><A NAME=tex2html658 HREF="node15.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif" ALT="next"></A><A NAME=tex2html660 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif" ALT="contents"></A><A HREF="../Guide6-18/node14.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM ALT="[BIG]"></A><A HREF="../Guide6-14/node14.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM ALT="[Normal]"></A><A HREF="../Guide6-10/node14.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM ALT="[small]"></A><BR><B> Next: </B> <A NAME=tex2html659 HREF="node15.html">2.4 Changing the type of a </A><B>Up: </B> <A NAME=tex2html657 HREF="node11.html">2 Manipulation of solutions</A><B> Prev: </B> <A NAME=tex2html651 HREF="node13.html">2.2 Symmetrization and "gluing together" of </A><B><A HREF="node2.html" >Contents</A></B><BR> <HR><P><ADDRESS></ADDRESS></BODY></HTML>⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -