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?
r.
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!
Ondrej
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 sfepyissues, I will summarize what was done recently on the bleeding edge:
 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).
 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.
participants (2)

Ondrej Certik

Robert Cimrman