/[escript]/trunk/escript/test/python/test_pdetools.py
ViewVC logotype

Diff of /trunk/escript/test/python/test_pdetools.py

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

revision 3891 by jfenwick, Wed Jan 18 02:30:48 2012 UTC revision 3892 by jfenwick, Tue Apr 10 08:57:23 2012 UTC
# Line 106  class Test_pdetools_noLumping(unittest.T Line 106  class Test_pdetools_noLumping(unittest.T
106  #        l=Locator(self.domain,numpy.ones((self.domain.getDim(),)))  #        l=Locator(self.domain,numpy.ones((self.domain.getDim(),)))
107  #        d=Data(0, Solution(self.domain))  #        d=Data(0, Solution(self.domain))
108  #        l.setValue(d, 7)  #        l.setValue(d, 7)
109  #        self.assertTrue(sup(d)>6, "value not set") # guarantees we have set something  #        self.assertTrue(sup(d)>6, "value not set")     # guarantees we have set something
110  #        self.assertTrue(Lsup(l.getValue(d)-7)<self.RES_TOL, "value not set in the correct place")          #        self.assertTrue(Lsup(l.getValue(d)-7)<self.RES_TOL, "value not set in the correct place")        
111    
112    
# Line 153  class Test_pdetools_noLumping(unittest.T Line 153  class Test_pdetools_noLumping(unittest.T
153  #        l=Locator(self.domain,numpy.ones((self.domain.getDim(),)))  #        l=Locator(self.domain,numpy.ones((self.domain.getDim(),)))
154  #        d=Data(0, Solution(self.domain))  #        d=Data(0, Solution(self.domain))
155  #        l.setValue(d, 7)  #        l.setValue(d, 7)
156  #        self.assertTrue(sup(d)>6, "value not set") # guarantees we have set something  #        self.assertTrue(sup(d)>6, "value not set")     # guarantees we have set something
157  #        self.assertTrue(Lsup(l.getValue(d)-7)<self.RES_TOL, "value not set in the correct place")  #        self.assertTrue(Lsup(l.getValue(d)-7)<self.RES_TOL, "value not set in the correct place")
158                        
159                    
# Line 585  class Test_pdetools_noLumping(unittest.T Line 585  class Test_pdetools_noLumping(unittest.T
585                  -8.7934289814322  ])                  -8.7934289814322  ])
586    
587        def Ap(x):        def Ap(x):
588        out=dot(A,x)            out=dot(A,x)
589            for i in range(size(x)):            for i in range(size(x)):
590              out[i]/=A[i,i]              out[i]/=A[i,i]
591            return out            return out

Legend:
Removed from v.3891  
changed lines
  Added in v.3892

  ViewVC Help
Powered by ViewVC 1.1.26