You are proposing
> return _var_first_local_df[2*var]
and
> return _var_first_local_df[2*var+1]
??
That is not consistent with the way the bounds are stored. As I mentioned,
consider the three-variable system with unknowns (u,v,w). For simplicity
say there are 150 dofs total, split 50-50-50
> BTW, there is a bug in dof_map.C, which I check out from svn,
> around line 915,
>
>>> if (var == n_vars-1)
> _var_first_local_df.push_back(next_free_dof);
>
> _var_first_local_df[nvar+1] record the offset of vars, like
> [u_begin, v_begin, w_begin, p_begin, p_end],
> but without "if (var == n
Hi,
I'm trying to implement the pressure convection-diffusion preconditioner
for
steady-state Navier-stokes equaions.
I am not sure if I am doing it in an approperate way... Here is the outline:
In general, use the pestc precondition shell matrix, specificly,
* build subblock Matrix manually
* e
