Hello Dealii mates, I am working on the schur complement for eigenvalue problems, following step - 55. My stiffness block matrix is like this:
*| A B.T |* *| B C |* and my mass block matrix is structured: *| M 0 |* *| 0 0 |* To obtain A_condensed, inversion of C is required. I have constructed the linear operators for C_inv and A_condensed using the inverse_operator and schur_complement functions based on the documentation here (https://dealii.org/developer/doxygen/deal.II/group__LAOperators.html#ga76acca911f21089cd3bb385d20ccc995). However, SLEPcWrappers eigensolvers, such as SLEPcWrappers::SolverKrylovSchur and SLEPcWrappers::::SolverGeneralizedDavidson, only accept PETScWrappers::MatrixBase rather than linear operators as said in their documentations. Is there an easy way to form a PETScWrappers::MPI::SparseMatrix from the linear operator? I am trying not to calculate the schur_complement explicitly, however, it seems i may have to. I would really appreciate it for any suggestions. Best regards, Yuchen Liu -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion visit https://groups.google.com/d/msgid/dealii/f3c3f566-65b3-4502-a81c-895e19e3fb64n%40googlegroups.com.
