극저온 방사광 X-선 나노회절 기반 마법각 꼬인 삼중층 그래핀(MATTG) 강상관 초전도 모아레 플랫밴드 능동 서보 아키텍처 (.NET 9 C# SIMD)

using System;
using System.Runtime.CompilerServices;
using System.Runtime.InteropServices;
using System.Runtime.Intrinsics;
using System.Runtime.Intrinsics.X86;

namespace ZeliDesk.Quantum.MoirSuperconductors
{
    /// <summary>
    /// 마법각 꼬인 삼중층 그래핀(MATTG) 강상관 초전도 모아레 플랫밴드 서보
    /// GC 무할당 ref struct 및 AVX-256 SIMD 가속화
    /// </summary>
    public ref struct MattgMoirSuperconductivityServo
    {
        private readonly ReadOnlySpan<float> _measuredAngles;     // 실측 국소 꼬임각 theta (deg)
        private readonly ReadOnlySpan<float> _carrierDensities;   // 국소 캐리어 밀도 n_e (10^12 cm^-2)
        private readonly Span<float> _topGateVoltages;           // 상부 게이트 전압 지령 (V)
        private readonly Span<float> _bottomGateVoltages;        // 하부 게이트 전압 지령 (V)

        private readonly float _magicAngle;     // 최적 마법각 (1.56 deg)
        private readonly float _targetDisp;     // 목표 수직 변위장 D (V/nm)
        private readonly float _angleGain;      // 각도 편차 보정 게인
        private readonly float _densityGain;    // 캐리어 밀도 보정 게인

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public MattgMoirSuperconductivityServo(
            ReadOnlySpan<float> measuredAngles,
            ReadOnlySpan<float> carrierDensities,
            Span<float> topGateVoltages,
            Span<float> bottomGateVoltages,
            float magicAngle = 1.56f,
            float targetDisp = 0.35f,
            float angleGain = 12.5f,
            float densityGain = 2.4f)
        {
            _measuredAngles = measuredAngles;
            _carrierDensities = carrierDensities;
            _topGateVoltages = topGateVoltages;
            _bottomGateVoltages = bottomGateVoltages;
            _magicAngle = magicAngle;
            _targetDisp = targetDisp;
            _angleGain = angleGain;
            _densityGain = densityGain;
        }

        /// <summary>
        /// 8개 모아레 도메인의 게이트 바이어스를 SIMD 병렬 연산
        /// </summary>
        [MethodImpl(MethodImplOptions.AggressiveOptimization)]
        public void ExecuteSuperconductivityLock()
        {
            int length = _measuredAngles.Length;
            int simdLength = length - (length % 8);

            Vector256<float> vMagic = Vector256.Create(_magicAngle);
            Vector256<float> vDisp = Vector256.Create(_targetDisp);
            Vector256<float> vAngleGain = Vector256.Create(_angleGain);
            Vector256<float> vDensGain = Vector256.Create(_densityGain);

            int i = 0;
            if (Avx2.IsSupported)
            {
                for (; i < simdLength; i += 8)
                {
                    Vector256<float> vTheta = Vector256.Create(MemoryMarshal.Cast<float, Vector256<float>>(_measuredAngles.Slice(i, 8))[0]);
                    Vector256<float> vNe = Vector256.Create(MemoryMarshal.Cast<float, Vector256<float>>(_carrierDensities.Slice(i, 8))[0]);

                    // Delta Theta = theta - magicAngle
                    Vector256<float> vDeltaTheta = Avx.Subtract(vTheta, vMagic);

                    // 각도 편차 피드백 항
                    Vector256<float> vAngleTerm = Avx.Multiply(vDeltaTheta, vAngleGain);
                    // 캐리어 밀도 피드백 항
                    Vector256<float> vDensTerm = Avx.Multiply(vNe, vDensGain);

                    // V_top = DensTerm + Disp + AngleTerm
                    Vector256<float> vTop = Avx.Add(Avx.Add(vDensTerm, vDisp), vAngleTerm);
                    // V_bottom = DensTerm - Disp - AngleTerm
                    Vector256<float> vBottom = Avx.Subtract(Avx.Subtract(vDensTerm, vDisp), vAngleTerm);

                    vTop.CopyTo(_topGateVoltages.Slice(i, 8));
                    vBottom.CopyTo(_bottomGateVoltages.Slice(i, 8));
                }
            }

            for (; i < length; i++)
            {
                float deltaTheta = _measuredAngles[i] - _magicAngle;
                float angleTerm = deltaTheta * _angleGain;
                float densTerm = _carrierDensities[i] * _densityGain;

                _topGateVoltages[i] = densTerm + _targetDisp + angleTerm;
                _bottomGateVoltages[i] = densTerm - _targetDisp - angleTerm;
            }
        }
    }
}