Re: [Libmesh-devel] constrain_element_matrix bug?

2016-02-08 Thread Robert Blake
Following up on this email. John Peterson, Roy Stogner, any ideas? On 02/04/2016 10:35 AM, Rob Blake wrote: tl;dr: I think the 4 argument version of DofMap::constrain_element_matrix is bugged, but I can't fix it without breaking the moose tests. Help? Is the 4 argument version of con

Re: [Libmesh-devel] Deadlock when using assign_global_indices, problems with hilbert indices

2014-10-30 Thread Robert Blake
Sorry, ignore the 5000^3 operation for a moment. I wrote up sample code without testing it. In my email to Roy, I included working code that manifests the bug in a < 1M node mesh. On 10/30/2014 05:06 PM, Kirk, Benjamin (JSC-EG311) wrote: > On Oct 29, 2014, at 10:10 PM, Roy Stogner wrote: > >

Re: [Libmesh-devel] Deadlock when using assign_global_indices, problems with hilbert indices

2014-10-30 Thread Robert Blake
seem to replicate the problem with 14.04. On 10/29/2014 11:10 PM, Roy Stogner wrote: On Wed, 29 Oct 2014, Robert Blake wrote: We're running into a intermittent deadlock when calling assign_global_indices(). The problem only appears on >= 5 processors, and even then it only appears intermi

[Libmesh-devel] Deadlock when using assign_global_indices, problems with hilbert indices

2014-10-29 Thread Robert Blake
Hi, We're running into a intermittent deadlock when calling assign_global_indices(). The problem only appears on >= 5 processors, and even then it only appears intermittently. Here's the error we see: Assertion `hi <= node_upper_bounds[communicator.rank()]' failed. hi = 887794606_2212644264

[Libmesh-devel] Custom Convergence Criterion for NewtonSolver

2014-07-14 Thread Robert
the functionality up to DiffSolver and implement it for PetscDiffSolver as well) Would this feature be useful for anyone else? Any further comments? Robert -- Want fast and easy access to all the code in your enterprise

Re: [Libmesh-devel] LibMesh FEMSystem Questions

2014-05-28 Thread Robert
Hello, On Wed, May 28, 2014 at 08:02:33AM -0500, Roy Stogner wrote: > > Copying to libmesh-devel in case anyone else wants to interject. > > On Wed, 28 May 2014, Robert Weidlich wrote: > > > My aim is to solve standard solid mechanic tasks with contact, so I > >

Re: [Libmesh-devel] Compilation Error with Intel Compiler

2013-03-22 Thread Robert
Hello, I did some additional builds and everything before f1718977c804424f0dd5befd9435f60b4e39b28d builds fine and with this revision (and probably with everything after that too) building fails. I will just use another compiler as advised by Paul. Robert On Fri, Mar 22, 2013 at 09:22:23AM

[Libmesh-devel] Compilation Error with Intel Compiler

2013-03-22 Thread Robert
lem->processor_id() = +std::min (libmesh_cast_int(i / target_size), + libmesh_cast_int(n-1)); But I cannot really see any error in there. Has anyone else an idea? Robert -- Everyone hates slow websites.

Re: [Libmesh-devel] [patch] VTK_IO: memleak, cleanup, parallel output

2013-01-24 Thread Robert
On Wed, Jan 23, 2013 at 12:43:17PM -0600, Roy Stogner wrote: > > On Thu, 10 Jan 2013, Robert wrote: > > > On Thu, Jan 10, 2013 at 07:27:43AM -0600, Roy Stogner wrote: > >> On Thu, 10 Jan 2013, Robert wrote: > >>> The compression feature depends on a somewhat r

Re: [Libmesh-devel] Monitor solution norm in NewtonSolvery

2013-01-11 Thread Robert
See > https://github.com/libMesh/libmesh/pull/29 Let me know, if there is anything to change. Robert On Thu, Jan 10, 2013 at 08:24:17AM -0600, Paul T. Bauman wrote: > On Thu, Jan 10, 2013 at 7:52 AM, Roy Stogner wrote: > > Would you also pass norm_total and any other

[Libmesh-devel] Monitor solution norm in NewtonSolver

2013-01-10 Thread Robert
function-pointer to NewtonSolver and the solve routine calls this function with the wanted value. Would this be a solution which could make it into the official libmesh? Does anyone has a better solution for this (beside parsing the output). Robert

Re: [Libmesh-devel] Automake merge complete

2012-11-20 Thread Robert
ile And one last thing: Does any know a possiblity to use the pkg-config version selector to switch between dbg and opt builds? Robert On Mon, Nov 19, 2012 at 10:43:20AM -0600, Kirk, Benjamin (JSC-EG311) wrote: > On Nov 19, 2012, at 10:07 AM, Robert wrote: > > > Hi Ben, > >

Re: [Libmesh-devel] Automake merge complete

2012-11-19 Thread Robert
Hi Ben, I am using the mpi installation provided by PetSc with --download-mpich=1. So the mpi used should be the one in ~/opt/build/petsc-3.2-p6/linux-gnu-opt/include Robert On Mon, Nov 19, 2012 at 09:28:22AM -0600, Kirk, Benjamin (JSC-EG311) wrote: > Very good, thanks. Where on your system

Re: [Libmesh-devel] Automake merge complete

2012-11-19 Thread Robert
llation of libmesh (r5498) works fine on the same computer with MPI. Even with "make V=0" I could not get more information about the compilerflags used, so I cannot really figure out, if it is a problem on my side or the build-system. Robert On Mon, Nov 19, 2012 at 08:19:15AM -0600, Kirk,

Re: [Libmesh-devel] Automake merge complete

2012-11-19 Thread Robert
match error. This is libtool 2.4.2, but the libtool: definition of this LT_INIT comes from libtool 2.2.6b. libtool: You should recreate aclocal.m4 with macros from libtool 2.4.2 libtool: and run autoconf again. make[1]: *** [gmvread.lo] Error 63 make[1]: Leaving directory `/home//opt/build/l

Re: [Libmesh-devel] [Libmesh-users] Averaging element fields to the nodes.

2012-10-29 Thread Robert
that Robert On Mon, Oct 29, 2012 at 06:06:05AM -0500, Kirk, Benjamin (JSC-EG311) wrote: > That method lets you conveniently access all the elements sharing a given > node, which is not what you want in this case. > > The usual way to do this is by solving an L2 projection onto

Re: [Libmesh-devel] Periodic BCs with Offset

2012-09-18 Thread Robert
of resolved according to the mailinglist, has anyone already coded this "easy" addition? Will this also work for the new vector_fe stuff? I am currently using the displacements (u,v,w) as primary variables and want to switch to using the nodes position (x,y,z) as primary variables

Re: [Libmesh-devel] FELagrangeVec

2012-06-21 Thread Robert
size. I tried 3-dimensional vectors for displacements and 6-dimensional vectors for postprocessed stresses. I currently stores these quantities in variables with suffixes "_x", "_y", "_z" and "_xx", "_xy", "_xz", "_yy", "_yx&qu

Re: [Libmesh-devel] [Libmesh-users] Elastic Solid Examples?

2012-01-12 Thread Robert
On Thu, Jan 12, 2012 at 10:13:03AM +0100, Robert wrote: > On Wed, Jan 11, 2012 at 02:10:30PM -0600, Roy Stogner wrote: > > On Tue, 10 Jan 2012, Robert wrote: > > > I have just tested this with the default configure options. Any > > > suggestions for improve

Re: [Libmesh-devel] [Libmesh-users] Elastic Solid Examples?

2012-01-12 Thread Robert
Hi, On Wed, Jan 11, 2012 at 02:10:30PM -0600, Roy Stogner wrote: > On Tue, 10 Jan 2012, Robert wrote: > > I have just tested this with the default configure options. Any > > suggestions for improvements are of course welcome. > > Would you try running it in parallel? Th

Re: [Libmesh-devel] [Libmesh-users] Elastic Solid Examples?

2012-01-10 Thread Robert
tensor classes, I discovered a small bug in TypeTensor::assign. See tensor.patch. Regards Robert On Mon, Jan 09, 2012 at 09:02:29AM +1300, Mark Davies wrote: > Hi Guys, > > I'm just getting started with the library and was wondering if there are > any examples dealing wi

Re: [Libmesh-devel] [patch] VTK_IO: memleak, cleanup, parallel output

2011-11-15 Thread Robert
esh and for the displacements and stresses at the nodes. I hope this are the informations you were looking for. Robert > > -- > John -- RSA(R) Conference 2012 Save $700 by Nov 18 Register now http://p.s

[Libmesh-devel] [patch] VTK_IO: memleak, cleanup, parallel output

2011-11-14 Thread Robert
ons, this is not tested very well, but hopefully useful anyway. Additionally I have added a flag and the correspoding api call for using the compressed VTK output. I attached the patch and uploaded it to this no-paste: > http://dpaste.de/4uEEA/ because at least the libmesh-users-list is st