Wideband autonomous SDR analysis engine forked from sdr-visual-suite
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

432 lines
10KB

  1. //go:build sdrplay
  2. package sdrplay
  3. /*
  4. #cgo windows LDFLAGS: -lsdrplay_api
  5. #cgo linux LDFLAGS: -lsdrplay_api
  6. #include "sdrplay_api.h"
  7. #include <stdlib.h>
  8. #include <string.h>
  9. extern void goStreamCallback(short *xi, short *xq, unsigned int numSamples, unsigned int reset, void *cbContext);
  10. static void StreamACallback(short *xi, short *xq, sdrplay_api_StreamCbParamsT *params, unsigned int numSamples, unsigned int reset, void *cbContext) {
  11. (void)params;
  12. goStreamCallback(xi, xq, numSamples, reset, cbContext);
  13. }
  14. static void EventCallback(sdrplay_api_EventT eventId, sdrplay_api_TunerSelectT tuner, sdrplay_api_EventParamsT *params, void *cbContext) {
  15. (void)eventId; (void)tuner; (void)params; (void)cbContext;
  16. }
  17. static sdrplay_api_CallbackFnsT sdrplay_get_callbacks() {
  18. sdrplay_api_CallbackFnsT cb;
  19. memset(&cb, 0, sizeof(cb));
  20. cb.StreamACbFn = StreamACallback;
  21. cb.StreamBCbFn = NULL;
  22. cb.EventCbFn = EventCallback;
  23. return cb;
  24. }
  25. static void sdrplay_set_fs(sdrplay_api_DeviceParamsT *p, double fsHz) {
  26. if (p && p->devParams) p->devParams->fsFreq.fsHz = fsHz;
  27. }
  28. static void sdrplay_set_rf(sdrplay_api_DeviceParamsT *p, double rfHz) {
  29. if (p && p->rxChannelA) p->rxChannelA->tunerParams.rfFreq.rfHz = rfHz;
  30. }
  31. static void sdrplay_set_gain(sdrplay_api_DeviceParamsT *p, unsigned int grDb) {
  32. if (p && p->rxChannelA) p->rxChannelA->tunerParams.gain.gRdB = grDb;
  33. }
  34. static void sdrplay_set_bw(sdrplay_api_DeviceParamsT *p, sdrplay_api_Bw_MHzT bw) {
  35. if (p && p->rxChannelA) p->rxChannelA->tunerParams.bwType = bw;
  36. }
  37. static void sdrplay_set_if_zero(sdrplay_api_DeviceParamsT *p) {
  38. if (p && p->rxChannelA) p->rxChannelA->tunerParams.ifType = sdrplay_api_IF_Zero;
  39. }
  40. static void sdrplay_disable_agc(sdrplay_api_DeviceParamsT *p) {
  41. if (p && p->rxChannelA) p->rxChannelA->ctrlParams.agc.enable = sdrplay_api_AGC_DISABLE;
  42. }
  43. static void sdrplay_set_agc(sdrplay_api_DeviceParamsT *p, int enable) {
  44. if (!p || !p->rxChannelA) return;
  45. if (enable) {
  46. p->rxChannelA->ctrlParams.agc.enable = sdrplay_api_AGC_100HZ;
  47. } else {
  48. p->rxChannelA->ctrlParams.agc.enable = sdrplay_api_AGC_DISABLE;
  49. }
  50. }
  51. static sdrplay_api_ErrT sdrplay_update(void *dev, int reason) {
  52. return sdrplay_api_Update(dev, sdrplay_api_Tuner_A, (sdrplay_api_ReasonForUpdateT)reason, sdrplay_api_Update_Ext1_None);
  53. }
  54. */
  55. import "C"
  56. import (
  57. "errors"
  58. "fmt"
  59. "runtime/cgo"
  60. "sync"
  61. "time"
  62. "unsafe"
  63. "sdr-visual-suite/internal/sdr"
  64. )
  65. type Source struct {
  66. mu sync.Mutex
  67. dev C.sdrplay_api_DeviceT
  68. params *C.sdrplay_api_DeviceParamsT
  69. devSelected bool
  70. ch chan []complex64
  71. handle cgo.Handle
  72. open bool
  73. sampleRate int
  74. centerHz float64
  75. gainDb float64
  76. agc bool
  77. buf []complex64
  78. capSamples int
  79. head int
  80. size int
  81. bwKHz int
  82. dropped uint64
  83. resets uint64
  84. lastSample time.Time
  85. cond *sync.Cond
  86. }
  87. func New(sampleRate int, centerHz float64, gainDb float64, bwKHz int) (sdr.Source, error) {
  88. s := &Source{
  89. ch: make(chan []complex64, 16),
  90. sampleRate: sampleRate,
  91. centerHz: centerHz,
  92. gainDb: gainDb,
  93. bwKHz: bwKHz,
  94. }
  95. s.cond = sync.NewCond(&s.mu)
  96. s.resizeBuffer(sampleRate, 0)
  97. s.handle = cgo.NewHandle(s)
  98. if err := s.configure(sampleRate, centerHz, gainDb, bwKHz); err != nil {
  99. _ = s.Stop()
  100. return nil, err
  101. }
  102. return s, nil
  103. }
  104. func (s *Source) configure(sampleRate int, centerHz float64, gainDb float64, bwKHz int) error {
  105. if err := cErr(C.sdrplay_api_Open()); err != nil {
  106. return fmt.Errorf("sdrplay_api_Open: %w", err)
  107. }
  108. s.open = true
  109. if err := cErr(C.sdrplay_api_LockDeviceApi()); err != nil {
  110. return fmt.Errorf("sdrplay_api_LockDeviceApi: %w", err)
  111. }
  112. var numDevs C.uint
  113. var devices [8]C.sdrplay_api_DeviceT
  114. if err := cErr(C.sdrplay_api_GetDevices(&devices[0], &numDevs, C.uint(len(devices)))); err != nil {
  115. _ = cErr(C.sdrplay_api_UnlockDeviceApi())
  116. return fmt.Errorf("sdrplay_api_GetDevices: %w", err)
  117. }
  118. if numDevs == 0 {
  119. _ = cErr(C.sdrplay_api_UnlockDeviceApi())
  120. return errors.New("no SDRplay devices found")
  121. }
  122. s.dev = devices[0]
  123. if err := cErr(C.sdrplay_api_SelectDevice(&s.dev)); err != nil {
  124. _ = cErr(C.sdrplay_api_UnlockDeviceApi())
  125. return fmt.Errorf("sdrplay_api_SelectDevice: %w", err)
  126. }
  127. s.devSelected = true
  128. _ = cErr(C.sdrplay_api_UnlockDeviceApi())
  129. var params *C.sdrplay_api_DeviceParamsT
  130. if err := cErr(C.sdrplay_api_GetDeviceParams(s.dev.dev, &params)); err != nil {
  131. return fmt.Errorf("sdrplay_api_GetDeviceParams: %w", err)
  132. }
  133. // Enable verbose debug from SDRplay service.
  134. _ = cErr(C.sdrplay_api_DebugEnable(s.dev.dev, C.sdrplay_api_DbgLvl_Verbose))
  135. s.params = params
  136. C.sdrplay_set_fs(s.params, C.double(sampleRate))
  137. C.sdrplay_set_rf(s.params, C.double(centerHz))
  138. C.sdrplay_set_gain(s.params, C.uint(gainDb))
  139. if bw := bwEnum(bwKHz); bw != 0 {
  140. C.sdrplay_set_bw(s.params, bw)
  141. if bwKHz > 0 {
  142. s.bwKHz = bwKHz
  143. }
  144. }
  145. C.sdrplay_set_if_zero(s.params)
  146. C.sdrplay_disable_agc(s.params)
  147. cb := C.sdrplay_get_callbacks()
  148. if err := cErr(C.sdrplay_api_Init(s.dev.dev, &cb, unsafe.Pointer(uintptr(s.handle)))); err != nil {
  149. return fmt.Errorf("sdrplay_api_Init: %w", err)
  150. }
  151. return nil
  152. }
  153. func (s *Source) Start() error { return nil }
  154. func (s *Source) UpdateConfig(sampleRate int, centerHz float64, gainDb float64, agc bool, bwKHz int) error {
  155. s.mu.Lock()
  156. defer s.mu.Unlock()
  157. if s.params == nil {
  158. return errors.New("sdrplay not initialized")
  159. }
  160. updateReasons := C.int(0)
  161. if sampleRate > 0 && sampleRate != s.sampleRate {
  162. C.sdrplay_set_fs(s.params, C.double(sampleRate))
  163. updateReasons |= C.int(C.sdrplay_api_Update_Dev_Fs)
  164. s.sampleRate = sampleRate
  165. s.resizeBuffer(sampleRate, 0)
  166. }
  167. if centerHz != 0 && centerHz != s.centerHz {
  168. C.sdrplay_set_rf(s.params, C.double(centerHz))
  169. updateReasons |= C.int(C.sdrplay_api_Update_Tuner_Frf)
  170. s.centerHz = centerHz
  171. }
  172. if gainDb != s.gainDb {
  173. C.sdrplay_set_gain(s.params, C.uint(gainDb))
  174. updateReasons |= C.int(C.sdrplay_api_Update_Tuner_Gr)
  175. s.gainDb = gainDb
  176. }
  177. if agc != s.agc {
  178. if agc {
  179. C.sdrplay_set_agc(s.params, 1)
  180. } else {
  181. C.sdrplay_set_agc(s.params, 0)
  182. }
  183. updateReasons |= C.int(C.sdrplay_api_Update_Ctrl_Agc)
  184. s.agc = agc
  185. }
  186. if bwKHz > 0 && bwKHz != s.bwKHz {
  187. if bw := bwEnum(bwKHz); bw != 0 {
  188. C.sdrplay_set_bw(s.params, bw)
  189. updateReasons |= C.int(C.sdrplay_api_Update_Tuner_BwType)
  190. s.bwKHz = bwKHz
  191. }
  192. }
  193. if updateReasons == 0 {
  194. return nil
  195. }
  196. if err := cErr(C.sdrplay_update(unsafe.Pointer(s.dev.dev), C.int(updateReasons))); err != nil {
  197. return err
  198. }
  199. return nil
  200. }
  201. func bwEnum(khz int) C.sdrplay_api_Bw_MHzT {
  202. switch khz {
  203. case 200:
  204. return C.sdrplay_api_BW_0_200
  205. case 300:
  206. return C.sdrplay_api_BW_0_300
  207. case 600:
  208. return C.sdrplay_api_BW_0_600
  209. case 1536:
  210. return C.sdrplay_api_BW_1_536
  211. case 5000:
  212. return C.sdrplay_api_BW_5_000
  213. case 6000:
  214. return C.sdrplay_api_BW_6_000
  215. case 7000:
  216. return C.sdrplay_api_BW_7_000
  217. case 8000:
  218. return C.sdrplay_api_BW_8_000
  219. default:
  220. return 0
  221. }
  222. }
  223. func (s *Source) resizeBuffer(sampleRate int, fftSize int) {
  224. capSamples := sampleRate
  225. if fftSize > 0 && fftSize*4 > capSamples {
  226. capSamples = fftSize * 4
  227. }
  228. if capSamples < 4096 {
  229. capSamples = 4096
  230. }
  231. if s.capSamples == capSamples && len(s.buf) == capSamples {
  232. return
  233. }
  234. newBuf := make([]complex64, capSamples)
  235. // copy existing data from ring
  236. toCopy := s.size
  237. if toCopy > capSamples {
  238. toCopy = capSamples
  239. }
  240. for i := 0; i < toCopy; i++ {
  241. newBuf[i] = s.buf[(s.head+i)%max(1, s.capSamples)]
  242. }
  243. s.buf = newBuf
  244. s.capSamples = capSamples
  245. s.head = 0
  246. s.size = toCopy
  247. }
  248. func (s *Source) appendRing(samples []complex64) {
  249. if len(samples) == 0 || s.capSamples == 0 {
  250. return
  251. }
  252. incoming := len(samples)
  253. over := s.size + incoming - s.capSamples
  254. if over > 0 {
  255. s.head = (s.head + over) % s.capSamples
  256. s.size -= over
  257. s.dropped += uint64(over)
  258. }
  259. start := (s.head + s.size) % s.capSamples
  260. first := min(incoming, s.capSamples-start)
  261. copy(s.buf[start:start+first], samples[:first])
  262. if first < incoming {
  263. copy(s.buf[0:incoming-first], samples[first:])
  264. }
  265. s.size += incoming
  266. if s.cond != nil {
  267. s.cond.Broadcast()
  268. }
  269. }
  270. func (s *Source) Stats() sdr.SourceStats {
  271. s.mu.Lock()
  272. defer s.mu.Unlock()
  273. ago := int64(-1)
  274. if !s.lastSample.IsZero() {
  275. ago = time.Since(s.lastSample).Milliseconds()
  276. }
  277. return sdr.SourceStats{BufferSamples: s.size, Dropped: s.dropped, Resets: s.resets, LastSampleAgoMs: ago}
  278. }
  279. func (s *Source) Flush() {
  280. s.mu.Lock()
  281. s.head = 0
  282. s.size = 0
  283. s.mu.Unlock()
  284. if s.cond != nil {
  285. s.cond.Broadcast()
  286. }
  287. }
  288. func min(a, b int) int {
  289. if a < b {
  290. return a
  291. }
  292. return b
  293. }
  294. func max(a, b int) int {
  295. if a > b {
  296. return a
  297. }
  298. return b
  299. }
  300. func (s *Source) Stop() error {
  301. s.mu.Lock()
  302. params := s.params
  303. s.params = nil
  304. s.mu.Unlock()
  305. if params != nil {
  306. _ = cErr(C.sdrplay_api_Uninit(s.dev.dev))
  307. }
  308. s.mu.Lock()
  309. defer s.mu.Unlock()
  310. if s.open {
  311. if s.devSelected {
  312. _ = cErr(C.sdrplay_api_ReleaseDevice(&s.dev))
  313. s.devSelected = false
  314. }
  315. _ = cErr(C.sdrplay_api_Close())
  316. s.open = false
  317. }
  318. if s.handle != 0 {
  319. s.handle.Delete()
  320. s.handle = 0
  321. }
  322. return nil
  323. }
  324. func (s *Source) ReadIQ(n int) ([]complex64, error) {
  325. deadline := time.Now().Add(5 * time.Second)
  326. s.mu.Lock()
  327. defer s.mu.Unlock()
  328. for s.size < n {
  329. remaining := time.Until(deadline)
  330. if remaining <= 0 {
  331. return nil, errors.New("timeout waiting for IQ samples")
  332. }
  333. if s.cond != nil {
  334. timer := time.AfterFunc(remaining, func() {
  335. s.cond.Broadcast()
  336. })
  337. s.cond.Wait()
  338. timer.Stop()
  339. } else {
  340. s.mu.Unlock()
  341. time.Sleep(20 * time.Millisecond)
  342. s.mu.Lock()
  343. }
  344. }
  345. out := make([]complex64, n)
  346. for i := 0; i < n; i++ {
  347. out[i] = s.buf[(s.head+i)%s.capSamples]
  348. }
  349. s.head = (s.head + n) % s.capSamples
  350. s.size -= n
  351. return out, nil
  352. }
  353. //export goStreamCallback
  354. func goStreamCallback(xi *C.short, xq *C.short, numSamples C.uint, reset C.uint, ctx unsafe.Pointer) {
  355. h := cgo.Handle(uintptr(ctx))
  356. src, ok := h.Value().(*Source)
  357. if !ok || src == nil {
  358. return
  359. }
  360. if reset != 0 {
  361. src.mu.Lock()
  362. src.head = 0
  363. src.size = 0
  364. src.resets++
  365. src.mu.Unlock()
  366. if src.cond != nil {
  367. src.cond.Broadcast()
  368. }
  369. }
  370. n := int(numSamples)
  371. if n <= 0 {
  372. return
  373. }
  374. iq := make([]complex64, n)
  375. xiSlice := unsafe.Slice((*int16)(unsafe.Pointer(xi)), n)
  376. xqSlice := unsafe.Slice((*int16)(unsafe.Pointer(xq)), n)
  377. const scale = 1.0 / 32768.0
  378. for i := 0; i < n; i++ {
  379. re := float32(float64(xiSlice[i]) * scale)
  380. im := float32(float64(xqSlice[i]) * scale)
  381. iq[i] = complex(re, im)
  382. }
  383. src.mu.Lock()
  384. src.lastSample = time.Now()
  385. src.appendRing(iq)
  386. src.mu.Unlock()
  387. }
  388. func cErr(err C.sdrplay_api_ErrT) error {
  389. if err == C.sdrplay_api_Success {
  390. return nil
  391. }
  392. return errors.New(C.GoString(C.sdrplay_api_GetErrorString(err)))
  393. }