On Fri, Jun 10, 2016 at 10:01 AM, <[email protected]> wrote:
> Dear developers,
>
> I enable the block storage for matrices, however, if I use P2P1 FE type
> for Navier-Stokes equatioins, the block type reduces to non-block type.
>
Hmm.. I don't think this is possible. I just looked at the documentation
for DofMap::has_blocked_representation(), and it says that you can only
have 1 variable group:
/**
* @returns true if the variables are capable of being stored in a blocked
* form. Presently, this means that there can only be one variable group,
* and that the group has more than one variable.
*/
bool has_blocked_representation() const
{
#ifdef LIBMESH_ENABLE_BLOCKED_STORAGE
return ((this->n_variable_groups() == 1) && (this->n_variables() > 1));
#else
return false;
#endif
}
If you want to use blocked storage, you'll possibly need to use some
stabilized equal-order formulation for incompressible flow?
The blocked storage feature was originally implemented for the compressible
Navier-Stokes equations, which use the same interpolation order for all
variables.
--
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
