On 12/01/2017 06:50 AM, Timo Heister wrote:
They're always going to be there because we keep constrained nodes in the
linear system.
If we modify the way ConstraintMatrix operates, we could work around
this though:
1. Without rescaling those equations (as we do inside ConstraintMatrix
right now
2017-12-01 8:50 GMT-05:00 Timo Heister :
> 3. Did we rip out the support for removing the constrained entries
> from the matrix completely?
That was my plan at first using ConstraintMatix::condense but according to
the documentation, the constrained entries are not removed.
Best,
Bruno
--
The
On Friday, December 1, 2017 at 2:50:38 PM UTC+1, Timo Heister wrote:
>
> > They're always going to be there because we keep constrained nodes in
> the
> > linear system.
>
> If we modify the way ConstraintMatrix operates, we could work around
> this though:
>
but there is really no need, is
> They're always going to be there because we keep constrained nodes in the
> linear system.
If we modify the way ConstraintMatrix operates, we could work around
this though:
1. Without rescaling those equations (as we do inside ConstraintMatrix
right now) the spurious EV would all be equal to 1 a
On 11/30/2017 05:47 PM, Bruno Turcksin wrote:
In step-36, there is an explanation on how Dirichlet boundary conditions
introduce spurious eigenvalues because some dofs are constrained. However,
there is no mention of hanging nodes. So I am wondering if I can treat them as
shown for the Dirich
Thanks Denis!
2017-11-30 17:31 GMT-05:00 Denis Davydov :
> Hi Bruno,
>
> AFAIK, there is a simple solution: make initial vector (or subspace)
> perpendicular to those constrained entries.
> That is, if you do Lancoz, set random initial vector and then zero out
> constrained DoFs.
> Then being Kry
Hi Bruno,
AFAIK, there is a simple solution: make initial vector (or subspace)
perpendicular to those constrained entries.
That is, if you do Lancoz, set random initial vector and then zero out
constrained DoFs.
Then being Krylov-based method it should form subspaces {x, Ax, A^2x,...}
orthogo
