[deal.II] get_function_values giving different results depending on the number of MPI ranks

Hi everyone, I've been trying to write a parallel distributed solution to a nonlinear PDE using Newton's method. It runs fine on one rank, but when I use two the solver fails to converge on the second iteration. I've tracked the difference (between one rank and two) down to the `get_function_va

[deal.II] Re: FullMatrix direct solver

Hello, You can use the invert function to get the inverse of the matrix and then you can use vmult with the new matrix. Best, Bruno On Wednesday, January 26, 2022 at 8:19:46 AM UTC-5 jossg...

[deal.II] FullMatrix direct solver

Good afternoon, I have a basic question. I can see that for sparse matrix the function SparseDirectUMFPACK.vmult() solves Ax=b. I am trying now to solve a system with full matrix. I can see that the function FullMatrix.vmult() does not apply the inverse as in the sparse analogy. What is the