Hi all,
I ran the continental extension prm file yesterday and I noticed that the randomized initial strain of the following function does not work (anymore). When looking at the paraview file I see that the plastic strain is constant and equal to 1.5 inside the designated area (see attached screen capture). Has the syntax changed? is this a (known) bug?
Cheers,
Cedric
subsection Initial composition model
set Model name = function
subsection Function
set Variable names = x,y
set Function expression = 0;
if(x>50.e3 && x<150.e3 && y>50.e3, 0.5 + rand_seed(1), 0);
if(y>=80.e3, 1, 0);
if(y<80.e3 && y>=60.e3, 1, 0);
if(y<60.e3, 1, 0);
end
end
– This is ASPECT, the Advanced Solver for Problems in Earth’s ConvecTion.
– . version 2.5.0-pre (main, 50d0c5740)
– . using deal.II 9.4.0
– . with 32 bit indices and vectorization level 1 (128 bits)
– . using Trilinos 12.18.1
– . using p4est 2.3.2
– . running in DEBUG mode
– . running with 1 MPI process
Thanks for reporting the issue, that is really odd. I just ran the cookbook with the most recent ASPECT version and got the expected random distribution of initial plastic strain. This was done on an ubuntu workstation.
I’m not sure what the issue could be here? Just to be sure, can you try reloading and rescaling the color bar?
Cheers,
John
– This is ASPECT, the Advanced Solver for Problems in Earth’s ConvecTion.
– . version 2.5.0-pre (main, bb7beb982)
– . using deal.II 9.4.0
– . with 32 bit indices and vectorization level 0 (64 bits)
– . using Trilinos 13.2.0
– . using p4est 2.2.0
– . running in DEBUG mode
– . running with 4 MPI processes
hi John,
the strain field is really ‘flat’ I am afraid. I am also running ubuntu. I will reinstall everything from scratch on my other workstation then.
Ce/
Oofta, that is really weird. There issue is likely with muparser, but that is built with deal.II and we are both on 9.4.0. What versions of gnu compilers are you on?
Cedric,
This is indeed weird. I have not seen anything like that before and I doubt recompiling everything would fix this. To be sure, this is the unmodified .prm file from the benchmarks folder?
(is always returning 1 not a decent random number generator?)
So I changed computer (still ubuntu machine), and installed candi from scratch today and re-compiled aspect. I ran the prm file just now without changing anything and still the same problem.
– This is ASPECT, the Advanced Solver for Problems in Earth’s ConvecTion.
– . version 2.5.0-pre (main, bb7beb982)
– . using deal.II 9.4.1
– . with 32 bit indices and vectorization level 1 (128 bits)
– . using Trilinos 13.2.0
– . using p4est 2.3.2
– . running in DEBUG mode
– . running with 1 MPI process
I am not able to reproduce so far. Can you try the following small deal.II program and report the output (if it returns only 1.0 then we see the same bug):