Search found 46 matches

by fgillet
04 Nov 2016, 15:52
Forum: Elmer/Ice
Topic: Transfer basal drag from one mesh to another mesh
Replies: 14
Views: 17944

Re: Transfer basal drag from one mesh to another mesh

Yes best with "savedata" "saveline" to save values on a boundary;
But it is always possible to save csv files with paraview if you only have vtus results.

then "awk" to produce "x,y,value"

fab
by fgillet
04 Nov 2016, 15:24
Forum: Elmer/Ice
Topic: Transfer basal drag from one mesh to another mesh
Replies: 14
Views: 17944

Re: Transfer basal drag from one mesh to another mesh

Hey, to do this I often save my results in an ascii file (x,y,value) and use the Scattered2DDataInterpolator ( http://elmerice.elmerfem.org/wiki/doku.php?id=solvers:scattered) to interpolate in the new mesh. This solver use external librairies that implement Natural neighbours and cubic spline inter...
by fgillet
04 Nov 2016, 14:02
Forum: Elmer/Ice
Topic: unique global node identifiers?
Replies: 2
Views: 5717

Re: unique global node identifiers?

Hi, look at the ParallelInfo_t type; the array GlobalDOFs (e.g. Mesh % ParallelInfo % GlobalDOFs) will give you the global node numbering. The other structures will give you informations about neighbours etc... Not sure if it is documented somewhere; better to directly look in source codes to see ho...
by fgillet
26 Apr 2016, 18:40
Forum: Elmer/Ice
Topic: AIFlow fabric directionality
Replies: 2
Views: 5843

Re: AIFlow fabric directionality

Hello, I have commited a fix for the computation of the euler angles for the fabric orientation in 3D. But I think, as I understand you cylinder set-up, that your results are correct, i.e. vertical one maximum fabrics are more favorable than horizontal 1 max fabrics. For perfectly aligned fabrics (i...
by fgillet
27 Aug 2015, 11:17
Forum: Elmer/Ice
Topic: Bedrock Sliding Boundary Condition
Replies: 4
Views: 9198

Re: Bedrock Sliding Boundary Condition

Hi, It's the iterative method for the linear system of the free surface (Solver 6) which is not converging. Increasing the verbosity level (Max Output Level = I in the Simulation section; max I=42) can sometimes help to see where the problem comes from. Have you checked your velocity field (you can ...
by fgillet
18 Aug 2015, 16:16
Forum: Elmer/Ice
Topic: Bedrock Sliding Boundary Condition
Replies: 4
Views: 9198

Re: Bedrock Sliding Boundary Condition

Hi,

Have you tried to use only Picard iterations (I see in your Navier-Stokes solver that you switch to Newton after 5 iterations)?
This is the first thing to try when the non linear iterations diverge.

cheers

Fabien
by fgillet
17 Jul 2015, 19:28
Forum: Elmer/Ice
Topic: Problems with large meshes
Replies: 3
Views: 7524

Re: Problems with large meshes

Hello; which version of MUMPS are you using? looking from the MUMPS versions changelog.: from 4.8.4 to 4.9: 18. Message "problem with NIV2_FLOPS message" should no more occur from 4.10.0 to 5.0.0: 24. Message "problem with NIV2_FLOPS message" should no more occur (there was still...
by fgillet
16 Jul 2015, 18:46
Forum: Elmer/Ice
Topic: Inverse Methods with SSA
Replies: 1
Views: 5541

Re: Inverse Methods with SSA

Hi David, you're right the Adjoint and Robin methods in the Elmer/Ice distribution are for Stokes only. I have developped the adjoint tools to invert viscosity or/and basal slipperiness with the SSA solver. You can find an application in this paper in discussion for The Cryosphere: http://www.the-cr...
by fgillet
29 Aug 2014, 14:11
Forum: Bug reports
Topic: SuperLU
Replies: 2
Views: 4300

Re: SuperLU

Hi,

tested on my Mac and a linux cluster.
It's looking good.

Thanks

Fabien
by fgillet
29 Aug 2014, 10:28
Forum: Bug reports
Topic: SuperLU
Replies: 2
Views: 4300

SuperLU

Hello, found the memory leak with SuperLU; for proper deallocation of all the objects in SolveSuperLU.c: REPLACE: Destroy_CompCol_Permuted(&AC); WITH: pxgstrf_finalize(&options, &AC); and REPLACE: Destroy_CompCol_Matrix(LUfactors->U); Destroy_SuperNode_Matrix(LUfactors->L); WITH: Destroy...