Hi, thank you for pointing this out what I obviously meant was "...phi_u[j[ will be of Tensor element and phi_p[i] a scalar..."
This comes from the discretization of a control problem where a Q1 control "meets" a Q2 test function. Sorry about the confusion. Best, Martin On Tue, Feb 1, 2011 at 2:13 PM, Markus Bürg <[email protected]> wrote: > Hello Martin, > > perhabs I missed a fact, but why should phi_u[i] and phi_u[j] have > different dimensions? What is the mathematical setting of your problem? > > Best Regards, > Markus > > > > Am 01.02.11 13:34, schrieb Martin Stoll: > > Hi all, > > thank you again for all the previous help. > > I am looking for a(n elegant) solution to construct a mass matrix for a > vector valued problem of the following form > [<phi_u,phi_u> <phi_u,phi_p>; <phi_p,phi_u> <phi_p,phi_p>] where the > diagonal block are just the velocity and pressure mass matrices. I also need > the mixed terms and > I realized that I couldn;t use the mass_matrix_creator. I naively simply > put > local_mass(i,j) += (phi_u[i] * phi_u[j] > + phi_u[i] * phi_p[j] > + phi_p[i] * phi_u[j] > + phi_p[i] * phi_p[j]) > * fe_values.JxW(q); > until I realized that phi_u[j[ will be of Tensor element and phi_u[i] a > scalar. > > I was hoping someone could suggest an elegant way (short?) for the above > assembly as I somewhat stuck here and cannot seem to come up with a sensible > way to do this. > > Thank you for your help. > > Best, > Martin > > -- > *Martin Stoll* > *Postdoctoral Research Fellow* > > Computational Methods in Systems and Control Theory Max Planck Institute for > Dynamics of Complex Technical Systems > > Sandtorstr. 1D-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 > > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii > > -- *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
