/[escript]/trunk/doc/user/finley.tex
ViewVC logotype

Diff of /trunk/doc/user/finley.tex

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/esys2/doc/user/finley.tex revision 107 by jgs, Thu Jan 27 06:21:48 2005 UTC trunk/doc/user/finley.tex revision 1318 by ksteube, Wed Sep 26 04:39:14 2007 UTC
# Line 1  Line 1 
1    %
2  % $Id$  % $Id$
3    %
4    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5    %
6    %           Copyright 2003-2007 by ACceSS MNRF
7    %       Copyright 2007 by University of Queensland
8    %
9    %                http://esscc.uq.edu.au
10    %        Primary Business: Queensland, Australia
11    %  Licensed under the Open Software License version 3.0
12    %     http://www.opensource.org/licenses/osl-3.0.php
13    %
14    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
15    %
16    
17    \chapter{ The Module \finley}
 \chapter{ The module \finley}  
18   \label{CHAPTER ON FINLEY}   \label{CHAPTER ON FINLEY}
19    
20  \begin{figure}  \begin{figure}
21  \centerline{\includegraphics[width=\figwidth]{FinleyMesh}}  \centerline{\includegraphics[width=\figwidth]{figures/FinleyMesh.eps}}
22  \caption{Subdivision of an Ellipse into triangles order 1 (\finleyelement{Tri3})}  \caption{Subdivision of an Ellipse into triangles order 1 (\finleyelement{Tri3})}
23  \label{FINLEY FIG 0}  \label{FINLEY FIG 0}
24  \end{figure}  \end{figure}
25    
26  \begin{figure}  \begin{figure}
27  \centerline{\includegraphics[width=\figwidth]{FinleyContact}}  \centerline{\includegraphics[width=\figwidth]{figures/FinleyContact.eps}}
28  \caption{Mesh around a contact region (\finleyelement{Rec4})}  \caption{Mesh around a contact region (\finleyelement{Rec4})}
29  \label{FINLEY FIG 01}  \label{FINLEY FIG 01}
30  \end{figure}  \end{figure}
# Line 26  It supports unstructured, 1D, 2D and 3D Line 39  It supports unstructured, 1D, 2D and 3D
39  library through the \LinearPDE class of \escript supporting its full functionality. {\it finley}  library through the \LinearPDE class of \escript supporting its full functionality. {\it finley}
40  is parallelized using the OpenMP \index{OpenMP} paradigm.  is parallelized using the OpenMP \index{OpenMP} paradigm.
41    
42  \subsection{Meshes}  \section{Formulation}
43    
44    For a single PDE with a solution with a single component the linear PDE is defined in the
45    following form:
46    \begin{equation}\label{FINLEY.SINGLE.1}
47    \begin{array}{cl} &
48    \displaystyle{
49    \int\hackscore{\Omega}
50    A\hackscore{jl} \cdot v\hackscore{,j}u\hackscore{,l}+ B\hackscore{j} \cdot v\hackscore{,j} u+ C\hackscore{l} \cdot v u\hackscore{,l}+D \cdot vu \; d\Omega }  \\
51    + & \displaystyle{\int\hackscore{\Gamma} d \cdot vu \; d{\Gamma} }
52    +  \displaystyle{\int\hackscore{\Gamma^{contact}} d^{contact} \cdot [v][u] \; d{\Gamma} } \\
53    = & \displaystyle{\int\hackscore{\Omega}  X\hackscore{j} \cdot v\hackscore{,j}+ Y \cdot v \; d\Omega }\\
54    + & \displaystyle{\int\hackscore{\Gamma} y \cdot v \; d{\Gamma}}  +
55    \displaystyle{\int\hackscore{\Gamma^{contact}} y^{contact}\cdot [v] \; d{\Gamma}} \\
56    \end{array}
57    \end{equation}
58    
59    \section{Meshes}
60  To understand the usage of \finley one needs to have an understanding of how the finite element meshes  To understand the usage of \finley one needs to have an understanding of how the finite element meshes
61  \index{FEM!mesh} are defined. \fig{FINLEY FIG 0} shows an example of the  \index{FEM!mesh} are defined. \fig{FINLEY FIG 0} shows an example of the
62  subdivision of an ellipse into so called elements \index{FEM!elements} \index{element}.  subdivision of an ellipse into so called elements \index{FEM!elements} \index{element}.
# Line 56  to second node the domain has to lie on Line 86  to second node the domain has to lie on
86  the domain has to lie on the left hand side when moving counterclockwise). If the gradient on the  the domain has to lie on the left hand side when moving counterclockwise). If the gradient on the
87  surface of the domain is to be calculated rich face elements face to be used. Rich elements on a face  surface of the domain is to be calculated rich face elements face to be used. Rich elements on a face
88  are identical to interior elements but with a modified order of nodes such that the 'first' face of the element aligns  are identical to interior elements but with a modified order of nodes such that the 'first' face of the element aligns
89  with the surface of the domian. In \fig{FINLEY FIG 0}  with the surface of the domain. In \fig{FINLEY FIG 0}
90  elements of the type \finleyelement{Tri3Face} are used.  elements of the type \finleyelement{Tri3Face} are used.
91  The face element reference number $20$ as a rich face element is defined by the nodes  The face element reference number $20$ as a rich face element is defined by the nodes
92  with reference numbers $11$, $0$ and $9$. Notice that the face element $20$ is identical to the  with reference numbers $11$, $0$ and $9$. Notice that the face element $20$ is identical to the
# Line 107  the nodes within an element. Line 137  the nodes within an element.
137  \linev{\finleyelement{Hex20}}{\finleyelement{Rec8}}{\finleyelement{Hex20Face}}{\finleyelement{Rec8_Contact}}{\finleyelement{Hex20Face_Contact}}  \linev{\finleyelement{Hex20}}{\finleyelement{Rec8}}{\finleyelement{Hex20Face}}{\finleyelement{Rec8_Contact}}{\finleyelement{Hex20Face_Contact}}
138  \end{tablev}  \end{tablev}
139  \caption{Finley elements and corresponding elements to be used on domain faces and contacts.  \caption{Finley elements and corresponding elements to be used on domain faces and contacts.
140  The rich types have to be used if the gradient of function is to be calculated on faces and contacts, resepctively.}  The rich types have to be used if the gradient of function is to be calculated on faces and contacts, respectively.}
141  \label{FINLEY TAB 1}  \label{FINLEY TAB 1}
142  \end{table}  \end{table}
143    
# Line 151  for i in range(ContactElement_Num): Line 181  for i in range(ContactElement_Num):
181     for j in range(ContactElement_numNodes): print " %d"%ContactElement_Nodes[i][j]     for j in range(ContactElement_numNodes): print " %d"%ContactElement_Nodes[i][j]
182     print "\n"     print "\n"
183  # point sources (not supported yet)  # point sources (not supported yet)
184  write("Point1 0",face_element_typ,numFaceElements)  write("Point1 0",face_element_type,numFaceElements)
185  \end{python}  \end{python}
186    
187  The following example of a mesh file defines the mesh shown in \fig{FINLEY FIG 01}:  The following example of a mesh file defines the mesh shown in \fig{FINLEY FIG 01}:
# Line 227  $7$, $10$, $15$ and $20$, respectively. Line 257  $7$, $10$, $15$ and $20$, respectively.
257    
258  \subsection{Linear Solvers in \LinearPDE}  \subsection{Linear Solvers in \LinearPDE}
259  Currently \finley supports the linear solvers \PCG, \GMRES, \PRESTWENTY and \BiCGStab.  Currently \finley supports the linear solvers \PCG, \GMRES, \PRESTWENTY and \BiCGStab.
260  For \GMRES the options \var{trancation} and \var{restart} of the \method{getSolution} can be  For \GMRES the options \var{truncation} and \var{restart} of the \method{getSolution} can be
261  used to control the trunction and restart during iteration. Default values are  used to control the truncation and restart during iteration. Default values are
262  \var{truncation}=5 and \var{restart}=20.  \var{truncation}=5 and \var{restart}=20.
263  The default solver is \BiCGStab  but if the symmetry flag is set \PCG is the default solver.  The default solver is \BiCGStab  but if the symmetry flag is set \PCG is the default solver.
264  \finley supports the solver options \var{iter_max} which specifies the maximum number of iterations steps,  \finley supports the solver options \var{iter_max} which specifies the maximum number of iterations steps,
# Line 237  In some installations \finley supports t Line 267  In some installations \finley supports t
267  solver options \var{reordering}=\constant{util.NO_REORDERING},  solver options \var{reordering}=\constant{util.NO_REORDERING},
268  \constant{util.MINIMUM_FILL_IN} or \constant{util.NESTED_DISSECTION} (default is \constant{util.NO_REORDERING}),  \constant{util.MINIMUM_FILL_IN} or \constant{util.NESTED_DISSECTION} (default is \constant{util.NO_REORDERING}),
269  \var{drop_tolerance} specifying the threshold for values to be dropped in the  \var{drop_tolerance} specifying the threshold for values to be dropped in the
270  incomplete elimation process (default is 0.01) and \var{drop_storage} specifying the maximum increase  incomplete elimination process (default is 0.01) and \var{drop_storage} specifying the maximum increase
271  in storage allowed in the  in storage allowed in the
272  incomplete elimation process (default is 1.20).  incomplete elimination process (default is 1.20).
273    
274  \subsection{Functions}  \subsection{Functions}
275  \begin{funcdesc}{Mesh}{fileName,integrationOrder=-1}  \begin{funcdesc}{Mesh}{fileName,integrationOrder=-1}
# Line 330  in $x_1$-direction and $x_2$-direction, Line 360  in $x_1$-direction and $x_2$-direction,
360  \end{funcdesc}  \end{funcdesc}
361    
362  \begin{funcdesc}{GlueFaces}{meshList,safetyFactor=0.2,tolerance=1.e-13}  \begin{funcdesc}{GlueFaces}{meshList,safetyFactor=0.2,tolerance=1.e-13}
363  Generates a new \Domain object from the list \var{mehList} of \finley meshes.  Generates a new \Domain object from the list \var{meshList} of \finley meshes.
364  Nodes in face elements whose difference of coordinates is less then \var{tolerance} times the  Nodes in face elements whose difference of coordinates is less then \var{tolerance} times the
365  diameter of the domain are merged. The corresponding face elements are removed from the mesh.    diameter of the domain are merged. The corresponding face elements are removed from the mesh.  
366    
# Line 338  TODO: explain \var{safetyFactor} and sho Line 368  TODO: explain \var{safetyFactor} and sho
368  \end{funcdesc}  \end{funcdesc}
369    
370  \begin{funcdesc}{JoinFaces}{meshList,safetyFactor=0.2,tolerance=1.e-13}  \begin{funcdesc}{JoinFaces}{meshList,safetyFactor=0.2,tolerance=1.e-13}
371  Generates a new \Domain object from the list \var{mehList} of \finley meshes.  Generates a new \Domain object from the list \var{meshList} of \finley meshes.
372  Face elements whose nodes coordinates have difference is less then \var{tolerance} times the  Face elements whose nodes coordinates have difference is less then \var{tolerance} times the
373  diameter of the domain are combined to form a contact element \index{element!contact}  diameter of the domain are combined to form a contact element \index{element!contact}
374  The corresponding face elements are removed from the mesh.    The corresponding face elements are removed from the mesh.  

Legend:
Removed from v.107  
changed lines
  Added in v.1318

  ViewVC Help
Powered by ViewVC 1.1.26