Ondrej Certik wrote:

On Wed, Oct 8, 2008 at 11:46 AM, Robert Cimrman cimr...@ntc.zcu.cz wrote:

...

Hi again,

the code is ready to unify the terms corresponding to the same integral (weak) forms. The unified term names are obvious to choose except one: how to call the unification of d_div, dw_div, d_grad, dw_grad? they are all defined as $\int_{\Omega} q\ \nabla \cdot \ul{v}$ but differ in type of variables (virtual vs. unknown). I propose dw_stokes, what do you think?

Sounds good!

As not all members of this list subscribe also sfepy-issues, I will summarize what was done recently on the bleeding edge:

1. term unification: new terms dw_stokes, dw_biot, dw_biot_th replaced

groups of related coupling terms, see input/navier_stokes.py, input/biot.py files.

Example - Stokes problem:

Instead of: "dw_div_grad.i2.Omega( fluid.viscosity, v, u )- dw_grad.i1.Omega( v, p ) = 0" "dw_div.i1.Omega( q, u ) = 0"""

we have now

"dw_div_grad.i2.Omega( fluid.viscosity, v, u )- dw_stokes.i1.Omega( v, p ) = 0" "dw_stokes.i1.Omega( u, q ) = 0"""

i.e. dw_stokes can behave as dw_grad, dw_div, depending on its arguments (which is an unknown and which is a test function in the FE sense).

2. time dependent problems: time derivatives of variables in term

arguments, see input/time_poisson.py

Example: transient heat conduction (no sources) "dw_mass_scalar.i1.Omega( s, dT/dt ) + dw_laplace.i1.Omega( coef.val, s, T ) = 0"

dT/dt is automatically replaced by the backward difference.

I have updated sfepy_manual.pdf on the google site to reflect those changes (though more decription is still needed).

cheers, r.