On Mon, Jun 13, 2016 at 11:02 AM, li.luo <[email protected]> wrote: > Thanks. > I thought of another way may work, that is, to use p2p2 for velocity and > pressure, then "disregard" those dofs belong to 2nd order elements, and > set the related rows in matrix to be identities. That I use a larger > problem, but v and p can be one variable group. But I don't known which > dofs / basis functions belong to 2nd order type. >
This isn't a good solution for a number of reasons. Aside from using extra memory, I'm not sure if it even makes sense algebraically to do this. It is also going to produce a non-sensical solution in any of the normal visualization formats (are you simply setting these pressure dofs to zero?) It is also going to be inefficient to zero the rows corresponding to second-order DOFs, which you will have to do in a separate step after assembly. Better to investigate the possibility of using a stabilized P1P1 scheme... what is your actual goal though -- just to use blocked storage? -- John ------------------------------------------------------------------------------ What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic patterns at an interface-level. Reveals which users, apps, and protocols are consuming the most bandwidth. Provides multi-vendor support for NetFlow, J-Flow, sFlow and other flows. Make informed decisions using capacity planning reports. https://ad.doubleclick.net/ddm/clk/305295220;132659582;e _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
