test/python/tp.py

from esys.escript import *

class TransportPDE(object):
    def __init__(domain,num_equations=1,theta=0.,dt_max=-1.,trace=True):
        self.__domain=domain
        self.__num_equations=num_equations
        self.__theta=theta
        self.__dt_max=dt_max
        self.__transport_problem=None

     def getDomain(self):
        return self.__domain
     def getTheta(self):
        return self.__theta
     def getDt_max(self):
        return self.__dt_max
     def getNumEquations(self):
        return self.__getNumEquations
     def getFunctionSpace(self):
        if self.redueced():
           return ReducedSolution(self.getDomain())
        else:
           return Solution(self.getDomain())


     def __getNewTransportProblem(self):
       """
       returns an instance of a new operator
       """
       self.trace("New Transport problem is allocated.")
       return self.getDomain().newTransportProblem( \
                               self.getTheta(),
                               self.getDt_max(),
                               self.getNumEquations(), \
                               self.getFunctionSpace(), \
                               self.__matrix_type)

     def setValue(M=Data(),A=Data(),B=Data(),C=Data(),D=Data(),X=Data(),Y=Data()):
         self.__transport_problem=self.__getNewTransportProblem()
         self.getDomain().

     def solve(X=Data(),Y=Data()):
1	from esys.escript import *
2
3	class TransportPDE(object):
4	def __init__(domain,num_equations=1,theta=0.,dt_max=-1.,trace=True):
5	self.__domain=domain
6	self.__num_equations=num_equations
7	self.__theta=theta
8	self.__dt_max=dt_max
9	self.__transport_problem=None
10
11	def getDomain(self):
12	return self.__domain
13	def getTheta(self):
14	return self.__theta
15	def getDt_max(self):
16	return self.__dt_max
17	def getNumEquations(self):
18	return self.__getNumEquations
19	def getFunctionSpace(self):
20	if self.redueced():
21	return ReducedSolution(self.getDomain())
22	else:
23	return Solution(self.getDomain())
24
25
26	def __getNewTransportProblem(self):
27	"""
28	returns an instance of a new operator
29	"""
30	self.trace("New Transport problem is allocated.")
31	return self.getDomain().newTransportProblem( \
32	self.getTheta(),
33	self.getDt_max(),
34	self.getNumEquations(), \
35	self.getFunctionSpace(), \
36	self.__matrix_type)
37
38	def setValue(M=Data(),A=Data(),B=Data(),C=Data(),D=Data(),X=Data(),Y=Data()):
39	self.__transport_problem=self.__getNewTransportProblem()
40	self.getDomain().
41
42	def solve(X=Data(),Y=Data()):