Yes, that part now works. Thank you, Wolfgang.
Anton.
On Wednesday, March 17, 2021 at 3:13:17 PM UTC-7 Wolfgang Bangerth wrote:
> On 3/17/21 10:41 PM, Anton Ermakov wrote:
> > Thanks for the reply. I think I got it. Basically, I had to manually
> place the
> > local matrices into the global ma
On 3/17/21 10:41 PM, Anton Ermakov wrote:
Thanks for the reply. I think I got it. Basically, I had to manually place the
local matrices into the global matrix, without taking into account the hanging
node constraints, which are later taken care of by
/constraints.condense(system_matrix, system
Thanks for the reply. I think I got it. Basically, I had to manually place
the local matrices into the global matrix, without taking into account the
hanging node constraints, which are later taken care of by
*constraints.condense(system_matrix, system_rhs);*
In the case of the mass matrix, I ha
On 3/16/21 6:12 AM, Anton Ermakov wrote:
I attach the minimally changed Step-26 code to illustrate the problem. It is
different from the original step-26 in three important places:
1) in the run function, the mesh is refined in a corner
2) in the setup_system the mass and Laplace matrices ar
