fm-rds-tx/internal/rds/encoder.go

package rds

import (
	"math"
	"sync/atomic"
	"time"
)

// RDS encoder — port of PiFmRds, adapted for arbitrary sample rates.
// At 228 kHz: uses exact {0,+1,0,-1} carrier (identical to PiFmRds).
// At other rates: uses sin() carrier, with integer samples-per-bit when possible.
// The biphase waveform from PiFmRds is resampled to the target rate.

const (
	refRate       = 228000.0  // PiFmRds native rate
	rdsBitRate    = 1187.5
	refSPB        = 192       // refRate / rdsBitRate
	bitsPerBlock  = 26
	bitsPerGroup  = 4 * bitsPerBlock
	refFilterSize = 576       // PiFmRds waveform length at 228k
)

const crcPoly = 0x1B9
var offsetWords = map[byte]uint16{'A': 0x0FC, 'B': 0x198, 'C': 0x168, 'c': 0x350, 'D': 0x1B4}

func crc10(data uint16) uint16 {
	var reg uint32 = uint32(data) << 10
	for i := 15; i >= 0; i-- {
		if reg&(1<<(uint(i)+10)) != 0 { reg ^= uint32(crcPoly) << uint(i) }
	}
	return uint16(reg & 0x3FF)
}
func encodeBlock(data uint16, offset byte) uint32 {
	return (uint32(data) << 10) | uint32(crc10(data)^offsetWords[offset])
}
// GroupScheduler manages priority-based scheduling of all RDS group types.
// Priority 1: 0A (PS/AF) and 2A (RT) — every cycle
// Priority 2: 4A (CT), 11A (RT+), 3A (ODA), 10A (PTYN) — every N cycles
// Priority 3: 14A (EON) — round-robin when slots available
type GroupScheduler struct {
	cfg RDSConfig

	// Group 0A state
	psIdx  int
	afPairs [][2]uint8 // precomputed AF code pairs

	// Group 2A state
	rtIdx    int
	rtABFlag bool

	// Group 10A state
	ptynIdx    int
	ptynABFlag bool

	// Group 11A state (RT+)
	rtPlusTag1 RTPlusTag
	rtPlusTag2 RTPlusTag
	rtPlusHas2 bool
	rtPlusToggle bool

	// Group 14A state (EON)
	eonStationIdx int
	eonVariantIdx int

	// Scheduling state
	cycle     int // overall group counter
	phase     int // position within current scheduling cycle
	lastCT    time.Time // last time CT was sent
}

func newGroupScheduler(cfg RDSConfig) *GroupScheduler {
	gs := &GroupScheduler{cfg: cfg}
	gs.afPairs = buildAFList(cfg.AF)
	if cfg.RTPlusEnabled && cfg.RT != "" {
		gs.rtPlusTag1, gs.rtPlusTag2, gs.rtPlusHas2 = ParseRTPlus(cfg.RT, cfg.RTPlusSeparator)
	}
	return gs
}

// refreshRTPlus re-parses RT+ tags after an RT text change.
func (gs *GroupScheduler) refreshRTPlus() {
	if gs.cfg.RTPlusEnabled && gs.cfg.RT != "" {
		gs.rtPlusTag1, gs.rtPlusTag2, gs.rtPlusHas2 = ParseRTPlus(gs.cfg.RT, gs.cfg.RTPlusSeparator)
		gs.rtPlusToggle = !gs.rtPlusToggle
	}
}

func (gs *GroupScheduler) NextGroup() [4]uint16 {
	gs.cycle++

	// --- Priority 1: 0A and 2A alternate ---
	// Pattern per cycle: 0A, 2A, 0A, 2A, 0A, 2A, 0A, 2A, [priority 2/3 slot]
	// This gives ~5.7 0A groups/sec and ~5.7 2A groups/sec at 11.4 groups/sec.

	// Every 8th group: try a priority 2/3 group
	if gs.cycle%8 == 0 {
		if g, ok := gs.nextPriority2Group(); ok {
			return g
		}
		// No priority 2/3 available — fall through to normal 0A/2A
	}

	// Alternate 0A and 2A
	if gs.cycle%2 == 0 {
		return gs.nextGroup0A()
	}
	return gs.nextGroup2A()
}

func (gs *GroupScheduler) nextGroup0A() [4]uint16 {
	var afPair [2]uint8
	if gs.psIdx < len(gs.afPairs) {
		afPair = gs.afPairs[gs.psIdx]
	}
	g := buildGroup0A(&gs.cfg, gs.psIdx, afPair)
	gs.psIdx = (gs.psIdx + 1) % 4
	return g
}

func (gs *GroupScheduler) nextGroup2A() [4]uint16 {
	g := buildGroup2A(gs.cfg.PI, gs.cfg.PTY, gs.cfg.TP, gs.rtABFlag, gs.rtIdx, gs.cfg.RT)
	gs.rtIdx++
	rtSegs := rtSegmentCount(gs.cfg.RT)
	if gs.rtIdx >= rtSegs {
		gs.rtIdx = 0
	}
	return g
}

// nextPriority2Group returns a lower-priority group if one is due.
// Round-robins through: CT → RT+ (3A+11A) → PTYN → EON
func (gs *GroupScheduler) nextPriority2Group() ([4]uint16, bool) {
	slot := (gs.cycle / 8) % 8 // 8 priority-2 slots, cycling

	switch {
	case slot == 0 && gs.cfg.CTEnabled:
		return gs.nextGroupCT(), true

	case slot == 1 && gs.cfg.RTPlusEnabled && gs.cfg.RT != "":
		// 3A: ODA announcement for RT+
		return buildGroup3A(gs.cfg.PI, gs.cfg.PTY, gs.cfg.TP, rtPlusODA, groupType11A), true

	case slot == 2 && gs.cfg.RTPlusEnabled && gs.cfg.RT != "":
		// 11A: RT+ tags
		return gs.nextGroupRTPlus(), true

	case slot == 3 && gs.cfg.PTYN != "":
		return gs.nextGroupPTYN(), true

	case slot >= 4 && len(gs.cfg.EON) > 0:
		return gs.nextGroupEON(), true
	}

	return [4]uint16{}, false
}

func (gs *GroupScheduler) nextGroupCT() [4]uint16 {
	now := time.Now().UTC()
	offset := gs.cfg.CTOffsetHalfHours
	if offset == 0 {
		// Auto-detect from OS timezone
		_, tzOffset := time.Now().Zone()
		offset = int8(tzOffset / 1800) // seconds → half-hours
	}
	gs.lastCT = now
	return buildGroup4A(gs.cfg.PI, gs.cfg.PTY, gs.cfg.TP, now, offset)
}

func (gs *GroupScheduler) nextGroupRTPlus() [4]uint16 {
	t1 := gs.rtPlusTag1
	t2 := gs.rtPlusTag2
	if !gs.rtPlusHas2 {
		t2 = RTPlusTag{} // dummy tag
	}
	// Length encoding: RT+ uses length-1 (0 = 1 char)
	t1Len := t1.Length
	if t1Len > 0 {
		t1Len--
	}
	t2Len := t2.Length
	if t2Len > 0 {
		t2Len--
	}
	return buildGroup11A(gs.cfg.PI, gs.cfg.PTY, gs.cfg.TP, true,
		t1.ContentType, t1.Start, t1Len,
		t2.ContentType, t2.Start, t2Len)
}

func (gs *GroupScheduler) nextGroupPTYN() [4]uint16 {
	g := buildGroup10A(gs.cfg.PI, gs.cfg.PTY, gs.cfg.TP, gs.ptynABFlag, gs.ptynIdx, gs.cfg.PTYN)
	gs.ptynIdx = (gs.ptynIdx + 1) % 2
	return g
}

func (gs *GroupScheduler) nextGroupEON() [4]uint16 {
	if len(gs.cfg.EON) == 0 {
		return gs.nextGroup0A() // fallback
	}
	eon := &gs.cfg.EON[gs.eonStationIdx]
	g := buildGroup14A(gs.cfg.PI, gs.cfg.PTY, gs.cfg.TP, gs.eonVariantIdx, eon)

	gs.eonVariantIdx++
	maxVariant := 4 // PS segments 0-3
	if len(eon.AF) >= 2 {
		maxVariant = 5 // +AF pair
	}
	if gs.eonVariantIdx >= maxVariant {
		gs.eonVariantIdx = 0
		gs.eonStationIdx = (gs.eonStationIdx + 1) % len(gs.cfg.EON)
	}
	return g
}

// Encoder generates RDS subcarrier samples at any sample rate.
// Uses the PiFmRds waveform resampled to the target rate.
type Encoder struct {
	config    RDSConfig
	scheduler *GroupScheduler
	sampleRate float64
	spb        int        // samples per bit at target rate
	waveform   []float64  // resampled PiFmRds waveform
	wfLen      int        // length of resampled waveform
	ring       []float64  // overlap-add ring buffer
	ringSize   int
	
	bitBuffer  [bitsPerGroup]int
	bitPos     int
	prevOutput int
	curOutput  int
	
	sampleCount  int
	inSampleIdx  int
	outSampleIdx int
	carrierPhase float64 // for sin() carrier at non-228k rates
	carrierStep  float64

	SampleRate float64

	// Live-updatable text — written by control API, read at group boundaries.
	// Zero-contention: atomic swap, checked once per RDS group (~88ms at 228kHz).
	// pendingText.set distinguishes "no pending update" from "update to empty string"
	// so that PS/RT can be explicitly cleared via UpdateText.
	livePS atomic.Value // pendingText
	liveRT atomic.Value // pendingText
	liveTA atomic.Int32 // -1=no change, 0=false, 1=true
	liveTP atomic.Int32 // -1=no change, 0=false, 1=true
}

// pendingText carries a pending text update for PS or RT.
// set=false means no update is pending; set=true means apply val (even if empty).
type pendingText struct {
	val string
	set bool
}

func NewEncoder(cfg RDSConfig) (*Encoder, error) {
	if cfg.SampleRate <= 0 { cfg.SampleRate = refRate }
	cfg.PS = normalizePS(cfg.PS); cfg.RT = normalizeRT(cfg.RT)
	
	rate := cfg.SampleRate
	spb := int(math.Round(rate / rdsBitRate))
	ratio := rate / refRate
	
	// Resample PiFmRds waveform to target rate
	wfLen := int(math.Round(float64(refFilterSize) * ratio))
	waveform := make([]float64, wfLen)
	for i := range waveform {
		srcPos := float64(i) / ratio
		idx := int(srcPos)
		frac := srcPos - float64(idx)
		if idx+1 < refFilterSize {
			waveform[i] = refWaveform[idx]*(1-frac) + refWaveform[idx+1]*frac
		} else if idx < refFilterSize {
			waveform[i] = refWaveform[idx]
		}
	}
	// Normalize to peak=1.0 so rdsInjection directly maps to injection %.
	// The raw PiFmRds waveform peaks at ~0.543, which would make config
	// values misleading (0.05 would give 2.7% instead of 5%).
	var peak float64
	for _, v := range waveform {
		if a := math.Abs(v); a > peak { peak = a }
	}
	if peak > 0 {
		for i := range waveform {
			waveform[i] /= peak
		}
	}
	
	ringSize := spb + wfLen
	
	enc := &Encoder{
		config:       cfg,
		scheduler:    newGroupScheduler(cfg),
		sampleRate:   rate,
		spb:          spb,
		waveform:     waveform,
		wfLen:        wfLen,
		ring:         make([]float64, ringSize),
		ringSize:     ringSize,
		bitPos:       bitsPerGroup,
		sampleCount:  spb,
		outSampleIdx: ringSize - 1,
		carrierStep:  57000.0 / rate,
		SampleRate:   rate,
	}
	enc.liveTA.Store(-1)
	enc.liveTP.Store(-1)
	return enc, nil
}

func (e *Encoder) Reset() {
	e.scheduler = newGroupScheduler(e.config)
	e.bitPos = bitsPerGroup; e.sampleCount = e.spb
	e.prevOutput = 0; e.curOutput = 0
	e.inSampleIdx = 0; e.outSampleIdx = e.ringSize - 1
	e.carrierPhase = 0
	for i := range e.ring { e.ring[i] = 0 }
}

// UpdateText hot-swaps PS and/or RT. Thread-safe — called from HTTP handlers,
// applied at the next RDS group boundary by the DSP goroutine.
//
// Pass empty string to leave a field unchanged. To explicitly clear a field
// (set PS to 8 spaces, or RT to empty), use ClearPS/ClearRT instead.
func (e *Encoder) UpdateText(ps, rt string) {
	if ps != "" {
		e.livePS.Store(pendingText{val: normalizePS(ps), set: true})
	}
	if rt != "" {
		e.liveRT.Store(pendingText{val: normalizeRT(rt), set: true})
	}
}

// ClearPS resets the Program Service name to 8 spaces at the next group boundary.
func (e *Encoder) ClearPS() {
	e.livePS.Store(pendingText{val: normalizePS(""), set: true})
}

// ClearRT resets RadioText to an empty string at the next group boundary.
// Per RDS spec, an empty RT causes receivers to clear their display.
func (e *Encoder) ClearRT() {
	e.liveRT.Store(pendingText{val: "", set: true})
}

// UpdateTA sets the Traffic Announcement flag live. Thread-safe.
// When TA goes true, receivers with TP-seek will interrupt CD playback.
func (e *Encoder) UpdateTA(ta bool) {
	if ta {
		e.liveTA.Store(1)
	} else {
		e.liveTA.Store(0)
	}
}

// UpdateTP sets the Traffic Program flag live. Thread-safe.
func (e *Encoder) UpdateTP(tp bool) {
	if tp {
		e.liveTP.Store(1)
	} else {
		e.liveTP.Store(0)
	}
}

// CurrentText returns the currently active PS and RT from the encoder scheduler.
// It reflects the last text applied at an RDS group boundary.
func (e *Encoder) CurrentText() (ps, rt string) {
	return e.scheduler.cfg.PS, e.scheduler.cfg.RT
}

// NextSample returns the next RDS subcarrier sample at the configured rate.
// Uses the internal free-running 57 kHz carrier. Prefer NextSampleWithCarrier
// for phase-locked operation in a stereo MPX chain.
func (e *Encoder) NextSample() float64 {
	carrier := math.Sin(2 * math.Pi * e.carrierPhase)
	e.carrierPhase += e.carrierStep
	if e.carrierPhase >= 1.0 { e.carrierPhase -= 1.0 }
	return e.NextSampleWithCarrier(carrier)
}

// NextSampleWithCarrier returns the next RDS sample modulated onto the
// supplied carrier value. The caller must provide sin(3 × pilotPhase × 2π)
// so that the 57 kHz RDS carrier is phase-locked to the 19 kHz pilot.
func (e *Encoder) NextSampleWithCarrier(carrier float64) float64 {
	if e.sampleCount >= e.spb {
		if e.bitPos >= bitsPerGroup {
			// Apply live text updates at group boundaries (~88ms at 228kHz).
			// Atomics are consumed (cleared) after reading to prevent
			// re-applying the same text every group and toggling A/B flag.
			if pt, ok := e.livePS.Load().(pendingText); ok && pt.set {
				e.scheduler.cfg.PS = pt.val
				e.livePS.Store(pendingText{}) // consumed
			}
			if pt, ok := e.liveRT.Load().(pendingText); ok && pt.set {
				e.scheduler.cfg.RT = pt.val
				e.scheduler.rtIdx = 0 // restart RT transmission for new text
				e.scheduler.rtABFlag = !e.scheduler.rtABFlag // toggle A/B per RDS spec
				e.scheduler.refreshRTPlus() // re-parse RT+ tags
				e.liveRT.Store(pendingText{}) // consumed
			}
			// Live TA/TP flags (no restart needed — applied immediately)
			if v := e.liveTA.Swap(-1); v >= 0 {
				e.scheduler.cfg.TA = v == 1
			}
			if v := e.liveTP.Swap(-1); v >= 0 {
				e.scheduler.cfg.TP = v == 1
			}
			e.getRDSGroup()
			e.bitPos = 0
		}
		// Differential encoding (PiFmRds-exact)
		curBit := e.bitBuffer[e.bitPos]
		e.prevOutput = e.curOutput
		e.curOutput = e.prevOutput ^ curBit
		inverting := (e.curOutput == 1)
		
		// Overlap-add waveform
		idx := e.inSampleIdx
		for j := 0; j < e.wfLen; j++ {
			val := e.waveform[j]
			if inverting { val = -val }
			e.ring[idx] += val
			idx++; if idx >= e.ringSize { idx = 0 }
		}
		e.inSampleIdx += e.spb
		if e.inSampleIdx >= e.ringSize { e.inSampleIdx -= e.ringSize }
		e.bitPos++
		e.sampleCount = 0
	}
	
	// Read envelope from ring buffer
	envelope := e.ring[e.outSampleIdx]
	e.ring[e.outSampleIdx] = 0
	e.outSampleIdx++; if e.outSampleIdx >= e.ringSize { e.outSampleIdx = 0 }
	
	e.sampleCount++
	return envelope * carrier
}

// NextSample228k is an alias for backward compat
func (e *Encoder) NextSample228k() float64 { return e.NextSample() }
func (e *Encoder) Generate(n int) []float64 {
	out := make([]float64, n); for i := range out { out[i] = e.NextSample() }; return out
}
func (e *Encoder) Symbol() float64 {
	// Populate the bit buffer on first call (bitPos starts at bitsPerGroup
	// after NewEncoder/Reset, so the guard below would return -1 immediately
	// without this bootstrap step).
	if e.bitPos >= bitsPerGroup {
		e.getRDSGroup()
		e.bitPos = 0
	}
	sym := 1.0
	if e.bitBuffer[e.bitPos] == 0 {
		sym = -1.0
	}
	e.sampleCount++
	if e.sampleCount >= e.spb {
		e.sampleCount = 0
		e.bitPos++
		if e.bitPos >= bitsPerGroup {
			e.getRDSGroup()
			e.bitPos = 0
		}
	}
	return sym
}

func (e *Encoder) getRDSGroup() {
	group := e.scheduler.NextGroup(); pos := 0
	for blk, off := range [4]byte{'A', 'B', 'C', 'D'} {
		encoded := encodeBlock(group[blk], off)
		for bit := bitsPerBlock - 1; bit >= 0; bit-- {
			e.bitBuffer[pos] = int((encoded >> uint(bit)) & 1); pos++
		}
	}
}


// refWaveform — exact PiFmRds waveform_biphase (576 samples at 228 kHz)
var refWaveform = [576]float64{
	2.532651330219518743e-03, 2.555044910373570170e-03, 2.566671021257086772e-03, 2.567238549701184643e-03,
	2.556496746672830885e-03, 2.534237165729987841e-03, 2.500295472534516759e-03, 2.454553115509106563e-03,
	2.396938848057333059e-03, 2.327430093141364103e-03, 2.246054141429317842e-03, 2.152889174683249481e-03,
	2.048065106558019741e-03, 1.931764233519830953e-03, 1.804221689169753116e-03, 1.665725695871127275e-03,
	1.516617608227482320e-03, 1.357291743639731877e-03, 1.188194995883951532e-03, 1.009826228393920439e-03,
	8.227354447016149448e-04, 6.275227342843667493e-04, 4.248369928834821877e-04, 2.153744171968115500e-04,
	-1.232252981576898894e-07, -2.208705497650993482e-04, -4.460407281680368592e-04, -6.747678807736505729e-04,
	-9.061496807071534867e-04, -1.139250166374730904e-03, -1.373102755669571148e-03, -1.606713454985107814e-03,
	-1.839064255174424622e-03, -2.069116705719896768e-03, -2.295815657515722953e-03, -2.518093163822523149e-03,
	-2.734872528130168675e-03, -2.945072486864876073e-03, -3.147611514104867864e-03, -3.341412234726898710e-03,
	-3.525405931697330516e-03, -3.698537132549940647e-03, -3.859768259460998226e-03, -4.008084326742814528e-03,
	-4.142497669033530158e-03, -4.262052682966127777e-03, -4.365830564655405997e-03, -4.452954024951411884e-03,
	-4.522591964073326663e-03, -4.573964086961633882e-03, -4.606345440470235343e-03, -4.619070853366329284e-03,
	-4.611539260015793187e-03, -4.583217888607013166e-03, -4.533646294807957243e-03, -4.462440221861094410e-03,
	-4.369295268298792508e-03, -4.253990344709691375e-03, -4.116390901303888447e-03, -3.956451908412099601e-03,
	-3.774220572511344761e-03, -3.569838770897127041e-03, -3.343545188717922216e-03, -3.095677142753055708e-03,
	-2.826672077047463587e-03, -2.537068716315400187e-03, -2.227507863888385838e-03, -1.898732831910291503e-03,
	-1.551589492469313079e-03, -1.187025939403796478e-03, -8.060917516227311674e-04, -4.099368499382796552e-04,
	1.900593832358028364e-07, 4.229434598897815903e-04, 8.568834232677739271e-04, 1.300478433616307234e-03,
	1.752108639186909398e-03, 2.210069531690153112e-03, 2.672576051828062533e-03, 3.137767118244720512e-03,
	3.603710575596331869e-03, 4.068408555938784511e-03, 4.529803246113349099e-03, 4.985783052286746689e-03,
	5.434189151275972511e-03, 5.872822416766181226e-03, 6.299450707012693891e-03, 6.711816499117106809e-03,
	7.107644853482626354e-03, 7.484651690592414135e-03, 7.840552360821295350e-03, 8.173070486592225875e-03,
	8.479947054827158270e-03, 8.758949736324453048e-03, 9.007882407426009638e-03, 9.224594848121719579e-03,
	9.406992589581098657e-03, 9.553046883008396370e-03, 9.660804760689966145e-03, 9.728399159148424027e-03,
	9.754059073439845171e-03, 9.736119710831808022e-03, 9.673032611387005070e-03, 9.563375702351037400e-03,
	9.405863252709357331e-03, 9.199355693838612638e-03, 8.942869271836061465e-03, 8.635585496870153838e-03,
	8.276860354756325477e-03, 7.866233245929400361e-03, 7.403435617058454557e-03, 6.888399250731210011e-03,
	6.321264178928198696e-03, 5.702386186409919011e-03, 5.032343870653735625e-03, 4.311945225603512621e-03,
	3.542233717233159249e-03, 2.724493819771388933e-03, 1.860255982396226328e-03, 9.513009972738021430e-04,
	-3.362590087872507610e-07, -9.923637372837323424e-04, -2.022229806201048148e-03, -3.087121544982383836e-03,
	-4.183963760744797630e-03, -5.309418750929982035e-03, -6.459886829603163697e-03, -7.631507634530486361e-03,
	-8.820162230001177273e-03, -1.002147601834123269e-02, -1.123082247097366773e-02, -1.244332768771695338e-02,
	-1.365387579078745396e-02, -1.485711515769077770e-02, -1.604746549484995649e-02, -1.721912575143924595e-02,
	-1.836608287047301891e-02, -1.948212137174597930e-02, -2.056083375874448449e-02, -2.159563173915068912e-02,
	-2.257975824605283968e-02, -2.350630024446240945e-02, -2.436820230522675906e-02, -2.515828092592677437e-02,
	-2.586923957586502801e-02, -2.649368443979239693e-02, -2.702414083259563338e-02, -2.745307025478293042e-02,
	-2.777288805626721910e-02, -2.797598167366384739e-02, -2.805472940409954943e-02, -2.800151967637884431e-02,
	-2.780877077828078359e-02, -2.746895099676937430e-02, -2.697459912600155829e-02, -2.631834529621684959e-02,
	-2.549293207489373297e-02, -2.449123578997222314e-02, -2.330628802347053594e-02, -2.193129722247495056e-02,
	-2.035967037326144244e-02, -1.858503468321714980e-02, -1.660125921427981280e-02, -1.440247641080375819e-02,
	-1.198310346409471074e-02, -9.337863455345842695e-03, -6.461806218342813941e-03, -3.350328863101498331e-03,
	8.040984433930472399e-07, 3.595441083468180719e-03, 7.437024284502863521e-03, 1.152857108893870150e-02,
	1.587265612925233688e-02, 2.047139402062919666e-02, 2.532642284173815955e-02, 3.043888841331177791e-02,
	3.580942942793266526e-02, 4.143816348300550373e-02, 4.732467406731786369e-02, 5.346799854986208217e-02,
	5.986661721770201311e-02, 6.651844340763700403e-02, 7.342081477423043068e-02, 8.057048573445227402e-02,
	8.796362112672609368e-02, 9.559579111958180220e-02, 1.034619674024025021e-01, 1.115565206879204074e-01,
	1.198732195531727329e-01, 1.284052306425759737e-01, 1.371451202536290992e-01, 1.460848573225170255e-01,
	1.552158178235605868e-01, 1.645287905930657713e-01, 1.740139845846588318e-01, 1.836610375596755829e-01,
	1.934590262126028026e-01, 2.033964777279658187e-01, 2.134613827614141035e-01, 2.236412098341045707e-01,
	2.339229211258132823e-01, 2.442929896485369901e-01, 2.547374177786905780e-01, 2.652417571223665282e-01,
	2.757911296845135252e-01, 2.863702503093192853e-01, 2.969634503555686478e-01, 3.075547025672780710e-01,
	3.181276470965210823e-01, 3.286656186320493500e-01, 3.391516745840755798e-01, 3.495686242724773130e-01,
	3.598990590626389707e-01, 3.701253833902611312e-01, 3.802298466136770361e-01, 3.901945756295767120e-01,
	4.000016081911345056e-01, 4.096329268202659746e-01, 4.190704933607786176e-01, 4.282962839026540069e-01,
	4.372923241989370990e-01, 4.460407253805190875e-01, 4.545237199298762798e-01, 4.627236978285061975e-01,
	4.706232427973597865e-01, 4.782051685485750880e-01, 4.854525549661181105e-01, 4.923487841323866965e-01,
	4.988775761175122669e-01, 5.050230244479789743e-01, 5.107696311712632831e-01, 5.161023414335079718e-01,
	5.210065774877549183e-01, 5.254682720509856741e-01, 5.294739009291574705e-01, 5.330105148305670504e-01,
	5.360657702892279719e-01, 5.386279596215074461e-01, 5.406860398410284763e-01, 5.422296604588032753e-01,
	5.432491900977250987e-01, 5.437357418528805386e-01, 5.436811973316860724e-01, 5.430782293105560488e-01,
	5.419203229476984296e-01, 5.402017954946758405e-01, 5.379178144525937899e-01, 5.350644141221254646e-01,
	5.316385105000952516e-01, 5.276379144789866693e-01, 5.230613433094983833e-01, 5.179084302901690862e-01,
	5.121797326520726168e-01, 5.058767376106826363e-01, 4.990018665611814508e-01, 4.915584773977406674e-01,
	4.835508649416218607e-01, 4.749842594673260865e-01, 4.658648233204350508e-01, 4.561996456252422338e-01,
	4.459967350847497403e-01, 4.352650108800791839e-01, 4.240142916808554152e-01, 4.122552827825399224e-01,
	3.999995613968218566e-01, 3.872595600805870397e-01, 3.740485484506272384e-01, 3.603806130210642222e-01,
	3.462706353977916263e-01, 3.317342687539049928e-01, 3.167879126713797899e-01, 3.014486863933045213e-01,
	2.857344005404113818e-01, 2.696635273493745433e-01, 2.532551694939167986e-01, 2.365290275532226649e-01,
	2.195053661954201318e-01, 2.022049791470396651e-01, 1.846491530223047517e-01, 1.668596300888892936e-01,
	1.488585700493792463e-01, 1.306685109200720063e-01, 1.123123290909592703e-01, 9.381319865274144465e-02,
	7.519455007851447159e-02, 5.648002834935082067e-02, 3.769345061436135680e-02, 1.885876347696378158e-02,
	0.000000000000000000e+00, -1.885876347696378158e-02, -3.769345061436135680e-02, -5.648002834935082067e-02,
	-7.519455007851447159e-02, -9.381319865274144465e-02, -1.123123290909592703e-01, -1.306685109200720063e-01,
	-1.488585700493792463e-01, -1.668596300888892936e-01, -1.846491530223047517e-01, -2.022049791470396651e-01,
	-2.195053661954201318e-01, -2.365290275532226649e-01, -2.532551694939167986e-01, -2.696635273493745433e-01,
	-2.857344005404113818e-01, -3.014486863933045213e-01, -3.167879126713797899e-01, -3.317342687539049928e-01,
	-3.462706353977916263e-01, -3.603806130210642222e-01, -3.740485484506272384e-01, -3.872595600805870397e-01,
	-3.999995613968218566e-01, -4.122552827825399224e-01, -4.240142916808554152e-01, -4.352650108800791839e-01,
	-4.459967350847497403e-01, -4.561996456252422338e-01, -4.658648233204350508e-01, -4.749842594673260865e-01,
	-4.835508649416218607e-01, -4.915584773977406674e-01, -4.990018665611814508e-01, -5.058767376106826363e-01,
	-5.121797326520726168e-01, -5.179084302901690862e-01, -5.230613433094983833e-01, -5.276379144789866693e-01,
	-5.316385105000952516e-01, -5.350644141221254646e-01, -5.379178144525937899e-01, -5.402017954946758405e-01,
	-5.419203229476984296e-01, -5.430782293105560488e-01, -5.436811973316860724e-01, -5.437357418528805386e-01,
	-5.432491900977250987e-01, -5.422296604588032753e-01, -5.406860398410284763e-01, -5.386279596215074461e-01,
	-5.360657702892279719e-01, -5.330105148305670504e-01, -5.294739009291574705e-01, -5.254682720509856741e-01,
	-5.210065774877549183e-01, -5.161023414335079718e-01, -5.107696311712632831e-01, -5.050230244479789743e-01,
	-4.988775761175122669e-01, -4.923487841323866965e-01, -4.854525549661181105e-01, -4.782051685485750880e-01,
	-4.706232427973597865e-01, -4.627236978285061975e-01, -4.545237199298762798e-01, -4.460407253805190875e-01,
	-4.372923241989370990e-01, -4.282962839026540069e-01, -4.190704933607786176e-01, -4.096329268202659746e-01,
	-4.000016081911345056e-01, -3.901945756295767120e-01, -3.802298466136770361e-01, -3.701253833902611312e-01,
	-3.598990590626389707e-01, -3.495686242724773130e-01, -3.391516745840755798e-01, -3.286656186320493500e-01,
	-3.181276470965210823e-01, -3.075547025672780710e-01, -2.969634503555686478e-01, -2.863702503093192853e-01,
	-2.757911296845135252e-01, -2.652417571223665282e-01, -2.547374177786905780e-01, -2.442929896485369901e-01,
	-2.339229211258132823e-01, -2.236412098341045707e-01, -2.134613827614141035e-01, -2.033964777279658187e-01,
	-1.934590262126028026e-01, -1.836610375596755829e-01, -1.740139845846588318e-01, -1.645287905930657713e-01,
	-1.552158178235605868e-01, -1.460848573225170255e-01, -1.371451202536290992e-01, -1.284052306425759737e-01,
	-1.198732195531727329e-01, -1.115565206879204074e-01, -1.034619674024025021e-01, -9.559579111958180220e-02,
	-8.796362112672609368e-02, -8.057048573445227402e-02, -7.342081477423043068e-02, -6.651844340763700403e-02,
	-5.986661721770201311e-02, -5.346799854986208217e-02, -4.732467406731786369e-02, -4.143816348300550373e-02,
	-3.580942942793266526e-02, -3.043888841331177791e-02, -2.532642284173815955e-02, -2.047139402062919666e-02,
	-1.587265612925233688e-02, -1.152857108893870150e-02, -7.437024284502863521e-03, -3.595441083468180719e-03,
	-8.040984433930472399e-07, 3.350328863101498331e-03, 6.461806218342813941e-03, 9.337863455345842695e-03,
	1.198310346409471074e-02, 1.440247641080375819e-02, 1.660125921427981280e-02, 1.858503468321714980e-02,
	2.035967037326144244e-02, 2.193129722247495056e-02, 2.330628802347053594e-02, 2.449123578997222314e-02,
	2.549293207489373297e-02, 2.631834529621684959e-02, 2.697459912600155829e-02, 2.746895099676937430e-02,
	2.780877077828078359e-02, 2.800151967637884431e-02, 2.805472940409954943e-02, 2.797598167366384739e-02,
	2.777288805626721910e-02, 2.745307025478293042e-02, 2.702414083259563338e-02, 2.649368443979239693e-02,
	2.586923957586502801e-02, 2.515828092592677437e-02, 2.436820230522675906e-02, 2.350630024446240945e-02,
	2.257975824605283968e-02, 2.159563173915068912e-02, 2.056083375874448449e-02, 1.948212137174597930e-02,
	1.836608287047301891e-02, 1.721912575143924595e-02, 1.604746549484995649e-02, 1.485711515769077770e-02,
	1.365387579078745396e-02, 1.244332768771695338e-02, 1.123082247097366773e-02, 1.002147601834123269e-02,
	8.820162230001177273e-03, 7.631507634530486361e-03, 6.459886829603163697e-03, 5.309418750929982035e-03,
	4.183963760744797630e-03, 3.087121544982383836e-03, 2.022229806201048148e-03, 9.923637372837323424e-04,
	3.362590087872507610e-07, -9.513009972738021430e-04, -1.860255982396226328e-03, -2.724493819771388933e-03,
	-3.542233717233159249e-03, -4.311945225603512621e-03, -5.032343870653735625e-03, -5.702386186409919011e-03,
	-6.321264178928198696e-03, -6.888399250731210011e-03, -7.403435617058454557e-03, -7.866233245929400361e-03,
	-8.276860354756325477e-03, -8.635585496870153838e-03, -8.942869271836061465e-03, -9.199355693838612638e-03,
	-9.405863252709357331e-03, -9.563375702351037400e-03, -9.673032611387005070e-03, -9.736119710831808022e-03,
	-9.754059073439845171e-03, -9.728399159148424027e-03, -9.660804760689966145e-03, -9.553046883008396370e-03,
	-9.406992589581098657e-03, -9.224594848121719579e-03, -9.007882407426009638e-03, -8.758949736324453048e-03,
	-8.479947054827158270e-03, -8.173070486592225875e-03, -7.840552360821295350e-03, -7.484651690592414135e-03,
	-7.107644853482626354e-03, -6.711816499117106809e-03, -6.299450707012693891e-03, -5.872822416766181226e-03,
	-5.434189151275972511e-03, -4.985783052286746689e-03, -4.529803246113349099e-03, -4.068408555938784511e-03,
	-3.603710575596331869e-03, -3.137767118244720512e-03, -2.672576051828062533e-03, -2.210069531690153112e-03,
	-1.752108639186909398e-03, -1.300478433616307234e-03, -8.568834232677739271e-04, -4.229434598897815903e-04,
	-1.900593832358028364e-07, 4.099368499382796552e-04, 8.060917516227311674e-04, 1.187025939403796478e-03,
	1.551589492469313079e-03, 1.898732831910291503e-03, 2.227507863888385838e-03, 2.537068716315400187e-03,
	2.826672077047463587e-03, 3.095677142753055708e-03, 3.343545188717922216e-03, 3.569838770897127041e-03,
	3.774220572511344761e-03, 3.956451908412099601e-03, 4.116390901303888447e-03, 4.253990344709691375e-03,
	4.369295268298792508e-03, 4.462440221861094410e-03, 4.533646294807957243e-03, 4.583217888607013166e-03,
	4.611539260015793187e-03, 4.619070853366329284e-03, 4.606345440470235343e-03, 4.573964086961633882e-03,
	4.522591964073326663e-03, 4.452954024951411884e-03, 4.365830564655405997e-03, 4.262052682966127777e-03,
	4.142497669033530158e-03, 4.008084326742814528e-03, 3.859768259460998226e-03, 3.698537132549940647e-03,
	3.525405931697330516e-03, 3.341412234726898710e-03, 3.147611514104867864e-03, 2.945072486864876073e-03,
	2.734872528130168675e-03, 2.518093163822523149e-03, 2.295815657515722953e-03, 2.069116705719896768e-03,
	1.839064255174424622e-03, 1.606713454985107814e-03, 1.373102755669571148e-03, 1.139250166374730904e-03,
	9.061496807071534867e-04, 6.747678807736505729e-04, 4.460407281680368592e-04, 2.208705497650993482e-04,
	1.232252981576898894e-07, -2.153744171968115500e-04, -4.248369928834821877e-04, -6.275227342843667493e-04,
	-8.227354447016149448e-04, -1.009826228393920439e-03, -1.188194995883951532e-03, -1.357291743639731877e-03,
	-1.516617608227482320e-03, -1.665725695871127275e-03, -1.804221689169753116e-03, -1.931764233519830953e-03,
	-2.048065106558019741e-03, -2.152889174683249481e-03, -2.246054141429317842e-03, -2.327430093141364103e-03,
	-2.396938848057333059e-03, -2.454553115509106563e-03, -2.500295472534516759e-03, -2.534237165729987841e-03,
	-2.556496746672830885e-03, -2.567238549701184643e-03, -2.566671021257086772e-03, -2.555044910373570170e-03,
}