I have written the script and it runs for a smaller RVE. But when I try to increase the RVE size and hence increase the number of spherical inclusions, it fails to run. I have attached the .py file herewith for your reference. It would be really helpful if someone could look into it and point out what I'm doing wrong.