/[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 4603 by gross, Tue Dec 17 23:37:08 2013 UTC revision 4604 by sshaw, Wed Dec 18 03:29:26 2013 UTC
# Line 562  class HTIWave(WaveBase): Line 562  class HTIWave(WaveBase):
562                        
563    
564             if DIM ==2 :             if DIM ==2 :
565            source_vector= [source_vector[0],source_vector[2]]                source_vector= [source_vector[0],source_vector[2]]
566                        
567             self.__r=Vector(0, DiracDeltaFunctions(self.__mypde.getDomain()))             self.__r=Vector(0, DiracDeltaFunctions(self.__mypde.getDomain()))
568             self.__r.setTaggedValue(self.__source_tag, source_vector)             self.__r.setTaggedValue(self.__source_tag, source_vector)
# Line 604  class HTIWave(WaveBase): Line 604  class HTIWave(WaveBase):
604                  sigma[1,0]=s                  sigma[1,0]=s
605    
606               else:               else:
607          e11=du[0,0]                  e11=du[0,0]
608          e33=du[1,1]                  e33=du[1,1]
609          sigma[0,0]=self.c11*e11+self.c13*e33                  sigma[0,0]=self.c11*e11+self.c13*e33
610          sigma[1,1]=self.c13*e11+self.c33*e33                  sigma[1,1]=self.c13*e11+self.c33*e33
611    
612          s=self.c66*(du[1,0]+du[0,1])                  s=self.c66*(du[1,0]+du[0,1])
613          sigma[0,1]=s                  sigma[0,1]=s
614          sigma[1,0]=s                  sigma[1,0]=s
615                            
616               self.__mypde.setValue(X=-sigma, y_dirac= self.__r * self.__wavelet.getAcceleration(t))               self.__mypde.setValue(X=-sigma, y_dirac= self.__r * self.__wavelet.getAcceleration(t))
617               return self.__mypde.getSolution()                           return self.__mypde.getSolution()            

Legend:
Removed from v.4603  
changed lines
  Added in v.4604

  ViewVC Help
Powered by ViewVC 1.1.26