Re: [deal.II] Direct solver matrix.free

2021-10-21 Thread Daniel Arndt
Joss, The MatrixFree framework relies on the operator to be only defined by its action. Accessing individual elements (as direct solvers typically require) falls out of this concept. That is to say that the MatrixFree framework only works well with an iterative solver. Best, Daniel Am Do., 21.

[deal.II] Direct solver matrix.free

2021-10-21 Thread Joss G.
Hi all, I was using the solver PETScWrappers::SparseDirectMUMPS for a sparse matrix approach. Now I am changing the code to a matrix-free implementation and I am not using PETSc anymore, instead, I use LinearAlegebra::distributed::Vector. I am not able to find an equivalent direct solver

[deal.II] deal.II Newsletter #187

2021-10-21 Thread 'Rene Gassmoeller' via deal.II User Group
Hello everyone! This is deal.II newsletter #187. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #12856: Fix typo (proposed by gfcas)

[deal.II] Preconditioner for asymmetric matrix-free style matrix

2021-10-21 Thread 殷承江
Hi all, As described in the title, I am looking for an appropriate preconditioner for an asymmetric matrix implemented in matrix-free style. I've gone though all the matrix-free tutorials and none of them talks about the asymmetric case. (Almost all of them use the Chebyshev smoother with

Re: [deal.II] Dirichlet bc via apply_boundary_values function

2021-10-21 Thread Мария Бронзова
Dear Wolfgang, I believe I answered you privately, but I cannot find my email anywhere. In case you didn"t receive it, I am sending such a test case here again. The testcase is a box porous material model with 8 cells to be able to interpolate BC's and to see the problematics clearlier. The

[deal.II] n components with matrix-free and MPI

2021-10-21 Thread Hermes Sampedro
Dear all, I am having some trouble understanding how to have access to the different components (2 for complex values) in a matrix-free implementation and MPI. I need to have access to each component to perform the different operations to them in the *local_apply()* function. I saw this