> On Sep 3, 2015, at 5:12 AM, David Knezevic <david.kneze...@akselos.com> wrote: > >> 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.
+1 ------------------------------------------------------------------------------ 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
