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 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 iterations take about 18
seconds with a single processor.
With AMR, using a 2 x 2 x 2 grid and 3 levels of refinement, 100 iterations
take about 800 seconds.
I’m attaching the code I’m using.
Without AMR, I build the matrix ( mass + dt * stiffness ) once and I update the
rhs at every timestep.
Conversely, with AMR I am building the matrix and the rhs at every timestep for
all the refinement levels.
Do you have any suggestions?
Thanks a lot for your help,
All the best,
Simone
elX = 16
elY = 16
elZ = 16
num_steps = 101
dt = 0.01
D = 3e-3
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
