On Thursday, 7 March 2019 23:35:15 UTC-5, Wolfgang Bangerth wrote:
>
> On 3/7/19 7:45 AM, Bruno Blais wrote:
> >
> > However, my GMRES stops very quickly after a certain number of newton
> > iteration with the following : AztecOO::Iterate error code -4: GMRES
> > Hessenberg ill-conditione
On 3/7/19 7:45 AM, Bruno Blais wrote:
>
> However, my GMRES stops very quickly after a certain number of newton
> iteration with the following : AztecOO::Iterate error code -4: GMRES
> Hessenberg ill-conditioned
This is surprising and suggests that something is wrong with the matrix, not
t
Dear Wolfgang,
Sorry for an additional question.
I have tried to implement by myself something similar to what is done in
ASPECT (I guess?, it is a bit harder for me to grasp ASPECT concepts
because of the introspection and etc.).
However, I seem to have failed miserably at one point.
My understa
> A quick question. I think I understand what is done in
> https://github.com/geodynamics/aspect/blob/master/source/simulator/helper_functions.cc#L1041
> Weirdfully, I found the "pickaxe" version easier to understand (Love the
> comments by the way, this is awesome)
:-) I suspect that Timo gets
Dear Wolfgang,
A quick question. I think I understand what is done in
https://github.com/geodynamics/aspect/blob/master/source/simulator/helper_functions.cc#L1041
Weirdfully, I found the "pickaxe" version easier to understand (Love the
comments by the way, this is awesome)
A question I have is r
Dear Bruno and Wolfgang,
Thank you for your answers.
I believe Wolfgang's answer is exactly what i had in mind (but said in
clear words...). I will look at the Aspect code and try to port that to
mine.
Thank you for the very detailed answer.
Best
Bruno
On Tuesday, 5 March 2019 10:08:29 UTC-5,
> One of the issue of my system is that the pressure is defined up to a
> constant. On coarse mesh this does not affect the GMRES solver. However, on
> finer mesh, it seems that the GMRES Solver is greatly affected by this
> near-singularity of the matrix system.
> I have often read in the lit
Hello everyone,
I am currently solving a GLS stabilized form of the Navier-Stokes equation
using DEALII.
The residual of the system looks similar to the regular Incompressible
Navier-Stokes, except that a stiffness matrix that is dependent on the
element size is added to the P-P block.
I have a