-
Notifications
You must be signed in to change notification settings - Fork 225
Implement second ipd branch from scfly #5543
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
Implement second ipd branch from scfly #5543
Conversation
925926e to
e9c1ad4
Compare
|
@BrianMarre I made this a draft since there are still open PRs to be merged beforehand and I would like to avoid an "accidental" merge. |
|
@BrianMarre #5542 has been merged and you could start resolving the merge conflicts. |
1449b27 to
da59b45
Compare
PrometheusPi
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Various parts are unclearly documented or differently implemented as the cited reference.
Also some minor code changes have been introduced that appear unnecessary.
include/picongpu/particles/atomicPhysics/ionizationPotentialDepression/StewartPyattIPD.hpp
Outdated
Show resolved
Hide resolved
...e/picongpu/particles/atomicPhysics/ionizationPotentialDepression/stage/CalculateIPDInput.hpp
Show resolved
Hide resolved
include/picongpu/particles/atomicPhysics/debug/TestIonizationPotentialDepression.hpp
Outdated
Show resolved
Hide resolved
include/picongpu/particles/atomicPhysics/ionizationPotentialDepression/StewartPyattIPD.hpp
Outdated
Show resolved
Hide resolved
include/picongpu/particles/atomicPhysics/ionizationPotentialDepression/StewartPyattIPD.hpp
Show resolved
Hide resolved
include/picongpu/particles/atomicPhysics/ionizationPotentialDepression/StewartPyattIPD.hpp
Outdated
Show resolved
Hide resolved
include/picongpu/particles/atomicPhysics/ionizationPotentialDepression/StewartPyattIPD.hpp
Outdated
Show resolved
Hide resolved
include/picongpu/particles/atomicPhysics/ionizationPotentialDepression/StewartPyattIPD.hpp
Show resolved
Hide resolved
include/picongpu/particles/atomicPhysics/ionizationPotentialDepression/StewartPyattIPD.hpp
Show resolved
Hide resolved
...congpu/particles/atomicPhysics/ionizationPotentialDepression/kernel/CalculateIPDInput.kernel
Outdated
Show resolved
Hide resolved
|
@BrianMarre Could it be that you marked some issues raised as fixed but did not push these fixes yet? |
|
yes, sorry, give me a bit of time |
|
@PrometheusPi the new version including your comments is now ready for review. |
ci:picongpu
f25fcec to
f9835b8
Compare
ci:picongpu
f9835b8 to
255b167
Compare
PrometheusPi
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
just a typo - looks good.
| // (1/(1/unit_length^3))^(1/3) = unit_length | ||
| float_X const ionSphereRadius = math::pow( | ||
| 0.75_X * (chargeStateAsFloat + 1._X) / (picongpu::PI * superCellConstantInput.freeElectronDensity), | ||
| 1. / 3.); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
agrees with code
| float_X const stewartPyattIPD | ||
| = superCellConstantInput.temperatureTimesk_Boltzman | ||
| * (math::pow((3._X * (superCellConstantInput.zStar + 1._X) * K + 1._X), 2._X / 3._X) - 1._X) | ||
| / (2._X * (superCellConstantInput.zStar + 1._X)); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
agrees with equation
include/picongpu/particles/atomicPhysics/ionizationPotentialDepression/StewartPyattIPD.hpp
Outdated
Show resolved
Hide resolved
Co-authored-by: Richard Pausch <[email protected]>
a774def
into
ComputationalRadiationPhysics:dev
implements the second branch used in the SCFLY-IPD calculation