avf0_prc0.c
498 lines
| 21.3 KiB
| text/x-c
|
CLexer
paul
|
r126 | /** 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 | ||||
|
paul
|
r172 | 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 ]; | ||||
|
paul
|
r173 | int buffer_asm_f0 [ NB_RING_NODES_ASM_F0 * TOTAL_SIZE_SM ]; | ||
|
paul
|
r126 | |||
| 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 ]; | ||||
|
paul
|
r179 | |||
| 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 | ||||
|
paul
|
r126 | |||
| //************ | ||||
| // 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; | ||||
|
paul
|
r179 | ring_node *nodeForAveraging; | ||
|
paul
|
r173 | ring_node *ring_node_tab[8]; | ||
|
paul
|
r126 | 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 | ||||
|
paul
|
r150 | |||
| //**************************************** | ||||
| // 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; | ||||
|
paul
|
r173 | msgForMATR.event = 0x00; // this composite event will be sent to the PRC0 task | ||
|
paul
|
r150 | // | ||
| //**************************************** | ||||
|
paul
|
r179 | nodeForAveraging = getRingNodeForAveraging( 0 ); | ||
| ring_node_tab[NB_SM_BEFORE_AVF0-1] = nodeForAveraging; | ||||
|
paul
|
r126 | for ( i = 2; i < (NB_SM_BEFORE_AVF0+1); i++ ) | ||
| { | ||||
|
paul
|
r179 | nodeForAveraging = nodeForAveraging->previous; | ||
| ring_node_tab[NB_SM_BEFORE_AVF0-i] = nodeForAveraging; | ||||
|
paul
|
r126 | } | ||
| // 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, | ||||
|
paul
|
r179 | nb_norm_bp1, nb_sbm_bp1, | ||
| &msgForMATR ); | ||||
|
paul
|
r126 | |||
| // 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; | ||||
|
paul
|
r149 | if ( lfrCurrentMode == LFR_MODE_BURST ) | ||
|
paul
|
r126 | { | ||
|
paul
|
r149 | 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; | ||||
|
paul
|
r126 | } | ||
| } | ||||
| if (nb_sbm_bp2 == nb_sm_before_f0.burst_sbm_bp2) | ||||
| { | ||||
| nb_sbm_bp2 = 0; | ||||
|
paul
|
r149 | if ( lfrCurrentMode == LFR_MODE_BURST ) | ||
|
paul
|
r126 | { | ||
|
paul
|
r149 | 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; | ||||
|
paul
|
r126 | } | ||
| } | ||||
| 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; | ||||
| // | ||||
|
paul
|
r149 | unsigned char sid; | ||
|
paul
|
r126 | rtems_status_code status; | ||
| rtems_id queue_id; | ||||
| rtems_id queue_id_q_p0; | ||||
|
paul
|
r179 | bp_packet_with_spare packet_norm_bp1; | ||
| bp_packet packet_norm_bp2; | ||||
| bp_packet packet_sbm_bp1; | ||||
| bp_packet packet_sbm_bp2; | ||||
|
paul
|
r173 | ring_node *current_ring_node_to_send_asm_f0; | ||
|
paul
|
r126 | |||
| unsigned long long int localTime; | ||||
|
paul
|
r172 | // 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; | ||||
|
paul
|
r126 | //************* | ||
| // NORM headers | ||||
|
paul
|
r179 | BP_init_header_with_spare( &packet_norm_bp1.header, | ||
|
paul
|
r126 | APID_TM_SCIENCE_NORMAL_BURST, SID_NORM_BP1_F0, | ||
| PACKET_LENGTH_TM_LFR_SCIENCE_NORM_BP1_F0, NB_BINS_COMPRESSED_SM_F0 ); | ||||
|
paul
|
r179 | BP_init_header( &packet_norm_bp2, | ||
|
paul
|
r126 | 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 ) | ||||
| { | ||||
|
paul
|
r179 | BP_init_header( &packet_sbm_bp1, | ||
|
paul
|
r126 | APID_TM_SCIENCE_NORMAL_BURST, SID_BURST_BP1_F0, | ||
| PACKET_LENGTH_TM_LFR_SCIENCE_SBM_BP1_F0, NB_BINS_COMPRESSED_SM_SBM_F0); | ||||
|
paul
|
r179 | BP_init_header( &packet_sbm_bp2, | ||
|
paul
|
r126 | 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 ) | ||||
| { | ||||
|
paul
|
r179 | BP_init_header( &packet_sbm_bp1, | ||
|
paul
|
r126 | APID_TM_SCIENCE_SBM1_SBM2, SID_SBM1_BP1_F0, | ||
| PACKET_LENGTH_TM_LFR_SCIENCE_SBM_BP1_F0, NB_BINS_COMPRESSED_SM_SBM_F0); | ||||
|
paul
|
r179 | BP_init_header( &packet_sbm_bp2, | ||
|
paul
|
r126 | 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 ) | ||||
| { | ||||
|
paul
|
r179 | BP_init_header( &packet_sbm_bp1, | ||
|
paul
|
r126 | APID_TM_SCIENCE_SBM1_SBM2, SID_SBM2_BP1_F0, | ||
| PACKET_LENGTH_TM_LFR_SCIENCE_SBM_BP1_F0, NB_BINS_COMPRESSED_SM_SBM_F0); | ||||
|
paul
|
r179 | BP_init_header( &packet_sbm_bp2, | ||
|
paul
|
r126 | 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( ); | ||||
|
paul
|
r135 | |||
|
paul
|
r126 | //**************** | ||
| //**************** | ||||
| // BURST SBM1 SBM2 | ||||
| //**************** | ||||
| //**************** | ||||
|
paul
|
r149 | if ( (incomingMsg->event & RTEMS_EVENT_BURST_BP1_F0 ) || (incomingMsg->event & RTEMS_EVENT_SBM_BP1_F0 ) ) | ||
|
paul
|
r126 | { | ||
|
paul
|
r149 | sid = getSID( incomingMsg->event ); | ||
|
paul
|
r126 | // 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 | ||||
|
paul
|
r179 | BP1_set( compressed_sm_sbm_f0, k_coeff_intercalib_f0_sbm, NB_BINS_COMPRESSED_SM_SBM_F0, packet_sbm_bp1.data ); | ||
|
paul
|
r126 | // 3) send the BP1 set | ||
|
paul
|
r179 | 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, | ||||
|
paul
|
r149 | PACKET_LENGTH_TM_LFR_SCIENCE_SBM_BP1_F0 + PACKET_LENGTH_DELTA, | ||
| sid); | ||||
|
paul
|
r126 | // 4) compute the BP2 set if needed | ||
|
paul
|
r149 | if ( (incomingMsg->event & RTEMS_EVENT_BURST_BP2_F0) || (incomingMsg->event & RTEMS_EVENT_SBM_BP2_F0) ) | ||
|
paul
|
r126 | { | ||
| // 1) compute the BP2 set | ||||
|
paul
|
r179 | BP2_set( compressed_sm_sbm_f0, NB_BINS_COMPRESSED_SM_SBM_F0, packet_sbm_bp2.data ); | ||
|
paul
|
r126 | // 2) send the BP2 set | ||
|
paul
|
r179 | 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, | ||||
|
paul
|
r149 | PACKET_LENGTH_TM_LFR_SCIENCE_SBM_BP2_F0 + PACKET_LENGTH_DELTA, | ||
| sid); | ||||
|
paul
|
r126 | } | ||
| } | ||||
| //***** | ||||
| //***** | ||||
| // 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 | ||||
|
paul
|
r179 | BP1_set( compressed_sm_norm_f0, k_coeff_intercalib_f0_norm, NB_BINS_COMPRESSED_SM_F0, packet_norm_bp1.data ); | ||
|
paul
|
r126 | // 3) send the BP1 set | ||
|
paul
|
r179 | 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, | ||||
|
paul
|
r133 | PACKET_LENGTH_TM_LFR_SCIENCE_NORM_BP1_F0 + PACKET_LENGTH_DELTA, | ||
| SID_NORM_BP1_F0 ); | ||||
|
paul
|
r126 | if (incomingMsg->event & RTEMS_EVENT_NORM_BP2_F0) | ||
| { | ||||
| // 1) compute the BP2 set using the same ASM as the one used for BP1 | ||||
|
paul
|
r179 | BP2_set( compressed_sm_norm_f0, NB_BINS_COMPRESSED_SM_F0, packet_norm_bp2.data ); | ||
|
paul
|
r126 | // 2) send the BP2 set | ||
|
paul
|
r179 | 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, | ||||
|
paul
|
r133 | PACKET_LENGTH_TM_LFR_SCIENCE_NORM_BP2_F0 + PACKET_LENGTH_DELTA, | ||
| SID_NORM_BP2_F0); | ||||
|
paul
|
r179 | |||
| // < 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> | ||||
|
paul
|
r126 | } | ||
| } | ||||
| if (incomingMsg->event & RTEMS_EVENT_NORM_ASM_F0) | ||||
| { | ||||
| // 1) reorganize the ASM and divide | ||||
| ASM_reorganize_and_divide( incomingMsg->norm->matrix, | ||||
|
paul
|
r173 | asm_f0_reorganized, | ||
| nb_sm_before_f0.norm_bp1 ); | ||||
|
paul
|
r126 | // 2) convert the float array in a char array | ||
|
paul
|
r173 | ASM_convert( asm_f0_reorganized, (char*) current_ring_node_to_send_asm_f0->buffer_address ); | ||
|
paul
|
r179 | current_ring_node_to_send_asm_f0->coarseTime = incomingMsg->coarseTimeNORM; | ||
| current_ring_node_to_send_asm_f0->fineTime = incomingMsg->fineTimeNORM; | ||||
|
paul
|
r173 | current_ring_node_to_send_asm_f0->sid = SID_NORM_ASM_F0; | ||
|
paul
|
r179 | |||
| // < 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> | ||||
|
paul
|
r126 | // 3) send the spectral matrix packets | ||
|
paul
|
r173 | 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; | ||||
|
paul
|
r126 | } | ||
| } | ||||
| } | ||||
| //********** | ||||
| // 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; | ||||
|
paul
|
r150 | nb_sm_before_f0.sbm1_bp1 = parameter_dump_packet.sy_lfr_s1_bp_p0 * 24; // 0.25 s per digit | ||
|
paul
|
r126 | 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; | ||||
| } | ||||
| } | ||||
|
paul
|
r179 | |||
| 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); | ||||
| } | ||||