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

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

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

revision 3502 by gross, Thu Apr 28 05:06:24 2011 UTC revision 3510 by gross, Fri May 13 06:09:49 2011 UTC
# Line 265  class DarcyFlow(object): Line 265  class DarcyFlow(object):
265        if self.solver  == self.SIMPLE or self.solver  == self.POST :        if self.solver  == self.SIMPLE or self.solver  == self.POST :
266          self.__pde_p.setValue(X=self.__g ,          self.__pde_p.setValue(X=self.__g ,
267                                Y=self.__f,                                Y=self.__f,
268                                y=-util.inner(self.domain.getNormal(),u0 * self.location_of_fixed_flux),                                y= - util.inner(self.domain.getNormal(),u0 * self.location_of_fixed_flux),
269                                r=p0)                                r=p0)
270          p=self.__pde_p.getSolution()          p=self.__pde_p.getSolution()
271          u = self.getFlux(p, u0)          u = self.getFlux(p, u0)
# Line 343  class DarcyFlow(object): Line 343  class DarcyFlow(object):
343                                   X=0.5*(self.__g - u - util.tensor_mult(self.__permeability,util.grad(p0)) ),                                   X=0.5*(self.__g - u - util.tensor_mult(self.__permeability,util.grad(p0)) ),
344                                   y= escript.Data(),                                   y= escript.Data(),
345                                   r=escript.Data())                                   r=escript.Data())
346    
347        dp=self.__pde_p.getSolution()        dp=self.__pde_p.getSolution()
348        p=GMRES(dp,        p=GMRES(dp,
349                self.__STAB_Aprod,                self.__STAB_Aprod,
# Line 362  class DarcyFlow(object): Line 363  class DarcyFlow(object):
363        u=self.getFlux(p0, u0)          u=self.getFlux(p0, u0)  
364            self.__pde_p.setValue( Y= self.__f,            self.__pde_p.setValue( Y= self.__f,
365                                   X=  0.5*(self.__g + u - util.tensor_mult(self.__permeability,util.grad(p0)) ),                                   X=  0.5*(self.__g + u - util.tensor_mult(self.__permeability,util.grad(p0)) ),
366                                   y= -   util.inner(self.domain.getNormal(), u),                                   y=  -  util.inner(self.domain.getNormal(), u),
367                                   r=escript.Data())                                   r=escript.Data())
368        dp=self.__pde_p.getSolution()        dp=self.__pde_p.getSolution()
369          
370          print dp
371              print p0+dp
372              
373        p=GMRES(dp,        p=GMRES(dp,
374                self.__SYMSTAB_Aprod,                self.__SYMSTAB_Aprod,
375            p0,            p0,
# Line 410  class DarcyFlow(object): Line 415  class DarcyFlow(object):
415        u = -self.__pde_v.getSolution()        u = -self.__pde_v.getSolution()
416        self.__pde_p.setValue(Y=escript.Data(),        self.__pde_p.setValue(Y=escript.Data(),
417                              X=0.5*(-u+util.tensor_mult(self.__permeability,gp)),                              X=0.5*(-u+util.tensor_mult(self.__permeability,gp)),
418                              y=   util.inner(self.domain.getNormal(), u),                              y=escript.Data(),
419                              r=escript.Data())                              r=escript.Data())
420            
421        return  self.__pde_p.getSolution()        return  self.__pde_p.getSolution()

Legend:
Removed from v.3502  
changed lines
  Added in v.3510

  ViewVC Help
Powered by ViewVC 1.1.26