> I think it is rather an issue of implementation to get the terms > contributing to the off-diagonal blocks in the mass matrix right. > Essentially its the velocity basis functions in an inner product with the > bases coming from the pressure space.
What Markus is saying that the correct term can not have the form (phi_u, phi_p) since the first factor is a vector and the second is a scalar. The result of the integration is a vector, but in a bilinear form you need a scalar as a result. It could be that you need to take a *particular component* of phi_u, but that there is no telling without knowing where the problem came from. As it stands, the bilinear form you should doesn't make any sense because it adds scalars and vectors together. Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
