Tearing.Dispersion - Require a Hermitian Δ′ where the SLAYER residual conjugates Δ̂

Đang mở
  • Có gì không ổn?
#606 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
35/100
Loại issue
Lỗi
Độ rõ ràng
Khá rõ ràng
Mức độ hoạt động
Sôi nổi
Công nghệ
julia
Lĩnh vực
backend, testing

Hướng nghiên cứu

Start with InnerLayer/SLAYER/Riccati.jl and its solve_inner comment, then trace how the residual uses Δ′ in the uncoupled and coupled cases, including SurfaceCoupling.jl. Run the shipped cases to measure the proposed imaginary and anti-Hermitian components, and inspect the existing test setup for an anti-Cauchy–Riemann finite-difference check at one Q. Done means the condition or residual handling is settled, the comment states the condition, and the test verifies it.

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

Mô tả

SLAYER's solve_inner returns conj(Δ̂(ĝ)) (InnerLayer/SLAYER/Riccati.jl). The comment there says "conjugation maps zeros to zeros", so the roots are unchanged. That holds only if the outer data are effectively real.

  • Uncoupled: the residual is dp_diag − S^(1/3)·conj(Δ̂) − Δ_crit, so its zeros solve Δ̂ = conj((Δ′ − Δ_crit)/S^(1/3)). dp_diag is ComplexF64 (SurfaceCoupling.jl). If Im Δ′ ≠ 0, the roots found are those of the conjugated outer stability, not the intended one. Δ_crit is real, so it is not a concern.
  • Coupled: conj(det(D′ − diag Δ̂)) = det(conj D′ − diag conj Δ̂), so the coupled roots are right when conj(D′) is D′ up to transpose, i.e. when D′ is Hermitian. A non-Hermitian Δ′, for example from a kinetic or resistive-wall outer region, would shift them.

Proposal:

  1. Measure Im Δ′ (diagonal) and the anti-Hermitian part of the coupled Δ′ on the shipped cases.
  2. Either assert or warn when they exceed a stated fraction, or conjugate Δ′ consistently in the residual so the roots are exact for any Δ′.
  3. Fix the Riccati.jl comment to state the condition.
  4. Add a test that checks the premise: an anti-Cauchy–Riemann finite-difference check that SLAYER's residual is anti-holomorphic at one Q.

Note that the conjugation is load-bearing for consumers that use the sign of Im Δ (the critical-resonant-field torque balance), so it cannot simply be removed.

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

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.