finley/py_src/readers.py


########################################################
#
# Copyright (c) 2003-2008 by University of Queensland
# Earth Systems Science Computational Center (ESSCC)
# http://www.uq.edu.au/esscc
#
# Primary Business: Queensland, Australia
# Licensed under the Open Software License version 3.0
# http://www.opensource.org/licenses/osl-3.0.php
#
########################################################

__copyright__="""Copyright (c) 2003-2008 by University of Queensland
Earth Systems Science Computational Center (ESSCC)
http://www.uq.edu.au/esscc
Primary Business: Queensland, Australia"""
__license__="""Licensed under the Open Software License version 3.0
http://www.opensource.org/licenses/osl-3.0.php"""
__url__="http://www.uq.edu.au/esscc/escript-finley"

"""
some mesh handling 

@var __author__: name of author
@var __licence__: licence agreement
@var __url__: url entry point on documentation
@var __version__: version
@var __date__: date of the version
"""

__author__="Lutz Gross, l.gross@uq.edu.au"

from esys.escript import *
from esys.pycad.gmsh import Design as GMSHDesign
from finleycpp import ReadGmsh

def MakeDomain(design,integrationOrder=-1, reducedIntegrationOrder=-1, optimizeLabeling=True):
        """
        creates a Finley L{Domain} from a L{esys.pycad.design.Design} object. Currently only
        gmsh is supported.

        @param design: the geometry
        @type design: L{esys.pycad.design.Design}
        @param integrationOrder: integration order. If -1 the default is used.
        @type integrationOrder: C{int}
        @param reducedIntegrationOrder: reduced integration order. If -1 the default is used.
        @type reducedIntegrationOrder: C{int}
        @param optimizeLabeling: if set the labeling of the mesh nodes is optimized
        @type optimizeLabeling: C{bool}
        @return: the Finley domain defined by the designs 
        @rtype: L{Domain}
        """
        if isinstance(design, GMSHDesign):
            mshname=design.getMeshHandler()
            dom = ReadGmsh(mshname,
                           design.getDim(),
                           integrationOrder,
                           reducedIntegrationOrder,
                           optimizeLabeling)
        else: 
            raise TypeError("Finley does not support %s designs."%design.__class__.__name__)
        # fill in the tag map
        design.getTagMap().passToDomain(dom)
        return dom
1	ksteube	1809
2			########################################################
3	ksteube	1312	#
4	ksteube	1809	# Copyright (c) 2003-2008 by University of Queensland
5			# Earth Systems Science Computational Center (ESSCC)
6			# http://www.uq.edu.au/esscc
7	ksteube	1312	#
8	ksteube	1809	# Primary Business: Queensland, Australia
9			# Licensed under the Open Software License version 3.0
10			# http://www.opensource.org/licenses/osl-3.0.php
11	ksteube	1312	#
12	ksteube	1809	########################################################
13	ksteube	1312
14	ksteube	1809	__copyright__="""Copyright (c) 2003-2008 by University of Queensland
15			Earth Systems Science Computational Center (ESSCC)
16			http://www.uq.edu.au/esscc
17			Primary Business: Queensland, Australia"""
18			__license__="""Licensed under the Open Software License version 3.0
19			http://www.opensource.org/licenses/osl-3.0.php"""
20			__url__="http://www.uq.edu.au/esscc/escript-finley"
21
22	gross	934	"""
23			some mesh handling
24
25			@var __author__: name of author
26			@var __licence__: licence agreement
27			@var __url__: url entry point on documentation
28			@var __version__: version
29			@var __date__: date of the version
30			"""
31
32			__author__="Lutz Gross, l.gross@uq.edu.au"
33
34			from esys.escript import *
35			from esys.pycad.gmsh import Design as GMSHDesign
36			from finleycpp import ReadGmsh
37
38			def MakeDomain(design,integrationOrder=-1, reducedIntegrationOrder=-1, optimizeLabeling=True):
39			"""
40			creates a Finley L{Domain} from a L{esys.pycad.design.Design} object. Currently only
41			gmsh is supported.
42
43			@param design: the geometry
44			@type design: L{esys.pycad.design.Design}
45			@param integrationOrder: integration order. If -1 the default is used.
46			@type integrationOrder: C{int}
47			@param reducedIntegrationOrder: reduced integration order. If -1 the default is used.
48			@type reducedIntegrationOrder: C{int}
49			@param optimizeLabeling: if set the labeling of the mesh nodes is optimized
50	ksteube	990	@type optimizeLabeling: C{bool}
51	gross	934	@return: the Finley domain defined by the designs
52			@rtype: L{Domain}
53			"""
54			if isinstance(design, GMSHDesign):
55			mshname=design.getMeshHandler()
56			dom = ReadGmsh(mshname,
57			design.getDim(),
58			integrationOrder,
59			reducedIntegrationOrder,
60			optimizeLabeling)
61			else:
62			raise TypeError("Finley does not support %s designs."%design.__class__.__name__)
63	gross	1044	# fill in the tag map
64			design.getTagMap().passToDomain(dom)
65	gross	934	return dom