On Thu, Sep 3, 2015 at 4:18 AM, Julian Andrej <j...@tf.uni-kiel.de> wrote:
> I tried Davids tip with the CondensedEigenSystem. For the scalar > valued problem of a laplace operator it works very good, but i still > have no correct solution for the Stokes equations. > > It seems i'm assembling the Mass matrix in a wrong way. Using the > notation of the example, is it correct, that the Mass matrix has > contributions from the coupling of uu,vv,uv,vu? So i tried using > > Me(i,j) += JxW[qp]*(phi[i][qp]*phi[j][qp]); > > in the Kuu, Kvv, Kuv, Kvu loops, without success. > > Any help appreciated :) > It sounds like you aren't sure about the weak formulation of your PDE in this case, so you need to figure that out first... Having said that, normally the mass matrix would not have off-diagonal blocks, so I suggest you try removing the Kuv and Kvu blocks. David ------------------------------------------------------------------------------ Monitor Your Dynamic Infrastructure at Any Scale With Datadog! Get real-time metrics from all of your servers, apps and tools in one place. SourceForge users - Click here to start your Free Trial of Datadog now! http://pubads.g.doubleclick.net/gampad/clk?id=241902991&iu=/4140 _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
