ERROR: Photon outside of NEXT cell #13
-
|
Hello, I am Rua Sulaiman, a grad Physics student. I have been trying to produce some simulations using Radmc3d for a disk with multiple rings to fit a certain SED. I was testing the model for higher mass for the inner ring. When I am using only 10^6 photons in the simulation everything goes well, however, for some values of the inner mass, when the photons number increases to 10^7 I start getting the following error: ERROR: Photon outside of NEXT cell I do not know exactly how to fix this. Can I get some idea of the main issue that I am facing so that I know exactly what to modify? Appreciated. Best, |
Beta Was this translation helpful? Give feedback.
Replies: 7 comments 5 replies
-
|
Dear Rua, |
Beta Was this translation helpful? Give feedback.
-
|
Hello again, Thank you for the quick answer! I will double-check with my parameters again. Thank you for pointing this out. Best, |
Beta Was this translation helpful? Give feedback.
-
|
Hi Rua,
could you send me the full error message? I would like to figure out what is going on here.
Cheers,
Kees
… On 18. Oct 2021, at 14:29, ruaSulaiman ***@***.***> wrote:
Hello,
I am Rua Sulaiman, a grad Physics student. I have been trying to produce some simulations using Radmc3d for a disk with multiple rings to fit a certain SED. I was testing the model for higher mass for the inner ring. When I am using only 10^6 photons in the simulation everything goes well, however, for some values of the inner mass, when the photons number increases to 10^7 I start getting the following error:
ERROR: Photon outside of NEXT cell
....
amrray_icross = 5 , cross_ds = 53392651762.109512 , idir = 3
Started inside of cell
....
sin2 cos2 = 1.0000000000000000 0.0000000000000000
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_OVERFLOW_FLAG IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP 1132
I do not know exactly how to fix this. Can I get some idea of the main issue that I am facing so that I know exactly what to modify?
Appreciated.
Best,
Rua
—
You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub <#13>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AFNQ7EUUNM4OJQ4QFWJEWMTUHQHMNANCNFSM5GGQZE2A>.
Triage notifications on the go with GitHub Mobile for iOS <https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android <https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>.
|
Beta Was this translation helpful? Give feedback.
-
|
Hi Rua, |
Beta Was this translation helpful? Give feedback.
-
|
Hi! The complete error message that I am getting is attached.
|
Beta Was this translation helpful? Give feedback.
-
|
Did you use ‘git checkout fixoutofcell’ ?
…Sent from my mobile phone
On 23. Oct 2021, at 18:37, ruaSulaiman ***@***.***> wrote:
Hi!
Thank you for following up on that! I have just tried it again and I am still getting the same error message (as long as switching the branch to fixoutofcell is all that is needed to use the new code files that you have written).
The complete error message that I am getting is attached.
—
You are receiving this because you commented.
Reply to this email directly, view it on GitHub, or unsubscribe.
Triage notifications on the go with GitHub Mobile for iOS or Android.
|
Beta Was this translation helpful? Give feedback.
-
|
Yes I did that |
Beta Was this translation helpful? Give feedback.
-
|
And recompiled?
…Sent from my mobile phone
On 23. Oct 2021, at 18:43, ruaSulaiman ***@***.***> wrote:
Yes I did that
—
You are receiving this because you commented.
Reply to this email directly, view it on GitHub, or unsubscribe.
Triage notifications on the go with GitHub Mobile for iOS or Android.
|
Beta Was this translation helpful? Give feedback.
-
|
I am sorry I forget to compile it in the beginning :) now I recompiled. Another message appeared and the code was aborted. Here is another screenshot. There was an attempt to correct the position of the photon, however, it did not work. Is this only because the mass of the inner ring is too high and so is still very unrealistic?
|
Beta Was this translation helpful? Give feedback.
-
|
I’ll have a look.
…Sent from my mobile phone
On 23. Oct 2021, at 19:29, ruaSulaiman ***@***.***> wrote:
I am sorry I forget to compile it in the beginning :) now I recompiled. Another message appeared and the code was aborted. Here is another screenshot. There was an attempt to correct the position of the photon, however, it did not work. Is this only because the mass of the inner ring is too high and so is still very unrealistic?
—
You are receiving this because you commented.
Reply to this email directly, view it on GitHub, or unsubscribe.
Triage notifications on the go with GitHub Mobile for iOS or Android.
|
Beta Was this translation helpful? Give feedback.
-
|
Sorry for having lost sight of this issue. Meanwhile we think the issue has been fixed in the branch |
Beta Was this translation helpful? Give feedback.
-
|
Hello! I know that this is a very old discussion but I wanted to bump this error up in priority! I am going to try to pull from 
|
Beta Was this translation helpful? Give feedback.
Dear Rua,
yes, sometimes that happens. It is related to the quadratic equations RADMC-3D has to solve for the spherical grid, and is very rare, but if you have very high optical depth and high number of photons, then even events that are very rare do happen. One day I should fix this once and for all, but usually one can avoid it if you have realistic optical depths and/or use the modified random walk method.
It is on my todo list, but because it does not happen often (and you can use a different random seed), it has not yet made the top of the list. But if more people complain, it might ;-).
Best regards,
Kees