← Back to team overview

dolfin team mailing list archive

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

 

On Saturday 20 December 2008 22:59:48 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.

I done some work on the IsOutFlowFacet in the direction previously discussed. 
The function is now mapped to python. Please have a look on the 
implementation. It could probably had been made easier.

The demo runs but the result differ from the cpp. Please have a look.

It also raises a lot of warings during assemble :

  *** Warning: Point belongs to more than one cell, picking first.

Johan

>
> changeset:   5427:36c52ce3d9c0ae6b57642915b52fbb503a2b795f
> user:        Anders Logg <logg@xxxxxxxxx>
> date:        Sat Dec 20 22:39:22 2008 +0100
> files:       demo/pde/poisson1D/python/demo.py
> description:
> Remove rescale=False in 1D demo after bug fix in Viper. Only one Viper bug
> remains.
>
>
> changeset:   5426:b5e8f9ca5863725071a68d4c1bcbfd825def81d0
> user:        Anders Logg <logg@xxxxxxxxx>
> date:        Sat Dec 20 08:58:26 2008 +0100
> files:       site-packages/dolfin/__init__.py
> description:
> Don't import ComplexODE in Python
>
> ----------------------------------------------------------------------
> For more details, visit http://www.fenics.org/hg/dolfin
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev




Follow ups

References