[Libmesh-devel] VTK output fails in parallel

2011-02-09 Thread Vetter Roman
Dear libmesh developers VTK output crashes when run in parallel (mpirun -n 2 ...). The following assertion fails: Assertion `it!=_global_to_local_map.end()' failed. [...]/libmesh/include/numerics/petsc_vector.h, line 1073, [...] The problem is reproducible with the following minimal working exa

[Libmesh-devel] convergence reason in PETSc nonlinear solver

2011-02-22 Thread Vetter Roman
Hi guys I propose to add the attached patch to trunk. It adds an integer storing the convergence reason returned by the PETSc nonlinear solver, closely analogous to the existing convergence boolean flag, which is somewhat meager. Without this, the convergence reason is lost right after solving, wi

Re: [Libmesh-devel] Recent library changes causing issues

2011-03-17 Thread Vetter Roman
>> I've just checked in a change which will hopefully please both the >> people who need to call this->clear() at the end of >> PetscNonlinearSolver::solve and the people who want >> PetscNonlinearSolver::print_converged_reason() to work. The basic >> idea is just to store the SNESConvergedReason

Re: [Libmesh-devel] Recent library changes causing issues

2011-03-19 Thread Vetter Roman
On Tue, Mar 8, 2011 at 5:01 PM, Cody Permann wrote: > I'm pretty sure "if (_snes)" is guaranteed to be false after a call to > SNESDestroy, but I could be wrong about that... Vetter since you are > the one actually using this functionality could you please let me know > if it's still working OK fo

[Libmesh-devel] Random seed for MeshTools:Modification::distort

2011-12-22 Thread Vetter Roman
Dear libMesh devs please consider adding the attached patch (for r5017) to the trunk. It adds the possibility to specify a custom random seed for the mesh distortion in MeshTools:Modification::distort in a backwards-compatible way. This is useful for solving the same problem with differently di

Re: [Libmesh-devel] Random seed for MeshTools:Modification::distort

2011-12-25 Thread Vetter Roman
> Would it make sense to parallelize it now > (breaking the consistency with previous distort() > results) *before* anyone has a bunch of > differently-seeded distorted meshes depending > on that consistency? In case you're asking me: I don't feel like I can answer that question. I guess I won't

Re: [Libmesh-devel] Random seed for MeshTools:Modification::distort

2012-01-09 Thread Vetter Roman
*bump* May I ask if something of the proposed kind will be implemented anytime soon? If not, I will have to continue using my own hacked "branch" of libmesh instead of the trunk... ;-) Roman Von: Vetter Roman [vette...@ethz.ch] Gesendet: S

[Libmesh-devel] OFFIO for 1D mesh

2012-09-14 Thread Vetter Roman
Hello everybody, I would like to propose the attached small patch for the trunk. It adds support for 1D meshes to the simple OFF mesh reader class, which is currently limited to 2D meshes. Best regards, Roman Index: src/mesh/off_io.C =

[Libmesh-devel] Rotating a 1D mesh

2012-09-14 Thread Vetter Roman
Hi! What's the reason for asserting mesh.mesh_dimension() != 1 in MeshTools::Modification::rotate()? Rotating a 1D mesh seems to work for me in opt mode, but the assertion obviously fails in dbg mode. Best regards, Roman

Re: [Libmesh-devel] OFFIO for 1D mesh

2012-09-20 Thread Vetter Roman
*bump* Von: Vetter Roman [vette...@ethz.ch] Gesendet: Freitag, 14. September 2012 18:28 An: libmesh-devel@lists.sourceforge.net Betreff: [Libmesh-devel] OFFIO for 1D mesh Hello everybody, I would like to propose the attached small patch for the trunk