← Back to team overview

dolfin team mailing list archive

Re: [HG DOLFIN] Work on getting dg advection-diffusion demo to work in python

 



Anders Logg wrote:
On Mon, Dec 22, 2008 at 05:52:01PM +0100, Anders Logg wrote:
On Mon, Dec 22, 2008 at 03:56:10PM +0000, Garth N. Wells wrote:

DOLFIN wrote:
One or more new changesets pushed to the primary dolfin repository.
A short summary of the last three changesets is included below.

changeset:   5428:5ee944b2216752dd3078f0906c9500837a65097a
tag:         tip
user:        "Johan Hake <hake@xxxxxxxxx>"
date:        Sat Dec 20 22:59:46 2008 +0100
files:       demo/pde/dg/advection-diffusion/cpp/main.cpp demo/pde/dg/advection-diffusion/python/demo.py dolfin/function/SpecialFunctions.cpp dolfin/function/SpecialFunctions.h site-packages/dolfin/specialfunctions.py
description:
Work on getting dg advection-diffusion demo to work in python
 - Changed name of OutFlowFacet function to Is...
 - IsOutFlowFacet is instantiated using a FunctionSpace and a field funtion
 - the evaluation of IsOutFlowFacet is done by evaluating the field in the midpoint
   of the facet
 - Does not work properly yet.

What's the status of this? Is anyone working on it and what needs to be done?

Garth
We need to move part of FunctionSpace::eval() to another function so
we can reuse it from IsOutFlowFacet::eval(). I'll have a look as soon
as I get a chance, maybe even tonight.

Once this is in place, the two advection-diffusion demos need to be
looked at, and then the two nonlinear demos.

I've made the changes to Function, FunctionSpace, and IsOutflowFacet.


The C++ and Python results look different.

Garth

Wouldn't it be natural to use IsOutflowFacet in both the C++ and
Python demos? It seems unnecessary to reimplement IsOutflowFacet in
the demo.



------------------------------------------------------------------------

_______________________________________________
DOLFIN-dev mailing list
DOLFIN-dev@xxxxxxxxxx
http://www.fenics.org/mailman/listinfo/dolfin-dev




Follow ups

References