Search found 47 matches
- 09 Apr 2024, 10:33
- Forum: Bug reports
- Topic: Bug in SaveScalars in parallel?
- Replies: 2
- Views: 291
Re: Bug in SaveScalars in parallel?
Hello, yes as today (rev: 71005ec) your test returns wrong results in parallel; issue is related to changes introduced in rev:e396adc on Feb 21rst 2023; The test in the BoundaryStatics routine in SaveScalars.F90 to skip nodes not owned by the current PE is : IF( Var % Solver % Matrix % ParallelInfo ...
- 08 Mar 2021, 14:58
- Forum: Elmer/Ice
- Topic: [SOLVED] AIFlow - Segmentation Fault with Calculate Loads
- Replies: 3
- Views: 5149
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.
- 08 Mar 2021, 14:55
- Forum: Elmer/Ice
- Topic: [SOLVED] AIFlow - Segmentation Fault with Calculate Loads
- Replies: 3
- Views: 5149
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...
- 21 Jan 2021, 12:04
- Forum: Elmer/Ice
- Topic: Elmer/Ice can't pass 'make'
- Replies: 3
- Views: 5864
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
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
- 02 Dec 2020, 10:00
- Forum: Elmer/Ice
- Topic: Where to start on a transient grounding line model
- Replies: 2
- Views: 5468
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...
- 10 Nov 2020, 19:41
- Forum: Elmer/Ice
- Topic: Problems with basal sliding when applying Porous Solver
- Replies: 7
- Views: 8418
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.
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.
- 02 Jul 2020, 16:23
- Forum: Elmer/Ice
- Topic: getFrictionLoads User Function resulting in segmentation fault error
- Replies: 8
- Views: 11094
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 ...
- 02 Jul 2020, 15:09
- Forum: Elmer/Ice
- Topic: getFrictionLoads User Function resulting in segmentation fault error
- Replies: 8
- Views: 11094
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
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
- 04 Jun 2020, 10:53
- Forum: Elmer/Ice
- Topic: Spiky Ice Sheet Surface in transient flow line model
- Replies: 6
- Views: 8195
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...
- 10 Apr 2018, 11:04
- Forum: Elmer/Ice
- Topic: GL dynamics and Stokes in Elmer v8.3
- Replies: 2
- Views: 6250
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...