Equilibrium - Make the magnetic-axis search's singular-Jacobian threshold relative

Đang mở
  • Có gì không ổn?
#608 0 bình luận 0 reaction 0 người được giao Xem trên GitHub

Maintainer thường phản hồi trong vòng 3 ngày

Chưa có ai nhận issue này.

Đánh giá

Độ khó
5/5
Thời gian dự kiến
Hơn một tuần
Mức phù hợp với người mới
38/100
Loại issue
Tính năng
Độ rõ ràng
Khá rõ ràng
Mức độ hoạt động
Sôi nổi
Công nghệ
julia

Hướng nghiên cứu

Start with Equilibrium/Axisymmetric/Direct.jl to inspect the Newton search and its singular_det check. Then check FastInterpolations' implementation or documentation for whether 2-D cubic_interp is C² or C¹. Done means the singular-Jacobian threshold scales with the equilibrium, and the fallback accounts for interpolation behavior that can make Newton cycle.

Do mô hình lập chỉ mục viết ra từ nội dung của issue.

Mô tả

The Newton search for the magnetic axis (Equilibrium/Axisymmetric/Direct.jl) treats the Hessian of ψ as singular below singular_det = 1e-20, an absolute threshold in the units of ψ_RR·ψ_ZZ. Its meaning therefore depends on the device size and the field strength.

Proposal: compare the determinant with a scale built from the equilibrium, for example (ψ_axis − ψ_boundary)²/L⁴ with L the minor radius, so the same relative threshold applies to a compact device and to an ITER-sized one. Also confirm whether FastInterpolations' 2-D cubic_interp is C² or only C¹. If only C¹, the Hessian jumps at grid nodes and an axis near a node can make Newton cycle, which the fallback should then expect.

Ngôn ngữ chính
Julia
Star
11
Fork
7
Merge trung bình
4 ngày 1 giờ
Pull request đã merge (30 ngày)
68

Chuẩn bị môi trường

  • Không có Dockerfile hay tệp Docker Compose
  • Có mẫu pull request
  • Không có hướng dẫn đóng góp

Bắt đầu từ đâu

  1. Đọc hết issue, rồi đọc hướng dẫn đóng góp của dự án.
  2. Bình luận trên issue rằng bạn sẽ nhận — tránh hai người làm cùng một việc.
  3. Fork repository và làm thay đổi trên một nhánh.
  4. Mở pull request có tham chiếu số hiệu của issue.

Issue khác của OpenFUSIONToolkit/GPEC

Tất cả issue của OpenFUSIONToolkit/GPEC

Issue tương tự

Thêm issue về Julia

Nhận issue mới trong hộp thư của bạn

Bản tóm tắt ngắn những issue GitHub phù hợp với người mới.