dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

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

Attachment: signature.asc
Description: Digital signature