Search found 12 matches

by Evan
08 Mar 2016, 07:54
Forum: Elmer/Ice
Topic: Pressure at the base in Full Stokes
Replies: 2
Views: 5671

Re: Pressure at the base in Full Stokes

Hi Peter, thanks for the answer. I looked into it some more, and it seems that the output I am getting has messed up pressure values in even numbered rows of nodes, while the odd numbered rows have correct values that are close to hydrostatic. This likely indicates there is an integer division error...
by Evan
07 Mar 2016, 06:26
Forum: Elmer/Ice
Topic: Pressure at the base in Full Stokes
Replies: 2
Views: 5671

Pressure at the base in Full Stokes

Hi, I've been looking at the pressure values at the bottom row of nodes when solving full Stokes, and noticed that the values are different (and usually lower) than purely hydrostatic pressure. I figured out that it is because the pressure being output is for the normal-tangential coordinate system,...
by Evan
20 Jul 2015, 07:26
Forum: Elmer/Ice
Topic: Problems with large meshes
Replies: 3
Views: 7551

Re: Problems with large meshes

BTW, if anyone else runs into this problem, upgrading to Mumps 5.0 fixed it.
by Evan
18 Jul 2015, 03:44
Forum: Elmer/Ice
Topic: Problems with large meshes
Replies: 3
Views: 7551

Re: Problems with large meshes

Hi Fabien, thanks for the response. I believe the supercomputer I am using only has Mumps 4.9. My laptop has 4.10 and has never run into this problem. I'll insist more that they upgrade Mumps, and hopefully that will solve the problem. Unfortunately using Parstokes is not an option, since we are usi...
by Evan
17 Jul 2015, 11:14
Forum: Elmer/Ice
Topic: Problems with large meshes
Replies: 3
Views: 7551

Problems with large meshes

Hello, I am currently trying to run a large (read: ice sheet sized) simulation, which demands a fairly fine resolution mesh. I am having problems with Mumps when using the Navier Stokes solver: ListToCRSMatrix: Matrix format changed from CRS to List problem with NIV2_FLOPS message -1.192092895507812...
by Evan
29 May 2015, 16:08
Forum: Elmer/Ice
Topic: Rotated coordinate system for basal boundary condition
Replies: 3
Views: 6927

Re: Rotated coordinate system for basal boundary condition

Hi Peter, Thank you for the explanation. It turns out the problem I am having is not a problem with the rotation itself, but that we are using a reduced A matrix for DefaultSolve after the first time step. Somewhere in our code we have missed changing a pointer, so it is rotating too many nodes. Bec...
by Evan
28 May 2015, 18:12
Forum: Elmer/Ice
Topic: Rotated coordinate system for basal boundary condition
Replies: 3
Views: 6927

Re: Rotated coordinate system for basal boundary condition

Just to follow up, after investigating the problem I am having with our code, it seems that the Back-Rotate routine is rotating not only the bottom layer of nodes, but also the layer directly above it. However, after exit from DefaultSolve, the second layer is not rotated back. Below is a stack of n...
by Evan
26 May 2015, 14:28
Forum: Elmer/Ice
Topic: Rotated coordinate system for basal boundary condition
Replies: 3
Views: 6927

Rotated coordinate system for basal boundary condition

Hi, I am tracking down some issues in our code, and I am seeking some clarification on the way that Elmer handles the non-penetrating condition at the basal boundary condition (i.e. the ice cannot penetrate through the bedrock) in the Navier-Stokes solver. As I understand it, Elmer rotates the coord...
by Evan
02 Apr 2015, 15:44
Forum: Software development
Topic: Memory leak problem with ParallelInitMatrix call
Replies: 0
Views: 3374

Memory leak problem with ParallelInitMatrix call

Hi there, First off, I am using Elmer revision 7111 on Linux Mint 17. I'm working on Josefin's modification to the Stokes solver in Elmer/Ice. I have been running into problems with memory leaks when running in parallel, which are significant (my computer runs out of memory on a 11,000 element simul...
by Evan
13 Jan 2015, 15:48
Forum: Elmer/Ice
Topic: Elmer/Ice test cases
Replies: 1
Views: 6813

Re: Elmer/Ice test cases

Hi, Cheng Gong and I have been trying to run the Elmer/Ice test cases as well. Pretty much failing in the same spots, with the same errors. I have version 6874 on my computer, while Cheng Gong has something closer to the latest version (7098). We both had failures on the same tests (I assume the las...