합성생물학 유전자 리프레실레이터(Repressilator) 발진 & 1,000Hz 0-Alloc C# 전사-번역 지연 서보 엔진

using System;
using System.Runtime.CompilerServices;
using System.Runtime.InteropServices;

namespace ZeliDesk.Stem.SyntheticBiology
{
    /// <summary>
    /// 3채널 생체 형광 센서 신호 (GFP, YFP, RFP 계측값)
    /// </summary>
    [StructLayout(LayoutKind.Sequential, Pack = 8)]
    public readonly struct RepressilatorFluorescence
    {
        public readonly double GfpLacI;
        public readonly double YfpTetR;
        public readonly double RfpCI;
        public readonly double CultureTemperatureC;

        public RepressilatorFluorescence(double gfp, double yfp, double rfp, double tempC)
        {
            GfpLacI = gfp;
            YfpTetR = yfp;
            RfpCI = rfp;
            CultureTemperatureC = tempC;
        }
    }

    /// <summary>
    /// 미세유체 칩 유도물질 주입 밸브 제어 명령
    /// </summary>
    public ref struct BioInducerCommand
    {
        public double IptgPumpFlowNlLMin; // IPTG 분당 나노리터 주입량
        public double AtcPumpFlowNlLMin;  // aTc 분당 나노리터 주입량
        public double CurrentOscillationPeriodHours;
        public bool SynchronizationLocked;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public BioInducerCommand(double iptg, double atc, double period, bool locked)
        {
            IptgPumpFlowNlLMin = iptg;
            AtcPumpFlowNlLMin = atc;
            CurrentOscillationPeriodHours = period;
            SynchronizationLocked = locked;
        }
    }

    /// <summary>
    /// 1,000Hz 무할당 유전자 리프레실레이터 서보 엔진
    /// </summary>
    public ref struct RepressilatorServoEngine
    {
        private const double TargetPeriodHours = 2.8; // 목표 발진 주기 (2.8시간)
        private double _prevGfp;
        private double _peakTimeHours;
        private double _lastPeriodHours;
        private double _elapsedTimeHours;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void ExecuteControlCycle(
            in RepressilatorFluorescence bio,
            double dtSeconds,
            out BioInducerCommand command)
        {
            _elapsedTimeHours += dtSeconds / 3600.0;

            // 1. GFP 형광 피크 검출을 통한 발진 주기 계측 (극대점 판별)
            bool isPeak = (_prevGfp > bio.GfpLacI) && (_prevGfp > 120.0);
            if (isPeak && (_elapsedTimeHours - _peakTimeHours > 1.2)) // 최소 불응기 1.2시간
            {
                _lastPeriodHours = _elapsedTimeHours - _peakTimeHours;
                _peakTimeHours = _elapsedTimeHours;
            }
            _prevGfp = bio.GfpLacI;

            // 2. 주기 오차 계산 및 유도물질 투입 피드백 제어
            double periodError = TargetPeriodHours - _lastPeriodHours;
            double kp = 15.0;

            // 주기가 너무 길면 IPTG를 주입하여 억제 완화 (전사 가속)
            double iptgFlow = 0.0;
            double atcFlow = 0.0;

            if (periodError < -0.15) // 주기가 길어짐
            {
                iptgFlow = Math.Clamp(-periodError * kp, 0.0, 50.0);
            }
            else if (periodError > 0.15) // 주기가 너무 짧음
            {
                atcFlow = Math.Clamp(periodError * kp, 0.0, 50.0);
            }

            bool locked = Math.Abs(periodError) <= 0.15;
            command = new BioInducerCommand(iptgFlow, atcFlow, _lastPeriodHours, locked);
        }
    }
}