Skip to content
Projects
Groups
Snippets
Help
Loading...
Help
Support
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
O
OpenXG-RAN
Project overview
Project overview
Details
Activity
Releases
Repository
Repository
Files
Commits
Branches
Tags
Contributors
Graph
Compare
Issues
0
Issues
0
List
Boards
Labels
Milestones
Merge Requests
0
Merge Requests
0
CI / CD
CI / CD
Pipelines
Jobs
Schedules
Analytics
Analytics
CI / CD
Repository
Value Stream
Wiki
Wiki
Snippets
Snippets
Members
Members
Collapse sidebar
Close sidebar
Activity
Graph
Create a new issue
Jobs
Commits
Issue Boards
Open sidebar
canghaiwuhen
OpenXG-RAN
Commits
b18bfe24
Commit
b18bfe24
authored
May 13, 2020
by
OpInConnect_NCTU Terngyin/NY/GK/KM
Browse files
Options
Browse Files
Download
Email Patches
Plain Diff
NCTU OpInConnect
parent
e0712dd8
Changes
9
Show whitespace changes
Inline
Side-by-side
Showing
9 changed files
with
409 additions
and
53 deletions
+409
-53
cmake_targets/build_oai
cmake_targets/build_oai
+10
-1
common/utils/LOG/vcd_signal_dumper.c
common/utils/LOG/vcd_signal_dumper.c
+20
-0
common/utils/LOG/vcd_signal_dumper.h
common/utils/LOG/vcd_signal_dumper.h
+20
-0
common/utils/T/T_defs.h
common/utils/T/T_defs.h
+9
-10
common/utils/T/T_messages.txt
common/utils/T/T_messages.txt
+51
-0
openair1/PHY/CODING/TESTBENCH/ldpctest.c
openair1/PHY/CODING/TESTBENCH/ldpctest.c
+267
-11
openair1/PHY/CODING/nrLDPC_encoder/defs.h
openair1/PHY/CODING/nrLDPC_encoder/defs.h
+32
-0
openair1/PHY/NR_TRANSPORT/nr_dlsch.c
openair1/PHY/NR_TRANSPORT/nr_dlsch.c
+0
-10
openair1/SIMULATION/NR_PHY/dlsim.c
openair1/SIMULATION/NR_PHY/dlsim.c
+0
-21
No files found.
cmake_targets/build_oai
View file @
b18bfe24
...
...
@@ -24,6 +24,15 @@
# brief OAI automated build tool that can be used to install, compile, run OAI.
# author Navid Nikaein, Lionel GAUTHIER, Laurent Thomas
# file build_oai
# brief Adjust simlist
# author Terngyin, NY, GK, KM (ISIP)
# email tyhsu@cs.nctu.edu.tw
# date 07-04-2020
# version 1.0
# note
# warning
set
-e
# Include helper functions
...
...
@@ -661,7 +670,7 @@ function main() {
if
[
"
$SIMUS_PHY
"
=
"1"
]
;
then
echo_info
"Compiling physical unitary tests simulators"
# TODO: fix: dlsim_tm4 pucchsim prachsim pdcchsim pbchsim mbmssim
simlist
=
"nr_dlsim "
simlist
=
"nr_dlsim
ldpctest
"
#simlist="dlsim ulsim ldpctest polartest smallblocktest nr_pbchsim nr_dlschsim nr_ulschsim nr_dlsim nr_ulsim nr_pucchsim"
for
f
in
$simlist
;
do
compilations
\
...
...
common/utils/LOG/vcd_signal_dumper.c
View file @
b18bfe24
...
...
@@ -31,6 +31,16 @@
* \warning
*/
/*!\file vcd_signal_dumper.c
* \brief Add some VCD signals of multi_ldpc_encoder
* \author Terngyin, NY, GK, KM (ISIP)
* \email tyhsu@cs.nctu.edu.tw
* \date 07-04-2020
* \version 1.0
* \note
* \warning
*/
#include <pthread.h>
#include <stdio.h>
#include <stdlib.h>
...
...
@@ -63,6 +73,7 @@ struct vcd_module_s {
const
char
*
eurecomVariablesNames
[]
=
{
"frame_number_TX0_eNB"
,
"multi_enc_finish"
,
"mask_ru"
,
"mask_tx_ru"
,
"frame_number_TX1_eNB"
,
...
...
@@ -254,6 +265,15 @@ const char* eurecomVariablesNames[] = {
const
char
*
eurecomFunctionsNames
[]
=
{
/* softmodem signals */
"rt_sleep"
,
"multi_enc"
,
"multi_enc_0"
,
"multi_enc_1"
,
"multi_enc_2"
,
"multi_enc_3"
,
"multi_enc_4"
,
"multi_enc_5"
,
"multi_enc_6"
,
"multi_enc_7"
,
"trx_read"
,
"trx_write"
,
"trx_read_ue"
,
...
...
common/utils/LOG/vcd_signal_dumper.h
View file @
b18bfe24
...
...
@@ -32,6 +32,16 @@
* \warning
*/
/*!\file vcd_signal_dumper.h
* \brief Add some VCD signals of multi_ldpc_encoder
* \author Terngyin, NY, GK, KM (ISIP)
* \email tyhsu@cs.nctu.edu.tw
* \date 07-04-2020
* \version 1.0
* \note
* \warning
*/
#ifndef VCD_SIGNAL_DUMPER_H_
#define VCD_SIGNAL_DUMPER_H_
...
...
@@ -40,6 +50,7 @@
/* WARNING: if you edit the enums below, update also string definitions in vcd_signal_dumper.c */
typedef
enum
{
VCD_SIGNAL_DUMPER_VARIABLES_FRAME_NUMBER_TX0_ENB
=
0
,
VCD_SIGNAL_DUMPER_VARIABLES_MULTI_ENC_FINISH
,
VCD_SIGNAL_DUMPER_VARIABLES_MASK_RU
,
VCD_SIGNAL_DUMPER_VARIABLES_MASK_TX_RU
,
VCD_SIGNAL_DUMPER_VARIABLES_FRAME_NUMBER_TX1_ENB
,
...
...
@@ -234,6 +245,15 @@ typedef enum {
typedef
enum
{
/* softmodem signals */
VCD_SIGNAL_DUMPER_FUNCTIONS_RT_SLEEP
=
0
,
VCD_SIGNAL_DUMPER_FUNCTIONS_MULTI_ENC
,
VCD_SIGNAL_DUMPER_FUNCTIONS_MULTI_ENC_0
,
VCD_SIGNAL_DUMPER_FUNCTIONS_MULTI_ENC_1
,
VCD_SIGNAL_DUMPER_FUNCTIONS_MULTI_ENC_2
,
VCD_SIGNAL_DUMPER_FUNCTIONS_MULTI_ENC_3
,
VCD_SIGNAL_DUMPER_FUNCTIONS_MULTI_ENC_4
,
VCD_SIGNAL_DUMPER_FUNCTIONS_MULTI_ENC_5
,
VCD_SIGNAL_DUMPER_FUNCTIONS_MULTI_ENC_6
,
VCD_SIGNAL_DUMPER_FUNCTIONS_MULTI_ENC_7
,
VCD_SIGNAL_DUMPER_FUNCTIONS_TRX_READ
,
VCD_SIGNAL_DUMPER_FUNCTIONS_TRX_WRITE
,
VCD_SIGNAL_DUMPER_FUNCTIONS_TRX_READ_UE
,
...
...
common/utils/T/T_defs.h
View file @
b18bfe24
...
...
@@ -19,16 +19,15 @@
* contact@openairinterface.org
*/
/*!\file common/utils/T/T_defs.h
* \brief
* \author
* \date
* \version
* \company EURECOM
* \email
/*!\file T_defs.h
* \brief Update VCD_NUM_FUNCTIONS and VCD_NUM_VARIABLES
* \author Terngyin, NY, GK, KM (ISIP)
* \email tyhsu@cs.nctu.edu.tw
* \date 07-04-2020
* \version 1.0
* \note
* \warning
*/
*/
#ifndef _T_defs_H_
#define _T_defs_H_
...
...
@@ -73,10 +72,10 @@ typedef struct {
}
T_cache_t
;
/* number of VCD functions (to be kept up to date! see in T_messages.txt) */
#define VCD_NUM_FUNCTIONS (2
47)
#define VCD_NUM_FUNCTIONS (2
56) //origin : 247
/* number of VCD variables (to be kept up to date! see in T_messages.txt) */
#define VCD_NUM_VARIABLES (18
6)
#define VCD_NUM_VARIABLES (18
7) //origin : 186
/* first VCD function (to be kept up to date! see in T_messages.txt) */
#define VCD_FIRST_FUNCTION ((uintptr_t)T_VCD_FUNCTION_RT_SLEEP)
...
...
common/utils/T/T_messages.txt
View file @
b18bfe24
#general logs
ID = ENB_MASTER_TICK
DESC = eNodeB master tick - one tick per ms, to be used as "reference clock", mostly for ticktime view
...
...
@@ -1211,6 +1212,11 @@ ID = VCD_VARIABLE_FRAME_NUMBER_TX0_ENB
GROUP = ALL:VCD:ENB:VCD_VARIABLE
FORMAT = ulong,value
VCD_NAME = frame_number_TX0_eNB
ID = VCD_VARIABLE_MULTI_ENC_FINISH
DESC = VCD variable MULTI_ENC_FINISH
GROUP = ALL:VCD:ENB:VCD_VARIABLE
FORMAT = ulong,value
VCD_NAME = multi_enc_finish
ID = VCD_VARIABLE_MASK_RU
DESC = VCD variable MASK_RU
GROUP = ALL:VCD:ENB:VCD_VARIABLE
...
...
@@ -2146,6 +2152,51 @@ ID = VCD_FUNCTION_RT_SLEEP
GROUP = ALL:VCD:ENB:VCD_FUNCTION
FORMAT = int,value
VCD_NAME = rt_sleep
ID = VCD_FUNCTION_MULTI_ENC
DESC = VCD function MULTI_ENC
GROUP = ALL:VCD:ENB:VCD_FUNCTION
FORMAT = int,value
VCD_NAME = multi_enc
ID = VCD_FUNCTION_MULTI_ENC_0
DESC = VCD function MULTI_ENC_0
GROUP = ALL:VCD:ENB:VCD_FUNCTION
FORMAT = int,value
VCD_NAME = multi_enc_0
ID = VCD_FUNCTION_MULTI_ENC_1
DESC = VCD function MULTI_ENC_1
GROUP = ALL:VCD:ENB:VCD_FUNCTION
FORMAT = int,value
VCD_NAME = multi_enc_1
ID = VCD_FUNCTION_MULTI_ENC_2
DESC = VCD function MULTI_ENC_2
GROUP = ALL:VCD:ENB:VCD_FUNCTION
FORMAT = int,value
VCD_NAME = multi_enc_2
ID = VCD_FUNCTION_MULTI_ENC_3
DESC = VCD function MULTI_ENC_3
GROUP = ALL:VCD:ENB:VCD_FUNCTION
FORMAT = int,value
VCD_NAME = multi_enc_3
ID = VCD_FUNCTION_MULTI_ENC_4
DESC = VCD function MULTI_ENC_4
GROUP = ALL:VCD:ENB:VCD_FUNCTION
FORMAT = int,value
VCD_NAME = multi_enc_4
ID = VCD_FUNCTION_MULTI_ENC_5
DESC = VCD function MULTI_ENC_5
GROUP = ALL:VCD:ENB:VCD_FUNCTION
FORMAT = int,value
VCD_NAME = multi_enc_5
ID = VCD_FUNCTION_MULTI_ENC_6
DESC = VCD function MULTI_ENC_6
GROUP = ALL:VCD:ENB:VCD_FUNCTION
FORMAT = int,value
VCD_NAME = multi_enc_6
ID = VCD_FUNCTION_MULTI_ENC_7
DESC = VCD function MULTI_ENC_7
GROUP = ALL:VCD:ENB:VCD_FUNCTION
FORMAT = int,value
VCD_NAME = multi_enc_7
ID = VCD_FUNCTION_TRX_READ
DESC = VCD function TRX_READ
GROUP = ALL:VCD:ENB:VCD_FUNCTION
...
...
openair1/PHY/CODING/TESTBENCH/ldpctest.c
100755 → 100644
View file @
b18bfe24
...
...
@@ -19,6 +19,16 @@
* contact@openairinterface.org
*/
/*!\file ldpctest.c
* \brief Implementation of multi_ldpc_encoder
* \author Terngyin, NY, GK, KM (ISIP)
* \email tyhsu@cs.nctu.edu.tw
* \date 07-04-2020
* \version 1.0
* \note
* \warning
*/
#include <stdlib.h>
#include <math.h>
#include <stdio.h>
...
...
@@ -28,6 +38,11 @@
#include "PHY/CODING/nrLDPC_encoder/defs.h"
#include "PHY/CODING/nrLDPC_decoder/nrLDPC_decoder.h"
#include "openair1/SIMULATION/NR_PHY/nr_unitary_defs.h"
#include <time.h>
//#include "T.h"
//#include "common/utils/LOG/log.h"
#include "common/utils/LOG/vcd_signal_dumper.h"
#define MAX_NUM_DLSCH_SEGMENTS 16
#define MAX_BLOCK_LENGTH 8448
...
...
@@ -88,6 +103,80 @@ uint16_t NB_UE_INST = 1;
short
lift_size
[
51
]
=
{
2
,
3
,
4
,
5
,
6
,
7
,
8
,
9
,
10
,
11
,
12
,
13
,
14
,
15
,
16
,
18
,
20
,
22
,
24
,
26
,
28
,
30
,
32
,
36
,
40
,
44
,
48
,
52
,
56
,
60
,
64
,
72
,
80
,
88
,
96
,
104
,
112
,
120
,
128
,
144
,
160
,
176
,
192
,
208
,
224
,
240
,
256
,
288
,
320
,
352
,
384
};
/*! \file openair1/PHY/CODING/TESTBENCH/ldpctest.c
* \brief NCTU OpInConnect
* \author Terngyin Hsu, Sendren Xu, Nungyi Kuo, Kuankai Hsiung, Kaimi Yang (OpInConnect_NCTU)
* \email tyhsu@cs.nctu.edu.tw
* \date 13-05-2020
* \version 2.0
* \note
* \warning
*/
// ==[START]multi_lpdc_encoder
struct
timespec
start_ts
,
end_ts
,
start_per_ts
,
end_per_ts
,
start_enc_ts
[
thread_num_max
],
end_enc_ts
[
thread_num_max
],
start_perenc_ts
[
thread_num_max
],
end_perenc_ts
[
thread_num_max
];
multi_ldpc_encoder
ldpc_enc
[
thread_num_max
];
int
thread_num
=
2
;
//Craete 2 threads
int
ifRand
=
0
;
int
vcd
=
0
;
// unsigned char *test_input0[MAX_NUM_DLSCH_SEGMENTS]={NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL};;
// unsigned char *test_input1[MAX_NUM_DLSCH_SEGMENTS]={NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL};;
// unsigned char *channel_input_optim0[MAX_NUM_DLSCH_SEGMENTS];
// unsigned char *channel_input_optim1[MAX_NUM_DLSCH_SEGMENTS];
static
void
*
multi_ldpc_encoder_proc
(
void
*
ptr
){
//ldpc_encoder_optim_8seg_multi(dlsch->harq_processes[harq_pid]->c,dlsch->harq_processes[harq_pid]->d,*Zc,Kb,Kr,BG,dlsch->harq_processes[harq_pid]->C,j,NULL,NULL,NULL,NULL);
//ldpc_encoder_optim_8seg_multi(unsigned char **test_input,unsigned char **channel_input,int Zc,int Kb,short block_length, short BG, int n_segments,unsigned int macro_num, time_stats_t *tinput,time_stats_t *tprep,time_stats_t *tparity,time_stats_t *toutput)
multi_ldpc_encoder
*
test
=
(
multi_ldpc_encoder
*
)
ptr
;
printf
(
"[READY] : %d
\n
"
,
test
->
id
);
while
(
pthread_cond_wait
(
&
ldpc_enc
[
test
->
id
].
cond
,
&
ldpc_enc
[
test
->
id
].
mutex
)
!=
0
);
//Print params of multi_ldpc_enc
// printf("[b]ldpc_encoder_optim_8seg: BG %d, Zc %d, Kb %d, block_length %d, segments %d\n",
// ldpc_enc[test->id].BG,
// ldpc_enc[test->id].Zc,
// ldpc_enc[test->id].Kb,
// ldpc_enc[test->id].block_length,
// ldpc_enc[test->id].n_segments);
int
j_start
,
j_end
;
//Check if our situation((n_segments == 15)&(thread_num == 2))
if
((
ldpc_enc
[
test
->
id
].
n_segments
==
15
)
&&
(
thread_num
==
2
)){
j_start
=
test
->
id
;
j_end
=
j_start
+
1
;
}
else
{
j_start
=
0
;
j_end
=
(
ldpc_enc
[
test
->
id
].
n_segments
/
8
+
1
);
}
int
offset
=
test
->
id
>
7
?
7
:
test
->
id
;
//printf("[OFFSET] : %d %d\n", offset, test->id);
VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME
(
VCD_SIGNAL_DUMPER_FUNCTIONS_MULTI_ENC_0
+
offset
,
1
);
clock_gettime
(
CLOCK_MONOTONIC
,
&
start_enc_ts
[
test
->
id
]);
//timing
for
(
int
j
=
j_start
;
j
<
j_end
;
j
++
){
//printf(" Movement No. Round Cost time \n");
printf
(
" Active %d %d
\n
"
,
test
->
id
,
j
);
clock_gettime
(
CLOCK_MONOTONIC
,
&
start_perenc_ts
[
test
->
id
]);
//timing
ldpc_encoder_optim_8seg_multi
(
ldpc_enc
[
test
->
id
].
test_input
,
ldpc_enc
[
test
->
id
].
channel_input_optim
,
ldpc_enc
[
test
->
id
].
Zc
,
ldpc_enc
[
test
->
id
].
Kb
,
ldpc_enc
[
test
->
id
].
block_length
,
ldpc_enc
[
test
->
id
].
BG
,
ldpc_enc
[
test
->
id
].
n_segments
,
j
,
NULL
,
NULL
,
NULL
,
NULL
);
clock_gettime
(
CLOCK_MONOTONIC
,
&
end_perenc_ts
[
test
->
id
]);
//timing
//printf(" Movement No. Round Cost time \n");
printf
(
" Done %d %d %.2f usec
\n
"
,
test
->
id
,
j
,
(
end_perenc_ts
[
test
->
id
].
tv_nsec
-
start_perenc_ts
[
test
->
id
].
tv_nsec
)
*
1
.
0
/
1000
);
}
clock_gettime
(
CLOCK_MONOTONIC
,
&
end_enc_ts
[
test
->
id
]);
//timing
VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME
(
VCD_SIGNAL_DUMPER_FUNCTIONS_MULTI_ENC_0
+
offset
,
0
);
//printf("[OK] : %d : %.2f usec\n", test->id, (end_perenc_ts[test->id].tv_nsec - start_perenc_ts[test->id].tv_nsec) *1.0 / 1000);
ldpc_enc
[
test
->
id
].
complete
=
1
;
return
0
;
}
// ==[END]multi_lpdc_encoder
int
test_ldpc
(
short
No_iteration
,
int
nom_rate
,
int
denom_rate
,
...
...
@@ -119,6 +208,7 @@ int test_ldpc(short No_iteration,
cpu_freq_GHz
=
get_cpu_freq_GHz
();
//short test_input[block_length];
unsigned
char
*
test_input
[
MAX_NUM_DLSCH_SEGMENTS
]
=
{
NULL
,
NULL
,
NULL
,
NULL
,
NULL
,
NULL
,
NULL
,
NULL
};;
unsigned
char
*
test_input_multi
[
thread_num_max
][
MAX_NUM_DLSCH_SEGMENTS
];
// ==For multi_ldpc_enc ==
//short *c; //padded codeword
unsigned
char
*
estimated_output
[
MAX_NUM_DLSCH_SEGMENTS
];
unsigned
char
*
estimated_output_bit
[
MAX_NUM_DLSCH_SEGMENTS
];
...
...
@@ -126,6 +216,7 @@ int test_ldpc(short No_iteration,
unsigned
char
*
channel_input
[
MAX_NUM_DLSCH_SEGMENTS
];
unsigned
char
*
channel_output_uncoded
[
MAX_NUM_DLSCH_SEGMENTS
];
unsigned
char
*
channel_input_optim
[
MAX_NUM_DLSCH_SEGMENTS
];
unsigned
char
*
channel_input_optim_multi
[
thread_num_max
][
MAX_NUM_DLSCH_SEGMENTS
];
double
*
channel_output
;
double
*
modulated_input
[
MAX_NUM_DLSCH_SEGMENTS
];
char
*
channel_output_fixed
[
MAX_NUM_DLSCH_SEGMENTS
];
...
...
@@ -159,6 +250,10 @@ int test_ldpc(short No_iteration,
test_input
[
j
]
=
(
unsigned
char
*
)
malloc16
(
sizeof
(
unsigned
char
)
*
block_length
/
8
);
channel_input
[
j
]
=
(
unsigned
char
*
)
malloc16
(
sizeof
(
unsigned
char
)
*
68
*
384
);
channel_input_optim
[
j
]
=
(
unsigned
char
*
)
malloc16
(
sizeof
(
unsigned
char
)
*
68
*
384
);
for
(
int
th
=
0
;
th
<
thread_num
;
th
++
){
// ==For multi_ldpc_enc ==
test_input_multi
[
th
][
j
]
=
(
unsigned
char
*
)
malloc16
(
sizeof
(
unsigned
char
)
*
block_length
/
8
);
channel_input_optim_multi
[
th
][
j
]
=
(
unsigned
char
*
)
malloc16
(
sizeof
(
unsigned
char
)
*
68
*
384
);
}
channel_output_uncoded
[
j
]
=
(
unsigned
char
*
)
malloc16
(
sizeof
(
unsigned
char
)
*
68
*
384
);
estimated_output
[
j
]
=
(
unsigned
char
*
)
malloc16
(
sizeof
(
unsigned
char
)
*
block_length
);
estimated_output_bit
[
j
]
=
(
unsigned
char
*
)
malloc16
(
sizeof
(
unsigned
char
)
*
block_length
);
...
...
@@ -202,8 +297,24 @@ int test_ldpc(short No_iteration,
for
(
j
=
0
;
j
<
MAX_NUM_DLSCH_SEGMENTS
;
j
++
)
{
for
(
i
=
0
;
i
<
block_length
/
8
;
i
++
)
{
test_input
[
j
][
i
]
=
(
unsigned
char
)
rand
();
//test_input[j][i]=j%256;
//test_input[j][i]=252;
for
(
int
th
=
0
;
th
<
thread_num
;
th
++
){
if
(
ifRand
){
test_input_multi
[
th
][
j
][
i
]
=
(
unsigned
char
)
rand
();
// use rand value
}
else
{
test_input_multi
[
th
][
j
][
i
]
=
test_input
[
j
][
i
];
// use the same value
}
}
}
}
//Print out test_input
printf
(
"===============[test_input]================
\n
"
);
for
(
j
=
0
;
j
<
1
;
j
++
)
{
for
(
i
=
0
;
i
<
10
;
i
++
)
{
printf
(
"%u"
,
test_input
[
j
][
i
]);
for
(
int
th
=
0
;
th
<
thread_num
;
th
++
){
printf
(
", %u"
,
test_input_multi
[
th
][
j
][
i
]);
}
printf
(
"
\n
"
);
}
}
...
...
@@ -281,7 +392,7 @@ int test_ldpc(short No_iteration,
break
;
}
}
printf
(
"============== Check params ===============
\n
"
);
printf
(
"ldpc_test: codeword_length %d, n_segments %d, block_length %d, BG %d, Zc %d, Kb %d
\n
"
,
n_segments
*
block_length
,
n_segments
,
block_length
,
BG
,
Zc
,
Kb
);
no_punctured_columns
=
(
int
)((
nrows
-
2
)
*
Zc
+
block_length
-
block_length
*
(
1
/
((
float
)
nom_rate
/
(
float
)
denom_rate
)))
/
Zc
;
// printf("puncture:%d\n",no_punctured_columns);
...
...
@@ -305,21 +416,103 @@ int test_ldpc(short No_iteration,
ldpc_encoder_optim(test_input[j],channel_input_optim[j],Zc,Kb,block_length,BG,&tinput,&tprep,&tparity,&toutput);
}
stop_meas(time_optim);*/
// unsigned char c_extension_try[2*22*384*32] __attribute__((aligned(32))); //double size matrix of c
// ==[START]multi_ldpc_enc
printf
(
"============== Multi threads ==============
\n
"
);
printf
(
" [Movement] [No.] [Part] [Cost time]
\n
"
);
for
(
int
th
=
0
;
th
<
thread_num
;
th
++
){
ldpc_enc
[
th
].
Zc
=
Zc
;
ldpc_enc
[
th
].
Kb
=
Kb
;
ldpc_enc
[
th
].
block_length
=
block_length
;
ldpc_enc
[
th
].
BG
=
BG
;
ldpc_enc
[
th
].
n_segments
=
n_segments
;
// ==Do all n_segments ==
//ldpc_enc[th].n_segments = th+1; // ==Do n_segments increasing ==
//ldpc_enc[th].n_segments = 1; // ==Just do first segment ==
ldpc_enc
[
th
].
test_input
=
test_input_multi
[
th
];
// It's a ptr @ the head of "th"th line
ldpc_enc
[
th
].
channel_input_optim
=
channel_input_optim_multi
[
th
];
// It's a ptr @ the head of "th"th line
}
VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME
(
VCD_SIGNAL_DUMPER_FUNCTIONS_MULTI_ENC
,
1
);
clock_gettime
(
CLOCK_MONOTONIC
,
&
start_ts
);
//timing
for
(
int
th
=
0
;
th
<
thread_num
;
th
++
){
pthread_cond_signal
(
&
ldpc_enc
[
th
].
cond
);
}
for
(
int
th
=
0
;
th
<
thread_num
;
th
++
){
while
(
ldpc_enc
[
th
].
complete
!=
1
);
// ==check if multi_ldpc_enc done ==
VCD_SIGNAL_DUMPER_DUMP_VARIABLE_BY_NAME
(
VCD_SIGNAL_DUMPER_VARIABLES_MULTI_ENC_FINISH
,
th
);
}
clock_gettime
(
CLOCK_MONOTONIC
,
&
end_ts
);
//timing
VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME
(
VCD_SIGNAL_DUMPER_FUNCTIONS_MULTI_ENC
,
0
);
for
(
int
th
=
0
;
th
<
thread_num
;
th
++
){
pthread_mutex_destroy
(
&
ldpc_enc
[
th
].
mutex
);
pthread_join
(
ldpc_enc
[
th
].
pthread
,
NULL
);
}
printf
(
"========== Time of Multi threads ==========
\n
"
);
printf
(
" [No.] [Total time]
\n
"
);
for
(
int
th
=
0
;
th
<
thread_num
;
th
++
){
printf
(
" %d %.2f usec
\n
"
,
th
,
(
end_enc_ts
[
th
].
tv_nsec
-
start_enc_ts
[
th
].
tv_nsec
)
*
1
.
0
/
1000
);
}
// printf(" t %.2f usec\n", (end_enc_ts[0].tv_nsec - start_enc_ts[0].tv_nsec) *1.0 / 1000);
// printf(" t %.2f usec\n", (end_enc_ts[1].tv_nsec - start_enc_ts[0].tv_nsec) *1.0 / 1000);
// printf(" t %.2f usec\n", (end_enc_ts[0].tv_nsec - start_enc_ts[1].tv_nsec) *1.0 / 1000);
// printf(" t %.2f usec\n", (end_enc_ts[1].tv_nsec - start_enc_ts[1].tv_nsec) *1.0 / 1000);
printf
(
" Total %.2f usec
\n
"
,
(
end_ts
.
tv_nsec
-
start_ts
.
tv_nsec
)
*
1
.
0
/
1000
);
// ==[END]multi_ldpc_enc
printf
(
"============== Single thread ==============
\n
"
);
printf
(
" [Movement] [No.] [Part] [Cost time]
\n
"
);
clock_gettime
(
CLOCK_MONOTONIC
,
&
start_ts
);
//timing
for
(
j
=
0
;
j
<
(
n_segments
/
8
+
1
);
j
++
)
{
printf
(
" Active Single %d
\n
"
,
j
);
start_meas
(
time_optim
);
clock_gettime
(
CLOCK_MONOTONIC
,
&
start_per_ts
);
//timing
ldpc_encoder_optim_8seg_multi
(
test_input
,
channel_input_optim
,
Zc
,
Kb
,
block_length
,
BG
,
n_segments
,
j
,
&
tinput
,
&
tprep
,
&
tparity
,
&
toutput
);
clock_gettime
(
CLOCK_MONOTONIC
,
&
end_per_ts
);
//timing
stop_meas
(
time_optim
);
printf
(
" Done Single %d %.2f usec
\n
"
,
j
,
(
end_per_ts
.
tv_nsec
-
start_per_ts
.
tv_nsec
)
*
1
.
0
/
1000
);
}
clock_gettime
(
CLOCK_MONOTONIC
,
&
end_ts
);
//timing
printf
(
"========== Time of Single thread ==========
\n
"
);
printf
(
" [No.] [Total time]
\n
"
);
printf
(
" Total %.2f usec
\n
"
,
(
end_ts
.
tv_nsec
-
start_ts
.
tv_nsec
)
*
1
.
0
/
1000
);
if
((
n_segments
!=
15
)
||
(
thread_num
!=
2
))
//if not our situation
printf
(
"Total*%d %.2f usec
\n
"
,
thread_num
,
(
end_ts
.
tv_nsec
-
start_ts
.
tv_nsec
)
*
1
.
0
*
thread_num
/
1000
);
//Print out channel_input
// printf("[channel_input]:\n");
// for (j=0;j<n_segments;j++){
// for (i = 0; i < 10; i++){
// printf("%u,%u\n", channel_input[j][i], channel_input_optim_multi[0][j][i]);
// }
// }
// ==dec here ==
if
(
ntrials
==
1
)
for
(
j
=
0
;
j
<
n_segments
;
j
++
)
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
]);
printf
(
"
[ERROR][Single]
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
);
}
//[START]Check the value of channel_input and channel_input_multi
if
((
n_segments
!=
15
)
||
(
thread_num
!=
2
)){
//if not our situation
if
(
ntrials
==
1
)
for
(
int
th
=
0
;
th
<
thread_num
;
th
++
)
for
(
j
=
0
;
j
<
n_segments
;
j
++
)
for
(
i
=
0
;
i
<
block_length
+
(
nrows
-
no_punctured_columns
)
*
Zc
-
removed_bit
;
i
++
)
if
(
channel_input
[
j
][
i
]
!=
channel_input_optim_multi
[
0
][
j
][
i
])
{
printf
(
"[ERROR][Multi(%d)]differ in seg %u pos %u (%u,%u)
\n
"
,
th
,
j
,
i
,
channel_input
[
j
][
i
],
channel_input_optim_multi
[
th
][
j
][
i
]);
free
(
channel_output
);
return
(
-
1
);
}
printf
(
"============ [OK]Result compare ===========
\n
"
);
}
else
{
printf
(
"=================== [OK] ==================
\n
"
);
}
//[END]Check the value of channel_input and channel_input_multi
//else{
// printf("NOT differ in seg %d pos %d (%d,%d)\n",j,i,channel_input[j][i],channel_input_optim[j][i]);
// }
...
...
@@ -345,12 +538,17 @@ int test_ldpc(short No_iteration,
if
((
i
&
0xf
)
==
0
)
printf
(
"
\n
e %u..%u: "
,
i
,
i
+
15
);
#endif
// ==dec here ==
if
(
channel_input_optim
[
j
][
i
-
2
*
Zc
]
==
0
)
modulated_input
[
j
][
i
]
=
1
.
0
;
///sqrt(2); //QPSK
else
modulated_input
[
j
][
i
]
=-
1
.
0
;
///sqrt(2);
// if (channel_input_optim_multi[0][i-2*Zc]==0)
// modulated_input[j][i]=1.0;///sqrt(2); //QPSK
// else
// modulated_input[j][i]=-1.0;///sqrt(2);
///channel_output[i] = modulated_input[i] + gaussdouble(0.0,1.0) * 1/sqrt(2*SNR);
//channel_output_fixed[i] = (char) ((channel_output[i]*128)<0?(channel_output[i]*128-0.5):(channel_output[i]*128+0.5)); //fixed point 9-7
//printf("llr[%d]=%d\n",i,channel_output_fixed[i]);
...
...
@@ -369,10 +567,11 @@ int test_ldpc(short No_iteration,
channel_output_uncoded
[
j
][
i
]
=
1
;
//QPSK demod
else
channel_output_uncoded
[
j
][
i
]
=
0
;
// ==dec here ==
if
(
channel_output_uncoded
[
j
][
i
]
!=
channel_input_optim
[
j
][
i
-
2
*
Zc
])
*
errors_bit_uncoded
=
(
*
errors_bit_uncoded
)
+
1
;
// if (channel_output_uncoded[j][i] != channel_input_optim_multi[0][j][i-2*Zc])
// *errors_bit_uncoded = (*errors_bit_uncoded) + 1;
}
}
// End segments
...
...
@@ -416,18 +615,27 @@ int test_ldpc(short No_iteration,
/// 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]);
// ==dec here ==
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]);
segment_bler
=
segment_bler
+
1
;
break
;
}
// if (estimated_output[j][i] != test_input_multi[0][j][i])
// {
// //////printf("error pos %d (%d, %d)\n\n",i,estimated_output[i],test_input[0][i]);
// segment_bler = segment_bler + 1;
// break;
// }
}
for
(
i
=
0
;
i
<
block_length
;
i
++
)
{
estimated_output_bit
[
j
][
i
]
=
(
estimated_output
[
j
][
i
/
8
]
&
(
1
<<
(
i
&
7
)))
>>
(
i
&
7
);
// ==dec here ==
test_input_bit
[
i
]
=
(
test_input
[
j
][
i
/
8
]
&
(
1
<<
(
i
&
7
)))
>>
(
i
&
7
);
// Further correct for multiple segments
// test_input_bit[i] = (test_input_multi[0][j][i/8]&(1<<(i&7)))>>(i&7); // Further correct for multiple segments
if
(
estimated_output_bit
[
j
][
i
]
!=
test_input_bit
[
i
])
{
*
errors_bit
=
(
*
errors_bit
)
+
1
;
...
...
@@ -465,6 +673,10 @@ int test_ldpc(short No_iteration,
free
(
channel_input
[
j
]);
free
(
channel_output_uncoded
[
j
]);
free
(
channel_input_optim
[
j
]);
for
(
int
th
=
0
;
th
<
thread_num
;
th
++
){
// ==For multi_ldpc_enc ==
free
(
test_input_multi
[
th
][
j
]);
free
(
channel_input_optim_multi
[
th
][
j
]);
}
free
(
modulated_input
[
j
]);
free
(
channel_output_fixed
[
j
]);
free
(
estimated_output
[
j
]);
...
...
@@ -507,7 +719,8 @@ int main(int argc, char *argv[])
short
block_length
=
8448
;
// decoder supports length: 1201 -> 1280, 2401 -> 2560
short
No_iteration
=
5
;
int
n_segments
=
1
;
//int n_segments=1;
int
n_segments
=
15
;
//Set n_segments for test
//double rate=0.333;
int
nom_rate
=
1
;
int
denom_rate
=
3
;
...
...
@@ -527,7 +740,7 @@ int main(int argc, char *argv[])
short
BG
=
0
,
Zc
,
Kb
;
while
((
c
=
getopt
(
argc
,
argv
,
"q:r:
s:S:l:n:d:i:t
:u:h"
))
!=
-
1
)
while
((
c
=
getopt
(
argc
,
argv
,
"q:r:
R:s:S:l:n:d:i:t:T:V
:u:h"
))
!=
-
1
)
switch
(
c
)
{
case
'q'
:
...
...
@@ -538,6 +751,12 @@ int main(int argc, char *argv[])
nom_rate
=
atoi
(
optarg
);
break
;
case
'R'
:
ifRand
=
atoi
(
optarg
);
if
(
ifRand
)
printf
(
"[WARNING]It will cause ERROR when comparing results
\n
"
);
break
;
case
'd'
:
denom_rate
=
atoi
(
optarg
);
break
;
...
...
@@ -562,6 +781,21 @@ int main(int argc, char *argv[])
SNR_step
=
atof
(
optarg
);
break
;
case
'T'
:
thread_num
=
atoi
(
optarg
);
if
(
thread_num
>
8
){
thread_num
=
8
;
printf
(
"[WARNING]Too many threads !!! Set thread_num 8
\n
"
);
}
else
if
(
thread_num
<
0
){
thread_num
=
0
;
printf
(
"[WARNING]Too less threads !!! Set thread_num 0
\n
"
);
}
break
;
case
'V'
:
vcd
=
atoi
(
optarg
);
break
;
case
'i'
:
No_iteration
=
atoi
(
optarg
);
break
;
...
...
@@ -578,6 +812,7 @@ int main(int argc, char *argv[])
printf
(
"-h This message
\n
"
);
printf
(
"-q Quantization bits, Default: 8
\n
"
);
printf
(
"-r Nominator rate, (1, 2, 22), Default: 1
\n
"
);
printf
(
"-R Will test_input be different for each thread, Default: 0
\n
"
);
printf
(
"-d Denominator rate, (3, 5, 25), Default: 1
\n
"
);
printf
(
"-l Block length (l > 3840 -> BG1, rest BG2 ), Default: 8448
\n
"
);
printf
(
"-n Number of simulation trials, Default: 1
\n
"
);
...
...
@@ -585,11 +820,32 @@ int main(int argc, char *argv[])
printf
(
"-s SNR per information bit (EbNo) in dB, Default: -2
\n
"
);
printf
(
"-S Number of segments (Maximum: 8), Default: 1
\n
"
);
printf
(
"-t SNR simulation step, Default: 0.1
\n
"
);
printf
(
"-T Number of threads, Default: 2
\n
"
);
printf
(
"-V VCD ON(1), OFF(0), Default: 0
\n
"
);
printf
(
"-i Max decoder iterations, Default: 5
\n
"
);
printf
(
"-u Set SNR per coded bit, Default: 0
\n
"
);
exit
(
1
);
break
;
}
// ==[START]VCD ==
if
(
vcd
==
1
){
T_init
(
2021
,
1
,
0
);
}
// ==[END]VCD ==
// ==[START]Create multi_ldpc_encoder ==
for
(
int
th
=
0
;
th
<
thread_num
;
th
++
){
pthread_attr_init
(
&
ldpc_enc
[
th
].
attr
);
pthread_mutex_init
(
&
ldpc_enc
[
th
].
mutex
,
NULL
);
pthread_cond_init
(
&
ldpc_enc
[
th
].
cond
,
NULL
);
ldpc_enc
[
th
].
id
=
th
;
ldpc_enc
[
th
].
flag_wait
=
1
;
ldpc_enc
[
th
].
complete
=
0
;
pthread_create
(
&
(
ldpc_enc
[
th
].
pthread
),
&
ldpc_enc
[
th
].
attr
,
multi_ldpc_encoder_proc
,
&
ldpc_enc
[
th
]);
printf
(
"[CREATE] LDPC encoder thread %d
\n
"
,
ldpc_enc
[
th
].
id
);
}
// ==[END]Create multi_ldpc_encoder ==
cpu_freq_GHz
=
get_cpu_freq_GHz
();
//printf("the decoder supports BG2, Kb=10, Z=128 & 256\n");
//printf(" range of blocklength: 1201 -> 1280, 2401 -> 2560\n");
...
...
openair1/PHY/CODING/nrLDPC_encoder/defs.h
100755 → 100644
View file @
b18bfe24
...
...
@@ -29,8 +29,40 @@
* \warning
*/
/*!\file defs.h
* \brief Declarations of parameters of multi_ldpc_encoder
* \author Terngyin, NY, GK, KM (ISIP)
* \email tyhsu@cs.nctu.edu.tw
* \date 07-04-2020
* \version 1.0
* \note
* \warning
*/
#include "PHY/TOOLS/time_meas.h"
// ==[START]multi_ldpc_encoder
//#define thread_num 2 // ==Change thread_num here ==
#define thread_num_max 8 // ==Change max tread_num here ==
typedef
struct
{
int
id
;
volatile
int
flag_wait
;
pthread_t
pthread
;
pthread_cond_t
cond
;
pthread_mutex_t
mutex
;
pthread_attr_t
attr
;
volatile
uint8_t
complete
;
unsigned
char
**
test_input
;
unsigned
char
**
channel_input_optim
;
int
Zc
;
int
Kb
;
short
block_length
;
short
BG
;
int
n_segments
;
//unsigned int macro_num; //Not necessary to do
}
multi_ldpc_encoder
;
// ==[END]multi_ldpc_encoder
/*ldpc_encoder.c*/
int
encode_parity_check_part_orig
(
unsigned
char
*
c
,
unsigned
char
*
d
,
short
BG
,
short
Zc
,
short
Kb
,
short
block_length
);
...
...
openair1/PHY/NR_TRANSPORT/nr_dlsch.c
View file @
b18bfe24
...
...
@@ -125,16 +125,6 @@ for (int i=0; i<encoded_length>>3; i++) {
printf
(
"
\t
"
);
}
printf
(
"
\n
"
);
/*! \file openair1/PHY/NR_TRANSPORT/nr_dlsch.c
* \brief parallel_architecture of 5G NR V2
* \author Terng-Yin Hsu, WEI-YING,LIN (OpInConnect_NCTU)
* \email tyhsu@cs.nctu.edu.tw
* \date 13-05-2020
* \version 1.3
* \note
* \warning
*/
#endif
long
sum
=
0
;
#ifdef thread_for_scrambling_modulation
...
...
openair1/SIMULATION/NR_PHY/dlsim.c
View file @
b18bfe24
...
...
@@ -284,18 +284,6 @@ static void *dlsch_encoding_proc(void *ptr){
static
void
*
scrambling_proc
(
void
*
ptr
){
/*! \file openair1/SIMULATION/NR_PHY/dlsim.c
* \brief parallel_architecture of 5G NR V2
* \author Terng-Yin Hsu, WEI-YING,LIN (OpInConnect_NCTU)
* \email tyhsu@cs.nctu.edu.tw
* \date 13-05-2020
* \version 1.3
* \note
* \warning
*/
scrambling_channel
*
test
=
(
scrambling_channel
*
)
ptr
;
int
q_id
=
test
->
q_id
;
struct
timespec
tt1
,
tt2
,
tt3
,
tt4
;
...
...
@@ -420,15 +408,6 @@ static void *modulation_proc(void *ptr){
struct
timespec
tt1
,
tt2
;
PHY_VARS_gNB
*
gNB
=
RC
.
gNB
[
0
][
0
];
NR_gNB_DLSCH_t
*
dlsch
=
gNB
->
dlsch
[
0
][
0
];
/*! \file openair1/SIMULATION/NR_PHY/dlsim.c
* \brief parallel_architecture of 5G NR V2
* \author Terng-Yin Hsu, WEI-YING,LIN (OpInConnect_NCTU)
* \email tyhsu@cs.nctu.edu.tw
* \date 13-05-2020
* \version 1.3
* \note
* \warning
*/
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment