[deal.II] Re: How to specify finite element component in MatrixFree?

Thank you peter. It looks awesome! I've got a few question about it. 1. In your github code, you actually assemble a block matrix system(in MatrixFree) but didn't renumber dofs as done in most Stokes examples? There's no multi dof_handlers. The solution in HeatProblem is a non-complex

[deal.II] Re: How to specify finite element component in MatrixFree?

Hi The block version of MGTransferGlobalCoarsening is MGTransferBlockGlobalCoarsening (see https://www.dealii.org/developer/doxygen/deal.II/classMGTransferBlockGlobalCoarsening.html). I don't know the details of step-29. But what I would recommend you to do is to take a look at

[deal.II] Re: Apply boundary integral to diagonal in MatrixFree

Hi, you are right such a function is unfortunately missing in deal.II. It has been on my TODO list for a long time but I didn't get around to implement it, since I haven't had to use DG for a while and in most of my applications I can move the BCs to the RHS. Setting up the matrix is not

Re: [deal.II] SLEPcWrappers::SolverKrylovSchur Zero pivot problem

On 10/31/22 22:12, Sebastian D.-R. wrote: I started to work a few days ago on an eigensolver for linear elasticity with the use of the SLEPc wrappers in deal.II  (based on step-8 and step-36), and I came across the same issue Leonhard described above. After some debugging I realised the

Re: [deal.II] error while running step_17

I have received following error while running step_17. How do I solve the problem? Himal, the error suggests that the version of gmsh you have on your machine is missing a function. Do any of the other (non-MPI) tutorials run? Best WB --

[deal.II] Re: How to specify finite element component in MatrixFree?

I tried to modify my code that FEEvaluation is constructed without matrix_free_data as in the above .cc in dealii/test, but get following error information: "FEEvaluation was initialized without a matrix-free object. Integer indexing is not possible". Guess the only method to do

[deal.II] How to specify finite element component in MatrixFree?

In test/matrix-free/assemble_matrix_02.cc a Stokes matrix is assembled in MatrixFree form, and DoFs correspond to velocity and pressure are called by FEEvaluation<...,dim>