-
Notifications
You must be signed in to change notification settings - Fork 23
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
How to solve the severe :"norm is drifting - with norm, normref=........."? #21
Comments
Dear user, When ray-tracing without multithreading (so that the messages arrive in order on the terminal), I see that this message appears only for a few pixels between row number 23 and row number 32. On rows 25 and 26, I count 14 such messages (each). The dark region in the center of the disk there is 20 pixels wide. I can make an educated guess that those messages arise from photons that come too close to the horizon so that time evolution stalls there. The vicinity of the horizon is also a place where the norm could legitimately diverge (numerically). To avoid such problems, we implement a sink radius in KerrBL::isStopCondition(). You could use the same mechanism. The "white edge" is another issue. Note that it is not filled with zeros, but with NaNs. This means that the integration worked correctly and that gyoto detects that the photon trajectory crosses the disk but that something fishy happens when computing the emission. I ran Gyoto in debug mode for a single pixel: My educated guess here is that you prescription for angular velocity is non physical at this distance (v>c?), more specifically this line must yield a division by zero or the square root of a negative number: |
Dear author, at Ww=0 Aa=-0.4, if you still get the warning : "SEVERE: Photon::hit: too many iterations (100001 vs. 100000), please do a simulation at Ww=0 Aa=0 and Ww=0 Aa=0.4 respectively Because the Thindisk need the parameter inner radius, I put the data for Ww, Aa and Inner Radius in this table, and This table may help you to reproduce my problem more clearly About the severe warning : "SEVERE: Photon::hit: too many iterations (100001 vs. 100000), break", I also encountered this warning, but only when I used my metric to calculate PageThorneDisk, and this question will not be discussed now And about the issue "white edge" , I checked my formula for angular velocity, and it doesn't seem to yield a division by zero or the square root of a negative number. I guess there's a limitation to my metric that this angular velocity is non-physical at Ww=0 Aa=-0.4, so this issue will not be discussed now. But the angular velocity seem to is physical at Ww=0 Aa=0,and I hope you to do a simulation at Ww=0 Aa=0. Due to my problem is very complicated, I am very sorry to bother you and hope you can forgive my disrespect and offense |
Dear Author, |
Dear user, it is defined in GyotoDefs.h:
Cheers |
Dear author, |
Dear author
I'm sorry to bother you, but I have a problem that is very confusing to me. When I used my metric to draw my black hole shadow, the black hole shadow is based on the thin disk model, I found that my picture was very strange, there was a white edge, it looked like it had not been calculated, And in the process of calculation, I got a large number of severe: "norm is drifting - with norm, normref=.........".
I looked at my metric a hundred times, and I still couldn't find any problem, and at the same time, I calculated the innermost stable orbital under the metric, and there is a innermost stable orbital under the parameter in my xml file. But from a theoretical analysis, it seems possible to draw a normal black hole shadow like a Schwarzschild black hole.
The specific details, such as the images I get, my plugins, and the xml files , I put in the zip package.
I would appreciate it if you could reply me in your busy schedule
my plugin.zip
The text was updated successfully, but these errors were encountered: