ldpc_on_GPU_isip

cmake_targets/CMakeLists.txt

@@ -22,7 +22,24 @@
 # Author: laurent THOMAS, Lionel GAUTHIER
 
 cmake_minimum_required (VERSION 3.0)
+##############################################
+# Base CUDA setting
+##############################################
+find_package(CUDA)
+message ("cuda include ${CUDA_INCLUDE_DIRS}")
+message ("cuda library ${CUDA_LIBRARY_DIRS}")
 
+add_definitions("-L/usr/local/cuda/lib64")
+
+SET(CUDA_NVCC_FLAGS
+  "${CUDA_NVCC_FLAGS};-arch=sm_75;")
+
+# Disable warnings for CUDA
+SET(CUDA_NVCC_FLAGS "${CUDA_NVCC_FLAGS};-lpthread;-w;-O3;--default-stream;per-thread;-I/usr/local/cuda/inc;-L/usr/local/cuda/lib -lcutil;-rdc=true;-lcudadevrt")
+SET(CUDA_VERBOSE_BUILD ON)
+SET(CUDA_HOST_COMPILER "/usr/bin/g++")
+
+SET(CUDA_SEPARABLE_COMPILATION ON)
 #########################################################
 # Base directories, compatible with legacy OAI building #
 #########################################################
@@ -174,7 +191,7 @@ set(CMAKE_C_FLAGS
 	"${CMAKE_C_FLAGS} -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_FCNTL_H=1 -DHAVE_ARPA_INET_H=1 -DHAVE_SYS_TIME_H=1 -DHAVE_SYS_SOCKET_H=1 -DHAVE_STRERROR=1 -DHAVE_SOCKET=1 -DHAVE_MEMSET=1 -DHAVE_GETTIMEOFDAY=1 -DHAVE_STDLIB_H=1 -DHAVE_MALLOC=1 -DHAVE_LIBSCTP -D${MKVER}"
 )
 set(CMAKE_CXX_FLAGS
-	"${CMAKE_CXX_FLAGS} ${C_FLAGS_PROCESSOR} -Wno-packed-bitfield-compat -fPIC -Wall -fno-strict-aliasing -rdynamic -std=c++11 -D${MKVER}"
+	"${CMAKE_CXX_FLAGS} ${C_FLAGS_PROCESSOR} -Wno-packed-bitfield-compat -fPIC -Wall -fno-strict-aliasing -rdynamic -std=c++11" #-D${MKVER}
 )
 
 add_definitions("-DASN_DISABLE_OER_SUPPORT")
@@ -2704,16 +2721,34 @@ target_link_libraries(smallblocktest
   -Wl,--start-group UTIL SIMU PHY_NR CONFIG_LIB -Wl,--end-group
   m pthread ${ATLAS_LIBRARIES} dl
   )
+###################################################
+# For CUDA library 
+###################################################
+CUDA_ADD_LIBRARY(LDPC_CU  
+  ${OPENAIR1_DIR}/PHY/CODING/nrLDPC_decoder_LYC/nrLDPC_decoder_LYC.cu
+)
 
-add_executable(ldpctest  
+cuda_add_executable(ldpctest
   ${OPENAIR1_DIR}/PHY/CODING/TESTBENCH/ldpctest.c
   ${T_SOURCE}
   ${SHLIB_LOADER_SOURCES}
-  )
-target_link_libraries(ldpctest
-  -Wl,--start-group UTIL SIMU PHY_NR CONFIG_LIB -Wl,--end-group
+)
+
+target_link_libraries(ldpctest -ldl
+  -Wl,--start-group 
+  LDPC_CU UTIL SIMU PHY_NR CONFIG_LIB 
+  -Wl,--end-group
   m pthread ${ATLAS_LIBRARIES} dl
-  )
+)
+#add_executable(ldpctest  
+#  ${OPENAIR1_DIR}/PHY/CODING/TESTBENCH/ldpctest.c
+#  ${T_SOURCE}
+#  ${SHLIB_LOADER_SOURCES}
+#  )
+#target_link_libraries(ldpctest
+#  -Wl,--start-group UTIL SIMU PHY_NR CONFIG_LIB -Wl,--end-group
+#  m pthread ${ATLAS_LIBRARIES} dl
+#  )
 
 add_executable(nr_dlschsim  
   ${OPENAIR1_DIR}/SIMULATION/NR_PHY/dlschsim.c 

cmake_targets/ldpctest_BG_1_Zc_384_rate_1-3_block_length_8448_maxit_5.txt

+SNR BLER BER UNCODED_BER ENCODER_MEAN ENCODER_STD ENCODER_MAX DECODER_TIME_MEAN DECODER_TIME_STD DECODER_TIME_MAX DECODER_ITER_MEAN DECODER_ITER_STD DECODER_ITER_MAX
+-2.000000 1.000000 0.346709 0.261600 386.906955 0.000000 386.906955 237.331605 0.000000 237.331605 4.000000 0.000000 4 
+-1.500000 1.000000 0.331321 0.245620 121.535501 -nan 121.535501 275.051825 -nan 275.051825 4.000000 0.000000 4 
+-1.000000 1.000000 0.324574 0.232521 110.700942 -nan 110.700942 273.851297 0.000004 273.851297 4.000000 0.000000 4 
+-0.500000 1.000000 0.314512 0.220723 106.537005 -nan 106.537005 264.439320 0.000004 264.439320 4.000000 0.000000 4 
+0.000000 1.000000 0.292022 0.206400 106.478256 0.000001 106.478256 265.911008 0.000005 265.911008 4.000000 0.000000 4 
+0.500000 1.000000 0.247751 0.193695 102.759587 0.000000 102.759587 263.615009 -nan 263.615009 4.000000 0.000000 4 
+1.000000 1.000000 0.219460 0.181818 103.813742 -nan 103.813742 260.650414 0.000000 260.650414 4.000000 0.000000 4 
+1.500000 1.000000 0.096354 0.165759 103.604808 0.000000 103.604808 261.180111 0.000000 261.180111 4.000000 0.000000 4 
+2.000000 1.000000 0.055990 0.155382 103.614034 0.000000 103.614034 259.823410 0.000000 259.823410 4.000000 0.000000 4 
+2.500000 1.000000 0.017401 0.140349 103.535678 0.000000 103.535678 260.160979 0.000000 260.160979 4.000000 0.000000 4 
+3.000000 1.000000 0.000118 0.124211 103.698956 -nan 103.698956 261.710344 0.000000 261.710344 4.000000 0.000000 4 
+3.500000 0.000000 0.000000 0.110835 102.942444 0.000000 102.942444 261.678181 0.000004 261.678181 4.000000 0.000000 4 

openair1/PHY/CODING/TESTBENCH/ldpctest.c

@@ -27,6 +27,7 @@
 #include "SIMULATION/TOOLS/sim.h"
 #include "PHY/CODING/nrLDPC_encoder/defs.h"
 #include "PHY/CODING/nrLDPC_decoder/nrLDPC_decoder.h"
+#include "PHY/CODING/nrLDPC_decoder_LYC/nrLDPC_decoder_LYC.h"
 #include "openair1/SIMULATION/NR_PHY/nr_unitary_defs.h"
 
 #define MAX_NUM_DLSCH_SEGMENTS 16
@@ -317,7 +318,6 @@ int test_ldpc(short No_iteration,
         for (i = 0; i < block_length+(nrows-no_punctured_columns) * Zc - removed_bit; i++)
           if (channel_input[j][i]!=channel_input_optim[j][i]) {
             printf("differ in seg %u pos %u (%u,%u)\n", j, i, channel_input[j][i], channel_input_optim[j][i]);
-            free(channel_output);
             return (-1);
           }
       //else{
@@ -392,19 +392,25 @@ int test_ldpc(short No_iteration,
       decParams.numMaxIter=No_iteration;
       decParams.outMode = nrLDPC_outMode_BIT;
       //decParams.outMode =nrLDPC_outMode_LLRINT8;
+      //decParams.outMode =nrLDPC_outMode_LLRINT8;
 
 
       for(j=0;j<n_segments;j++) {
-    	  start_meas(time_decoder);
       // decode the sequence
       // decoder supports BG2, Z=128 & 256
       //esimated_output=ldpc_decoder(channel_output_fixed, block_length, No_iteration, (double)((float)nom_rate/(float)denom_rate));
       ///nrLDPC_decoder(&decParams, channel_output_fixed, estimated_output, NULL);
-          n_iter = nrLDPC_decoder(&decParams, (int8_t*)channel_output_fixed[j], (int8_t*)estimated_output[j], p_nrLDPC_procBuf, p_decoder_profiler);
-      
-	stop_meas(time_decoder);
+#ifdef __NR_LDPC_DECODER_LYC__H__
+		 n_iter = nrLDPC_decoder_LYC(&decParams, (int8_t*)channel_output_fixed[j], (int8_t*)estimated_output[j], block_length, time_decoder);
+#else
+		 start_meas(time_decoder);
+		 n_iter = nrLDPC_decoder(&decParams, (int8_t*)channel_output_fixed[j], (int8_t*)estimated_output[j], p_nrLDPC_procBuf, p_decoder_profiler);
+		 stop_meas(time_decoder);
+#endif      
+	
       }
 
+
       //for (i=(Kb+nrows) * Zc-5;i<(Kb+nrows) * Zc;i++)
       //  printf("esimated_output[%d]=%d\n",i,esimated_output[i]);
 
@@ -413,9 +419,9 @@ int test_ldpc(short No_iteration,
       for (i=0; i<block_length>>3; i++)
       {
           //printf("block_length>>3: %d \n",block_length>>3);
-         /// printf("i: %d \n",i);
-          ///printf("estimated_output[%d]: %d \n",i,estimated_output[i]);
-          ///printf("test_input[0][%d]: %d \n",i,test_input[0][i]);
+          // printf("i: %d \n",i);
+          // printf("estimated_output[%d]: %d \n",i,estimated_output[i]);
+          // printf("test_input[0][%d]: %d \n",i,test_input[0][i]);
         if (estimated_output[j][i] != test_input[j][i])
         {
       //////printf("error pos %d (%d, %d)\n\n",i,estimated_output[i],test_input[0][i]);
@@ -460,6 +466,7 @@ int test_ldpc(short No_iteration,
 
   *errors_bit_uncoded = *errors_bit_uncoded / (double)((Kb+nrows-no_punctured_columns-2) * Zc-removed_bit);
 
+int k=0;
   for(j=0;j<MAX_NUM_DLSCH_SEGMENTS;j++) {
     free(test_input[j]);
     free(channel_input[j]);

openair1/PHY/CODING/nrLDPC_decoder_LYC/nrLDPC_decoder_LYC.h

+#ifndef __NR_LDPC_DECODER_LYC__H__
+#define __NR_LDPC_DECODER_LYC__H__
+
+#include "nrLDPC_types.h"
+#include "nrLDPC_init_mem.h"
+
+/**
+   \brief LDPC decoder
+   \param p_decParams LDPC decoder parameters
+   \param p_llr Input LLRs
+   \param p_llrOut Output vector
+   \param p_profiler LDPC profiler statistics
+*/
+
+int32_t nrLDPC_decoder_LYC(t_nrLDPC_dec_params* p_decParams, int8_t* p_llr, int8_t* p_out, int block_length, time_stats_t *time_decoder);
+
+#endif

openair1/PHY/CODING/nrLDPC_decoder_LYC/nrLDPC_init_mem.h

+/*
+ * Licensed to the OpenAirInterface (OAI) Software Alliance under one or more
+ * contributor license agreements.  See the NOTICE file distributed with
+ * this work for additional information regarding copyright ownership.
+ * The OpenAirInterface Software Alliance licenses this file to You under
+ * the OAI Public License, Version 1.1  (the "License"); you may not use this file
+ * except in compliance with the License.
+ * You may obtain a copy of the License at
+ *
+ *      http://www.openairinterface.org/?page_id=698
+ *
+ * Unless required by applicable law or agreed to in writing, software
+ * distributed under the License is distributed on an "AS IS" BASIS,
+ * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+ * See the License for the specific language governing permissions and
+ * limitations under the License.
+ *-------------------------------------------------------------------------------
+ * For more information about the OpenAirInterface (OAI) Software Alliance:
+ *      contact@openairinterface.org
+ */
+
+/*!\file nrLDPC_init_mem.h
+ * \brief Defines the function to initialize the LDPC decoder and sets correct LUTs.
+ * \author Sebastian Wagner (TCL Communications) Email: <mailto:sebastian.wagner@tcl.com>
+ * \date 07-12-2018
+ * \version 1.0
+ * \note
+ * \warning
+ */
+
+#ifndef __NR_LDPC_INIT_MEM__H__
+#define __NR_LDPC_INIT_MEM__H__
+
+#include <stdlib.h>
+#include "nrLDPC_defs.h"
+#include "nrLDPC_types.h"
+
+#ifndef malloc32_clear
+/**
+   \brief Allocates 32 byte aligned memory and initializes to zero
+   \param size Input size in bytes
+   \return Pointer to memory
+*/
+static inline void* malloc32_clear(size_t size)
+{
+    void* ptr = (void*) memalign(32, size+32);
+    memset(ptr, 0, size);
+    return ptr;
+}
+#endif
+
+/**
+   \brief Allocates and initializes the internal decoder processing buffers
+   \param p_decParams Pointer to decoder parameters
+   \param p_lut Pointer to decoder LUTs
+   \return Number of LLR values
+*/
+static inline t_nrLDPC_procBuf* nrLDPC_init_mem(void)
+{
+    t_nrLDPC_procBuf* p_procBuf = (t_nrLDPC_procBuf*) malloc32_clear(sizeof(t_nrLDPC_procBuf));
+
+    if (p_procBuf)
+    {
+        p_procBuf->cnProcBuf    = (int8_t*) malloc32_clear(NR_LDPC_SIZE_CN_PROC_BUF*sizeof(int8_t));
+        p_procBuf->cnProcBufRes = (int8_t*) malloc32_clear(NR_LDPC_SIZE_CN_PROC_BUF*sizeof(int8_t));
+        p_procBuf->bnProcBuf    = (int8_t*) malloc32_clear(NR_LDPC_SIZE_BN_PROC_BUF*sizeof(int8_t));
+        p_procBuf->bnProcBufRes = (int8_t*) malloc32_clear(NR_LDPC_SIZE_BN_PROC_BUF*sizeof(int8_t));
+        p_procBuf->llrRes       = (int8_t*) malloc32_clear(NR_LDPC_MAX_NUM_LLR     *sizeof(int8_t));
+        p_procBuf->llrProcBuf   = (int8_t*) malloc32_clear(NR_LDPC_MAX_NUM_LLR     *sizeof(int8_t));
+    }
+
+    return(p_procBuf);
+}
+
+static inline void nrLDPC_free_mem(t_nrLDPC_procBuf* p_procBuf)
+{
+    free(p_procBuf->cnProcBuf);
+    free(p_procBuf->cnProcBufRes);
+    free(p_procBuf->bnProcBuf);
+    free(p_procBuf->bnProcBufRes);
+    free(p_procBuf->llrRes);
+    free(p_procBuf->llrProcBuf);
+
+    free(p_procBuf);
+}
+#endif

openair1/PHY/CODING/nrLDPC_decoder_LYC/nrLDPC_types.h

+/*
+ * Licensed to the OpenAirInterface (OAI) Software Alliance under one or more
+ * contributor license agreements.  See the NOTICE file distributed with
+ * this work for additional information regarding copyright ownership.
+ * The OpenAirInterface Software Alliance licenses this file to You under
+ * the OAI Public License, Version 1.1  (the "License"); you may not use this file
+ * except in compliance with the License.
+ * You may obtain a copy of the License at
+ *
+ *      http://www.openairinterface.org/?page_id=698
+ *
+ * Unless required by applicable law or agreed to in writing, software
+ * distributed under the License is distributed on an "AS IS" BASIS,
+ * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+ * See the License for the specific language governing permissions and
+ * limitations under the License.
+ *-------------------------------------------------------------------------------
+ * For more information about the OpenAirInterface (OAI) Software Alliance:
+ *      contact@openairinterface.org
+ */
+
+/*!\file nrLDPC_types.h
+ * \brief Defines all types for the LDPC decoder
+ * \author Sebastian Wagner (TCL Communications) Email: <mailto:sebastian.wagner@tcl.com>
+ * \date 27-03-2018
+ * \version 1.0
+ * \note
+ * \warning
+ */
+
+#ifndef __NR_LDPC_TYPES__H__
+#define __NR_LDPC_TYPES__H__
+
+#include "PHY/TOOLS/time_meas.h"
+#include "nrLDPC_defs.h"
+
+// ==============================================================================
+// TYPES
+
+/**
+   Structure containing the pointers to the LUTs.
+ */
+typedef struct nrLDPC_lut {
+    const uint32_t* startAddrCnGroups; /**< Start addresses for CN groups in CN processing buffer */
+    const uint8_t*  numCnInCnGroups; /**< Number of CNs in every CN group */
+    const uint8_t*  numBnInBnGroups; /**< Number of CNs in every BN group */
+    const uint32_t* startAddrBnGroups; /**< Start addresses for BN groups in BN processing buffer  */
+    const uint16_t* startAddrBnGroupsLlr; /**< Start addresses for BN groups in LLR processing buffer  */
+    const uint16_t** circShift[NR_LDPC_NUM_CN_GROUPS_BG1]; /**< LUT for circular shift values for all CN groups and Zs */
+    const uint32_t** startAddrBnProcBuf[NR_LDPC_NUM_CN_GROUPS_BG1]; /**< LUT of start addresses of CN groups in BN proc buffer */
+    const uint8_t**  bnPosBnProcBuf[NR_LDPC_NUM_CN_GROUPS_BG1]; /**< LUT of BN positions in BG for CN groups */
+    const uint16_t* llr2llrProcBufAddr; /**< LUT for transferring input LLRs to LLR processing buffer */
+    const uint8_t*  llr2llrProcBufBnPos; /**< LUT BN position in BG */
+    const uint8_t** posBnInCnProcBuf[NR_LDPC_NUM_CN_GROUPS_BG1]; /**< LUT for llr2cnProcBuf */
+} t_nrLDPC_lut;
+
+/**
+   Enum with possible LDPC output formats.
+ */
+typedef enum nrLDPC_outMode {
+    nrLDPC_outMode_BIT, /**< 32 bits per uint32_t output */
+    nrLDPC_outMode_BITINT8, /**< 1 bit per int8_t output */
+    nrLDPC_outMode_LLRINT8 /**< Single LLR value per int8_t output */
+} e_nrLDPC_outMode;
+
+/**
+   Structure containing LDPC decoder parameters.
+ */
+typedef struct nrLDPC_dec_params {
+    uint8_t BG; /**< Base graph */
+    uint16_t Z; /**< Lifting size */
+    uint8_t R; /**< Decoding rate: Format 15,13,... for code rates 1/5, 1/3,... */
+    uint8_t numMaxIter; /**< Maximum number of iterations */
+    e_nrLDPC_outMode outMode; /**< Output format */
+} t_nrLDPC_dec_params;
+
+/**
+   Structure containing LDPC decoder processing time statistics.
+ */
+typedef struct nrLDPC_time_stats {
+    time_stats_t llr2llrProcBuf; /**< Statistics for function llr2llrProcBuf */
+    time_stats_t llr2CnProcBuf; /**< Statistics for function llr2CnProcBuf */
+    time_stats_t cnProc; /**< Statistics for function cnProc */
+    time_stats_t cnProcPc; /**< Statistics for function cnProcPc */
+    time_stats_t bnProcPc; /**< Statistics for function bnProcPc */
+    time_stats_t bnProc; /**< Statistics for function bnProc */
+    time_stats_t cn2bnProcBuf; /**< Statistics for function cn2bnProcBuf */
+    time_stats_t bn2cnProcBuf; /**< Statistics for function bn2cnProcBuf */
+    time_stats_t llrRes2llrOut; /**< Statistics for function llrRes2llrOut */
+    time_stats_t llr2bit; /**< Statistics for function llr2bit */
+    time_stats_t total; /**< Statistics for total processing time */
+} t_nrLDPC_time_stats;
+
+/**
+   Structure containing the processing buffers
+ */
+typedef struct nrLDPC_procBuf {
+    int8_t* cnProcBuf; /**< CN processing buffer */
+    int8_t* cnProcBufRes; /**< Buffer for CN processing results */
+    int8_t* bnProcBuf; /**< BN processing buffer */
+    int8_t* bnProcBufRes; /**< Buffer for BN processing results */
+    int8_t* llrRes; /**< Buffer for LLR results */
+    int8_t* llrProcBuf; /**< LLR processing buffer */
+} t_nrLDPC_procBuf;
+
+#endif