Skip to content
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

Uf23 #490

Merged
merged 6 commits into from
Jun 3, 2024
Merged

Uf23 #490

merged 6 commits into from
Jun 3, 2024

Conversation

MichaelUnger
Copy link
Contributor

added new models of the coherent magnetic field from 2311.12120

@lukasmerten
Copy link
Member

Hi @MichaelUnger
Thanks for providing your new magnetic field model. I will take a look at the code later this week and get back to you.

@MichaelUnger
Copy link
Contributor Author

MichaelUnger commented May 28, 2024

Hi @lukasmerten
thanks for looking into this! This week I am at a workshop with @asavelie @avvliet and @rafaelab so ideally it would be good to know if there any issues with the pull request, then I can solve them here with them here in person.

Cheers,
Michael

*/
UF23Field(const ModelType mt);
/// no default constructor
UF23Field() = delete;
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

You probably chose to not provide a default model because none of the models yielded a significantly better fit than the others, right?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yes, I think that not providing a default model makes the user aware that this is a collection of 8 models to choose from.

double& fPoloidalR = fParameters[ePoloidalR];
double& fPoloidalW = fParameters[ePoloidalW];
double& fPoloidalZ = fParameters[ePoloidalZ];
double& fPoloidalXi = fParameters[ePoloidalXi];
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

alignment

}

namespace crpropa {
UF23Field::UF23Field(const ModelType mt) :
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Could you add a comment that the UF23 model is internally using 1e-6*gauss=1, kpc=1, etc.

If you have a better place to make this statement I am also fine with it.

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Added a comment it at this point of cpp file, since it is internal to the code and not visible to the outside user.

@lukasmerten
Copy link
Member

Please add a comment to the CHANGELOG file under new features.

@lukasmerten
Copy link
Member

Apart from the small changes I suggested this PR looks good to me.

@MichaelUnger
Copy link
Contributor Author

Thanks for the review, I changed the fork accordingly.

@lukasmerten lukasmerten merged commit 95ff182 into CRPropa:master Jun 3, 2024
4 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants