Hi Timo,

thanks for the suggestion. One of the things I might not have pointed out is
that I want to assemble a Stokes system once and then get N_T contributions
to the right-hand-side depending on the different boundary condition for
each time step.

I realize that step-22 uses

constraints.distribute_local_to_global
<http://www.dealii.org/developer/doxygen/deal.II/classConstraintMatrix.html#ab1ae692d24f925661364552c9a4606ae>
(local_matrix, local_rhs,
                                               local_dof_indices,
                                               system_matrix, system_rhs);


but in my case I have one system_matrix and N_T system_rhs.

Any help is greatly appreciated.

Cheers,
Martin

On Wed, Jan 12, 2011 at 11:46 AM, Timo Heister <
[email protected]> wrote:

> Hi Martin,
>
> > For non-time-dependent problems I would simply use
> apply_boundary_conditions
> > to get the right-hand-side contributions. Is there any way to do this
> > efficiently for N_T (number of time-steps) boundary conditions?
>
> I would advice to use the ConstraintMatrix class for that, see step-22.
>
> --
> Timo Heister
> http://num.math.uni-goettingen.de/~heister<http://num.math.uni-goettingen.de/%7Eheister>
>



-- 
*Martin Stoll*
*Postdoctoral Research Fellow*

Computational Methods in Systems and Control Theory
Max Planck Institute for Dynamics of Complex Technical Systems
Sandtorstr. 1
D-39106 Magdeburg
Germany


Email: [email protected]
URL : http://www.mpi-magdeburg.mpg.de/people/stollm
Tel :+49 391 6110 384
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to