[Libmesh-devel] Mesh and DOF numbering

2014-05-21 Thread Dana Christen
Hi all, I would like the DOF numbering of a system to match that of the original mesh. Does setting allow_renumbering(false) in the Mesh instance guarantee that the global DOF returned by DofMap.dof_indices will match the original numbering from the mesh file? My code runs on one process. The reas

Re: [Libmesh-devel] Convergence reason flag in LinearSolver

2014-02-13 Thread Dana Christen
6, > // MOOSE_DIVERGED_NONSYMMETRIC= -7, > // MOOSE_DIVERGED_INDEFINITE_PC = -8, > // MOOSE_DIVERGED_NANORINF= -9, > // MOOSE_DIVERGED_INDEFINITE_MAT = -10 > }; > > > You can see that we're not really using all of those linear ones yet.

Re: [Libmesh-devel] Convergence reason flag in LinearSolver

2014-02-12 Thread Dana Christen
Hehe ok sounds good! I'll probably start with some PETSc-specific fix for the time being. I'll get back to you when I have a generic patch. Cheers, Dana On Wed, Feb 12, 2014 at 5:05 PM, Roy Stogner wrote: > > On Wed, 12 Feb 2014, Dana Christen wrote: > > I need

Re: [Libmesh-devel] Convergence reason flag in LinearSolver

2014-02-12 Thread Dana Christen
class... Dana On Wed, Feb 12, 2014 at 5:10 PM, John Peterson wrote: > > > > On Wed, Feb 12, 2014 at 9:05 AM, Roy Stogner wrote: > >> >> On Wed, 12 Feb 2014, Dana Christen wrote: >> >> > I need a way to retrieve the "convergence reason" flag

[Libmesh-devel] Convergence reason flag in LinearSolver

2014-02-12 Thread Dana Christen
Hello, I need a way to retrieve the "convergence reason" flag provided by the PETSc iterative solvers from libMesh. Currently, the only way to access this information is to print it with PetscLinearSolver::print_converged_reason(), which is not very convenient. I've been thinking about extending t

[Libmesh-devel] Mesh stitching and find_neighbors()

2014-01-31 Thread Dana Christen
Hello, I'm trying to optimize the stitching_helper() function in SerialMesh, which is currently a bottleneck in the work I'm doing (stitching together a large number of small components to assemble a full mesh). I've found out that over 85% of the running time of this function is actually spent i