Re: [deal.II] Multigrid Preconditioner for HDG

2023-05-04 Thread Wolfgang Bangerth
On 5/4/23 21:31, Aditya Karthik Saravanakumar wrote: Note: dofh_uhat is a DoFHandler object corresponding to the FE_FaceQ finite element space. _ERROR MESSAGE_ An error occurred in line <374> of file in function     const dealii::FullMatrix& dealii::FiniteElement<,

Re: [deal.II] Multigrid Preconditioner for HDG

2023-05-04 Thread Aditya Karthik Saravanakumar
That's fair. To be more specific, I'm essentially trying to incorporate the multigrid code from step 16 into an HDG solver that's very similar to step 51. My code fails when I try to call the build(DoFHandler) function on my MGTransferPrebuilt object as shown below (similar to what's done in

Re: [deal.II] Multigrid Preconditioner for HDG

2023-05-04 Thread Wolfgang Bangerth
On 5/4/23 19:26, Aditya Karthik Saravanakumar wrote: I was able to set up all the MGLevelObjects but the MGTransferPrebuilt class isn't able to generate the prolongation matrices. Aditya, I don't know enough about multigrid to help, but you might get better answers if you are specific about

[deal.II] Multigrid Preconditioner for HDG

2023-05-04 Thread Aditya Karthik Saravanakumar
Hello everyone, I've been trying to set up a geometric multigrid preconditioner (following step 16) for an HDG Poisson solver. I was able to set up all the MGLevelObjects but the MGTransferPrebuilt class isn't able to generate the prolongation matrices. Has anyone been successful in

[deal.II] deal.II Newsletter #252

2023-05-04 Thread 'Rene Gassmoeller' via deal.II User Group
Hello everyone! This is deal.II newsletter #252. 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: #15177: Make some variables in a few tests 'const'. (proposed by

Re: [deal.II] Understanding MeshWorker::mesh_loop order with adaptive refinement

2023-05-04 Thread Timo Heister
I think this is just a typo in the documentation exactly like you pointed out. It would be great to see that fixed. Are you willing to do that? I am not sure I understand your problem you had to work around. If you go into detail we might be able to figure out the issue (if there is one inside

Re: [deal.II] Finding AVX instructions always fails at cmake configure step

2023-05-04 Thread blais...@gmail.com
Dear Martin, Now I feel kind off dumb. This is what I was doing to compile an AVX test on the side, did not think to include it this way in the CMAKE procedure. Thanks! Worked very well and now it detects my AVX2 instructions :). Best Bruno On Thursday, May 4, 2023 at 8:57:43 a.m. UTC-4 Martin

Re: [deal.II] Finding AVX instructions always fails at cmake configure step

2023-05-04 Thread Martin Kronbichler
Dear Bruno, You need to specify appropriate CXX flags. For example, I always use   -D CMAKE_CXX_FLAGS="-march=native" on my local machines (that run the code I'm compiling) or   -D CMAKE_CXX_FLAGS="-march=icelake-server" on clusters to be sure the right code gets generated. For example GCC

[deal.II] Finding AVX instructions always fails at cmake configure step

2023-05-04 Thread blais...@gmail.com
Dear all, Hope you are well. I have tried a few ways to enable AVX instructions (which I think my 12th Gen Intel(R) Core(TM) i9-12900K) should support, yet I am always unable to do so. At the Cmake configuration step, the AVX steps all fails whereas the SSE ones succeed. How can I remedy this

[deal.II] Poor parallel efficiency of FEValues::reinit() with Workstream

2023-05-04 Thread Jonas Knoch
Dear deal.II community, I have implemented a multiscale finite element problem using deal.II. It includes PDEs on different levels, namely a macroscopic/effective PDE and so-called "cell problems" (not to be confused with cells of the grid), which have to be solved in each quadrature point of