> 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
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
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
> 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
Bruno,
On Tuesday, March 5, 2019 at 9:00:13 AM UTC-5, Bruno Blais wrote:
> Are there any examples in DEALII where a Poisson problem is solved but
> with strictly Neumann boundaries? That would be a very similar problem that
> could guide me in the right direction.
>
Yes, take a look at
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
As initial debug approach, I used the function local_assemble_system
instead of assemble_cell_terms, which worked (left the other functions
alone). Then I tried to assemble just the boundary terms, by adding
local_assemble_boundaries, uncommenting assemble_Dirichlet_boundary_terms()
and
Dear Bruno,
> Thank you for the rapid and very detailed answer. This makes this community
> so great.
Bruno T. has given you far more insight into the inner workings of the Trilinos
linear algebra packages than I’d have been able to (thanks :-) ), so I’m really
glad that we were able to help