[Libmesh-users] SparseMatrix Determinant

2018-08-08 Thread Barna Becsek
Hello, is there a way to compute a SparseMatrix’s determinant? I would like to solve a linear system but my final residual always equals exactly 0 after just one iteration (where the solve terminates). Hence, I would like to check the matrix for singularity. Thank you. Sincerely, Barna --

[Libmesh-users] Checking for NaNs in LinearImplicitSystems Solution

2018-08-07 Thread Barna Becsek
Hi, I was wondering what the easiest way were to check if any of the values in the solution of a linear system are NaNs. best, Barna -- Check out the vibrant tech community on one of the world's most engaging tech site

[Libmesh-users] LinearImplicitSystem and LinearSolver Methods

2018-08-06 Thread Barna Becsek
Hello, I have a problem, where I need to solve a linear system in a transient simulation. However, the system I solve has no transient character, i.e. in each time step I just need to solve a new system. Anyway, at some point of my simulation, it seems like the result of my linear system become

[Libmesh-users] Reference Elements with Boundary Nodes

2017-06-08 Thread Barna Becsek
Hello, What would be the easiest way to reference elements that contain a node on a boundary. Ideally, I would like to have a mesh object with these elements but unchanged node ids. Does the BoundaryMesh class handle this sort of thing or is that class for actual boundary element (e.g. 2D els i

Re: [Libmesh-users] Lean setup of FE

2017-06-08 Thread Barna Becsek
Thank you, that was very helpful. > On 6 Jun 2017, at 14:57, Paul T. Bauman wrote: > > > > On Tue, Jun 6, 2017 at 8:52 AM, Barna Becsek <mailto:barnabec...@gmail.com>> wrote: > Hello, > > I am in a situation where I need to setup test functions for an L2

[Libmesh-users] Lean setup of FE

2017-06-06 Thread Barna Becsek
Hello, I am in a situation where I need to setup test functions for an L2projection. Essentially, all I have to access in the end is get_phi() and get_JxW(). Now, I need to call fe->reinit() for every element that I add. Since I do not need anything else than get_phi() and get_JxW(), I would li

Re: [Libmesh-users] Parallel Mesh Generation in Libmesh

2017-03-13 Thread Barna Becsek
ote: > > > On Fri, 3 Mar 2017, Barna Becsek wrote: > >>> On 3 Mar 2017, at 17:47, Roy Stogner wrote: >>> >>> On Fri, 3 Mar 2017, Barna Becsek wrote: >>> >>>>> On 3 Mar 2017, at 17:28, Roy Stogner wrote: >>>>> >

Re: [Libmesh-users] Parallel Mesh Generation in Libmesh

2017-03-09 Thread Barna Becsek
Dear Roy, I was finally able to install in a location where I can both read from compute nodes and edit the files. In the attachment is what I found. Basically, it seems like get_weights() returns an incorrect address. –Barna > On 7 Mar 2017, at 11:19, Barna Becsek wrote: > >&

Re: [Libmesh-users] Parallel Mesh Generation in Libmesh

2017-03-07 Thread Barna Becsek
> On 3 Mar 2017, at 22:30, Roy Stogner wrote: > > > On Fri, 3 Mar 2017, Barna Becsek wrote: > >>> On 3 Mar 2017, at 17:47, Roy Stogner wrote: >>> >>> On Fri, 3 Mar 2017, Barna Becsek wrote: >>> >>>>> On 3 Mar 2017, at 1

Re: [Libmesh-users] Parallel Mesh Generation in Libmesh

2017-03-01 Thread Barna Becsek
Dear Roy, I found the place the assertion is tripped. In the attachment you can see the stack and location of the casting. Haven’t dug into why this happens yet… –Barna > On 1 Mar 2017, at 14:50, Barna Becsek wrote: > > Alright, thank you for giving it a look. As far as I can

Re: [Libmesh-users] Parallel Mesh Generation in Libmesh

2017-03-01 Thread Barna Becsek
Roy Stogner wrote: > > > On Tue, 28 Feb 2017, Barna Becsek wrote: > >> I actually attached the stack to the previous email but I am not >> sure whether I can actually send attachments to this mailing list. > > Depends on attachment size, but I simply overlooked it. I

Re: [Libmesh-users] Parallel Mesh Generation in Libmesh

2017-02-28 Thread Barna Becsek
ay to include pictures in here? Thank you for your patience. –Barna > On 28 Feb 2017, at 17:57, Roy Stogner wrote: > > > On Tue, 28 Feb 2017, Barna Becsek wrote: > >> To remedy this I called prepare_for_use once additionally before >> gather_neigboring_elements (to m

Re: [Libmesh-users] Parallel Mesh Generation in Libmesh

2017-02-28 Thread Barna Becsek
gt; > On Mon, 27 Feb 2017, Barna Becsek wrote: > >> I hope this does not come too late. > > That depends entirely on your deadlines, I fear. ;-) > >> But I finally got the debugger to work with some help of support on >> the Supercomputer that we are using. This

Re: [Libmesh-users] Parallel Mesh Generation in Libmesh

2017-02-27 Thread Barna Becsek
017, at 17:04, Roy Stogner wrote: > > > On Fri, 13 Jan 2017, Barna Becsek wrote: > >> Hmm, you are right but I cannot find a process that would correspond to the >> children of any of these. > > Is this a cluster, by any chance? Where slurm is running o

Re: [Libmesh-users] Parallel Mesh Generation in Libmesh

2017-01-13 Thread Barna Becsek
pipe_w 15:14 pts/107 00:00:00 /opt/slurm/16.05.8/bin/srun -n 8 whale-dbg -i IMP/RunImpact2D.i 0 S becsekba 77353 61771 0 80 0 - 2637 pipe_w 15:19 pts/117 00:00:00 grep --color=auto whale-dbg > On 13 Jan 2017, at 15:12, Roy Stogner wrote: > > > On Fri, 13 Jan 2017,

Re: [Libmesh-users] Parallel Mesh Generation in Libmesh

2017-01-13 Thread Barna Becsek
0x0041b7f1 in reader_loop () #6 0x0041b4db in main () > On 13 Jan 2017, at 12:44, Barna Becsek wrote: > > Ok, this is the backtrace of the running processes. There are two processes > running: > > 0 S becsekba 54451 54421 0 80 0 - 76108 futex_ 12:39 pts/92 00:0

Re: [Libmesh-users] Parallel Mesh Generation in Libmesh

2017-01-13 Thread Barna Becsek
–Barna > On 12 Jan 2017, at 17:51, Roy Stogner wrote: > > > On Thu, 12 Jan 2017, Barna Becsek wrote: > >> What I meant was the program will not exit gather_neighboring_elements. I >> think the processes are still running. > > Right. But you can e.g. attach

Re: [Libmesh-users] Parallel Mesh Generation in Libmesh

2017-01-12 Thread Barna Becsek
What I meant was the program will not exit gather_neighboring_elements. I think the processes are still running. > On Jan 12, 2017, at 4:26 PM, Roy Stogner wrote: > > >> On Thu, 12 Jan 2017, Barna Becsek wrote: >> >> thank you very much for these thorough elaborat

Re: [Libmesh-users] Parallel Mesh Generation in Libmesh

2017-01-12 Thread Barna Becsek
coarse meshes (9x9). Maybe I can try this manually. I added some more comments/questions below. > On 12 Jan 2017, at 00:06, Roy Stogner wrote: > > > On Tue, 10 Jan 2017, Barna Becsek wrote: > >> What we DON’T HAVE: >> - A clear understanding of how ghost elements in Lib

[Libmesh-users] Parallel Mesh Generation in Libmesh

2017-01-10 Thread Barna Becsek
Dear Users, We are writing a multiphysics tool where we need to replicate a FD mesh in a FE way. We have been almost able to translate one into the other but are facing some minor difficulties (hopefully). I would greatly appreciate if someone could help us find the missing steps for having a c