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
,
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
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
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
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
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
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
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
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
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
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
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
12 matches
Mail list logo