|
|
/** Functions related to data processing.
|
|
|
*
|
|
|
* @file
|
|
|
* @author P. LEROY
|
|
|
*
|
|
|
* These function are related to data processing, i.e. spectral matrices averaging and basic parameters computation.
|
|
|
*
|
|
|
*/
|
|
|
|
|
|
#include "avf0_prc0.h"
|
|
|
#include "fsw_processing.h"
|
|
|
|
|
|
nb_sm_before_bp_asm_f0 nb_sm_before_f0;
|
|
|
|
|
|
//***
|
|
|
// F0
|
|
|
ring_node_asm asm_ring_norm_f0 [ NB_RING_NODES_ASM_NORM_F0 ];
|
|
|
ring_node_asm asm_ring_burst_sbm_f0 [ NB_RING_NODES_ASM_BURST_SBM_F0 ];
|
|
|
|
|
|
ring_node ring_to_send_asm_f0 [ NB_RING_NODES_ASM_F0 ];
|
|
|
int buffer_asm_f0 [ NB_RING_NODES_ASM_F0 * TOTAL_SIZE_SM ];
|
|
|
|
|
|
float asm_f0_reorganized [ TOTAL_SIZE_SM ];
|
|
|
char asm_f0_char [ TIME_OFFSET_IN_BYTES + (TOTAL_SIZE_SM * 2) ];
|
|
|
float compressed_sm_norm_f0[ TOTAL_SIZE_COMPRESSED_ASM_NORM_F0];
|
|
|
float compressed_sm_sbm_f0 [ TOTAL_SIZE_COMPRESSED_ASM_SBM_F0 ];
|
|
|
|
|
|
float k_coeff_intercalib_f0_norm[ NB_BINS_COMPRESSED_SM_F0 * NB_K_COEFF_PER_BIN ]; // 11 * 32 = 352
|
|
|
float k_coeff_intercalib_f0_sbm[ NB_BINS_COMPRESSED_SM_SBM_F0 * NB_K_COEFF_PER_BIN ]; // 22 * 32 = 704
|
|
|
|
|
|
//************
|
|
|
// RTEMS TASKS
|
|
|
|
|
|
rtems_task avf0_task( rtems_task_argument lfrRequestedMode )
|
|
|
{
|
|
|
int i;
|
|
|
|
|
|
rtems_event_set event_out;
|
|
|
rtems_status_code status;
|
|
|
rtems_id queue_id_prc0;
|
|
|
asm_msg msgForMATR;
|
|
|
ring_node *nodeForAveraging;
|
|
|
ring_node *ring_node_tab[8];
|
|
|
ring_node_asm *current_ring_node_asm_burst_sbm_f0;
|
|
|
ring_node_asm *current_ring_node_asm_norm_f0;
|
|
|
|
|
|
unsigned int nb_norm_bp1;
|
|
|
unsigned int nb_norm_bp2;
|
|
|
unsigned int nb_norm_asm;
|
|
|
unsigned int nb_sbm_bp1;
|
|
|
unsigned int nb_sbm_bp2;
|
|
|
|
|
|
nb_norm_bp1 = 0;
|
|
|
nb_norm_bp2 = 0;
|
|
|
nb_norm_asm = 0;
|
|
|
nb_sbm_bp1 = 0;
|
|
|
nb_sbm_bp2 = 0;
|
|
|
|
|
|
reset_nb_sm_f0( lfrRequestedMode ); // reset the sm counters that drive the BP and ASM computations / transmissions
|
|
|
ASM_generic_init_ring( asm_ring_norm_f0, NB_RING_NODES_ASM_NORM_F0 );
|
|
|
ASM_generic_init_ring( asm_ring_burst_sbm_f0, NB_RING_NODES_ASM_BURST_SBM_F0 );
|
|
|
current_ring_node_asm_norm_f0 = asm_ring_norm_f0;
|
|
|
current_ring_node_asm_burst_sbm_f0 = asm_ring_burst_sbm_f0;
|
|
|
|
|
|
BOOT_PRINTF1("in AVFO *** lfrRequestedMode = %d\n", (int) lfrRequestedMode)
|
|
|
|
|
|
status = get_message_queue_id_prc0( &queue_id_prc0 );
|
|
|
if (status != RTEMS_SUCCESSFUL)
|
|
|
{
|
|
|
PRINTF1("in MATR *** ERR get_message_queue_id_prc0 %d\n", status)
|
|
|
}
|
|
|
|
|
|
while(1){
|
|
|
rtems_event_receive(RTEMS_EVENT_0, RTEMS_WAIT, RTEMS_NO_TIMEOUT, &event_out); // wait for an RTEMS_EVENT0
|
|
|
|
|
|
//****************************************
|
|
|
// initialize the mesage for the MATR task
|
|
|
msgForMATR.norm = current_ring_node_asm_norm_f0;
|
|
|
msgForMATR.burst_sbm = current_ring_node_asm_burst_sbm_f0;
|
|
|
msgForMATR.event = 0x00; // this composite event will be sent to the PRC0 task
|
|
|
//
|
|
|
//****************************************
|
|
|
|
|
|
nodeForAveraging = getRingNodeForAveraging( 0 );
|
|
|
|
|
|
ring_node_tab[NB_SM_BEFORE_AVF0-1] = nodeForAveraging;
|
|
|
for ( i = 2; i < (NB_SM_BEFORE_AVF0+1); i++ )
|
|
|
{
|
|
|
nodeForAveraging = nodeForAveraging->previous;
|
|
|
ring_node_tab[NB_SM_BEFORE_AVF0-i] = nodeForAveraging;
|
|
|
}
|
|
|
|
|
|
// compute the average and store it in the averaged_sm_f1 buffer
|
|
|
SM_average( current_ring_node_asm_norm_f0->matrix,
|
|
|
current_ring_node_asm_burst_sbm_f0->matrix,
|
|
|
ring_node_tab,
|
|
|
nb_norm_bp1, nb_sbm_bp1,
|
|
|
&msgForMATR );
|
|
|
|
|
|
// update nb_average
|
|
|
nb_norm_bp1 = nb_norm_bp1 + NB_SM_BEFORE_AVF0;
|
|
|
nb_norm_bp2 = nb_norm_bp2 + NB_SM_BEFORE_AVF0;
|
|
|
nb_norm_asm = nb_norm_asm + NB_SM_BEFORE_AVF0;
|
|
|
nb_sbm_bp1 = nb_sbm_bp1 + NB_SM_BEFORE_AVF0;
|
|
|
nb_sbm_bp2 = nb_sbm_bp2 + NB_SM_BEFORE_AVF0;
|
|
|
|
|
|
if (nb_sbm_bp1 == nb_sm_before_f0.burst_sbm_bp1)
|
|
|
{
|
|
|
nb_sbm_bp1 = 0;
|
|
|
// set another ring for the ASM storage
|
|
|
current_ring_node_asm_burst_sbm_f0 = current_ring_node_asm_burst_sbm_f0->next;
|
|
|
if ( lfrCurrentMode == LFR_MODE_BURST )
|
|
|
{
|
|
|
msgForMATR.event = msgForMATR.event | RTEMS_EVENT_BURST_BP1_F0;
|
|
|
}
|
|
|
else if ( (lfrCurrentMode == LFR_MODE_SBM1) || (lfrCurrentMode == LFR_MODE_SBM2) )
|
|
|
{
|
|
|
msgForMATR.event = msgForMATR.event | RTEMS_EVENT_SBM_BP1_F0;
|
|
|
}
|
|
|
}
|
|
|
|
|
|
if (nb_sbm_bp2 == nb_sm_before_f0.burst_sbm_bp2)
|
|
|
{
|
|
|
nb_sbm_bp2 = 0;
|
|
|
if ( lfrCurrentMode == LFR_MODE_BURST )
|
|
|
{
|
|
|
msgForMATR.event = msgForMATR.event | RTEMS_EVENT_BURST_BP2_F0;
|
|
|
}
|
|
|
else if ( (lfrCurrentMode == LFR_MODE_SBM1) || (lfrCurrentMode == LFR_MODE_SBM2) )
|
|
|
{
|
|
|
msgForMATR.event = msgForMATR.event | RTEMS_EVENT_SBM_BP2_F0;
|
|
|
}
|
|
|
}
|
|
|
|
|
|
if (nb_norm_bp1 == nb_sm_before_f0.norm_bp1)
|
|
|
{
|
|
|
nb_norm_bp1 = 0;
|
|
|
// set another ring for the ASM storage
|
|
|
current_ring_node_asm_norm_f0 = current_ring_node_asm_norm_f0->next;
|
|
|
if ( (lfrCurrentMode == LFR_MODE_NORMAL)
|
|
|
|| (lfrCurrentMode == LFR_MODE_SBM1) || (lfrCurrentMode == LFR_MODE_SBM2) )
|
|
|
{
|
|
|
msgForMATR.event = msgForMATR.event | RTEMS_EVENT_NORM_BP1_F0;
|
|
|
}
|
|
|
}
|
|
|
|
|
|
if (nb_norm_bp2 == nb_sm_before_f0.norm_bp2)
|
|
|
{
|
|
|
nb_norm_bp2 = 0;
|
|
|
if ( (lfrCurrentMode == LFR_MODE_NORMAL)
|
|
|
|| (lfrCurrentMode == LFR_MODE_SBM1) || (lfrCurrentMode == LFR_MODE_SBM2) )
|
|
|
{
|
|
|
msgForMATR.event = msgForMATR.event | RTEMS_EVENT_NORM_BP2_F0;
|
|
|
}
|
|
|
}
|
|
|
|
|
|
if (nb_norm_asm == nb_sm_before_f0.norm_asm)
|
|
|
{
|
|
|
nb_norm_asm = 0;
|
|
|
if ( (lfrCurrentMode == LFR_MODE_NORMAL)
|
|
|
|| (lfrCurrentMode == LFR_MODE_SBM1) || (lfrCurrentMode == LFR_MODE_SBM2) )
|
|
|
{
|
|
|
msgForMATR.event = msgForMATR.event | RTEMS_EVENT_NORM_ASM_F0;
|
|
|
}
|
|
|
}
|
|
|
|
|
|
//*************************
|
|
|
// send the message to MATR
|
|
|
if (msgForMATR.event != 0x00)
|
|
|
{
|
|
|
status = rtems_message_queue_send( queue_id_prc0, (char *) &msgForMATR, MSG_QUEUE_SIZE_PRC0);
|
|
|
}
|
|
|
|
|
|
if (status != RTEMS_SUCCESSFUL) {
|
|
|
printf("in AVF0 *** Error sending message to MATR, code %d\n", status);
|
|
|
}
|
|
|
}
|
|
|
}
|
|
|
|
|
|
rtems_task prc0_task( rtems_task_argument lfrRequestedMode )
|
|
|
{
|
|
|
char incomingData[MSG_QUEUE_SIZE_SEND]; // incoming data buffer
|
|
|
size_t size; // size of the incoming TC packet
|
|
|
asm_msg *incomingMsg;
|
|
|
//
|
|
|
unsigned char sid;
|
|
|
rtems_status_code status;
|
|
|
rtems_id queue_id;
|
|
|
rtems_id queue_id_q_p0;
|
|
|
bp_packet_with_spare packet_norm_bp1;
|
|
|
bp_packet packet_norm_bp2;
|
|
|
bp_packet packet_sbm_bp1;
|
|
|
bp_packet packet_sbm_bp2;
|
|
|
ring_node *current_ring_node_to_send_asm_f0;
|
|
|
|
|
|
unsigned long long int localTime;
|
|
|
|
|
|
// init the ring of the averaged spectral matrices which will be transmitted to the DPU
|
|
|
init_ring( ring_to_send_asm_f0, NB_RING_NODES_ASM_F0, (volatile int*) buffer_asm_f0, TOTAL_SIZE_SM );
|
|
|
current_ring_node_to_send_asm_f0 = ring_to_send_asm_f0;
|
|
|
|
|
|
//*************
|
|
|
// NORM headers
|
|
|
BP_init_header_with_spare( &packet_norm_bp1.header,
|
|
|
APID_TM_SCIENCE_NORMAL_BURST, SID_NORM_BP1_F0,
|
|
|
PACKET_LENGTH_TM_LFR_SCIENCE_NORM_BP1_F0, NB_BINS_COMPRESSED_SM_F0 );
|
|
|
BP_init_header( &packet_norm_bp2,
|
|
|
APID_TM_SCIENCE_NORMAL_BURST, SID_NORM_BP2_F0,
|
|
|
PACKET_LENGTH_TM_LFR_SCIENCE_NORM_BP2_F0, NB_BINS_COMPRESSED_SM_F0);
|
|
|
|
|
|
//****************************
|
|
|
// BURST SBM1 and SBM2 headers
|
|
|
if ( lfrRequestedMode == LFR_MODE_BURST )
|
|
|
{
|
|
|
BP_init_header( &packet_sbm_bp1,
|
|
|
APID_TM_SCIENCE_NORMAL_BURST, SID_BURST_BP1_F0,
|
|
|
PACKET_LENGTH_TM_LFR_SCIENCE_SBM_BP1_F0, NB_BINS_COMPRESSED_SM_SBM_F0);
|
|
|
BP_init_header( &packet_sbm_bp2,
|
|
|
APID_TM_SCIENCE_NORMAL_BURST, SID_BURST_BP2_F0,
|
|
|
PACKET_LENGTH_TM_LFR_SCIENCE_SBM_BP2_F0, NB_BINS_COMPRESSED_SM_SBM_F0);
|
|
|
}
|
|
|
else if ( lfrRequestedMode == LFR_MODE_SBM1 )
|
|
|
{
|
|
|
BP_init_header( &packet_sbm_bp1,
|
|
|
APID_TM_SCIENCE_SBM1_SBM2, SID_SBM1_BP1_F0,
|
|
|
PACKET_LENGTH_TM_LFR_SCIENCE_SBM_BP1_F0, NB_BINS_COMPRESSED_SM_SBM_F0);
|
|
|
BP_init_header( &packet_sbm_bp2,
|
|
|
APID_TM_SCIENCE_SBM1_SBM2, SID_SBM1_BP2_F0,
|
|
|
PACKET_LENGTH_TM_LFR_SCIENCE_SBM_BP2_F0, NB_BINS_COMPRESSED_SM_SBM_F0);
|
|
|
}
|
|
|
else if ( lfrRequestedMode == LFR_MODE_SBM2 )
|
|
|
{
|
|
|
BP_init_header( &packet_sbm_bp1,
|
|
|
APID_TM_SCIENCE_SBM1_SBM2, SID_SBM2_BP1_F0,
|
|
|
PACKET_LENGTH_TM_LFR_SCIENCE_SBM_BP1_F0, NB_BINS_COMPRESSED_SM_SBM_F0);
|
|
|
BP_init_header( &packet_sbm_bp2,
|
|
|
APID_TM_SCIENCE_SBM1_SBM2, SID_SBM2_BP2_F0,
|
|
|
PACKET_LENGTH_TM_LFR_SCIENCE_SBM_BP2_F0, NB_BINS_COMPRESSED_SM_SBM_F0);
|
|
|
}
|
|
|
else
|
|
|
{
|
|
|
PRINTF1("in PRC0 *** lfrRequestedMode is %d, several headers not initialized\n", (unsigned int) lfrRequestedMode)
|
|
|
}
|
|
|
|
|
|
status = get_message_queue_id_send( &queue_id );
|
|
|
if (status != RTEMS_SUCCESSFUL)
|
|
|
{
|
|
|
PRINTF1("in PRC0 *** ERR get_message_queue_id_send %d\n", status)
|
|
|
}
|
|
|
status = get_message_queue_id_prc0( &queue_id_q_p0);
|
|
|
if (status != RTEMS_SUCCESSFUL)
|
|
|
{
|
|
|
PRINTF1("in PRC0 *** ERR get_message_queue_id_prc0 %d\n", status)
|
|
|
}
|
|
|
|
|
|
BOOT_PRINTF1("in PRC0 *** lfrRequestedMode = %d\n", (int) lfrRequestedMode)
|
|
|
|
|
|
while(1){
|
|
|
status = rtems_message_queue_receive( queue_id_q_p0, incomingData, &size, //************************************
|
|
|
RTEMS_WAIT, RTEMS_NO_TIMEOUT ); // wait for a message coming from AVF0
|
|
|
|
|
|
incomingMsg = (asm_msg*) incomingData;
|
|
|
|
|
|
localTime = getTimeAsUnsignedLongLongInt( );
|
|
|
|
|
|
//****************
|
|
|
//****************
|
|
|
// BURST SBM1 SBM2
|
|
|
//****************
|
|
|
//****************
|
|
|
if ( (incomingMsg->event & RTEMS_EVENT_BURST_BP1_F0 ) || (incomingMsg->event & RTEMS_EVENT_SBM_BP1_F0 ) )
|
|
|
{
|
|
|
sid = getSID( incomingMsg->event );
|
|
|
// 1) compress the matrix for Basic Parameters calculation
|
|
|
ASM_compress_reorganize_and_divide( incomingMsg->burst_sbm->matrix, compressed_sm_sbm_f0,
|
|
|
nb_sm_before_f0.burst_sbm_bp1,
|
|
|
NB_BINS_COMPRESSED_SM_SBM_F0, NB_BINS_TO_AVERAGE_ASM_SBM_F0,
|
|
|
ASM_F0_INDICE_START);
|
|
|
// 2) compute the BP1 set
|
|
|
BP1_set( compressed_sm_sbm_f0, k_coeff_intercalib_f0_sbm, NB_BINS_COMPRESSED_SM_SBM_F0, packet_sbm_bp1.data );
|
|
|
// 3) send the BP1 set
|
|
|
set_time( packet_sbm_bp1.time, (unsigned char *) &incomingMsg->coarseTimeSBM );
|
|
|
set_time( packet_sbm_bp1.acquisitionTime, (unsigned char *) &incomingMsg->coarseTimeSBM );
|
|
|
BP_send( (char *) &packet_sbm_bp1, queue_id,
|
|
|
PACKET_LENGTH_TM_LFR_SCIENCE_SBM_BP1_F0 + PACKET_LENGTH_DELTA,
|
|
|
sid);
|
|
|
// 4) compute the BP2 set if needed
|
|
|
if ( (incomingMsg->event & RTEMS_EVENT_BURST_BP2_F0) || (incomingMsg->event & RTEMS_EVENT_SBM_BP2_F0) )
|
|
|
{
|
|
|
// 1) compute the BP2 set
|
|
|
BP2_set( compressed_sm_sbm_f0, NB_BINS_COMPRESSED_SM_SBM_F0, packet_sbm_bp2.data );
|
|
|
// 2) send the BP2 set
|
|
|
set_time( packet_sbm_bp2.time, (unsigned char *) &incomingMsg->coarseTimeSBM );
|
|
|
set_time( packet_sbm_bp2.acquisitionTime, (unsigned char *) &incomingMsg->coarseTimeSBM );
|
|
|
BP_send( (char *) &packet_sbm_bp2, queue_id,
|
|
|
PACKET_LENGTH_TM_LFR_SCIENCE_SBM_BP2_F0 + PACKET_LENGTH_DELTA,
|
|
|
sid);
|
|
|
}
|
|
|
}
|
|
|
|
|
|
//*****
|
|
|
//*****
|
|
|
// NORM
|
|
|
//*****
|
|
|
//*****
|
|
|
if (incomingMsg->event & RTEMS_EVENT_NORM_BP1_F0)
|
|
|
{
|
|
|
// 1) compress the matrix for Basic Parameters calculation
|
|
|
ASM_compress_reorganize_and_divide( incomingMsg->norm->matrix, compressed_sm_norm_f0,
|
|
|
nb_sm_before_f0.norm_bp1,
|
|
|
NB_BINS_COMPRESSED_SM_F0, NB_BINS_TO_AVERAGE_ASM_F0,
|
|
|
ASM_F0_INDICE_START );
|
|
|
// 2) compute the BP1 set
|
|
|
BP1_set( compressed_sm_norm_f0, k_coeff_intercalib_f0_norm, NB_BINS_COMPRESSED_SM_F0, packet_norm_bp1.data );
|
|
|
// 3) send the BP1 set
|
|
|
set_time( packet_norm_bp1.header.time, (unsigned char *) &incomingMsg->coarseTimeNORM );
|
|
|
set_time( packet_norm_bp1.header.acquisitionTime, (unsigned char *) &incomingMsg->coarseTimeNORM );
|
|
|
BP_send( (char *) &packet_norm_bp1, queue_id,
|
|
|
PACKET_LENGTH_TM_LFR_SCIENCE_NORM_BP1_F0 + PACKET_LENGTH_DELTA,
|
|
|
SID_NORM_BP1_F0 );
|
|
|
if (incomingMsg->event & RTEMS_EVENT_NORM_BP2_F0)
|
|
|
{
|
|
|
// 1) compute the BP2 set using the same ASM as the one used for BP1
|
|
|
BP2_set( compressed_sm_norm_f0, NB_BINS_COMPRESSED_SM_F0, packet_norm_bp2.data );
|
|
|
// 2) send the BP2 set
|
|
|
set_time( packet_norm_bp2.time, (unsigned char *) &incomingMsg->coarseTimeNORM );
|
|
|
set_time( packet_norm_bp2.acquisitionTime, (unsigned char *) &incomingMsg->coarseTimeNORM );
|
|
|
BP_send( (char *) &packet_norm_bp2, queue_id,
|
|
|
PACKET_LENGTH_TM_LFR_SCIENCE_NORM_BP2_F0 + PACKET_LENGTH_DELTA,
|
|
|
SID_NORM_BP2_F0);
|
|
|
|
|
|
// < TMP DATA>
|
|
|
#define INDEX_COMPRESSED 1
|
|
|
unsigned int signif;
|
|
|
float significand;
|
|
|
unsigned int nbitexp = 6;
|
|
|
unsigned int nbitsig = 16 - nbitexp; // number of bits for the significand
|
|
|
unsigned int rangesig = (1 << nbitsig)-1; // == 2^nbitsig - 1
|
|
|
int expmax = 32;
|
|
|
int expmin = expmax - ((int) (1 << nbitexp)) + 1;
|
|
|
int exponent;
|
|
|
float auto_a0;
|
|
|
exponent = ( (int) ( (packet_norm_bp2.data[INDEX_COMPRESSED * NB_BYTES_PER_BP2] & 0xfc) >> 2) ) + expmin; // [1111 1100]
|
|
|
printf("exponent = %x, computed with exp = %x, expmin = %d\n",
|
|
|
exponent,
|
|
|
(packet_norm_bp2.data[INDEX_COMPRESSED * NB_BYTES_PER_BP2] & 0xfc) >> 2,
|
|
|
expmin);
|
|
|
signif = ( (packet_norm_bp2.data[INDEX_COMPRESSED * NB_BYTES_PER_BP2] & 0x3) << 8 ) + packet_norm_bp2.data[INDEX_COMPRESSED * NB_BYTES_PER_BP2+1];
|
|
|
significand = ( ( (float) signif ) / ( (float) rangesig) + 1) / 2;
|
|
|
auto_a0 = significand * pow(2,exponent);
|
|
|
printf("(BP2) [%d] compressed = %f *** AUTO A0 = %x, %x, exponent = %x, significand = %f ===> %f\n",
|
|
|
INDEX_COMPRESSED,
|
|
|
compressed_sm_norm_f0[INDEX_COMPRESSED * NB_VALUES_PER_SM],
|
|
|
packet_norm_bp2.data[ INDEX_COMPRESSED * NB_BYTES_PER_BP2],
|
|
|
packet_norm_bp2.data[ INDEX_COMPRESSED * NB_BYTES_PER_BP2 + 1],
|
|
|
exponent, significand, auto_a0 );
|
|
|
// printf("(BP2) 0 = %f, 1 = %f, 2 = %f, 3 = %f, 4 = %f, 5 = %f, 6 = %f, 7 = %f, 8 = %f, 9 = %f, 10 = %f,\n",
|
|
|
// compressed_sm_norm_f0[0 * NB_VALUES_PER_SM],
|
|
|
// compressed_sm_norm_f0[1 * NB_VALUES_PER_SM],
|
|
|
// compressed_sm_norm_f0[2 * NB_VALUES_PER_SM],
|
|
|
// compressed_sm_norm_f0[3 * NB_VALUES_PER_SM],
|
|
|
// compressed_sm_norm_f0[4 * NB_VALUES_PER_SM],
|
|
|
// compressed_sm_norm_f0[5 * NB_VALUES_PER_SM],
|
|
|
// compressed_sm_norm_f0[6 * NB_VALUES_PER_SM],
|
|
|
// compressed_sm_norm_f0[7 * NB_VALUES_PER_SM],
|
|
|
// compressed_sm_norm_f0[8 * NB_VALUES_PER_SM],
|
|
|
// compressed_sm_norm_f0[9 * NB_VALUES_PER_SM],
|
|
|
// compressed_sm_norm_f0[10 * NB_VALUES_PER_SM]);
|
|
|
// </TMP DATA>
|
|
|
|
|
|
}
|
|
|
}
|
|
|
|
|
|
if (incomingMsg->event & RTEMS_EVENT_NORM_ASM_F0)
|
|
|
{
|
|
|
// 1) reorganize the ASM and divide
|
|
|
ASM_reorganize_and_divide( incomingMsg->norm->matrix,
|
|
|
asm_f0_reorganized,
|
|
|
nb_sm_before_f0.norm_bp1 );
|
|
|
// 2) convert the float array in a char array
|
|
|
ASM_convert( asm_f0_reorganized, (char*) current_ring_node_to_send_asm_f0->buffer_address );
|
|
|
current_ring_node_to_send_asm_f0->coarseTime = incomingMsg->coarseTimeNORM;
|
|
|
current_ring_node_to_send_asm_f0->fineTime = incomingMsg->fineTimeNORM;
|
|
|
current_ring_node_to_send_asm_f0->sid = SID_NORM_ASM_F0;
|
|
|
|
|
|
// < TMP DATA>
|
|
|
#define INDEX_TO_LOOK_AT 31
|
|
|
float b11;
|
|
|
unsigned char *b11_charPtr;
|
|
|
b11_charPtr = (unsigned char*) &b11;
|
|
|
b11_charPtr[0] = ((unsigned char *) current_ring_node_to_send_asm_f0->buffer_address)[(INDEX_TO_LOOK_AT * NB_VALUES_PER_SM) * 2];
|
|
|
b11_charPtr[1] = ((unsigned char *) current_ring_node_to_send_asm_f0->buffer_address)[(INDEX_TO_LOOK_AT * NB_VALUES_PER_SM) * 2 +1];
|
|
|
b11_charPtr[2] = 0x00;
|
|
|
b11_charPtr[3] = 0x00;
|
|
|
printf("(ASM) initial = %f, reorganized and divided = %f, converted = %f\n",
|
|
|
incomingMsg->norm->matrix[INDEX_TO_LOOK_AT], // 32 * 96 = 3072 Hz
|
|
|
asm_f0_reorganized[ INDEX_TO_LOOK_AT * NB_VALUES_PER_SM ],
|
|
|
b11);
|
|
|
// </TMP DATA>
|
|
|
|
|
|
// 3) send the spectral matrix packets
|
|
|
status = rtems_message_queue_send( queue_id, ¤t_ring_node_to_send_asm_f0, sizeof( ring_node* ) );
|
|
|
// change asm ring node
|
|
|
current_ring_node_to_send_asm_f0 = current_ring_node_to_send_asm_f0->next;
|
|
|
}
|
|
|
|
|
|
}
|
|
|
}
|
|
|
|
|
|
//**********
|
|
|
// FUNCTIONS
|
|
|
|
|
|
void reset_nb_sm_f0( unsigned char lfrMode )
|
|
|
{
|
|
|
nb_sm_before_f0.norm_bp1 = parameter_dump_packet.sy_lfr_n_bp_p0 * 96;
|
|
|
nb_sm_before_f0.norm_bp2 = parameter_dump_packet.sy_lfr_n_bp_p1 * 96;
|
|
|
nb_sm_before_f0.norm_asm = (parameter_dump_packet.sy_lfr_n_asm_p[0] * 256 + parameter_dump_packet.sy_lfr_n_asm_p[1]) * 96;
|
|
|
nb_sm_before_f0.sbm1_bp1 = parameter_dump_packet.sy_lfr_s1_bp_p0 * 24; // 0.25 s per digit
|
|
|
nb_sm_before_f0.sbm1_bp2 = parameter_dump_packet.sy_lfr_s1_bp_p1 * 96;
|
|
|
nb_sm_before_f0.sbm2_bp1 = parameter_dump_packet.sy_lfr_s2_bp_p0 * 96;
|
|
|
nb_sm_before_f0.sbm2_bp2 = parameter_dump_packet.sy_lfr_s2_bp_p1 * 96;
|
|
|
nb_sm_before_f0.burst_bp1 = parameter_dump_packet.sy_lfr_b_bp_p0 * 96;
|
|
|
nb_sm_before_f0.burst_bp2 = parameter_dump_packet.sy_lfr_b_bp_p1 * 96;
|
|
|
|
|
|
if (lfrMode == LFR_MODE_SBM1)
|
|
|
{
|
|
|
nb_sm_before_f0.burst_sbm_bp1 = nb_sm_before_f0.sbm1_bp1;
|
|
|
nb_sm_before_f0.burst_sbm_bp2 = nb_sm_before_f0.sbm1_bp2;
|
|
|
}
|
|
|
else if (lfrMode == LFR_MODE_SBM2)
|
|
|
{
|
|
|
nb_sm_before_f0.burst_sbm_bp1 = nb_sm_before_f0.sbm2_bp1;
|
|
|
nb_sm_before_f0.burst_sbm_bp2 = nb_sm_before_f0.sbm2_bp2;
|
|
|
}
|
|
|
else if (lfrMode == LFR_MODE_BURST)
|
|
|
{
|
|
|
nb_sm_before_f0.burst_sbm_bp1 = nb_sm_before_f0.burst_bp1;
|
|
|
nb_sm_before_f0.burst_sbm_bp2 = nb_sm_before_f0.burst_bp2;
|
|
|
}
|
|
|
else
|
|
|
{
|
|
|
nb_sm_before_f0.burst_sbm_bp1 = nb_sm_before_f0.burst_bp1;
|
|
|
nb_sm_before_f0.burst_sbm_bp2 = nb_sm_before_f0.burst_bp2;
|
|
|
}
|
|
|
}
|
|
|
|
|
|
void init_k_coefficients_f0( void )
|
|
|
{
|
|
|
init_k_coefficients( k_coeff_intercalib_f0_norm, NB_BINS_COMPRESSED_SM_F0 );
|
|
|
init_k_coefficients( k_coeff_intercalib_f0_sbm, NB_BINS_COMPRESSED_SM_SBM_F0);
|
|
|
}
|
|
|
|
|
|
void test_TCH( void )
|
|
|
{
|
|
|
#define NB_BINS_COMPRESSED_MATRIX_TCH 1
|
|
|
|
|
|
unsigned char LFR_BP1_f0[NB_BINS_COMPRESSED_MATRIX_TCH*NB_BYTES_BP1];
|
|
|
unsigned char LFR_BP2_f0[NB_BINS_COMPRESSED_MATRIX_TCH*NB_BYTES_BP2];
|
|
|
float k_coefficients[NB_BINS_COMPRESSED_MATRIX_TCH * NB_K_COEFF_PER_BIN];
|
|
|
|
|
|
float compressed_spectral_matrix_TCH[ NB_BINS_COMPRESSED_MATRIX_TCH * NB_VALUES_PER_SPECTRAL_MATRIX ] = {
|
|
|
1.02217712e+06,
|
|
|
-8.58216250e+04,
|
|
|
-3.22199043e+04,
|
|
|
1.01597820e+05,
|
|
|
8.10333875e+05,
|
|
|
1.19030141e+05,
|
|
|
-8.69636688e+05,
|
|
|
5.01504031e+05,
|
|
|
-1.01948547e+05,
|
|
|
1.35475020e+04,
|
|
|
-3.67825469e+04,
|
|
|
-1.10950273e+05,
|
|
|
2.10715000e+04,
|
|
|
4.49727383e+04,
|
|
|
-4.37282031e+04,
|
|
|
3.83337695e+03,
|
|
|
1.05317175e+06,
|
|
|
-4.04155312e+05,
|
|
|
-1.32987891e+05,
|
|
|
1.49277250e+05,
|
|
|
-4.39122625e+05,
|
|
|
9.46006250e+05,
|
|
|
2.64386625e+05,
|
|
|
3.71843125e+05,
|
|
|
3.39770000e+05
|
|
|
};
|
|
|
|
|
|
init_k_coefficients( k_coefficients, NB_BINS_COMPRESSED_MATRIX_TCH );
|
|
|
|
|
|
printf("\n");
|
|
|
|
|
|
BP1_set(compressed_spectral_matrix_TCH, k_coefficients, NB_BINS_COMPRESSED_MATRIX_TCH, LFR_BP1_f0);
|
|
|
|
|
|
printf("\n");
|
|
|
|
|
|
BP2_set(compressed_spectral_matrix_TCH, NB_BINS_COMPRESSED_MATRIX_TCH, LFR_BP2_f0);
|
|
|
}
|
|
|
|
