Search found 46 matches

by fgillet
08 Mar 2021, 14:58
Forum: Elmer/Ice
Topic: [SOLVED] AIFlow - Segmentation Fault with Calculate Loads
Replies: 3
Views: 4450

Re: AIFlow - Segmentation Fault with Calculate Loads

Also, the bulkValues contains the stiffness matrix before applying the Dirichlet conditions; while the Values is the stiffness matrix after applying Dirichlet; so that computing the residual tells you what you have to add to your system to comply with your Dirichlet conditions.
by fgillet
08 Mar 2021, 14:55
Forum: Elmer/Ice
Topic: [SOLVED] AIFlow - Segmentation Fault with Calculate Loads
Replies: 3
Views: 4450

Re: AIFlow - Segmentation Fault with Calculate Loads

Hi Kevin, following some discussions on the elmerice slack, in the channel beginner-course-nov2020; I have updated the codes so that they call all the required Elemer default routines to compute the residual. This is now fixed in the elmerice branch; https://github.com/ElmerCSC/elmerfem/commit/b35e0...
by fgillet
21 Jan 2021, 12:04
Forum: Elmer/Ice
Topic: Elmer/Ice can't pass 'make'
Replies: 3
Views: 5081

Re: Elmer/Ice can't pass 'make'

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
by fgillet
02 Dec 2020, 10:00
Forum: Elmer/Ice
Topic: Where to start on a transient grounding line model
Replies: 2
Views: 4770

Re: Where to start on a transient grounding line model

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...
by fgillet
10 Nov 2020, 19:41
Forum: Elmer/Ice
Topic: Problems with basal sliding when applying Porous Solver
Replies: 7
Views: 7472

Re: Problems with basal sliding when applying Porous Solver

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.
by fgillet
02 Jul 2020, 16:23
Forum: Elmer/Ice
Topic: getFrictionLoads User Function resulting in segmentation fault error
Replies: 8
Views: 10010

Re: getFrictionLoads User Function resulting in segmentation fault error

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 ...
by fgillet
02 Jul 2020, 15:09
Forum: Elmer/Ice
Topic: getFrictionLoads User Function resulting in segmentation fault error
Replies: 8
Views: 10010

Re: getFrictionLoads User Function resulting in segmentation fault error

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
by fgillet
04 Jun 2020, 10:53
Forum: Elmer/Ice
Topic: Spiky Ice Sheet Surface in transient flow line model
Replies: 6
Views: 7793

Re: Spiky Ice Sheet Surface in transient flow line model

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...
by fgillet
10 Apr 2018, 11:04
Forum: Elmer/Ice
Topic: GL dynamics and Stokes in Elmer v8.3
Replies: 2
Views: 5877

GL dynamics and Stokes in Elmer v8.3

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...
by fgillet
10 Apr 2018, 10:48
Forum: Elmer/Ice
Topic: Elmer/Ice tests
Replies: 2
Views: 6056

Re: Elmer/Ice tests

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 ...