Next Article in Journal
Joint Optimal Power Allocation and Relay Selection Scheme in Energy Harvesting Two-Way Relaying Network
Previous Article in Journal
Tooth-Marked Tongue Recognition Using Gradient-Weighted Class Activation Maps
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Efficient Tensor Sensing for RF Tomographic Imaging on GPUs

1
Department of Computer Engineering and Science, Shanghai University, Shanghai 200444, China
2
Shanghai Institute for Advanced Communication and Data Science, Shanghai 200444, China
*
Author to whom correspondence should be addressed.
Future Internet 2019, 11(2), 46; https://doi.org/10.3390/fi11020046
Submission received: 4 January 2019 / Revised: 31 January 2019 / Accepted: 11 February 2019 / Published: 15 February 2019
(This article belongs to the Section Big Data and Augmented Intelligence)

Abstract

:
Radio-frequency (RF) tomographic imaging is a promising technique for inferring multi-dimensional physical space by processing RF signals traversed across a region of interest. Tensor-based approaches for tomographic imaging are superior at detecting the objects within higher dimensional spaces. The recently-proposed tensor sensing approach based on the transform tensor model achieves a lower error rate and faster speed than the previous tensor-based compress sensing approach. However, the running time of the tensor sensing approach increases exponentially with the dimension of tensors, thus not being very practical for big tensors. In this paper, we address this problem by exploiting massively-parallel GPUs. We design, implement, and optimize the tensor sensing approach on an NVIDIA Tesla GPU and evaluate the performance in terms of the running time and recovery error rate. Experimental results show that our GPU tensor sensing is as accurate as the CPU counterpart with an average of 44.79 × and up to 84.70 × speedups for varying-sized synthetic tensor data. For IKEA Model 3D model data of a smaller size, our GPU algorithm achieved 15.374× speedup over the CPU tensor sensing. We further encapsulate the GPU algorithm into an open-source library, called cuTensorSensing (CUDA Tensor Sensing), which can be used for efficient RF tomographic imaging.

1. Introduction

Radio frequency (RF) tomographic imaging, also called wireless tomography, is a technique for detecting objects within a specific region by analyzing radio frequency signals transmitted between wireless nodes, as shown in Figure 1. This technique uses radio reflections that bounce off the object body and does not require the object to carry any wireless device. This advantage has made it ideal for security checking, rescue operations, space applications, and smart buildings [1].
In general, wireless signal propagating on a link loses power due to shadowing, distance, and multipath fading. The task of RF tomographic imaging is to estimate a spatial distribution of the shadowing loss in a detecting region from the measured power of wireless signals. Classic radio tomography and vector-based methods have been proposed. The work in [2] presented a linear model for using received signal strength (RSS) measurements to obtain images of moving objects. The work in [3] addressed substantial issues for the practical implementation. A tensor-based compressed sensing method [1] was proposed for RF tomographic imaging in a three and higher dimensionally-structured region, which estimates the distribution by utilizing its low-rank property. However, the method requires computing tensor singular value decomposition (t-SVD) in each algorithmic iteration, which leads to high computational complexity. To address this problem, Deng et al. [4] proposed to use the transform-based tensor model to formulate the RF tomographic imaging as a tensor sensing problem, then used a fast iterative algorithm Alt-Min to solve the tensor sensing problem. Their method fully utilizes the geometric structure of the three-dimensional loss field tensor. Compared to the tensor-based compressed sensing method, Deng’s method achieved a lower error rate and faster computation speed.
However, the iterative algorithm Alt-Min proposed by Deng et al. [4] iteratively estimates a pair of tensors, which is a computationally-intensive process, thus not being very practical for RF tomographic imaging for big objects or regions, or real-time applications. In this paper, we address this problem by exploiting graphics processing units (GPUs). Because of massive hardware parallelism and high memory bandwidth, GPUs have been widely used in diverse applications including machine learning [5,6,7], graph processing [8,9,10], big data analytics [11,12], image processing [13], and fluid dynamics [14]. In order to reap the power of GPUs, the algorithmic steps need to be mapped delicately onto the architecture of GPUs, especially the thread and memory hierarchy. In this work, we design, implement, and optimize the transform model-based tensor sensing method of Deng et al. [4] on an NVIDIA Tesla V100 GPU with CUDA (Compute Unified Device Architecture) and evaluate it in terms of the running time and recovery error. Experiment results show that the GPU algorithm achieves similar relative error as the CPU counterpart. Moreover, the GPU tensor sensing outperforms the CPU tensor sensing on all tensor sizes with 45.39 × speedup on average and up to 84.70× speedup on bigger tensor sizes such as 120 × 120 × 6 . We encapsulate this GPU tensor sensing algorithm into an open-source library called “cuTensorSensing” (CUDA Tensor Sensing), which is available at [15].
Our contributions are summarized as follows. First, we analyze the steps of the transform model-based tensor sensing using the Alt-Min algorithm and discuss how to map them onto the GPU architecture. Second, we design, implement, and optimize the tensor sensing on an NVIDIA Tesla V100 GPU. Third, we evaluate and compare the GPU tensor sensing and the CPU tensor sensing in terms of running time and relative error with both synthetic data and IKEA Model data. Fourth, we encapsulate the GPU tensor sensing implementation into an open-source library such that it can be used in diverse applications.
The remainder of the paper is organized as follows. In Section 2, we discuss the related works. Section 3 presents the notations and briefly summarizes the RF tomographic imaging task as a tensor sensing problem. Section 4 describes the design, implementation and optimizations of the tensor sensing on the GPU. In Section 5, we describe the experiment methodology. Section 6 evaluates the GPU tensor sensing with both synthetic data and IKEA Model data. The conclusions are given in Section 7.

2. Related Works

Existing works on RF tomographic imaging can be classified into vector-based and tensor-based approaches. The vector-based approaches [16,17] aim at estimating a spatial distribution of the shadowing loss in 2D regions of interest. They are not able to infer three-dimensional regions due to the fact that spatial structures of the signal data are ignored. Therefore, researchers have proposed the tensor-based approaches [1,4] to infer three-dimensional spaces. The tensor-based compressed sensing [1] uses the tensor nuclear norm (TNN) [18] to extend RF tomographic problems to the three-dimensional case. This approach has a high computational complexity and error rate. Deng et al. [4] exploited the transform-based tensor model [19] to explore three-dimensional spatial structures for higher accuracy and speed. Our work aims at accelerate Deng’s work on GPUs to make it practical for real-time scenarios and bigger tensors.
GPUs have massive parallelism and high memory bandwidth, and many existing research works [20,21,22,23,24] have demonstrated the benefit of utilizing GPUs to accelerate general purpose computing. Due to the high dimensions of tensors, tensor computations are often computation-intensive and time consuming. Recently, GPUs have been increasingly adopted to accelerate diverse tensor computations. Some works focused on accelerating specific tensor operations including tensor contraction [25,26], factorization [27], transpose [28,29], and tensor-matrix multiplication [30]. These works propose parallel tensor algorithms specifically optimized for the GPU architectures. GPUTENSOR [31] is a parallel tensor factorization algorithm that splits a tensor into smaller blocks and exploits the inherent parallelism and high-memory bandwidth of GPUs. To handle dynamic tensor data, GPUTENSOR updates its previously-factorized components instead of recomputing them from the raw data. Sparse tensor-times-dense matrix multiplication is a critical bottleneck in data analysis and mining applications, and [32] proposed an efficient primitive on CPU and GPU platforms. Different from these works, our work considers the tensor sensing based on the transform tensor model for RF tomographic imaging.

3. Description of the Tensor Sensing Problem

Deng et al. [4] proposed to use the transform-based tensor model to formulate the RF tomographic imaging as a tensor sensing problem. Then, they utilized a fast iterative algorithm Alt-Min to solve the tensor sensing problem. Here, we briefly summarize their approach including the Alt-Min algorithm in order to map it onto the GPU architecture.

3.1. Alt-Min Algorithm

The goal of tensor sensing is to recover the loss field tensor X from the linear map matrix A and the measurement vector y , which is formulated as follows:
X ^ = arg min X y A X F 2 , s . t . rank ( X ) r .
This method uses the Alt-Min algorithm (Algorithm 1) to estimate two low rank matrices iteratively whose matrix product is the squeezed matrix of the object tensor X .
Algorithm 1 Alt-Min algorithm of the tensor sensing.
Input: linear map matrix A , measurement vector y , iteration number L.
Output: squeezed X : X
1:
Initialize U 0 randomly;
2:
for = 1 L do
3:
     V least squares minimization ( A , U 1 , y )
4:
     U least squares minimization ( A , V , y )
5:
end for
Output: Pair of tensors ( U L , V L ).

3.2. Implementation of the Tensor Sensing on CPU

Algorithm 2 shows the implementation of tensor sensing on CPU.
Algorithm 2 Implementation of the tensor sensing on CPU.
Input: Randomly-initialized matrix U , measurement vector y , linear map matrix A
Output: X
1:
for i = 1 I t e r N u m do
2:
    use U to form a block diagonal matrix U b
3:
     W A U b
4:
    vec( V ) ← perform least square minimization on W and y
5:
     V transform vec( V ) back to matrix form
6:
     V t transpose V
7:
    use V to form a block diagonal matrix V t b
8:
     A t transpose A
9:
     W = A t V t b
10:
    vec( U t ) ← perform least squares minimization on W and y
11:
     U t transform vec( U t ) back to matrix form
12:
     U transpose U t
13:
end for
14:
return X = U V
First, U of the previous iteration (in the first iteration, U is initialized randomly) is used to form a block diagonal matrix U b (Algorithm 2, Line 2):
U b = U U U
Next, the least squares minimization problem is solved to estimate V of the next iteration (Algorithm 2, Line 3 and Line 4). The least minimization problem is formulated as follows:
vec ( V ) = arg min V y A U b V F 2
As the estimated V is in vector form vec( V ), vec( V ) is transformed back to matrix form V (Algorithm 2, Line 5). To estimate U , the least squares minimization problem in Equation (3) should be transposed:
vec ( U t ) = arg min U t y A t V t U t F 2
The process of estimating U is similar to estimating V , which is implemented with the corresponding transposed matrix (Algorithm 2 Lines 6–12). Next, the above process is repeated until the number of iterations reaches the set value. Lastly, the two estimated matrices are multiplied to get the final result matrix (Algorithm 2, Line 14).

4. The Implementation and Optimization of Efficient GPU Tensor Sensing

To achieve high performance on GPUs, we need to consider the data representation, the mappings from computations to GPU threads, and memory accesses. We first design a basic GPU tensor sensing implementation, then optimize the implementation to further improve performance.

4.1. Design and Implementation of the GPU Tensor Sensing

4.1.1. Data Structure

In Algorithm 2, after least squares minimization (Lines 4 and 10 in Algorithm 2), we get vectorized matrices. For a matrix A R m × n , the corresponding vectorized A is A v R m n × 1 . The vectorized matrices are converted back to the original matrices (Lines 5 and 11 in Algorithm 2). In some scientific computing programming languages, such as MATLAB, this conversion must be done with the appropriate conversion function. We adopt the column-first storage format to store matrices and vectors, which not only ensures read and write continuity, but also avoids explicit vector-to-matrix conversions since vector A v and matrix A in memory are the same in this format, as shown in Figure 2.

4.1.2. Multiplication of Block Diagonal Matrices

Using the operational properties of the block matrix, we get:
[ A 1 , A 2 , , A N 2 ] U c U c U c = [ A 1 U c , A 2 U c , , A N 2 U c ]
This shows that multiplication of block diagonal matrices can be transformed into a batch of small matrix multiplications. As we use column-first format to store A , the batch of A i is stored in constant stride. Let p indicate the location of the first element of A 0 , then the location of the first element of A i is p + i × N 1 N 3 . We utilize the gemmStridedBatchd() routine in the NVIDIA cuBLAS Library to compute a batch of matrix multiplications simultaneously to achieve better performance. Figure 3 shows how this process performs on GPU and Table 1 shows the parameters setting of this routine.

4.1.3. Eliminating Explicit Transpose Operations

After each least squares method, the transpose of the target matrix is obtained. The transpose operation of the matrix needs to be performed (Lines 6 and 12 in Algorithm 2). However, the transpose operation takes much computing time and resource. As the operation after transpose of the matrix is the multiplication of diagonal matrices, we eliminate the explicit matrix transpose by enabling the transpose option in the gemmStridedBatched() routine. In this way, the gemmStridedBatched() routine will perform matrix transpose implicitly and efficiently before the matrix multiplications.

4.1.4. Least Squares Minimization

As shown in Algorithm 2, least squares minimization (LSM) is the major step of the tensor sensing approach, which is the most time-consuming part of the entire approach. There are many approaches to perform least squares minimization. QR factorization is one of the most efficient approaches, which is well supported by CUDA.

4.2. Optimizations of the GPU Tensor Sensing

During the computation flow of the GPU tensor sensing, frequent data transfer between the CPU and GPU will significantly degrade system performance. Therefore, we design a data reusing strategy to reduce data transfer overhead and resource consumption.
In the entire tensor sensing flow, we invoke data transfer only twice at the beginning and at the end. At the beginning, the input data are transferred from the CPU to the GPU. At the end, the final result matrix is sent back to the CPU. We optimize the computations in the tensor sensing flow such that they all perform in-place calculations. For instance, in the QR decomposition to solve the least squares problem, the input vector y is overrode by the result vectors (vec( U ) and vec( V )). Therefore, we need to reassign the vector y at the beginning of each least squares minimization iteration. However, it is an expensive operation to load the original data of vector y from the CPU memory every time. Instead, we pre-allocate a memory space named d y L on the GPU that stores the original data of the vector y . Every time we need to reassign the d y , we use the GPU device-to-device transfer routine cudaMemcpyDeviceToDevice() to copy the data from d y L to d y . Since the device-to-device bandwidth of 380 GB/s is much higher than the bandwidth between the CPU and GPU of 12 GB/s, this strategy significantly reduces data transfer overhead. Algorithm 3 and Figure 4 describe the computation flow and data organization of the tensor sensing on the GPU.
Algorithm 3 Computation flow of the tensor sensing on the GPU.
Input: Data on CPU memory: randomly-initialized U 0 , measurement vector y R M , matrix A R M × N 1 N 2 N 3 converted from M sensing tensors A m
Output: X R N 1 N 3 × N 2
1:
apply for memory on GPU device: d y , d A , d W , d y L ( d A ( A ) means that the content in d A is A , the same below)
2:
data transfer: A , y , U 0 c u d a M e m c p y H o s t T o D e v i c e d A ( A ) , d y ( U 0 ) , d y L ( y )
3:
for i = 0 I t e r N u m do
4:
     d W ( uninitialized ) , d A ( A ) , d y ( U i ) c u b l a s < t > g e m m S t r i d e d B a t c h d d W ( A diag ( U i ) ) , d A ( A ) , d y ( U i )
5:
     d y ( U i ) , d y L ( y ) c u d a M e m c p y D e v i c e T o D e v i c e d y ( y ) , d y L ( y )
6:
     d W ( A diag ( U i ) ) , d y ( y ) c u s o l v e r < t > q r d y ( V ) , d W ( uninitialized )
7:
     d W ( uninitialized , d A ( A ) , d y ( V i ) ) e c u b l a s < t > g e m m S t r i d e d B a t c h d d W ( A diag ( V i ) ) , d A ( A ) , d y ( V i )
8:
     d y ( V i ) , d y L ( y ) c u d a M e m c p y D e v i c e T o D e v i c e d y ( y ) , d y L ( y )
9:
     d W ( A diag ( V i ) ) , d y ( y ) c u s o l v e r < t > q r d y ( U i ) , d W ( uninitialized )
10:
end for
11:
return X = U V

5. Experiment Methodology

In this section, we describe the experimental methodology including the hardware and software platform, testing data, testing process, and the comparison metrics.

5.1. Hardware and Software Platform

We use an NVIDIA Tesla V100 GPU to evaluate the performance of GPU tensor sensing. The GPU incorporates 5120 CUDA cores @ 1.53 GHz and 32 GB DDR5 memory. It is installed on a server with 128 GB DDR memory and an Intel i7-7820x CPU with 8 cores @ 3.60 GHz supporting 16 hardware threads with hyper-threading technology. The server is running Ubuntu 18.04 LTS with Kernel Version 4.15.0. The CPU and GPU tensor sensing is running with MATLAB R2017b and CUDA 10.0, respectively.

5.2. Testing Data

In the experiment, we used both IKEA model and synthetic data in the evaluation. For IKEA model data, we used the IKEA 3D data [1] that generated a ground truth tensors of size 60 × 60 × 15 with rank 6. Each 3D model was placed in the middle of the “tensor” and occupied a part of the space. In this task, we mainly focused on the location and outline information, while the texture and color information were ignored. The synthetic data were generated according to the compressed sensing model [1]. The synthetic data did not correspond to a specific physical model. We used them to demonstrate the performance of our approach for different data sizes.

5.3. Testing Process

The synthetic and real data were processed by the CPU tensor sensing implementation in MATLAB and GPU tensor sensing implementation in CUDA, respectively. We evaluated and compared two versions of GPU tensor sensing: the unoptimized one and the optimized one. The unoptimized GPU tensor sensing adopted none of the optimization techniques in Section 4.2. We repeated each experiment five times and report the average results.

5.4. Comparison Metrics

We adopted two metrics for comparison: running time and relative error rate.
  • Running time: Varying the tensor size and fixing other parameters, we measured the execution time of the CPU tensor sensing, unoptimized GPU tensor sensing, and optimized GPU tensor sensing. Finally, we calculated speedups as the running time of the CPU tensor sensing divided by the running time of GPU tensor sensing.
  • Error rate: We adopted the metric relative square error, defined as R S E = X ^ X F / X F .

6. Results and Analysis

6.1. Running Time of Synthetic Data

Figure 5a shows that the running time of the CPU tensor sensing and two GPU tensor sensing implementations (unoptimized and optimized ones) for X of size n × n × 6 of rank one, where n varies from 40–120 at a step of 20. The sampling rate was set to 20%, and both CPU and GPU tensor sensing performed five iterations for completion. The detailed time value is listed in Table 2. While M = N 1 × N 2 × N 3 × sampling rate and A R M × N 1 N 2 N 3 , the scale of the main operation matrix A increased at a rate of four times as the increase rate of n.
We can see that the running time of the CPU tensor sensing was polynomial with the increase of n, while the running time of the GPU tensor sensing was linearly growing. The unoptimized and optimized GPU tensor sensing achieved an average of 4.24 × and 44.79 × and up to 9.77 × and 84.70 × speedups, respectively. This illustrates the effectiveness of the optimization methods proposed in Section 4.2. When the tensor size was small, the data transfer occupied a major portion of the entire execution time of the unoptimized GPU tensor sensing. As a result, its performance was even lower than the CPU tensor sensing.

6.2. Error Rate and Running Time of IKEA Model Data

This experiment evaluated the error rate of the CPU and GPU tensor sensing under different iterations for X of size 60 × 60 × 15 with rank six. The sampling rate was set to 50%. Figure 6b shows the recovery results with 10 iterations. The running time of the CPU tensor sensing at five iterations was 14.91 s on average, while the running time of the GPU tensor sensing was 0.97 s on average; thus, the speedup is 15.37 × . As shown in Figure 7, under increased iterations from 1–30, the RSEs dropped significantly. More importantly, the CPU tensor sensing and GPU implementation achieved almost the same RSEs at all iterations (the two curves in Figure 7 overlap with each other), which means that they achieved similar error rate performance in tensor sensing.

7. Conclusions

In this work, we present an open-source library named cuTensorSensing for efficient RF tomographic imaging on GPUs. The experiment evaluations show that the proposed GPU tensor sensing works effectively and accurately. Compared with the counterpart on CPU, the GPU tensor sensing achieved a similar error rate, but much faster speed. For synthetic data, the GPU tensor sensing achieved an average of 44.79 × and up to 84.70 × speedups versus the CPU tensor sensing for bigger tensors. For IKEA Model 3D objects’ data of a smaller tensor size, the GPU tensor sensing achieved a 15.37 × speedup over the CPU tensor sensing. The cuTensorSensing library is useful for efficient RF tomographic imaging.

Author Contributions

Abstract, introduction, related works, conclusions, and the revising of the entire manuscript, T.Z.; algorithms, experiments, and results, D.X.

Funding

This research is partially supported by the Natural Science Foundation of Shanghai under Grant No. 17ZR1409800 and the Science and Technology Committee of Shanghai Municipality under Grant No. 16010500400.

Acknowledgments

The authors would like to thank the anonymous reviewers for their fruitful feedback and comments, which have helped them improve the quality of this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Matsuda, T.; Yokota, K.; Takemoto, K.; Hara, S.; Ono, F.; Takizawa, K.; Miura, R. Multi-dimensional wireless tomography using tensor-based compressed sensing. Wirel. Pers. Commun. 2017, 96, 3361–3384. [Google Scholar] [CrossRef]
  2. Wilson, J.; Patwari, N. Radio tomographic imaging with wireless networks. IEEE Trans. Mob. Comput. 2010, 9, 621–632. [Google Scholar] [CrossRef]
  3. Beck, B.; Ma, X.; Baxley, R. Ultrawideband Tomographic Imaging in Uncalibrated Networks. IEEE Trans. Wirel. Commun. 2016, 15, 6474–6486. [Google Scholar] [CrossRef]
  4. Deng, T.; Qian, F.; Liu, X.Y.; Zhang, M.; Walid, A. Tensor Sensing for Rf Tomographic Imaging. In Proceedings of the 2018 IEEE International Conference on Multimedia and Expo (ICME), Miami, FL, USA, 10–12 April 2018; pp. 1–6. [Google Scholar]
  5. Cui, H.; Zhang, H.; Ganger, G.R.; Gibbons, P.B.; Xing, E.P. GeePS: Scalable deep learning on distributed GPUs with a GPU-specialized parameter server. In Proceedings of the Eleventh European Conference on Computer Systems, London, UK, 18–21 April 2016; pp. 1–16. [Google Scholar]
  6. Brito, R.; Fong, S.; Song, W.; Cho, K.; Bhatt, C.; Korzun, D. Detecting Unusual Human Activities Using GPU-Enabled Neural Network and Kinect Sensors. In Internet of Things and Big Data Technologies for Next Generation Healthcare; Springer: Berlin/Heidelberg, Germany, 2017; pp. 359–388. [Google Scholar]
  7. Campos, V.; Sastre, F.; Yagües, M.; Torres, J.; Giró-i Nieto, X. Scaling a Convolutional Neural Network for Classification of Adjective Noun Pairs with TensorFlow on GPU Clusters. In Proceedings of the 17th IEEE/ACM International Symposium on Cluster, Cloud and Grid Computing, Madrid, Spain, 14–17 May 2017; pp. 677–682. [Google Scholar]
  8. Shi, X.; Luo, X.; Liang, J.; Zhao, P.; Di, S.; He, B.; Jin, H. Frog: Asynchronous graph processing on GPU with hybrid coloring model. IEEE Trans. Knowl. Data Eng. 2018, 30, 29–42. [Google Scholar] [CrossRef]
  9. Zhong, W.; Sun, J.; Chen, H.; Xiao, J.; Chen, Z.; Cheng, C.; Shi, X. Optimizing Graph Processing on GPUs. IEEE Trans. Parallel Distrib. Syst. 2017, 28, 1149–1162. [Google Scholar] [CrossRef]
  10. Pan, Y.; Wang, Y.; Wu, Y.; Yang, C.; Owens, J.D. Multi-GPU graph analytics. In Proceedings of the 2017 IEEE International Parallel and Distributed Processing Symposium (IPDPS), Orlando, FL, USA, 29 May–2 June 2017; pp. 479–490. [Google Scholar]
  11. Gutiérrez, P.D.; Lastra, M.; Benítez, J.M.; Herrera, F. SMOTE-GPU: Big Data preprocessing on commodity hardware for imbalanced classification. Prog. Artif. Intell. 2017, 6, 1–8. [Google Scholar] [CrossRef]
  12. Rathore, M.M.; Son, H.; Ahmad, A.; Paul, A.; Jeon, G. Real-time big data stream processing using GPU with spark over hadoop ecosystem. Int. J. Parallel Program. 2017, 46, 1–17. [Google Scholar] [CrossRef]
  13. Devadithya, S.; Pedross-Engel, A.; Watts, C.M.; Landy, N.I.; Driscoll, T.; Reynolds, M.S. GPU-Accelerated Enhanced Resolution 3-D SAR Imaging With Dynamic Metamaterial Antennas. IEEE Trans. Microw. Theory Tech. 2017, 65, 5096–5103. [Google Scholar] [CrossRef]
  14. Verma, K.; Szewc, K.; Wille, R. Advanced load balancing for SPH simulations on multi-GPU architectures. In Proceedings of the High Performance Extreme Computing Conference (HPEC), Waltham, MA, USA, 12–14 September 2017; pp. 1–7. [Google Scholar]
  15. Intelligent Information Processing (IIP) Lab. Available online: http://www.findai.com (accessed on 15 February 2019).
  16. Kanso, M.A.; Rabbat, M.G. Compressed RF tomography for wireless sensor networks: Centralized and decentralized approaches. In International Conference on Distributed Computing in Sensor Systems; Springer: Berlin/Heidelberg, Germany, 2009; pp. 173–186. [Google Scholar]
  17. Mostofi, Y. Compressive cooperative sensing and mapping in mobile networks. IEEE Trans. Mob. Comput. 2011, 10, 1769–1784. [Google Scholar] [CrossRef]
  18. Li, Q.; Schonfeld, D.; Friedland, S. Generalized tensor compressive sensing. In Proceedings of the 2013 IEEE International Conference on Multimedia and Expo (ICME), San Jose, CA, USA, 15–19 July 2013; pp. 1–6. [Google Scholar]
  19. Liu, X.Y.; Wang, X. Fourth-order tensors with multidimensional discrete transforms. arXiv, 2017; arXiv:1705.01576. [Google Scholar]
  20. Jing, N.; Jiang, L.; Zhang, T.; Li, C.; Fan, F.; Liang, X. Energy-efficient eDRAM-based on-chip storage architecture for GPGPUs. IEEE Trans. Comput. 2016, 65, 122–135. [Google Scholar] [CrossRef]
  21. Zhang, T.; Jing, N.; Jiang, K.; Shu, W.; Wu, M.Y.; Liang, X. Buddy SM: Sharing Pipeline Front-End for Improved Energy Efficiency in GPGPUs. ACM Trans. Archit. Code Optim. (TACO) 2015, 12, 1–23. [Google Scholar] [CrossRef]
  22. Zhang, T.; Zhang, J.; Shu, W.; Wu, M.Y.; Liang, X. Efficient graph computation on hybrid CPU and GPU systems. J. Supercomput. 2015, 71, 1563–1586. [Google Scholar] [CrossRef]
  23. Zhang, T.; Shu, W.; Wu, M.Y. CUIRRE: An open-source library for load balancing and characterizing irregular applications on GPUs. J. Parallel Distrib. Comput. 2014, 74, 2951–2966. [Google Scholar] [CrossRef]
  24. Zhang, T.; Tong, W.; Shen, W.; Peng, J.; Niu, Z. Efficient Graph Mining on Heterogeneous Platforms in the Cloud. In Cloud Computing, Security, Privacy in New Computing Environments; Springer: Berlin/Heidelberg, Germany, 2016; pp. 12–21. [Google Scholar]
  25. Nelson, T.; Rivera, A.; Balaprakash, P.; Hall, M.; Hovland, P.D.; Jessup, E.; Norris, B. Generating efficient tensor contractions for gpus. In Proceedings of the 44th International Conference on Parallel Processing (ICPP), Beijing, China, 1–4 September 2015; pp. 969–978. [Google Scholar]
  26. Shi, Y.; Niranjan, U.; Anandkumar, A.; Cecka, C. Tensor contractions with extended BLAS kernels on CPU and GPU. In Proceedings of the IEEE 23rd International Conference on High Performance Computing (HiPC), Kochi, India, 16–19 December 2016; pp. 193–202. [Google Scholar]
  27. Antikainen, J.; Havel, J.; Josth, R.; Herout, A.; Zemcik, P.; Hautakasari, M. Nonnegative tensor factorization accelerated using GPGPU. IEEE Trans. Parallel Distrib. Syst. 2011, 22, 1135–1141. [Google Scholar] [CrossRef]
  28. Lyakh, D.I. An efficient tensor transpose algorithm for multicore CPU, Intel Xeon Phi, and NVidia Tesla GPU. Comput. Phys. Commun. 2015, 189, 84–91. [Google Scholar] [CrossRef] [Green Version]
  29. Hynninen, A.P.; Lyakh, D.I. cuTT: A high-performance tensor transpose library for CUDA compatible GPUs. arXiv, 2017; arXiv:1705.01598. [Google Scholar]
  30. Rogers, D.M. Efficient primitives for standard tensor linear algebra. In Proceedings of the XSEDE16 Conference on Diversity, Big Data, and Science at Scale, ACM, Miami, FL, USA, 17–21 July 2016; p. 14. [Google Scholar]
  31. Zou, B.; Li, C.; Tan, L.; Chen, H. GPUTENSOR: Efficient tensor factorization for context-aware recommendations. Inf. Sci. 2015, 299, 159–177. [Google Scholar] [CrossRef]
  32. Li, J.; Ma, Y.; Yan, C.; Vuduc, R. Optimizing sparse tensor times matrix on multi-core and many-core architectures. In Proceedings of the IEEE Workshop on Irregular Applications: Architecture and Algorithms (IA3), Salt Lake City, UT, USA, 13–18 November 2016; pp. 26–33. [Google Scholar]
Figure 1. An illustration of the RF tomographic imaging network.
Figure 1. An illustration of the RF tomographic imaging network.
Futureinternet 11 00046 g001
Figure 2. Vectorizing a matrix A into vec( A ) in memory.
Figure 2. Vectorizing a matrix A into vec( A ) in memory.
Futureinternet 11 00046 g002
Figure 3. Multiplication of block diagonal matrices on GPU.
Figure 3. Multiplication of block diagonal matrices on GPU.
Futureinternet 11 00046 g003
Figure 4. Memory organization in the calculation process.
Figure 4. Memory organization in the calculation process.
Futureinternet 11 00046 g004
Figure 5. (a) The running time of the CPU tensor sensing and GPU tensor sensing and (b) the speedups of the unoptimized and optimized GPU tensor sensing.
Figure 5. (a) The running time of the CPU tensor sensing and GPU tensor sensing and (b) the speedups of the unoptimized and optimized GPU tensor sensing.
Futureinternet 11 00046 g005
Figure 6. (a) The 3D visualization of the IKEA models and (b) the corresponding recovery results.
Figure 6. (a) The 3D visualization of the IKEA models and (b) the corresponding recovery results.
Futureinternet 11 00046 g006
Figure 7. RSE of the CPU tensor sensing and GPU tensor sensing.
Figure 7. RSE of the CPU tensor sensing and GPU tensor sensing.
Futureinternet 11 00046 g007
Table 1. The parameters of the gemmStridedBatchd() routine in the cuBLAS library.
Table 1. The parameters of the gemmStridedBatchd() routine in the cuBLAS library.
ParametersMeaningValue
transAoperation op( A ) that is non- or transposenon-transpose
transUoperation op( U ) that is non- or transposetranspose
A pointer to the A matrix corresponding to the first instance of the batch d A
U pointer to the U matrix d y
W pointer to the W matrix d W
strideAthe address offset between A i and A i + 1 M × N 1 N 3
strideUthe address offset between U i and U i + 1 0
strideWthe address offset between W i and W i + 1 M × N 3
batchNumnumber of g e m m to perform in the batch N 2
Table 2. Running time under varying tensor sizes ( n × n × 6 ).
Table 2. Running time under varying tensor sizes ( n × n × 6 ).
n406080100120
CPU tensor sensing time (s)3.0713.6550.43118.84251.54
Unoptimized GPU tensor sensing time (s)9.5011.4012.1520.7025.74
Optimized GPU tensor sensing time (s)0.440.630.982.012.97
Speedups-unoptimized0.321.204.155.749.77
Speedups-optimized6.9821.6751.4659.1284.70

Share and Cite

MDPI and ACS Style

Xu, D.; Zhang, T. Efficient Tensor Sensing for RF Tomographic Imaging on GPUs. Future Internet 2019, 11, 46. https://doi.org/10.3390/fi11020046

AMA Style

Xu D, Zhang T. Efficient Tensor Sensing for RF Tomographic Imaging on GPUs. Future Internet. 2019; 11(2):46. https://doi.org/10.3390/fi11020046

Chicago/Turabian Style

Xu, Da, and Tao Zhang. 2019. "Efficient Tensor Sensing for RF Tomographic Imaging on GPUs" Future Internet 11, no. 2: 46. https://doi.org/10.3390/fi11020046

APA Style

Xu, D., & Zhang, T. (2019). Efficient Tensor Sensing for RF Tomographic Imaging on GPUs. Future Internet, 11(2), 46. https://doi.org/10.3390/fi11020046

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop