Hi Wolfgang,
I'm working on a hypothesis that some dG degrees of freedom can be
constrained without compromising the stability for e.g. advection
diffusion equations, so I want a mesh that has some dG and some cG
elements, which is quite unusual.
Hope this helps,
John
John Chapman
www.maths.dur.ac.uk/~tphj28
On 08/09/10 12:58, [email protected] wrote:
I'm adapting tutorial step-27 to use discontinuous finite elements.
When assembling the local matrices to the system matrix step-27 uses
constraints.distribute_local_to_global (ui_vi_matrix, cell_rhs,
local_dof_indices,
system_matrix,
system_rhs);
Beyond all the other questions that have been raised on this list,
let me ask what constraints you have at all? If you use discontinuous
elements, there shouldn't be any and in that case the
distribute_local_to_global function does a straight copy of the local
matrix and vector entries into the global objects in the same way
as, for example, steps 3...6. You won't need the constraints object
for anything at all in that case.
Best
W.
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
