[Libmesh-users] SNES variational inequalities are broken

Hi, I just pulled and built the HEAD (0cceba267df9e8f10a9c739a928c0b109039c888) this morning. The SNES variational inequality solvers are break on the headnode when assembling the bounds. For example misc/ex7 now breaks with: (base) node004:/shared/libmesh/Linux-gnu-O-new/examples/miscellaneous

Re: [Libmesh-users] new libmeshdm fails with newtonrsls

Hi Roy, You can simply regenerate the error if you pass: ./example-opt --verbose dim=1 N=1024 initial_state=strip initial_center=0.5 initial_width=0.1 dt=1e-10 max_time=.5e-8 -snes_type vinewtonssls --use-petsc-dm This generates a seg fault on my mac with the PETSc and libMesh revisions I prov

Re: [Libmesh-users] pre-partitioning a big mesh

Thank you Derek, This is really helpful. Ata > On Jan 11, 2017, at 6:15 PM, Derek Gaston wrote: > > splitter also has the ability to do M->N parittioning. What this means is > that splitter itself can be run using MPI... and can generate partitions for > any number of processors. In fact, y

[Libmesh-users] Vector containing a quantity over all quad points

Dear users/developers, I’m wondering if there is a way to create a numeric vector (attach it to a system or not) that has the right shape (length, distribution, ghost points) so I can store a quantity calculated from solution of a system for each quadrature point for each element in the mesh.