Parallel Prony's method with multivariate matrix pencil approach and its numerical aspect

by   Nela Bosner, et al.

Prony's method is a standard tool exploited for solving many imaging and data analysis problems that result in parameter identification in sparse exponential sums f(k)=∑_j=1^Tc_je^-2π i⟨ t_j,k⟩, k∈ℤ^d, where the parameters are pairwise different { t_j}_j=1^M⊂ [0,1)^d, and { c_j}_j=1^M⊂ℂ∖{ 0} are nonzero. The focus of our investigation is on a Prony's method variant based on a multivariate matrix pencil approach. The method constructs matrices S_1, …, S_d from the sampling values, and their simultaneous diagonalization yields the parameters { t_j}_j=1^M. The parameters { c_j}_j=1^M are computed as the solution of an linear least squares problem, where the matrix of the problem is determined by { t_j}_j=1^M. Since the method involves independent generation and manipulation of certain number of matrices, there is intrinsic capacity for parallelization of the whole computation process on several levels. Hence, we propose parallel version of the Prony's method in order to increase its efficiency. The tasks concerning generation of matrices is divided among GPU's block of threads and CPU, where heavier load is put on the GPU. From the algorithmic point of view, the CPU is dedicated to the more complex tasks. With careful choice of algorithms solving the subtasks, the load between CPU and GPU is balanced. Besides the parallelization techniques, we are also concerned with some numerical issues, and we provide detailed numerical analysis of the method in case of noisy input data. Finally, we performed a set of numerical tests which confirm superior efficiency of the parallel algorithm and consistency of the forward error with the results of numerical analysis.


page 1

page 2

page 3

page 4


Parallel solver for shifted systems in a hybrid CPU-GPU framework

This paper proposes a combination of a hybrid CPU--GPU and a pure GPU so...

Optimized Multivariate Polynomial Determinant on GPU

We present an optimized algorithm calculating determinant for multivaria...

The electromagnetic waves propagation in unmagnetized plasma media using parallelized finite-difference time-domain method

The finite-difference time-domain (FDTD) method has been commonly utiliz...

A Variant RSA Acceleration with Parallelization

The standard RSA relies on multiple big-number modular exponentiation op...

Efficient parallelization strategy for real-time FE simulations

This paper introduces an efficient and generic framework for finite-elem...

Toeplitz matrices in the Boundary Control method

Solving inverse problems by dynamical variant of the BC-method is basica...

On the efficient parallel computing of long term reliable trajectories for the Lorenz system

In this work we propose an efficient parallelization of multiple-precisi...

Please sign up or login with your details

Forgot password? Click here to reset