License: CC BY 4.0
arXiv:2606.14992v1 [cs.AR] 12 Jun 2026

KATANA: A Fast, Low-Power Mapping of Kalman Filters onto Edge NPUs for Real-Time Tracking

Bodhisatwa Kundu1, Anish Rooj1, Sumit Saha1, Abhradeep Sarkar1, Arghadip Das2, Arnab Raha3, Mrinal K. Naskar1
Abstract

State estimation is the closed-loop core of every real-time tracking system, from radar surveillance, missile guidance, and counter-unmanned aerial vehicle (UAV) defense to autonomous driving and robotics. All of these deployments run on edge platforms: defense systems mount on vehicles, drones, and interceptors far from fixed infrastructure, while civilian pipelines live on cars, drones, robots, and handheld devices, where every additional watt of compute erodes mission duration or operational range. Two hard constraints follow: each new measurement must be fused before the next control cycle (a few milliseconds, scaled by the number of tracked targets), and the total compute must fit within a strict battery and thermal power envelope. The Linear and Extended Kalman Filters (LKF, EKF) are the dominant estimators on this class of system, but today they execute almost exclusively on the CPU, which serializes multi-object tracking (MOT) updates, or on custom FPGA/ASIC accelerators that lengthen design cycles and add silicon area. Contemporary AI personal computer (AI-PC) system-on-chip (SoC) silicon, such as the Intel Core Ultra Series 1 and Series 2, already integrates a low-power, data-parallel Neural Processing Unit (NPU) alongside the CPU and GPU; we therefore ask whether the Kalman filter (KF) can be mapped onto this existing matrix engine to meet real-time and low-power budgets simultaneously. Such a mapping avoids a dedicated accelerator and keeps the CPU and GPU free for their primary workloads. We present KATANA, a novel NPU-aware optimization framework that delivers the first end-to-end mapping of the LKF and EKF onto a commercial NPU, together with the first cross-platform (CPU, GPU, NPU) characterization of these filters on shipping AI-PC silicon. KATANA applies three algebraic graph rewrites: subtract-to-add reformulation via a precomputed negative-projection matrix 𝐇neg\mathbf{H}_{\text{neg}}, static-shape tensor fusion, and block-diagonal batched parallelization, so that 100% of operations execute on the Data Processing Unit (DPU) matrix engine. On the Series 2, the optimized batched EKF reaches 223.35 frames per second (FPS) at 13.43 W active power and the LKF reaches 408.73 FPS at 14.05 W, delivering up to a 97.9% reduction in dynamic energy versus the CPU implementation.

I Introduction

Real-time state estimation is the closed-loop heartbeat of modern tracking systems. In defense, ground-based and airborne radars are mounted on vehicles, ships, and aircraft far from fixed infrastructure [1, 2], while counter-UAV and missile-guidance loops run on drones and interceptors that themselves move at high speed. In civilian deployments, the same tracking workloads run on cars (for advanced driver assistance), drones, robots, and handheld devices [3, 4], where the tracker has to live on-device for latency and connectivity reasons. All of these are edge platforms: battery-powered and either fanless or weight-constrained, so every additional watt of compute erodes mission duration or operational range, while the control loop must still close within milliseconds of each new measurement and scaled by the number of tracked targets. Across all of these settings, the LKF and its non-linear counterpart, the EKF, have been the dominant estimators for over six decades [5, 6]: optimal in the minimum-mean-square-error (MMSE) sense, recursive (hence streaming-friendly), and structured around dense linear algebra that any modern compute substrate should be able to execute.

The hardware to do so already exists. Contemporary client SoCs, such as the Intel Core Ultra Series 1 and Series 2, Apple M-series, and Qualcomm Snapdragon X, are now heterogeneous, integrating a CPU, a GPU, and a dedicated NPU optimized for low-power AI inference. Yet, despite this hardware diversity, the KF today runs almost exclusively on the CPU in commodity software stacks, or on custom FPGA/ASIC accelerators in dedicated tracking systems [7]. Three problems follow. First, MOT with NN independent filters serializes on the CPU, capping throughput and burning power. Second, designing a dedicated accelerator extends time-to-deployment and adds die area on top of an already-busy SoC. Third, the on-die NPU, a data-parallel matrix engine designed for sustained low-power operation, sits idle whenever no neural workload is feeding it. Mapping the KF to this otherwise-idle NPU yields a double dividend: it avoids a separate accelerator and keeps the CPU and GPU free for their primary general-compute and graphics workloads, so the tracker coexists with rather than monopolizes the SoC, making the end-to-end system more responsive.

This work therefore asks the natural question: can the KF be mapped onto the existing AI-PC NPU and made real-time at low active power, without custom silicon? Fig. 1 captures the overall premise: take a classical signal-processing algorithm such as KF-based tracking, target a heterogeneous AI-PC platform, and offload the recursive estimator to the on-die NPU. The challenge, however, is that NPUs are designed around dense multiply-accumulate (MAC) dataflow on the DPU; any operation that falls outside that pattern (Subtract, Reshape, Transpose, Gather) is routed to a scalar Digital Signal Processor (DSP) and forces costly DPU\leftrightarrowDSP context switches that can consume 10–30% of inference time on the small kernels typical of KF tracking.

Refer to caption
Figure 1: Overview. KATANA targets a classical signal-processing workload, Kalman-filter-based object tracking (left), on heterogeneous AI-PC SoCs (centre) and maps the recursive estimator onto the on-die NPU (right) for faster and more energy-efficient execution.

Building on the above motivation, we present KATANA, a novel NPU-aware optimization framework for mapping traditional state estimators onto AI-PC silicon. Our contributions are:

  • We present the first cross-platform (CPU, GPU, NPU) characterization of the LKF and EKF inference on two consecutive AI-PC generations (Series 1 and Series 2).

  • We introduce three NPU-aware algebraic graph rewrites: a precomputed negative-projection matrix 𝐇neg\mathbf{H}_{\text{neg}} that converts subtract-heavy innovations into DPU-native adds, static-shape tensor fusion that removes runtime Reshape/Transpose nodes, and block-diagonal batched parallelization that packs NN independent filters into one inference call. Together, these rewrites move 100% of operations onto the DPU.

  • We demonstrate that the optimized NPU pipeline meets real-time latency budgets within a sustained 13–14 W active envelope, sustains >>200 FPS multi-filter throughput on the Series 2, and reduces dynamic energy by up to 97.9% versus CPU execution.

  • We validate end-to-end on a live video stream, where the NPU-resident LKF and EKF consume <<1% of a 33 ms frame budget at 30 FPS, leaving the CPU and GPU free for detection and downstream analytics.

II Background

We briefly review the algorithmic primitive we wish to accelerate (Section II-A) and the NPU substrate on which we will run it (Section II-B). Together, these two pieces frame the optimization problem addressed in Section IV.

II-A Kalman Filtering for State Estimation

The KF is a recursive estimator for the hidden state of a noisy dynamical system [5, 6]. The LKF assumes 𝐱k=𝐅𝐱k1+𝐰k1\mathbf{x}_{k}=\mathbf{F}\mathbf{x}_{k-1}+\mathbf{w}_{k-1} and 𝐳k=𝐇𝐱k+𝐯k\mathbf{z}_{k}=\mathbf{H}\mathbf{x}_{k}+\mathbf{v}_{k} with Gaussian process and measurement noise 𝐰k,𝐯k\mathbf{w}_{k},\mathbf{v}_{k}, and alternates a prediction step (propagating the state and covariance through 𝐅\mathbf{F}) with an update that corrects via the innovation 𝐲k=𝐳k𝐇𝐱^k|k1\mathbf{y}_{k}=\mathbf{z}_{k}-\mathbf{H}\hat{\mathbf{x}}_{k|k-1} scaled by the Kalman gain 𝐊k=𝐏k|k1𝐇T(𝐇𝐏k|k1𝐇T+𝐑)1\mathbf{K}_{k}=\mathbf{P}_{k|k-1}\mathbf{H}^{T}(\mathbf{H}\mathbf{P}_{k|k-1}\mathbf{H}^{T}+\mathbf{R})^{-1}. The EKF keeps the same linear-gain structure but linearizes the dynamics and observation maps about the current estimate via the Jacobians 𝐅k,𝐇k\mathbf{F}_{k},\mathbf{H}_{k}. Each recursion therefore reduces to a chain of dense matrix multiplications plus a single inversion.

II-B NPU Microarchitecture and Execution Flow

Fig. 2 shows the Intel NPU on the Series 1 and Series 2 SoCs: a Command Interface, a managed on-chip SRAM with DMA, and a Compute Cluster of identical Compute Engines that each pair a Systolic DPU (dense MAC) with two Vector DSP units, a Post-Compute Unit, and Load/Store units [8, 9]; the same architecture and its five-stage tensor-execution pipeline are described in detail in [14]. Only one property matters for the rest of the paper: the DPU dominates throughput and energy efficiency on dense GEMM, while any op routed to the DSP serializes with the DPU pipeline and triggers a DPU\leftrightarrowDSP context switch whose fixed cost is significant for the small tensors typical of KF tracking. Maximizing DPU occupancy is therefore the central design problem for Section IV.

III Related Work

Our work sits at the intersection of two research threads: hardware acceleration of state estimators, and algorithm–NPU co-design for non-CNN workloads.

Hardware acceleration of KF tracking. FPGA implementations of multi-dimensional KFs for object tracking achieve deterministic latency but require lengthy RTL design cycles and carry higher static power than current SoC fabrics [7]. GPU implementations exploit batch parallelism in many-target settings but exceed the thermal envelopes of fanless and battery-powered edge devices, and lower-precision strategies such as stochastic computing [10] reduce per-operation energy in principle without addressing the DSP\leftrightarrowDPU context-switch cost on heterogeneous AI accelerators.

NPU co-design for non-CNN workloads. Modern NPUs increasingly host workloads beyond convolutional networks: FlexNPU offers a dataflow-flexible substrate for energy-efficient edge inference [9], and recent frameworks have mapped LLMs [11], GNNs [12], SSMs [13], and Hyena/Kolmogorov–Arnold Networks [14] onto resource-constrained NPUs. A consistent finding is that every successful NPU mapping requires algorithm-level restructuring to align the workload with the matrix-engine dataflow.

Gap. Classical signal-processing kernels (KFs in particular) have not, to our knowledge, been mapped to a commercial NPU, and no cross-platform (CPU, GPU, NPU) characterization exists for the LKF and EKF on shipping AI-PC silicon. KATANA fills this gap.

Refer to caption
Figure 2: NPU microarchitecture: Command Interface, on-chip SRAM with managed L2 cache, and a Compute Cluster of Compute Engines, each combining a systolic DPU, two DSP units, a PCU, and Load/Store units.

IV NPU-Aware Graph Optimization for Kalman Filters

Our design goal follows directly from Section II-B: every operation in the LKF and EKF prediction-update recursion must execute on the DPU, with zero fall-back to the DSP, so that the NPU’s GEMM pipeline is the only critical path. All filters are authored in PyTorch, exported to ONNX, and compiled to each target backend through Intel OpenVINO 2024.5 [8, 4]. Three graph rewrites, applied prior to compilation, transform a naive ONNX export into a pure-DPU graph. Fig. 3 traces the resulting Netron graphs for the LKF (top panel) and EKF (bottom panel) through four stages: Baseline, Optimization-1 (subtract elimination), Optimization-2 (static-shape tensor fusion), and Batched (block-diagonal parallelization).

IV-A Optimization Pipeline Overview

In the baseline graphs of Fig. 3, the OpenVINO compiler falls back to the DSP for the subtraction in the innovation 𝐲k=𝐳k𝐇𝐱^k|k1\mathbf{y}_{k}=\mathbf{z}_{k}-\mathbf{H}\hat{\mathbf{x}}_{k|k-1} and to DMA helpers for the dynamic Reshape, Unsqueeze, and Gather nodes introduced by the default batch axis. We measure that these scalar and memory-management nodes account for roughly 15% of execution time on the LKF baseline (Subtract) and a further 12% on the EKF baseline (FP32-to-FP16 Convert and other control ops). Concretely, Optimization-1 removes the explicit Subtract module by precomputing a negative-projection matrix; Optimization-2 collapses the remaining dynamic shapes and Transposes so that no DSP-bound node survives; and the Batched stage expands the system matrices into block-diagonal form so that a single inference processes N=200N=200 independent filters in parallel. We describe each rewrite in the next three subsections.

Refer to caption
Refer to caption
Figure 3: Netron graphs of the LKF (top) and EKF (bottom) inference pipelines across four optimization stages. Baseline retains DSP-bound Subtract and dynamic-shape nodes (highlighted); Optimization-1 replaces the Subtract module with an Add over a precomputed 𝐇neg\mathbf{H}_{\text{neg}}; Optimization-2 fuses static tensor shapes and removes runtime Transposes, leaving the graph entirely on the DPU; Batched expands the system-matrix dimensions to process N=200N=200 filters per inference call.

IV-B Algebraic Reformulation: Subtract Elimination

We address the three sources of DPU fall-back identified in Section IV-A in turn, starting with the dominant Subtract op. The innovation update can be reformulated by absorbing the sign of the observation projection into a constant tensor. Defining a precomputed matrix 𝐇neg=𝐇\mathbf{H}_{\text{neg}}=-\mathbf{H} at model initialization transforms the innovation from 𝐲k=𝐳k𝐇𝐱^k|k1\mathbf{y}_{k}=\mathbf{z}_{k}-\mathbf{H}\hat{\mathbf{x}}_{k|k-1} into 𝐲k=𝐳k+𝐇neg𝐱^k|k1\mathbf{y}_{k}=\mathbf{z}_{k}+\mathbf{H}_{\text{neg}}\hat{\mathbf{x}}_{k|k-1}, which is a single GEMM followed by an element-wise Add, both DPU-native. Crucially, 𝐇neg\mathbf{H}_{\text{neg}} is folded into the ONNX graph as a constant and carries no runtime cost. The same rewrite applies to every other subtraction in the recursion (state-error and covariance update), so the entire critical path of the LKF and EKF reduces to GEMMs and Adds, as visible in the Optimization-1 column of Fig. 3.

IV-C Static Tensor Fusion for Pure-DPU Execution

Two further rewrites remove the residual DSP and DMA operations visible in the third column of Fig. 3. First, we lower the dynamic batch axis [B,][B,\cdot] that the ONNX exporter inserts by default to static 1-D and 2-D shapes [dim][\text{dim}]; this eliminates the runtime Unsqueeze, Squeeze, and Reshape nodes that OpenVINO had been dispatching to the DSP for shape bookkeeping. Second, we precompute every transposed matrix that appears in the recursion (𝐅T\mathbf{F}^{T}, 𝐇T\mathbf{H}^{T}, and 𝐇negT\mathbf{H}_{\text{neg}}^{T}) and fold them into the graph as constants, so no runtime Transpose op survives. After this stage, the LKF and EKF inference graphs consist exclusively of DPU-native MatMul, Add, and a small number of Concat primitives, as confirmed by the Perfetto traces of Fig. 4.

IV-D Block-Diagonal Batched Parallelization

While Sections IV-B and IV-C eliminate DSP fall-back for a single filter, MOT workloads require many filters to run concurrently. For NN independent filters, we therefore expand each per-filter matrix into a block-diagonal (Nn)×(Nn)(Nn)\times(Nn) system matrix, where nn is the per-filter state dimension (n=6n=6 for our LKF and n=8n=8 for our EKF). The measurement, noise, and gain matrices are expanded identically. As a result, this restructuring packs NN uncoupled filters into a single NPU inference call: the DPU’s MAC pipeline is saturated because the GEMM operands are now wide enough to amortize per-call dispatch overhead, while the block-diagonal sparsity guarantees that filters do not cross-couple numerically. All batched results in Table I use N=200N=200.

V Experimental Setup

We measure each filter variant of Section IV on representative AI-PC hardware, profiling latency, power, and dynamic energy across all three compute targets in the SoC.

Hardware platforms. We evaluate KATANA on two AI-PC class platforms: the Intel Core Ultra Series 1 and Series 2 reference systems. The CPU, GPU, and NPU are measured on the same chassis under the default thermal policy, so that platform variation is the only difference across compute targets.

Toolchain. All filters are authored in PyTorch, exported to ONNX, and compiled to each target backend through Intel OpenVINO 2024.5 [8]. We use FP16 precision uniformly across CPU, GPU, and NPU to isolate architectural effects from numerical-format effects.

Profiling. Per-operator latency and the DPU/DSP/DMA breakdown shown in Fig. 4 are captured with the OpenVINO Perfetto trace export. Whole-SoC active power is sampled with the Intel Power Gadget at 100 Hz over the measurement window. Each latency value in Table I is the mean of 1000 inference iterations after a 100-iteration warm-up.

Workloads. The LKF uses a state of dimension n=6n=6 (3-D position and velocity); the EKF uses n=8n=8 (constant-turn-rate with acceleration). Single-instance configurations process one filter per inference call; batched configurations process N=200N=200 filters in parallel through the block-diagonal expansion of Section IV-D.

End-to-end tracking pipeline. To validate KATANA on a realistic workload, we deploy it inside a live tracking pipeline derived from the OpenVINO tracking notebooks [4]. A Haar-cascade detector runs on the CPU and supplies bounding-box centroids as measurements; the NPU-resident LKF and EKF maintain independent state estimates in parallel and feed predicted positions back to the renderer overlay on the next frame. The resulting tracks over a 23-s sequence are shown in Fig. 5.

TABLE I: Latency, throughput, power, and dynamic energy of LKF and EKF variants across CPU (C), GPU (G), and NPU (N) on Intel Core Ultra Series 1 and Series 2. Bold marks the best NPU-side cell per (workload, metric).
Series 1 Series 2
Lat (ms) Thr (FPS) Pwr (W) Eng (mJ) Lat (ms) Thr (FPS) Pwr (W) Eng (mJ)
Config C G N C G N C G N C G N C G N C G N C G N C G N
LKF 0.05 0.17 0.30 15834.40 5238.68 3096.10 28.00 28.00 25.92 1.40 4.76 7.77 0.02 0.30 0.25 71753.38 96448.02 6268.29 19.98 17.11 10.42 0.40 5.13 2.61
LKF OPT 1 0.04 0.21 0.30 16648.75 4299.45 3110.31 28.00 28.00 26.66 1.12 5.88 8.00 0.02 0.31 0.19 75432.49 96555.99 8105.27 20.01 16.99 10.93 0.40 5.27 2.08
LKF OPT 2 0.05 0.17 0.29 14367.23 5345.74 3229.84 28.00 28.00 27.68 1.40 4.76 8.03 0.02 0.29 0.19 76877.72 5467.01 8099.65 19.93 11.88 9.71 0.40 3.45 1.84
LKF BATCHED 79.40 5.15 5.86 12.38 154.19 145.18 28.00 28.01 20.06 2223.20 144.23 117.55 69.05 1.77 2.96 23.75 640.10 408.73 17.41 25.61 14.05 1202.16 45.33 41.59
EKF 0.06 0.43 0.48 12819.82 2189.93 1964.70 28.00 28.00 23.00 1.68 12.04 11.04 0.03 0.57 0.41 66664.12 84068.90 3913.76 20.91 18.89 11.31 0.63 10.77 4.64
EKF OPT 1 0.06 0.45 0.47 12740.57 2087.28 1987.45 28.00 28.00 24.84 1.68 12.60 11.67 0.04 1.04 0.34 22850.51 94521.34 4650.89 22.00 15.75 7.47 0.88 16.38 2.54
EKF OPT 2 0.05 0.25 0.29 13999.15 3670.79 3278.13 28.00 28.00 25.28 1.40 7.00 7.33 0.03 0.38 0.18 58684.68 3078.13 7994.73 22.10 19.15 12.77 0.66 7.28 2.30
EKF BATCHED 122.27 15.82 10.84 8.06 55.09 78.51 28.00 28.00 19.54 3423.56 442.96 211.81 146.28 2.68 4.98 10.93 403.91 223.35 22.26 23.43 13.43 3256.46 62.79 66.88

VI Results and Discussion

We characterize KATANA along four axes: the per-stage compute breakdown that confirms DPU saturation (Fig. 4); the latency and throughput scaling with workload density (Table I); the sustained power and dynamic energy budgets; and the end-to-end tracking demonstration on live video (Fig. 5). Each axis verifies one design hypothesis from Section IV in turn.

VI-A Compute-Stage Breakdown

Fig. 4 reports the per-operator NPU compute breakdown across the four stages of Section IV. For the LKF, successive stages eliminate the baseline DPU Subtract bar and the residual non-MAC tail, leaving the Batched configuration with \approx95% of execution in DPU MatMul. The EKF baseline scatters work across several DSP and DMA tail ops that together account for roughly half of inference time; after the three rewrites, the batched EKF reaches \approx80% DPU MatMul with the rest in DPU Add and DMA Concat. The Netron-level transformations of Fig. 3 thus translate directly into measurable DPU occupancy.

Refer to caption
Figure 4: Per-operator compute-latency breakdown reported by the OpenVINO Perfetto trace for (a) LKF and (b) EKF across the four optimization stages of Section IV. Each successive stage migrates work off DSP and DMA paths and onto the DPU; the optimized batched configurations are dominated by DPU-MatMul, confirming the saturation of the matrix engine.

VI-B Latency and Throughput Scaling

Table I lists latency, throughput, power, and dynamic energy for every (filter, platform, configuration) combination; two regimes emerge.

Single-instance regime. At one filter per inference call, the CPU wins on raw latency (down to 0.02 ms on the Series 2) because the workload is too small to amortize NPU dispatch overhead; the NPU still meets any real-time budget at 0.18–0.30 ms but is not yet exercising its advantage.

Batched regime (N=200N{=}200). Once the workload is dense enough to keep the matrix engine busy, the NPU is the throughput-best target on the Series 1 SoC for the EKF (78.5 FPS vs. 8.1 CPU / 55.1 GPU). The Series 2 nearly triples NPU throughput to 223.35 FPS (EKF) and pushes the LKF to 408.73 FPS, validating the block-diagonal batching design of Section IV-D. The Series 2 GPU reaches higher peak FPS through race-to-idle execution but at substantially higher sustained power, as quantified next.

VI-C Power and Energy Efficiency

For the EKF N=200N{=}200 workload, the Series 2 NPU sustains only 13.43 W against 22–25 W for the CPU/GPU on the same chassis, and the corresponding 66.88 mJ of dynamic energy is a 97.9% reduction relative to the Series 2 CPU on the same workload (Table I). The GPU reaches competitive throughput through race-to-idle execution but at a peak power that fanless, battery-bound edge platforms cannot sustain. The NPU’s advantage is therefore sustained low power, not peak FLOPs: precisely the regime in which defense and mobile tracking platforms operate.

Refer to caption Refer to caption Refer to caption Refer to caption Refer to caption Refer to caption
Refer to caption Refer to caption Refer to caption Refer to caption Refer to caption Refer to caption
Figure 5: Frames extracted at odd timestamps from the live KATANA tracking sequence with the LKF and EKF running on the Intel Core Ultra Series 2 NPU. Both filters maintain a stable track through scale and motion changes across the 23-s sequence.

VI-D Real-Time Tracking on Live Video

Fig. 5 shows snapshots from the NPU-accelerated tracking pipeline of Section V, sampled at odd timestamps over a 23-s sequence. Both filters lock onto the target on the first frame and follow it through scale and motion changes without visible drift. On the Series 2, the optimized single-instance LKF and EKF take 0.19 ms and 0.18 ms per update, occupying <<1% of the 33 ms frame budget at 30 FPS and leaving the rest available to the Haar-cascade detector, a future re-identification head, or other concurrent SoC workloads. The end-to-end demonstration confirms that an NPU-resident KF can serve as the always-on tracking engine without monopolizing the platform.

VII Conclusion

KATANA delivers the first end-to-end mapping and cross-platform (CPU, GPU, NPU) characterization of the LKF and EKF on commercial AI-PC silicon. Three algebraic rewrites (subtract-to-add reformulation, static-shape fusion, and block-diagonal batching) move 100% of the recursion onto the DPU, so the optimized pipeline sustains >>200 FPS multi-filter throughput at sub-15 W and cuts dynamic energy by up to 97.9% versus the CPU. Classical signal-processing pipelines can therefore be re-targeted to existing AI-PC NPUs without custom accelerators.

References

  • [1] Y. Bar-Shalom et al., “Estimation with Applications to Tracking and Navigation: Theory, Algorithms, and Software,” Wiley, 2001.
  • [2] S. S. Blackman et al., “Design and Analysis of Modern Tracking Systems,” Artech House, 1999.
  • [3] Y. Zhang et al., “Application of Deep Learning Techniques in UAV Image Recognition and Tracking,” ResearchGate Preprint, 2026.
  • [4] OpenVINO Toolkit, “Person Tracking with OpenVINO,” GitHub Repo., 2024.
  • [5] R. E. Kalman, “A New Approach to Linear Filtering and Prediction Problems,” ASME J. Basic Eng., 1960.
  • [6] G. Welch et al., “An Introduction to the Kalman Filter,” Univ. North Carolina Chapel Hill, Tech. Rep. TR 95-041, 1995.
  • [7] P. Babu et al., “FPGA Implementation of Multi-Dimensional Kalman Filter for Object Tracking and Motion Detection,” Microprocess. Microsyst., 2021.
  • [8] Intel Corp., “Intel NPU Architecture and OpenVINO Optimization Guide,” Intel Tech. Doc., Rev. 2024.5, 2024.
  • [9] A. Raha et al., “FlexNPU: A Dataflow-Aware Flexible Deep Learning Accelerator for Energy-Efficient Edge Devices,” Front. High Perform. Comput., 2025.
  • [10] A. Alaghi et al., “Survey of Stochastic Computing,” ACM Trans. Embedded Comput. Syst., 2013.
  • [11] A. Raha et al., “LLM-NPU: Towards Efficient Foundation Model Inference on Low-Power Neural Processing Units,” IEEE COINS, 2025.
  • [12] A. Das et al., “GraNNite: Enabling High-Performance Execution of Graph Neural Networks on Resource-Constrained Neural Processing Units,” arXiv:2502.06921, 2025.
  • [13] A. Das et al., “XAMBA: Enabling Efficient State Space Models on Resource-Constrained Neural Processing Units,” ICLR Workshop, 2025.
  • [14] A. Das et al., “Towards Efficient Acceleration of Hyena and Kolmogorov-Arnold Networks on NPUs,” arXiv Preprint, 2025.