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. 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
