On Tue, Mar 03, 2009 at 08:08:56PM -0500, Hatef Monajemi wrote: > Hi > > I was just wondering if it is possible to implement the Petrov-Galerkin type > of > weighting functions in dolfin. I am trying to solve a convection dominated > problem in which the use of standard Galerkin weighting generates the > oscillations which can be removed by upwinding. Is there any demo which shows > the implementation of Petrov-Galerkin type of testfunctions. > > Thanks, > > Hatef
Yes, look at demo/pde/stokes/stabilized It's very straightforward, just write v = TestFunction(V) v = v + delta*A(v) where delta is a stabilization parameter and A is a suitable operator. -- Anders
signature.asc
Description: Digital signature
_______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
