On 12/20/23 04:35, Simon wrote:
But it is not clear to me how to compute the norms of velocity and pressure in
that case.
So far, I used a Block indices object as returned by
DoFTools::count_dofs_per_fe_block()
since I compute the l2-norm only of the unconstrained DoFs.
That said, how do you a
Ernesto,
Yes, something is wrong with documentation. I don't know why the functions
don't show up.
I am not sure what happened with double_contract being deprecated. Maybe it
was a different version maybe the function was reinstated.
Best,
Bruno
On Tuesday, December 19, 2023 at 1:29:54 AM U
Thank you so much! I did all of this and it works now. I had tried the dof
index before, but I think the output to a vtu file was wrong, so now I have
fixed that too.
Much appreciated,
Giselle
On Monday, December 18, 2023 at 4:55:25 p.m. UTC-5 Wolfgang Bangerth wrote:
>
> On 12/18/23 13:23, Gi
Timo,
" What I have done in several codes (and what is used in ASPECT), is to
always create block matrices and vectors but to put all DoFs into a single
block when a direct solver is used. This decision can be made at rintime.
Inside the solve routine you can just access block_matrix.block(0,0)
