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. 1
D-39106 Magdeburg
Germany
Email: [email protected] <mailto:[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
