test/python/debugIntegrate2.py


__copyright__="""  Copyright (c) 2006 by ACcESS MNRF
                    http://www.access.edu.au
                Primary Business: Queensland, Australia"""
__license__="""Licensed under the Open Software License version 3.0
             http://www.opensource.org/licenses/osl-3.0.php"""
from esys.escript import *
from esys.escript.linearPDEs import *
from esysi import finley

def makeParallelFaultSystem2D(xne=1,yne=1,xl=1.,yfaults=[0.,1.],integrationOrder=0,p1=0,p2=0):
 """
 generates a 2D rectangular domain and a system of parallel faults at yfaults 
   
 input data :
 xne - the number of elements over the horizontal edge of the rectangular mesh
 xl  - length of the horizontal edge of the domain
 yfault[y0,...yi,...y1] where y0 and y1 are the extremity of the block and yi the fault position in the (Oy) direction
   output data :
        a mesh with horizontal discontinuities at yi positions
 """

 #yne=xne*1./(xl+1.e-15)
 
 mshs=[]
 for i in range(len(yfaults)-1):
     yl=yfaults[i+1]-yfaults[i]
     msh=finley.Rectangle(xne,yne,1,xl,yl,p1,p2,useElementsOnFace=True)
     msh1=msh
     #print "mesh done"
     msh1.setX(msh.getX()+[0.,yfaults[i]])
     #print "Xnew"
     mshs.append(msh1)
     #print "apres append"
     #results=finley.joinFaces(mshs)
     mesh_joint=finley.JoinFaces(mshs)
     #print "apres faces jointes"
 return mesh_joint

def main():

    xne=10
    yne=20
    xl=10
    yl=10
    yfaults=[0.,5.,10.]
    numberfaults=1

    mesh=makeParallelFaultSystem2D(xne,yne,xl,yfaults,0,0,0)
    
    v=Vector(10.,ContinuousFunction(mesh))
    
    on_elem_v=interpolate(v,Function(mesh))
   
    value_integratedEL_vec=integrate(on_elem_v)

    print "vector--------"
    print "value interpolation",value_integratedEL_vec
    print "--------"
 
main()
1	jgs	147
2	elspeth	617	__copyright__=""" Copyright (c) 2006 by ACcESS MNRF
3			http://www.access.edu.au
4			Primary Business: Queensland, Australia"""
5			__license__="""Licensed under the Open Software License version 3.0
6			http://www.opensource.org/licenses/osl-3.0.php"""
7	jgs	149	from esys.escript import *
8			from esys.escript.linearPDEs import *
9			from esysi import finley
10	jgs	147
11			def makeParallelFaultSystem2D(xne=1,yne=1,xl=1.,yfaults=[0.,1.],integrationOrder=0,p1=0,p2=0):
12			"""
13			generates a 2D rectangular domain and a system of parallel faults at yfaults
14
15			input data :
16			xne - the number of elements over the horizontal edge of the rectangular mesh
17			xl - length of the horizontal edge of the domain
18			yfault[y0,...yi,...y1] where y0 and y1 are the extremity of the block and yi the fault position in the (Oy) direction
19			output data :
20			a mesh with horizontal discontinuities at yi positions
21			"""
22
23			#yne=xne*1./(xl+1.e-15)
24
25			mshs=[]
26			for i in range(len(yfaults)-1):
27			yl=yfaults[i+1]-yfaults[i]
28			msh=finley.Rectangle(xne,yne,1,xl,yl,p1,p2,useElementsOnFace=True)
29			msh1=msh
30			#print "mesh done"
31			msh1.setX(msh.getX()+[0.,yfaults[i]])
32			#print "Xnew"
33			mshs.append(msh1)
34			#print "apres append"
35			#results=finley.joinFaces(mshs)
36			mesh_joint=finley.JoinFaces(mshs)
37			#print "apres faces jointes"
38			return mesh_joint
39
40			def main():
41
42			xne=10
43			yne=20
44			xl=10
45			yl=10
46			yfaults=[0.,5.,10.]
47			numberfaults=1
48
49			mesh=makeParallelFaultSystem2D(xne,yne,xl,yfaults,0,0,0)
50
51			v=Vector(10.,ContinuousFunction(mesh))
52
53			on_elem_v=interpolate(v,Function(mesh))
54
55			value_integratedEL_vec=integrate(on_elem_v)
56
57			print "vector--------"
58			print "value interpolation",value_integratedEL_vec
59			print "--------"
60
61			main()
Name	Value
svn:eol-style	native
svn:keywords	Author Date Id Revision