On 4/20/23 17:27, Tao Jin wrote:
My question is, is there a way to involve two different dof_handlers when
using the WorkStream class, something like the following:
What you are looking for is an object that looks like a *pair* of iterators
and when you call ++it on it, then it increments
On 4/20/23 02:57, Vidyesh Dapse wrote:
I want a method that can give me a flag indicating "whether a given face of a
particular cell requires ghost data for computation of face integral?"
This can only happen if you are on a locally owned cell and the neighboring
cell (across from the face
Wasim,
we do not know what modifications you made to your program, and so there is
very little anyone can help you with the error. If I had to guess, you are
forgetting to call dof_handler.distribute_dofs() before you start to access
vertex dof indices.
Best
W.
On 4/19/23 12:14, Wasim
Dear Wolfgang,
Thank you so much for your inputs. Indeed, Step-31 is a good example for
what I am looking for. However, I have a follow-up question regarding using
the *WorkStream* class based on TBB during the assembly stage when two
different dof_handlers are involved.
To be specific, in
Hello,
You are using a non-ghosted vector. You need to created a ghosted vector
and they assign the non-ghosted vector to the ghosted one.
Best,
Bruno
On Thursday, April 20, 2023 at 7:07:09 AM UTC-4 Abbas wrote:
> I am trying to querry for the jump in my solution with
>
Hello everyone!
This is deal.II newsletter #250.
It automatically reports recently merged features and discussions about the
deal.II finite element library.
## Below you find a list of recently proposed or merged features:
#15118: PETScWrappers: update KSP support (proposed by stefanozampini)
I am trying to querry for the jump in my solution with
"get_jump_in_function_values"
. My solution vector is PETSc distributed. When I run a code that's like
the one below:
template
void output_jumps::interface_jump()
{
QGauss face_quadrature_formula(degree + 1);
FEInterfaceValues
Hello everyone.
I am working on the discontinuous-Galerkin solver, which is described in
tutorial step-76. In solver, time stepping is performed using the Low
Storage Runge Kutta method. The stage-wise computation is performed using
loop_cell_centric(). The solver uses MatrixFree class. I