[deal.II] residual for a static condensation problem

2024-10-14 Thread Yann Jobic
Dear Dealii group, I am currently learning Dealii (and i really like this lib). I would like to bench it in order to compare its performance with an in house FEM code at my lab. I want to solve the stokes equation with a static condensation of the pressure over the velocity dof. I started fr

Re: [deal.II] residual for a static condensation problem

2024-10-15 Thread Yann Jobic
, Yann Le 10/15/2024 à 5:02 AM, Wolfgang Bangerth a écrit : On 10/14/24 14:03, Yann Jobic wrote: However, i don't know how to compute the residual ! It should be easy, but i still haven't succeeded. If I extract the different submatrices of the system matrix at the element level, i

Re: [deal.II] residual for a static condensation problem

2024-10-16 Thread Yann Jobic
sense ? It could still be confusing. Many thanks, Best regards, Yann Le 15/10/2024 à 19:04, Wolfgang Bangerth a écrit : On 10/15/24 01:19, Yann Jobic wrote: My next step is to compute that way (static condensation of the pressure) the navier-stokes equations with a Newton non linear solver, as

Re: [deal.II] Printing PETScWrappers::MPI::SparseMatrix to a single file

2024-11-05 Thread Yann Jobic
Hello Sparsh, I use the function "block_write" for sequential matrix/vectors, which writes a binary file. It's way better than text files. I used the function this way : std::ofstream out5("MatrixSystem.bin"); system_matrix.block_write(out5); std::ofstream out6("RHS.bin"); sy

[deal.II] access coordinates of quadrature points

2024-10-28 Thread Yann Jobic
Hello, I'm trying to access the coordinates of the quadrature points. I'm doing : QGauss quadrature_formula(degree + 1); FEValues fe_values(fe, quadrature_formula, update_values | update_quadrature_points | update_JxW_values | update

[deal.II] Re: access coordinates of quadrature points

2024-10-28 Thread Yann Jobic
I found my mistake. The written code works fine, but i cannot inspect the values directly in the debugger. That confused me. Sorry for polluting the mailing list. Shall both email erased ? Regards, Yann Le 10/28/2024 à 6:36 PM, Yann Jobic a écrit : Hello, I'm trying to access the coordi

Re: [deal.II] Periodic boundaries

2024-12-02 Thread Yann Jobic
tanding how an error comes about gets you 90% the way towards understanding how to fix the code (in your own application, or perhaps in deal.II itself). Best  Wolfgang On 11/28/24 11:37, Yann Jobic wrote: ** Caution: EXTERNAL Sender ** Hello, I'm trying to define periodic constrain

[deal.II] parallel periodic boundary conditions

2024-12-06 Thread Yann Jobic
Hello everybody, I have a working NS parallel solver, with AMR. It's working fine (i.e. good convergence rate of the error in L2 norm) with imposed velocity and outflow. I now need to use parallel periodic boundary conditions, with volumetric forces. I found that: 1) It's working perfectly

[deal.II] PETSC SparseDirectMUMPS and factorization reuse

2025-02-12 Thread Yann Jobic
Hello all, I want to reuse the factorization of a PETSC/MUMPS solve for different right-hand side vectors. I found the following discussion: https://groups.google.com/g/dealii/c/Sw0wHKEFuiA and the issue: https://github.com/dealii/dealii/issues/2869 However, it's quite old. Is there a way to

Re: [deal.II] Regarding the run time of my code

2025-04-15 Thread Yann Jobic
Hello Md Mahmudul, I'm running also NSE, and with higher degree of freedom. It's not taking that much time. In order to understand what is going on, i've got some questions : 1) are you in debug compilation mode ? 2) it seems that you are solving a linear system with a direct method. Is there

Re: [deal.II] Regarding the run time of my code

2025-04-15 Thread Yann Jobic
ation number its fast. In that case, should I use an iterative solver? Md Mahmudul Islam On Tue, Apr 15, 2025, 2:30 AM Yann Jobic <mailto:yjo...@gmail.com>> wrote: Hello Md Mahmudul, I'm running also NSE, and with higher degree of freedom. It's not taking tha

Re: [deal.II] Installing dealii with METIS (for compiling a code gallery program)

2025-03-06 Thread Yann Jobic
Hello Sudip, If i may answer, i'm compiling petsc apart of dealii, and when compiling dealii, i'm using those kind of cmake options : -DHDF5_DIR=$PETSC_DIR -DPETSC_DIR=$PETSC_DIR -DPETSC_ARCH="" -DSLEPC_DIR=$PETSC_DIR -DP4EST_DIR=$PETSC_DIR -DDEAL_II_WITH_P4EST=ON -DMETIS_DIR=$PETSC_DIR -DDEAL