/[escript]/trunk/downunder/py_src/seismic.py
ViewVC logotype

Diff of /trunk/downunder/py_src/seismic.py

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

revision 4691 by sshaw, Wed Feb 19 05:23:17 2014 UTC revision 4821 by sshaw, Tue Apr 1 04:58:33 2014 UTC
# Line 30  import sys Line 30  import sys
30  import time  import time
31  from esys.escript import *  from esys.escript import *
32  import esys.escript.unitsSI as U  import esys.escript.unitsSI as U
33  from esys.escript.linearPDEs import LinearSinglePDE, LinearPDESystem, WavePDE  from esys.escript.linearPDEs import LinearSinglePDE, LinearPDESystem, WavePDE, SolverOptions
34    
35  class Wavelet(object):  class Wavelet(object):
36          """          """
# Line 356  class SonicWave(WaveBase): Line 356  class SonicWave(WaveBase):
356    
357             self.__wavelet=wavelet             self.__wavelet=wavelet
358             self.__mypde=LinearSinglePDE(domain)             self.__mypde=LinearSinglePDE(domain)
359             if lumping: self.__mypde.getSolverOptions().setSolverMethod(self.__mypde.getSolverOptions().HRZ_LUMPING)             if lumping: self.__mypde.getSolverOptions().setSolverMethod(SolverOptions.HRZ_LUMPING)
360             self.__mypde.setSymmetryOn()             self.__mypde.setSymmetryOn()
361             self.__mypde.setValue(D=1./v_p**2)             self.__mypde.setValue(D=1./v_p**2)
362             self.__source_tag=source_tag             self.__source_tag=source_tag
# Line 449  class VTIWave(WaveBase): Line 449  class VTIWave(WaveBase):
449              self.__mypde.setValue(X=self.__mypde.createCoefficient('X'))              self.__mypde.setValue(X=self.__mypde.createCoefficient('X'))
450    
451          if lumping:          if lumping:
452              self.__mypde.getSolverOptions().setSolverMethod(              self.__mypde.getSolverOptions().setSolverMethod(SolverOptions.HRZ_LUMPING)
                     self.__mypde.getSolverOptions().HRZ_LUMPING)  
453          self.__mypde.setSymmetryOn()          self.__mypde.setSymmetryOn()
454          self.__mypde.setValue(D=rho*kronecker(DIM))          self.__mypde.setValue(D=rho*kronecker(DIM))
455          self.__source_tag=source_tag          self.__source_tag=source_tag
# Line 581  class HTIWave(WaveBase): Line 580  class HTIWave(WaveBase):
580                  self.__mypde.setValue(X=self.__mypde.createCoefficient('X'))                  self.__mypde.setValue(X=self.__mypde.createCoefficient('X'))
581                        
582             if lumping:             if lumping:
583                  self.__mypde.getSolverOptions().setSolverMethod(self.__mypde.getSolverOptions().HRZ_LUMPING)                  self.__mypde.getSolverOptions().setSolverMethod(SolverOptions.HRZ_LUMPING)
584             self.__mypde.setSymmetryOn()             self.__mypde.setSymmetryOn()
585             self.__mypde.setValue(D=rho*kronecker(DIM))             self.__mypde.setValue(D=rho*kronecker(DIM))
586             self.__source_tag=source_tag             self.__source_tag=source_tag
# Line 706  class TTIWave(WaveBase): Line 705  class TTIWave(WaveBase):
705             self.__wavelet=wavelet             self.__wavelet=wavelet
706    
707             self.__mypde=LinearPDESystem(domain)             self.__mypde=LinearPDESystem(domain)
708             if lumping: self.__mypde.getSolverOptions().setSolverMethod(self.__mypde.getSolverOptions().HRZ_LUMPING)             if lumping: self.__mypde.getSolverOptions().setSolverMethod(SolverOptions.HRZ_LUMPING)
709             self.__mypde.setSymmetryOn()             self.__mypde.setSymmetryOn()
710             self.__mypde.setValue(D=rho*kronecker(DIM), X=self.__mypde.createCoefficient('X'))             self.__mypde.setValue(D=rho*kronecker(DIM), X=self.__mypde.createCoefficient('X'))
711             self.__source_tag=source_tag             self.__source_tag=source_tag
# Line 808  class SonicHTIWave(WaveBase): Line 807  class SonicHTIWave(WaveBase):
807             self.__wavelet=wavelet             self.__wavelet=wavelet
808                        
809             self.__mypde=LinearPDESystem(domain)             self.__mypde=LinearPDESystem(domain)
810             if lumping: self.__mypde.getSolverOptions().setSolverMethod(self.__mypde.getSolverOptions().HRZ_LUMPING)             if lumping: self.__mypde.getSolverOptions().setSolverMethod(SolverOptions.HRZ_LUMPING)
811             self.__mypde.setSymmetryOn()             self.__mypde.setSymmetryOn()
812             self.__mypde.setValue(D=kronecker(2), X=self.__mypde.createCoefficient('X'))             self.__mypde.setValue(D=kronecker(2), X=self.__mypde.createCoefficient('X'))
813             self.__source_tag=source_tag             self.__source_tag=source_tag

Legend:
Removed from v.4691  
changed lines
  Added in v.4821

  ViewVC Help
Powered by ViewVC 1.1.26