[Libmesh-devel] problem with introducing new element type to libmesh

2011-07-01 Thread Jan Biermann
Hi, I have a problem introducing a new element type (GEGENBAUER_2_5 -> in enums defined as 7) in libmesh. actually I had it done already in an older version, but getting it into a newer version creates problems. I put the compile output below and as you can see, everything compiles fine but the

Re: [Libmesh-devel] Possible race condition in threaded periodic constraints

2011-07-01 Thread Derek Gaston
On Jul 1, 2011, at 1:58 PM, Roy Stogner wrote: > > On Fri, 1 Jul 2011, Derek Gaston wrote: > >> Yes periodic BCs have never worked with threading... they deadlock. > > The disgusting thing about threading bugs is that there's no such > thing as as "never worked". We do test one periodic BC

Re: [Libmesh-devel] Possible race condition in threaded periodic constraints

2011-07-01 Thread Roy Stogner
On Fri, 1 Jul 2011, Derek Gaston wrote: > Yes periodic BCs have never worked with threading... they deadlock. The disgusting thing about threading bugs is that there's no such thing as as "never worked". We do test one periodic BC example, and it gets run with 2 threads at least 10 times a

Re: [Libmesh-devel] Possible race condition in threaded periodic constraints

2011-07-01 Thread Derek Gaston
Yes periodic BCs have never worked with threading... they deadlock. We've just never got around to fixing it. Glad you tracked it down. If you're not going to have time to get to it just let us know. Derek On Jul 1, 2011, at 12:40 PM, Roy Stogner wrote: > > I just noticed this while fu

[Libmesh-devel] Possible race condition in threaded periodic constraints

2011-07-01 Thread Roy Stogner
I just noticed this while futzing with some related code: In FEBase::compute_periodic_constraints, we lock DofConstraints writes with a spin_mutex, but those writes may be happening concurrently with the unlocked reads elsewhere in the same function; and in that case on our underlying STL map the

[Libmesh-devel] problem with introducing new element type to libmesh

2011-07-01 Thread Jan Biermann
Hi, I have a problem introducing a new element type (GEGENBAUER_2_5 -> in enums defined as 7) in libmesh. actually I had it done already in an older version, but getting it into a newer version creates problems. I put the compile output below and as you can see, everything compiles fine but the