ParStokes example .sif?

Extension of Elmer in computational glaciology

ParStokes example .sif?

Postby mewo2 » 05 Mar 2015, 17:09

Hi all,

I'm trying to get up to speed with the ParStokes solver using the information on the wiki. However, I'm having trouble getting a simple inclined slab test-case to run. When I try, I get the error "ERROR:: LUDecomp: Matrix is singluar.", followed by a segfault from within ParStokes. I suspect that I have something misconfigured (initial conditions?) but I'm not sure exactly what.

Does anybody have a working example of a simple test-case for ParStokes?

Cheers,
Martin
mewo2
 
Posts: 6
Joined: 26 Sep 2014, 15:30

Re: ParStokes example .sif?

Postby mika » 20 Mar 2015, 10:58

Hi,

The following sif-file should give a starting point for writing a working case definition.

Best regards,
Mika
Attachments
simplecase.sif
(2.9 KiB) Downloaded 1316 times
mika
 
Posts: 39
Joined: 15 Sep 2009, 07:44

Re: ParStokes example .sif?

Postby mewo2 » 20 Mar 2015, 17:33

Thanks, Mika - I managed to figure it out on my own, but hopefully this will help others in the future.

Also, am I right in thinking that the ParStokes solver is 3d-only? I haven't had any success with 2d cases.
mewo2
 
Posts: 6
Joined: 26 Sep 2014, 15:30

Re: ParStokes example .sif?

Postby tzwinger » 21 Mar 2015, 18:42

Checking the code, it seems that the dimension of the problem is being taken care of - i.e., I found no reason why it shouldn't work in 2D. Question is - do you really need it in a 2D problem or are you equally happy with the standard Navier-Stokes Solver and a direct solver?

Regards,

Thomas
tzwinger
Site Admin
 
Posts: 79
Joined: 24 Aug 2009, 12:20

Re: ParStokes example .sif?

Postby mewo2 » 23 Mar 2015, 13:38

No, I'm fine with just 3d. I was only trying 2d for a quick test comparison between Navier-Stokes and ParStokes.
mewo2
 
Posts: 6
Joined: 26 Sep 2014, 15:30


Return to Elmer/Ice

Who is online

Users browsing this forum: No registered users and 1 guest