/[escript]/trunk/escript/py_src/mountains.py
ViewVC logotype

Diff of /trunk/escript/py_src/mountains.py

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

revision 2626 by jfenwick, Fri Aug 21 06:30:25 2009 UTC revision 2627 by gross, Mon Aug 24 01:57:34 2009 UTC
# Line 61  class Mountains: Line 61  class Mountains:
61      self.__PDE_W = LinearPDE(domain)      self.__PDE_W = LinearPDE(domain)
62      self.__PDE_W.setSymmetryOn()      self.__PDE_W.setSymmetryOn()
63      A=kronecker(domain)*eps*0      A=kronecker(domain)*eps*0
64      A[self.__DIM-1,self.__DIM-1]=(sup(FunctionOnBoundary(self.__domain).getSize())/log(2.))**2      A[self.__DIM-1,self.__DIM-1]=(0.3*(sup(z)-inf(z))/log(2.))**2
65        # A[self.__DIM-1,self.__DIM-1]=(sup(FunctionOnBoundary(self.__domain).getSize())/log(2.))**2
66      self.__PDE_W.setValue(D=1, A=A, q=whereZero(sup(z)-z))      self.__PDE_W.setValue(D=1, A=A, q=whereZero(sup(z)-z))
67    
68      self.__PDE_H = LinearPDE(domain)      self.__PDE_H = LinearPDE(domain)
69      self.__PDE_H.setSymmetryOn()      self.__PDE_H.setSymmetryOn()
70      if reduced: self.__PDE_H.setReducedOrderOn()      if reduced: self.__PDE_H.setReducedOrderOn()
71      A=kronecker(domain)*0      # A=kronecker(domain)*0
72      A[self.__DIM-1,self.__DIM-1]=0.1      # A[self.__DIM-1,self.__DIM-1]=0.1
73      self.__PDE_H.setValue(D=1.0)      self.__PDE_H.setValue(D=1.0)
74      # self.__PDE_H.getSolverOptions().setSolverMethod(SolverOptions.LUMPING)      # self.__PDE_H.getSolverOptions().setSolverMethod(SolverOptions.LUMPING)
75    

Legend:
Removed from v.2626  
changed lines
  Added in v.2627

  ViewVC Help
Powered by ViewVC 1.1.26