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

More changes #33

Merged
merged 36 commits into from
Nov 20, 2024
Merged

More changes #33

merged 36 commits into from
Nov 20, 2024

Conversation

dschwoerer
Copy link
Collaborator

No description provided.

dschwoerer and others added 30 commits August 16, 2023 14:34
It seems in previous numpy versions the index was allowed to be of type
float. This is not the case anymore and we need to explicitly cast to
int.
Float has only a limited accuracy, and this can lead to ind being n in
some cases.  Doing the module later in integer space avoids this issue,
although theoretically an overflow could happen. However, in that case
the accuracy is probably very poor.
They are, sorted fastes to slow:
 * _fast: extremly fast, but works only for sufficiently circular shapes
* _convex_hull: still very fast, but can fail on beans, if the bean is
very thin. Requires shapely
* _alphashape: reasonably fast. So far worked for all tested cases.
Should be much faster then lines_from_points(). Requires shapely and
alphashape
* _two_opt: Should find very good solutions, but is very slow. Mostly
here so that I do not add it again, because I forgot that it is very
slow.
Mostly for generating plots for a paper. It should probably not be used
in practice ...
J is likely (a bit) wrong, as it misses the By/Bxyz correction for g22
Use smooth volumes, based on the curvilinear grid.
This should improve agreement of the volumes with the normal
interpretation of the cells, and thus improve confinement in the normal
interpretation.
Replaced geqdsk file with one from freeqdsk
Not sure when this changed exactly
bshanahan and others added 6 commits November 8, 2024 10:15
Used wrong number of indices to rzcoord array returned by trace_poincare.

Now saves 3D metrics rather than 2D.
g22 was constant, but should have been a function of R
Using a different magnetic field fixes this
@bendudson bendudson merged commit 15d692f into master Nov 20, 2024
11 checks passed
@bendudson bendudson deleted the more-6 branch November 20, 2024 17:01
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.

4 participants