On Mon, Jun 13, 2016 at 12:00 AM, <li....@siat.ac.cn> wrote:
> Dear developers,
>
> I created a 2nd order finite element, TET or CUBE, however for some
> reasons, I want to identify the 2nd order of dofs in this element.
> For example, the continuity equaiton in N-S equations, if I set 2nd order
> for pressure, but I want to disable the 2nd order dofs, which "i" is belong
> to 2or dofs? Thank you!
>
> for (unsigned int i=0; i<khi.size(); i++)
> {
> for (unsigned int j=0; j<phi.size(); j++)
> {
> Kpu(i,j) += -JxW[qp]*khi[i][qp]*dphi[j][qp](0);
> Kpv(i,j) += -JxW[qp]*khi[i][qp]*dphi[j][qp](1);
> }
> }
>
There's definitely a way to simply "identify" the 2nd-order dofs, but I
don't think this is going to actually do what you want, i.e. "disable" them
in the system of equations.
Maybe you can let us know what you are actually trying to do, and we could
give some suggestions.
--
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
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
