Re: [Libmesh-users] Assertion `mesh.comm().verify(mesh.n_elem())' failed.

2017-05-10 Thread Roy Stogner
On Wed, 10 May 2017, Harshad Sahasrabudhe wrote: > What libMesh release or git hash is this? > > This is LibMesh 0.9.5. More than a year old? Can you replicate the problem either with libMesh 1.2 or with the git master/HEAD? --- Roy --

Re: [Libmesh-users] Reusing the preconditioner or matrices onto a new refined grid

2017-05-10 Thread Paul T. Bauman
Boris, please point us to the branch, thanks. On Wed, May 10, 2017 at 2:10 PM, Roy Stogner wrote: > > On Tue, 9 May 2017, Paul T. Bauman wrote: > > On Tue, May 9, 2017 at 9:36 PM, Roy Stogner >> wrote: >> >> If so, then you might want to help us get the projection matrix code >> wor

Re: [Libmesh-users] Assertion `mesh.comm().verify(mesh.n_elem())' failed.

2017-05-10 Thread Harshad Sahasrabudhe
> > What libMesh release or git hash is this? This is LibMesh 0.9.5. Harshad -- Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot

Re: [Libmesh-users] Assertion `mesh.comm().verify(mesh.n_elem())' failed.

2017-05-10 Thread Roy Stogner
On Wed, 10 May 2017, Harshad Sahasrabudhe wrote: > I'm trying to import a mesh using the GmshIO class. However, I get the > following error while broadcasting the mesh when I run with 2 processes: > > Assertion `mesh.comm().verify(mesh.n_elem())' failed. > > [0] src/mesh/mesh_communication.C, lin

[Libmesh-users] Assertion `mesh.comm().verify(mesh.n_elem())' failed.

2017-05-10 Thread Harshad Sahasrabudhe
Hi, I'm trying to import a mesh using the GmshIO class. However, I get the following error while broadcasting the mesh when I run with 2 processes: Assertion `mesh.comm().verify(mesh.n_elem())' failed. [0] src/mesh/mesh_communication.C, line 814, compiled [1] src/mesh/mesh_communication.C, line

Re: [Libmesh-users] Reusing the preconditioner or matrices onto a new refined grid

2017-05-10 Thread Roy Stogner
On Tue, 9 May 2017, Paul T. Bauman wrote: On Tue, May 9, 2017 at 9:36 PM, Roy Stogner wrote: If so, then you might want to help us get the projection matrix code working.  I've got incomplete code in the roystgnr/projection_matrix branch; the matrix assignment code is writte

[Libmesh-users] assemble_before_solve not useful in FEMSystem ?

2017-05-10 Thread Salazar De Troya, Miguel
Hello, It seems that by default FEMSystem always assembles the matrix and ignores the flag System::assemble_before_solve. The DiffSolvers, NewtonSolver and PetscDiffSolver always call System::assembly(false,true). The workarounds I see would be to create my own TimeSolver to override TimeSolver

Re: [Libmesh-users] Reusing the preconditioner or matrices onto a new refined grid

2017-05-10 Thread Salazar De Troya, Miguel
Yes, it would be like that. I will look at that branch. Thanks. Miguel -- On 5/9/17, 6:36 PM, "Roy Stogner" wrote: On Tue, 9 May 2017, Salazar De Troya, Miguel wrote: > There is a function System::project_vector() to project the vectors > defined in one mesh onto a new