On 8/25/21 9:37 PM, Anton Ermakov wrote:
>>> You can do that, but I believe that it's also possible to work with the
original system where one of the two matrices is only semidefinite.
I have redone it with the original system, but it seems that PETSc is
complaining about missing diagonal
>>> *That might not easily be possible with our PETSc interfaces.*
*PETSc has the concept of MatShell, which is in essence what
deal.II'sLinearOperator is: A class that implements a matrix-vector
operation,without giving access to individual matrix elements. You could
trywhether you can
On 8/25/21 11:53 AM, Anton Ermakov wrote:
I I¯ A 0 ¯I *eigenvalue + I¯ 0 B ¯ I I * eigenvector =0
I I_ 0 0 _I I_ C D _I I
¯ ¯
where eigenvector consists of a vector
Dear Deal.II community,
I am working on implementing a vector eigenproblem in deal.ii, starting
from the step-36 example. The question I am studying is the normal
oscillation of a fluid planet. The eigenproblem is of this type:
_