[aspect-devel] cryptic error in particle generation

[Bart Niday]

Hi all, 

I’m trying to use the particle probability density function to generate a random uniform distribution over part of the model domain, but I’ve been running into this error:

Number of active cells: 16 (on 3 levels)
Number of degrees of freedom: 268 (162+25+81)

Assertion failed: (min_arg < max_arg), function uniform_real_distribution, file /Applications/deal.II-8.5-brew.app/Contents/Resources/brew/include/boost/random/uniform_real_distribution.hpp, line 159.

The attached parameter file replicates the problem, at least for me. If I change the probability density function from if(y>0.5,1,0) to if(x>0.5,1,0), it works normally. Changing y>0.5 to a different number, sometimes it crashes and sometimes it works. The problem also happens on my desktop, but not on my university cluster (running a slightly different version of Aspect, but the particle code is unchanged).

Any clues as to what this might be? Does anyone else see the same problem running the attached file?

Thanks,



- Bart Niday
Ph.D. Student
Department of Earth Sciences, University of Oregon
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20180611/b67ecdcf/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: particle-error.prm
Type: application/octet-stream
Size: 792 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20180611/b67ecdcf/attachment.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20180611/b67ecdcf/attachment-0001.html>