File detector2d_core.h¶

#ifndef DP_DETECTOR2D_CORE_H
#define DP_DETECTOR2D_CORE_H

#include "doppler/corr2d/corr2d_core.h"
#include "doppler/dp_state.h"
#include "doppler/f32_buffer/f32_buffer_core.h"

#ifdef __cplusplus
extern "C" {
#endif

/* ── Noise aggregation mode (shared definition with detector_core.h) ────── */

#ifndef DET_NOISE_MODE_T_DEFINED
#define DET_NOISE_MODE_T_DEFINED
typedef enum
{
  DET_NOISE_MEAN = 0,
  DET_NOISE_MEDIAN = 1,
  DET_NOISE_MIN = 2,
  DET_NOISE_MAX = 3,
} det_noise_mode_t;
#endif /* DET_NOISE_MODE_T_DEFINED */

/* ── Per-detection result ───────────────────────────────────────────────── */

typedef struct
{
  size_t row;       
  size_t col;       
  float peak_mag;   
  float noise_est;  
  float test_stat;  
} det_result2d_t;

#ifndef DET_PEAK_T_DEFINED
#define DET_PEAK_T_DEFINED
typedef struct
{
  size_t row;   
  size_t col;   
  float  value; 
} det_peak_t;
#endif

/* ── Detector2D state ───────────────────────────────────────────────────── */

typedef struct
{
  dp_corr2d_state_t *corr;     
  dp_f32_t *ring;             
  dp_f32_framer_t framer;     
  float _Complex *out_buf;   
  float *mag_buf;           
  float *noise_scratch;     
  size_t ny;                
  size_t nx;                
  size_t n;                 
  size_t ring_cap;          
  size_t noise_lo;          
  size_t noise_hi;          
  det_noise_mode_t noise_mode;
  float threshold;          
  /* Last dump results — updated on every dump regardless of threshold. */
  size_t peak_row;
  size_t peak_col;
  float peak_mag;
  float noise_est;
  float test_stat;
  int _last_corr_valid;     
  size_t consumed;          
} dp_detector2d_state_t;

/* ── Lifecycle ──────────────────────────────────────────────────────────── */

dp_detector2d_state_t *dp_detector2d_create (const float _Complex *ref, size_t ny,
                                       size_t nx, size_t dwell,
                                       size_t noise_lo, size_t noise_hi,
                                       det_noise_mode_t noise_mode,
                                       float threshold, int nthreads);

void dp_detector2d_destroy (dp_detector2d_state_t *state);

void dp_detector2d_reset (dp_detector2d_state_t *state);

int dp_detector2d_set_ref (dp_detector2d_state_t *state, const float _Complex *ref);

void dp_detector2d_set_threshold (dp_detector2d_state_t *state, float threshold);

/* ── Stream push ────────────────────────────────────────────────────────── */

size_t dp_detector2d_push (dp_detector2d_state_t *state, const float _Complex *in,
                        size_t n_in, det_result2d_t *result,
                        size_t max_results);

size_t dp_detector2d_consumed (const dp_detector2d_state_t *state);

/* ── Serializable state (standard bytes interface; see dp_state.h) ──────────
 * corr2d child + the framer's carry (its own child blob, fewer than ny*nx
 * samples) + the last-dump result fields; scratch is config (rebuilt by
 * create).  Version 2 is: the corr2d child, the framer child, peak_row and
 * peak_col (u64), then peak_mag, noise_est and test_stat (f32).  Version 1
 * carried the ring's raw contents zero-padded to ring_cap and a last_corr
 * flag; out_buf was never in the blob, so last_corr is None after
 * set_state. */
#define DETECTOR2D_STATE_MAGIC DP_FOURCC ('D','E','T','2')
#define DETECTOR2D_STATE_VERSION 2u
size_t dp_detector2d_state_bytes (const dp_detector2d_state_t *state);
void dp_detector2d_get_state (const dp_detector2d_state_t *state, void *blob);
int dp_detector2d_set_state (dp_detector2d_state_t *state, const void *blob);

#ifdef __cplusplus
}
#endif

#endif /* DETECTOR2D_CORE_H */