Elmer Discussion Forum

Hello;
mmg has too be compiled with the -fPIC flags;

You have to set
CMAKE_CXX_FLAGS=-fPIC
CMAKE_C_FLAGS=-fPIC
when using cmake to configure the compilation of mmg.

In general you will have to set this flag is most external librairies that you want to link with Elmer

best regards
Fabien

Hello, The Tete Rousse example during the beginner course is a good example; the cavity problem is the same as the grounding line (an ice shelf is a big open cavity ;) ). You also have a 2D example in the elmerice tests (under the elmer sources elmerice/Tests/GL_MISMIP). With Stokes you have to solv...

Hello;
if the iterative solver is not converging you can try to use a direct solver (umfpack in serial or Mumps in parallel) this should always work (at least if you have enough memory).

If the direct method is too slow you may still try to use another iterative solver and/or preconditionner.

Hi Thomas, Yes I agree that core issue is with the GetEquation in the USF and this should be changed. However, this does not crash for me even if I don't have a body for the bed; and I think this is why it has been un-noticed so far, espacially when the Aruco test case was designed. So to play safe ...

Hello,
code is ok in getFrictionLoads if Equation is not assigned, which should be the case if element do not belong to a body with a given equation?

I think this could be fixed in GetEquationId (DefUtils.F90) which should return with Found=.False. if body_id < 1 ?

Cheers
Fab

Hello; I think it is also sensitive to the stabilisation method. By default the Free surafce solver uses the Bubbles; SUPG stabilisation is also possible; However in this case the stabilisation parameter is optimal in steady-state. There is a lot of litterature on this subject and possible stabilisa...

Dear ElmerIce users, Following important changes on how Dirichlet Conditions are handle in Elmer v8.3, we found that the conditional dirichlet condition for the normal velocity in the grounded part of the ice sheet was not properly released during grounding line retreat. This has now be fixed in the...

Hello Rupert, thanks for the report and sorry for the late reply. At Grenoble we are running all the Elmer/Ice tests, after each code revision. I try to keep an eye on the results as much as I can. I have updated most of the tests with failure. On the 52 elmerice tests (ctest -L elmerice); There is ...

ok.
Thanks!

Fabien

Hi Peter; Yes a call to DefaultInitialize() will also solve our problem. Currently we have this problem with the Navier-Stokes solver, I see that it is calling InitializeToZero(). Replacing the InitializeToZero() line by CALL DefaultInitialize() will sove our problem. I see that there is few other s...