user/examples/mytools.py

# $Id$
from esys.linearPDEs import LinearPDE
import numarray
class Helmholtz(LinearPDE):
   def setValue(self,kappa=0,omega=1,f=0,eta=0,g=0)
        # get spatial dimension 
        ndim=self.getDim()
        # get kronecker symbol from numarray: kronecker[i,j]=1 for i=j and =0 else
        kronecker=numarray.identity(ndim)
        # map kappa, omega, f, eta, g to the coefficients of the general PDE
        self._setValue(A=kappa*kronecker,D=omega,Y=f,d=eta,y=g)
1	jgs	102	# $Id$
2			from esys.linearPDEs import LinearPDE
3			import numarray
4			class Helmholtz(LinearPDE):
5	jgs	108	def setValue(self,kappa=0,omega=1,f=0,eta=0,g=0)
6			# get spatial dimension
7			ndim=self.getDim()
8			# get kronecker symbol from numarray: kronecker[i,j]=1 for i=j and =0 else
9			kronecker=numarray.identity(ndim)
10			# map kappa, omega, f, eta, g to the coefficients of the general PDE
11			self._setValue(A=kappa*kronecker,D=omega,Y=f,d=eta,y=g)
Name	Value
svn:eol-style	native
svn:keywords	Author Date Id Revision