Re: [Libmesh-users] AMR speed

2017-04-28 Thread Rossi, Simone
gner <royst...@ices.utexas.edu<mailto:royst...@ices.utexas.edu>> wrote: On Thu, 27 Apr 2017, Rossi, Simone wrote: The run times for 100 timesteps using AMR can be more than 10 times slower than when using a fine uniform grid. For example, with a 16 x 16 x 16 uniform grid, 100 iteration

Re: [Libmesh-users] AMR speed

2017-04-28 Thread Rossi, Simone
| - On Apr 27, 2017, at 14:29, Rossi, Simone <sro...@email.unc.edu<mailto:sro...@email.unc.edu>> wrote: Ok, I ran again the tests with different max_h_levels with the perflog enabled.

Re: [Libmesh-users] AMR speed

2017-04-27 Thread Rossi, Simone
give you all the details as in the example. On Thu, Apr 27, 2017 at 10:54 AM, Rossi, Simone <sro...@email.unc.edu<mailto:sro...@email.unc.edu>> wrote: Dear Vikram, as in the examples, I am using the libmesh::KellyErrorEstimator. I’m compiling libmesh with the --ena

Re: [Libmesh-users] AMR speed

2017-04-27 Thread Rossi, Simone
ationSystems | | build_solution_vector() 1 0.0001 0.560.0001 0.640.79 0.90 | Thanks. On Wed, Apr 26, 2017 at 10:09 PM, Rossi, Simone <sro...@email.unc.edu<mailto:sro...@email.unc.edu

[Libmesh-users] AMR speed

2017-04-26 Thread Rossi, Simone
Dear Roy, dear Paul, dear all, I am testing AMR in libmesh using simple linear elements. My test case is a propagating front described by a reaction-diffusion equation with a cubic bistable reaction term. I followed the adaptivity examples to create this test case. The run times for 100

[Libmesh-users] DG implementation

2017-02-09 Thread Rossi, Simone
Dear all, I’m trying to implement a DG formulation for solid mechanics using linear elements. I’m currently using the following strategy: I define my variables as MONOMIALS but in the assembly I pretend my FEType is LAGRANGE. This gives me the correct values of the shape functions on the edges.

Re: [Libmesh-users] compilation error

2016-11-07 Thread Rossi, Simone
I have the --disable-unique-ptr option in my configuration script since it is required by another library I’m using on top of libMesh. Thanks, Simone On Nov 7, 2016, at 11:49, John Peterson <jwpeter...@gmail.com<mailto:jwpeter...@gmail.com>> wrote: On Mon, Nov 7, 2016 at 9:

[Libmesh-users] compilation error

2016-11-07 Thread Rossi, Simone
Dear all, I’m upgrading to PETSc 3.7.4 and I am now compiling the git version of libMesh. On my workstation (using gcc 5.4 and openmpi 2.0.0) I get a compilation error /not_backed_up/srossi/TPL/libmesh/libmesh-git/src/reduced_basis/rb_eim_assembly.C: In member function ‘virtual void

Re: [Libmesh-users] ksp monitor

2016-11-03 Thread Rossi, Simone
Hey John, I’m using PETSc 3.7.2. I’ll install a another version to check if that’s the problem. Thanks, Simone On Nov 3, 2016, at 17:36, John Peterson <jwpeter...@gmail.com<mailto:jwpeter...@gmail.com>> wrote: On Thu, Nov 3, 2016 at 2:45 PM, Rossi, Simone <sro...@email.unc

Re: [Libmesh-users] Systems of Equations, Ex4 -- Ex5

2016-11-03 Thread Rossi, Simone
elasticity problem, getting finally some instabilities in the pressure. Thanks for the help, Simone On Oct 13, 2016, at 13:12, David Knezevic <david.kneze...@akselos.com<mailto:david.kneze...@akselos.com>> wrote: On Thu, Oct 13, 2016 at 12:34 PM, Rossi, Simone <sro...@email.unc

[Libmesh-users] ksp monitor

2016-11-03 Thread Rossi, Simone
Dear all, I have been trying to run some of the time dependent examples. Every time I use the -ksp_monitor command line option (with PETSc) to see the output of the linear system, after 4 iterations the simulation crashes saying that there are too many kip_monitors. The same thing happen when

Re: [Libmesh-users] Systems of Equations, Ex4 -- Ex5

2016-10-13 Thread Rossi, Simone
e On Oct 13, 2016, at 12:00, David Knezevic <david.kneze...@akselos.com<mailto:david.kneze...@akselos.com>> wrote: On Thu, Oct 13, 2016 at 11:57 AM, Rossi, Simone <sro...@email.unc.edu<mailto:sro...@email.unc.edu>> wrote: Dear all, I’m playing around with the elasticit

[Libmesh-users] Systems of Equations, Ex4 -- Ex5

2016-10-13 Thread Rossi, Simone
Dear all, I’m playing around with the elasticity tests in the system of equations examples (more specifically Ex4 and Ex5). In particular I’m trying to set the poisson ratio to nu = 0.49995. With this choice the solution I get using second order lagrangian elements does not make any sense. For

Re: [Libmesh-users] Surface PDEs

2016-10-07 Thread Rossi, Simone
Thanks Paul, so it is already implemented. It sounds awesome. Thanks again, Best Simone On Oct 7, 2016, at 13:34, Paul T. Bauman <ptbau...@gmail.com<mailto:ptbau...@gmail.com>> wrote: On Fri, Oct 7, 2016 at 12:36 PM, Rossi, Simone <sro...@email.unc.edu<mailto:sro...@ema

[Libmesh-users] Surface PDEs

2016-10-07 Thread Rossi, Simone
Dear all, I’d like to extend my libmesh-based reaction-diffusion solver to a surface solver, taking the trace of the 3D shape functions on the surface mesh. Before I start implementing a new element, do you have any suggestion on how to approach it? Thanks, Simone

[Libmesh-users] IO runtime error

2016-08-19 Thread Rossi, Simone
Hello, I’m experiencing a runtime error when export the equation systems if I use more than one EquationSystems with different variables in each system. In particular, I created 2 EquationSystems: systems1 with one P1 variable and systems2 with two P1 variables. After initializing, I write the