Jose,
First, I have never used KINSOL so I don't know what it requires. With that
being said, this x.reinit(locally_relevant_dofs, MPI_COMM_WORLD); is
wrong. Either you only need the locally_owned_dofs and you write x.reinit
(locally_owned_dofs, MPI_COMM_WORLD); or you need the
Hello everyone!
This is deal.II newsletter #217.
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:
#14078: Remove remaining deprecations (proposed by masterleinad)
Dear dealii community,
I am working on a gradient enhanced crystal plasticity code where the
displacement field is continuous but the slip fields can be discontinuous.
To this end, I am using a hp::FECollection. For example, in the case
of two crystals with one slip system I have a
Hi Bruno,
sorry for the late response. I get the following error:
An error occurred in line <1469> of file
in function
dealii::IndexSet
dealii::TrilinosWrappers::MPI::Vector::locally_owned_elements() const
The violated condition was:
owned_elements.size() == size()
Additional