/** 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 "fsw_processing.h"
#include "fsw_processing_globals.c"

unsigned int nb_sm_f0;
unsigned int nb_sm_f0_aux_f1;
unsigned int nb_sm_f1;
unsigned int nb_sm_f0_aux_f2;

//************************
// spectral matrices rings
ring_node_sm sm_ring_f0[ NB_RING_NODES_SM_F0 ];
ring_node_sm sm_ring_f1[ NB_RING_NODES_SM_F1 ];
ring_node_sm sm_ring_f2[ NB_RING_NODES_SM_F2 ];
ring_node_sm *current_ring_node_sm_f0;
ring_node_sm *current_ring_node_sm_f1;
ring_node_sm *current_ring_node_sm_f2;
ring_node_sm *ring_node_for_averaging_sm_f0;
ring_node_sm *ring_node_for_averaging_sm_f1;
ring_node_sm *ring_node_for_averaging_sm_f2;

//***********************************************************
// Interrupt Service Routine for spectral matrices processing

void spectral_matrices_isr_f0( void )
{
    unsigned char status;
    unsigned long long int time_0;
    unsigned long long int time_1;
    unsigned long long int syncBit0;
    unsigned long long int syncBit1;

    status = spectral_matrix_regs->status & 0x03;   // [0011] get the status_ready_matrix_f0_x bits

    time_0 = get_acquisition_time( (unsigned char *) &spectral_matrix_regs->f0_0_coarse_time );
    time_1 = get_acquisition_time( (unsigned char *) &spectral_matrix_regs->f0_1_coarse_time );
    syncBit0 = ( (unsigned long long int) (spectral_matrix_regs->f0_0_coarse_time & 0x80000000) ) << 16;
    syncBit1 = ( (unsigned long long int) (spectral_matrix_regs->f0_1_coarse_time & 0x80000000) ) << 16;

    switch(status)
    {
    case 0:
        break;
    case 3:
        if ( time_0 < time_1 )
        {
            close_matrix_actions( &nb_sm_f0, NB_SM_BEFORE_AVF0, Task_id[TASKID_AVF0],
                                  ring_node_for_averaging_sm_f0, current_ring_node_sm_f0, time_0 | syncBit0);
            current_ring_node_sm_f0 = current_ring_node_sm_f0->next;
            spectral_matrix_regs->f0_0_address = current_ring_node_sm_f0->buffer_address;
            close_matrix_actions( &nb_sm_f0, NB_SM_BEFORE_AVF0, Task_id[TASKID_AVF0],
                                  ring_node_for_averaging_sm_f0, current_ring_node_sm_f0, time_1 | syncBit1);
            current_ring_node_sm_f0 = current_ring_node_sm_f0->next;
            spectral_matrix_regs->f0_1_address = current_ring_node_sm_f0->buffer_address;
        }
        else
        {
            close_matrix_actions( &nb_sm_f0, NB_SM_BEFORE_AVF0, Task_id[TASKID_AVF0],
                                  ring_node_for_averaging_sm_f0, current_ring_node_sm_f0, time_1 | syncBit1);
            current_ring_node_sm_f0 = current_ring_node_sm_f0->next;
            spectral_matrix_regs->f0_1_address = current_ring_node_sm_f0->buffer_address;
            close_matrix_actions( &nb_sm_f0, NB_SM_BEFORE_AVF0, Task_id[TASKID_AVF0],
                                  ring_node_for_averaging_sm_f0, current_ring_node_sm_f0, time_0 | syncBit0);
            current_ring_node_sm_f0 = current_ring_node_sm_f0->next;
            spectral_matrix_regs->f0_0_address = current_ring_node_sm_f0->buffer_address;
        }
        spectral_matrix_regs->status = 0x03;   // [0011]
        break;
    case 1:
        close_matrix_actions( &nb_sm_f0, NB_SM_BEFORE_AVF0, Task_id[TASKID_AVF0],
                              ring_node_for_averaging_sm_f0, current_ring_node_sm_f0, time_0 | syncBit0);
        current_ring_node_sm_f0 = current_ring_node_sm_f0->next;
        spectral_matrix_regs->f0_0_address = current_ring_node_sm_f0->buffer_address;
        spectral_matrix_regs->status = 0x01;   // [0001]
        break;
    case 2:
        close_matrix_actions( &nb_sm_f0, NB_SM_BEFORE_AVF0, Task_id[TASKID_AVF0],
                              ring_node_for_averaging_sm_f0, current_ring_node_sm_f0, time_1 | syncBit1);
        current_ring_node_sm_f0 = current_ring_node_sm_f0->next;
        spectral_matrix_regs->f0_1_address = current_ring_node_sm_f0->buffer_address;
        spectral_matrix_regs->status = 0x02;   // [0010]
        break;
    }
}

void spectral_matrices_isr_f1( void )
{
    unsigned char status;
    unsigned long long int time;
    unsigned long long int syncBit;
    rtems_status_code status_code;

    status = (spectral_matrix_regs->status & 0x0c) >> 2;   // [1100] get the status_ready_matrix_f0_x bits

    switch(status)
    {
    case 0:
        break;
    case 3:
        // UNEXPECTED VALUE
        spectral_matrix_regs->status = 0xc0;   // [1100]
        status_code = rtems_event_send( Task_id[TASKID_DUMB], RTEMS_EVENT_11 );
        break;
    case 1:
        time = get_acquisition_time( (unsigned char *) &spectral_matrix_regs->f1_0_coarse_time );
        syncBit = ( (unsigned long long int) (spectral_matrix_regs->f1_0_coarse_time & 0x80000000) ) << 16;
        close_matrix_actions( &nb_sm_f1, NB_SM_BEFORE_AVF1, Task_id[TASKID_AVF1],
                              ring_node_for_averaging_sm_f1, current_ring_node_sm_f1, time | syncBit);
        current_ring_node_sm_f1 = current_ring_node_sm_f1->next;
        spectral_matrix_regs->f1_0_address = current_ring_node_sm_f1->buffer_address;
        spectral_matrix_regs->status = 0x04;   // [0100]
        break;
    case 2:
        time = get_acquisition_time( (unsigned char *) &spectral_matrix_regs->f1_1_coarse_time );
        syncBit = ( (unsigned long long int) (spectral_matrix_regs->f1_1_coarse_time & 0x80000000) ) << 16;
        close_matrix_actions( &nb_sm_f1, NB_SM_BEFORE_AVF1, Task_id[TASKID_AVF1],
                              ring_node_for_averaging_sm_f1, current_ring_node_sm_f1, time | syncBit);
        current_ring_node_sm_f1 = current_ring_node_sm_f1->next;
        spectral_matrix_regs->f1_1_address = current_ring_node_sm_f1->buffer_address;
        spectral_matrix_regs->status = 0x08;   // [1000]
        break;
    }
}

void spectral_matrices_isr_f2( void )
{
    unsigned char status;
    rtems_status_code status_code;

    status = (spectral_matrix_regs->status & 0x30) >> 4;   // [0011 0000] get the status_ready_matrix_f0_x bits

    ring_node_for_averaging_sm_f2 = current_ring_node_sm_f2;

    current_ring_node_sm_f2 = current_ring_node_sm_f2->next;

    switch(status)
    {
    case 0:
        break;
    case 3:
        // UNEXPECTED VALUE
        spectral_matrix_regs->status = 0x30;   // [0011 0000]
        status_code = rtems_event_send( Task_id[TASKID_DUMB], RTEMS_EVENT_11 );
        break;
    case 1:
        ring_node_for_averaging_sm_f2->coarseTime = spectral_matrix_regs->f2_0_coarse_time;
        ring_node_for_averaging_sm_f2->fineTime = spectral_matrix_regs->f2_0_fine_time;
        spectral_matrix_regs->f2_0_address = current_ring_node_sm_f2->buffer_address;
        spectral_matrix_regs->status = 0x10;   // [0001 0000]
        if (rtems_event_send( Task_id[TASKID_AVF2], RTEMS_EVENT_0 ) != RTEMS_SUCCESSFUL)
        {
            status_code = rtems_event_send( Task_id[TASKID_DUMB], RTEMS_EVENT_3 );
        }
        break;
    case 2:
        ring_node_for_averaging_sm_f2->coarseTime = spectral_matrix_regs->f2_1_coarse_time;
        ring_node_for_averaging_sm_f2->fineTime = spectral_matrix_regs->f2_1_fine_time;
        spectral_matrix_regs->f2_1_address = current_ring_node_sm_f2->buffer_address;
        spectral_matrix_regs->status = 0x20;   // [0010 0000]
        if (rtems_event_send( Task_id[TASKID_AVF2], RTEMS_EVENT_0 ) != RTEMS_SUCCESSFUL)
        {
            status_code = rtems_event_send( Task_id[TASKID_DUMB], RTEMS_EVENT_3 );
        }
        break;
    }
}

void spectral_matrix_isr_error_handler( void )
{
    rtems_status_code status_code;

    if (spectral_matrix_regs->status & 0x7c0)    // [0111 1100 0000]
    {
        status_code = rtems_event_send( Task_id[TASKID_DUMB], RTEMS_EVENT_8 );
    }
}

rtems_isr spectral_matrices_isr( rtems_vector_number vector )
{    
    // STATUS REGISTER
    // input_fifo_write(2) *** input_fifo_write(1) *** input_fifo_write(0)
    //           10                    9                       8
    // buffer_full ** bad_component_err ** f2_1 ** f2_0 ** f1_1 ** f1_0 ** f0_1 ** f0_0
    //      7                  6             5       4       3       2       1       0

    spectral_matrices_isr_f0();

    spectral_matrices_isr_f1();

    spectral_matrices_isr_f2();

//    spectral_matrix_isr_error_handler();
}

rtems_isr spectral_matrices_isr_simu( rtems_vector_number vector )
{
    rtems_status_code status_code;

    //***
    // F0
    nb_sm_f0 = nb_sm_f0 + 1;
    if (nb_sm_f0 == NB_SM_BEFORE_AVF0 )
    {
        ring_node_for_averaging_sm_f0 = current_ring_node_sm_f0;
        if (rtems_event_send( Task_id[TASKID_AVF0], RTEMS_EVENT_0 ) != RTEMS_SUCCESSFUL)
        {
            status_code = rtems_event_send( Task_id[TASKID_DUMB], RTEMS_EVENT_3 );
        }
        nb_sm_f0 = 0;
    }

    //***
    // F1
    nb_sm_f0_aux_f1 = nb_sm_f0_aux_f1 + 1;
    if (nb_sm_f0_aux_f1 == 6)
    {
        nb_sm_f0_aux_f1 = 0;
        nb_sm_f1 = nb_sm_f1 + 1;
    }
    if (nb_sm_f1 == NB_SM_BEFORE_AVF1 )
    {
        ring_node_for_averaging_sm_f1 = current_ring_node_sm_f1;
        if (rtems_event_send( Task_id[TASKID_AVF1], RTEMS_EVENT_0 ) != RTEMS_SUCCESSFUL)
        {
            status_code = rtems_event_send( Task_id[TASKID_DUMB], RTEMS_EVENT_3 );
        }
        nb_sm_f1 = 0;
    }

    //***
    // F2
    nb_sm_f0_aux_f2 = nb_sm_f0_aux_f2 + 1;
    if (nb_sm_f0_aux_f2 == 96)
    {
        nb_sm_f0_aux_f2 = 0;
        ring_node_for_averaging_sm_f2 = current_ring_node_sm_f2;
        if (rtems_event_send( Task_id[TASKID_AVF2], RTEMS_EVENT_0 ) != RTEMS_SUCCESSFUL)
        {
            status_code = rtems_event_send( Task_id[TASKID_DUMB], RTEMS_EVENT_3 );
        }
    }
}

//******************
// Spectral Matrices

void reset_nb_sm( void )
{
    nb_sm_f0 = 0;
    nb_sm_f0_aux_f1 = 0;
    nb_sm_f0_aux_f2 = 0;

    nb_sm_f1 = 0;
}

void SM_init_rings( void )
{
    unsigned char i;

    // F0 RING
    sm_ring_f0[0].next            = (ring_node_sm*) &sm_ring_f0[1];
    sm_ring_f0[0].previous        = (ring_node_sm*) &sm_ring_f0[NB_RING_NODES_SM_F0-1];
    sm_ring_f0[0].buffer_address  =
            (int) &sm_f0[ 0 ];

    sm_ring_f0[NB_RING_NODES_SM_F0-1].next           = (ring_node_sm*) &sm_ring_f0[0];
    sm_ring_f0[NB_RING_NODES_SM_F0-1].previous       = (ring_node_sm*) &sm_ring_f0[NB_RING_NODES_SM_F0-2];
    sm_ring_f0[NB_RING_NODES_SM_F0-1].buffer_address =
            (int) &sm_f0[ (NB_RING_NODES_SM_F0-1) * TOTAL_SIZE_SM ];

    for(i=1; i<NB_RING_NODES_SM_F0-1; i++)
    {
        sm_ring_f0[i].next            = (ring_node_sm*) &sm_ring_f0[i+1];
        sm_ring_f0[i].previous        = (ring_node_sm*) &sm_ring_f0[i-1];
        sm_ring_f0[i].buffer_address  =
                (int) &sm_f0[ i * TOTAL_SIZE_SM ];
    }

    // F1 RING
    sm_ring_f1[0].next            = (ring_node_sm*) &sm_ring_f1[1];
    sm_ring_f1[0].previous        = (ring_node_sm*) &sm_ring_f1[NB_RING_NODES_SM_F1-1];
    sm_ring_f1[0].buffer_address  =
            (int) &sm_f1[ 0 ];

    sm_ring_f1[NB_RING_NODES_SM_F1-1].next           = (ring_node_sm*) &sm_ring_f1[0];
    sm_ring_f1[NB_RING_NODES_SM_F1-1].previous       = (ring_node_sm*) &sm_ring_f1[NB_RING_NODES_SM_F1-2];
    sm_ring_f1[NB_RING_NODES_SM_F1-1].buffer_address =
            (int) &sm_f1[ (NB_RING_NODES_SM_F1-1) * TOTAL_SIZE_SM ];

    for(i=1; i<NB_RING_NODES_SM_F1-1; i++)
    {
        sm_ring_f1[i].next            = (ring_node_sm*) &sm_ring_f1[i+1];
        sm_ring_f1[i].previous        = (ring_node_sm*) &sm_ring_f1[i-1];
        sm_ring_f1[i].buffer_address  =
                (int) &sm_f1[ i * TOTAL_SIZE_SM ];
    }

    // F2 RING
    sm_ring_f2[0].next            = (ring_node_sm*) &sm_ring_f2[1];
    sm_ring_f2[0].previous        = (ring_node_sm*) &sm_ring_f2[NB_RING_NODES_SM_F2-1];
    sm_ring_f2[0].buffer_address  =
            (int) &sm_f2[ 0 ];

    sm_ring_f2[NB_RING_NODES_SM_F2-1].next           = (ring_node_sm*) &sm_ring_f2[0];
    sm_ring_f2[NB_RING_NODES_SM_F2-1].previous       = (ring_node_sm*) &sm_ring_f2[NB_RING_NODES_SM_F2-2];
    sm_ring_f2[NB_RING_NODES_SM_F2-1].buffer_address =
            (int) &sm_f2[ (NB_RING_NODES_SM_F2-1) * TOTAL_SIZE_SM ];

    for(i=1; i<NB_RING_NODES_SM_F2-1; i++)
    {
        sm_ring_f2[i].next            = (ring_node_sm*) &sm_ring_f2[i+1];
        sm_ring_f2[i].previous        = (ring_node_sm*) &sm_ring_f2[i-1];
        sm_ring_f2[i].buffer_address  =
                (int) &sm_f2[ i * TOTAL_SIZE_SM ];
    }

    DEBUG_PRINTF1("asm_ring_f0 @%x\n", (unsigned int) sm_ring_f0)
    DEBUG_PRINTF1("asm_ring_f1 @%x\n", (unsigned int) sm_ring_f1)
    DEBUG_PRINTF1("asm_ring_f2 @%x\n", (unsigned int) sm_ring_f2)

    spectral_matrix_regs->f0_0_address = sm_ring_f0[0].buffer_address;
    DEBUG_PRINTF1("spectral_matrix_regs->matrixF0_Address0 @%x\n", spectral_matrix_regs->f0_0_address)
}

void SM_generic_init_ring( ring_node_sm *ring, unsigned char nbNodes, volatile int sm_f[] )
{
    unsigned char i;

    //***************
    // BUFFER ADDRESS
    for(i=0; i<nbNodes; i++)
    {
        ring[ i ].buffer_address = (int) &sm_f[ i * TOTAL_SIZE_SM ];
    }

    //*****
    // NEXT
    ring[ nbNodes - 1 ].next = (ring_node_sm*) &ring[ 0 ];
    for(i=0; i<nbNodes-1; i++)
    {
        ring[ i ].next = (ring_node_sm*) &ring[ i + 1 ];
    }

    //*********
    // PREVIOUS
    ring[ 0 ].previous = (ring_node_sm*) &ring[ nbNodes -1 ];
    for(i=1; i<nbNodes; i++)
    {
        ring[ i ].previous = (ring_node_sm*) &ring[ i - 1 ];
    }
}

void ASM_generic_init_ring( ring_node_asm *ring, unsigned char nbNodes )
{
    unsigned char i;

    ring[ nbNodes - 1 ].next
            = (ring_node_asm*) &ring[ 0 ];

    for(i=0; i<nbNodes-1; i++)
    {
        ring[ i ].next            = (ring_node_asm*) &ring[ i + 1 ];
    }
}

void SM_reset_current_ring_nodes( void )
{
    current_ring_node_sm_f0 = sm_ring_f0[0].next;
    current_ring_node_sm_f1 = sm_ring_f1[0].next;
    current_ring_node_sm_f2 = sm_ring_f2[0].next;

    ring_node_for_averaging_sm_f0 = sm_ring_f0;
    ring_node_for_averaging_sm_f1 = sm_ring_f1;
    ring_node_for_averaging_sm_f2 = sm_ring_f2;
}

void ASM_init_header( Header_TM_LFR_SCIENCE_ASM_t *header)
{
    header->targetLogicalAddress = CCSDS_DESTINATION_ID;
    header->protocolIdentifier = CCSDS_PROTOCOLE_ID;
    header->reserved = 0x00;
    header->userApplication = CCSDS_USER_APP;
    header->packetID[0] = (unsigned char) (APID_TM_SCIENCE_NORMAL_BURST >> 8);
    header->packetID[1] = (unsigned char) (APID_TM_SCIENCE_NORMAL_BURST);
    header->packetSequenceControl[0] = 0xc0;
    header->packetSequenceControl[1] = 0x00;
    header->packetLength[0] = 0x00;
    header->packetLength[1] = 0x00;
    // DATA FIELD HEADER
    header->spare1_pusVersion_spare2 = 0x10;
    header->serviceType = TM_TYPE_LFR_SCIENCE; // service type
    header->serviceSubType = TM_SUBTYPE_LFR_SCIENCE; // service subtype
    header->destinationID = TM_DESTINATION_ID_GROUND;
    // AUXILIARY DATA HEADER
    header->sid = 0x00;
    header->biaStatusInfo = 0x00;
    header->pa_lfr_pkt_cnt_asm = 0x00;
    header->pa_lfr_pkt_nr_asm = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->pa_lfr_asm_blk_nr[0] = 0x00;  // BLK_NR MSB
    header->pa_lfr_asm_blk_nr[1] = 0x00;  // BLK_NR LSB
}

void ASM_send(Header_TM_LFR_SCIENCE_ASM_t *header, char *spectral_matrix,
                    unsigned int sid, spw_ioctl_pkt_send *spw_ioctl_send, rtems_id queue_id)
{
    unsigned int i;
    unsigned int length = 0;
    rtems_status_code status;

    for (i=0; i<2; i++)
    {
        // (1) BUILD THE DATA
        switch(sid)
        {
        case SID_NORM_ASM_F0:
            spw_ioctl_send->dlen = TOTAL_SIZE_ASM_F0_IN_BYTES / 2;  // 2 packets will be sent
            spw_ioctl_send->data = &spectral_matrix[
                    ( (ASM_F0_INDICE_START + (i*NB_BINS_PER_PKT_ASM_F0) ) * NB_VALUES_PER_SM ) * 2
                    ];
            length = PACKET_LENGTH_TM_LFR_SCIENCE_ASM_F0;
            header->pa_lfr_asm_blk_nr[0] = (unsigned char)  ( (NB_BINS_PER_PKT_ASM_F0) >> 8 ); // BLK_NR MSB
            header->pa_lfr_asm_blk_nr[1] = (unsigned char)    (NB_BINS_PER_PKT_ASM_F0);        // BLK_NR LSB
            break;
        case SID_NORM_ASM_F1:
            spw_ioctl_send->dlen = TOTAL_SIZE_ASM_F1_IN_BYTES / 2;  // 2 packets will be sent
            spw_ioctl_send->data = &spectral_matrix[
                    ( (ASM_F1_INDICE_START + (i*NB_BINS_PER_PKT_ASM_F1) ) * NB_VALUES_PER_SM ) * 2
                    ];
            length = PACKET_LENGTH_TM_LFR_SCIENCE_ASM_F1;
            header->pa_lfr_asm_blk_nr[0] = (unsigned char)  ( (NB_BINS_PER_PKT_ASM_F1) >> 8 ); // BLK_NR MSB
            header->pa_lfr_asm_blk_nr[1] = (unsigned char)    (NB_BINS_PER_PKT_ASM_F1);        // BLK_NR LSB
            break;
        case SID_NORM_ASM_F2:
            spw_ioctl_send->dlen = TOTAL_SIZE_ASM_F2_IN_BYTES / 2;  // 2 packets will be sent
            spw_ioctl_send->data = &spectral_matrix[
                    ( (ASM_F2_INDICE_START + (i*NB_BINS_PER_PKT_ASM_F2) ) * NB_VALUES_PER_SM ) * 2
                    ];
            length = PACKET_LENGTH_TM_LFR_SCIENCE_ASM_F2;
            header->pa_lfr_asm_blk_nr[0] = (unsigned char)  ( (NB_BINS_PER_PKT_ASM_F2) >> 8 ); // BLK_NR MSB
            header->pa_lfr_asm_blk_nr[1] = (unsigned char)    (NB_BINS_PER_PKT_ASM_F2);        // BLK_NR LSB
            break;
        default:
            PRINTF1("ERR *** in ASM_send *** unexpected sid %d\n", sid)
            break;
        }
        spw_ioctl_send->hlen = HEADER_LENGTH_TM_LFR_SCIENCE_ASM + CCSDS_PROTOCOLE_EXTRA_BYTES;
        spw_ioctl_send->hdr = (char *) header;
        spw_ioctl_send->options = 0;

        // (2) BUILD THE HEADER
        increment_seq_counter_source_id( header->packetSequenceControl, sid );
        header->packetLength[0] = (unsigned char) (length>>8);
        header->packetLength[1] = (unsigned char) (length);
        header->sid = (unsigned char) sid;   // SID
        header->pa_lfr_pkt_cnt_asm = 2;
        header->pa_lfr_pkt_nr_asm = (unsigned char) (i+1);

        // (3) SET PACKET TIME
        header->time[0] = (unsigned char) (time_management_regs->coarse_time>>24);
        header->time[1] = (unsigned char) (time_management_regs->coarse_time>>16);
        header->time[2] = (unsigned char) (time_management_regs->coarse_time>>8);
        header->time[3] = (unsigned char) (time_management_regs->coarse_time);
        header->time[4] = (unsigned char) (time_management_regs->fine_time>>8);
        header->time[5] = (unsigned char) (time_management_regs->fine_time);
        //
        header->acquisitionTime[0] = header->time[0];
        header->acquisitionTime[1] = header->time[1];
        header->acquisitionTime[2] = header->time[2];
        header->acquisitionTime[3] = header->time[3];
        header->acquisitionTime[4] = header->time[4];
        header->acquisitionTime[5] = header->time[5];

        // (4) SEND PACKET
        status =  rtems_message_queue_send( queue_id, spw_ioctl_send, ACTION_MSG_SPW_IOCTL_SEND_SIZE);
        if (status != RTEMS_SUCCESSFUL) {
            printf("in ASM_send *** ERR %d\n", (int) status);
        }
    }
}

//*****************
// Basic Parameters

void BP_init_header( Header_TM_LFR_SCIENCE_BP_t *header,
                     unsigned int apid, unsigned char sid,
                     unsigned int packetLength, unsigned char blkNr )
{
    header->targetLogicalAddress = CCSDS_DESTINATION_ID;
    header->protocolIdentifier = CCSDS_PROTOCOLE_ID;
    header->reserved = 0x00;
    header->userApplication = CCSDS_USER_APP;
    header->packetID[0] = (unsigned char) (apid >> 8);
    header->packetID[1] = (unsigned char) (apid);
    header->packetSequenceControl[0] = TM_PACKET_SEQ_CTRL_STANDALONE;
    header->packetSequenceControl[1] = 0x00;
    header->packetLength[0] = (unsigned char) (packetLength >> 8);
    header->packetLength[1] = (unsigned char) (packetLength);
    // DATA FIELD HEADER
    header->spare1_pusVersion_spare2 = 0x10;
    header->serviceType = TM_TYPE_LFR_SCIENCE; // service type
    header->serviceSubType = TM_SUBTYPE_LFR_SCIENCE; // service subtype
    header->destinationID = TM_DESTINATION_ID_GROUND;
    // AUXILIARY DATA HEADER
    header->sid = sid;
    header->biaStatusInfo = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->pa_lfr_bp_blk_nr[0] = 0x00;  // BLK_NR MSB
    header->pa_lfr_bp_blk_nr[1] = blkNr;  // BLK_NR LSB
}

void BP_init_header_with_spare(Header_TM_LFR_SCIENCE_BP_with_spare_t *header,
                                unsigned int apid, unsigned char sid,
                                unsigned int packetLength , unsigned char blkNr)
{
    header->targetLogicalAddress = CCSDS_DESTINATION_ID;
    header->protocolIdentifier = CCSDS_PROTOCOLE_ID;
    header->reserved = 0x00;
    header->userApplication = CCSDS_USER_APP;
    header->packetID[0] = (unsigned char) (apid >> 8);
    header->packetID[1] = (unsigned char) (apid);
    header->packetSequenceControl[0] = TM_PACKET_SEQ_CTRL_STANDALONE;
    header->packetSequenceControl[1] = 0x00;
    header->packetLength[0] = (unsigned char) (packetLength >> 8);
    header->packetLength[1] = (unsigned char) (packetLength);
    // DATA FIELD HEADER
    header->spare1_pusVersion_spare2 = 0x10;
    header->serviceType = TM_TYPE_LFR_SCIENCE; // service type
    header->serviceSubType = TM_SUBTYPE_LFR_SCIENCE; // service subtype
    header->destinationID = TM_DESTINATION_ID_GROUND;
    // AUXILIARY DATA HEADER
    header->sid = sid;
    header->biaStatusInfo = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->time[0] = 0x00;
    header->source_data_spare = 0x00;
    header->pa_lfr_bp_blk_nr[0] = 0x00;  // BLK_NR MSB
    header->pa_lfr_bp_blk_nr[1] = blkNr;  // BLK_NR LSB
}

void BP_send(char *data, rtems_id queue_id, unsigned int nbBytesToSend, unsigned int sid )
{
    rtems_status_code status;

    // SET THE SEQUENCE_CNT PARAMETER
    increment_seq_counter_source_id( (unsigned char*) &data[ PACKET_POS_SEQUENCE_CNT ], sid );
    // SEND PACKET
    status =  rtems_message_queue_send( queue_id, data, nbBytesToSend);
    if (status != RTEMS_SUCCESSFUL)
    {
        printf("ERR *** in BP_send *** ERR %d\n", (int) status);
    }
}

//******************
// general functions

void reset_spectral_matrix_regs( void )
{
    /** This function resets the spectral matrices module registers.
     *
     * The registers affected by this function are located at the following offset addresses:
     *
     * - 0x00 config
     * - 0x04 status
     * - 0x08 matrixF0_Address0
     * - 0x10 matrixFO_Address1
     * - 0x14 matrixF1_Address
     * - 0x18 matrixF2_Address
     *
     */

    spectral_matrix_regs->config = 0x00;
    spectral_matrix_regs->status = 0x00;

    spectral_matrix_regs->f0_0_address = current_ring_node_sm_f0->previous->buffer_address;
    spectral_matrix_regs->f0_1_address = current_ring_node_sm_f0->buffer_address;
    spectral_matrix_regs->f1_0_address = current_ring_node_sm_f1->previous->buffer_address;
    spectral_matrix_regs->f1_1_address = current_ring_node_sm_f1->buffer_address;
    spectral_matrix_regs->f2_0_address = current_ring_node_sm_f2->previous->buffer_address;
    spectral_matrix_regs->f2_1_address = current_ring_node_sm_f2->buffer_address;
}

void set_time( unsigned char *time, unsigned char * timeInBuffer )
{
    time[0] = timeInBuffer[0];
    time[1] = timeInBuffer[1];
    time[2] = timeInBuffer[2];
    time[3] = timeInBuffer[3];
    time[4] = timeInBuffer[6];
    time[5] = timeInBuffer[7];
}

unsigned long long int get_acquisition_time( unsigned char *timePtr )
{
    unsigned long long int acquisitionTimeAslong;
    acquisitionTimeAslong = 0x00;
    acquisitionTimeAslong = ( (unsigned long long int) (timePtr[0] & 0x7f) << 40 ) // [0111 1111] mask the synchronization bit
            + ( (unsigned long long int) timePtr[1] << 32 )
            + ( (unsigned long long int) timePtr[2] << 24 )
            + ( (unsigned long long int) timePtr[3] << 16 )
            + ( (unsigned long long int) timePtr[6] << 8  )
            + ( (unsigned long long int) timePtr[7]       );
    return acquisitionTimeAslong;
}

void close_matrix_actions(unsigned int *nb_sm, unsigned int nb_sm_before_avf, rtems_id avf_task_id,
                           ring_node_sm *node_for_averaging, ring_node_sm *ringNode,
                           unsigned long long int time )
{
    unsigned char *timePtr;
    unsigned char *coarseTimePtr;
    unsigned char *fineTimePtr;
    rtems_status_code status_code;

    timePtr = (unsigned char *) &time;
    coarseTimePtr = (unsigned char *) &node_for_averaging->coarseTime;
    fineTimePtr = (unsigned char *) &node_for_averaging->fineTime;

    *nb_sm = *nb_sm + 1;
    if (*nb_sm == nb_sm_before_avf)
    {
        node_for_averaging = ringNode;
        coarseTimePtr[0] = timePtr[2];
        coarseTimePtr[1] = timePtr[3];
        coarseTimePtr[2] = timePtr[4];
        coarseTimePtr[3] = timePtr[5];
        fineTimePtr[2] = timePtr[6];
        fineTimePtr[3] = timePtr[7];
        if (rtems_event_send( avf_task_id, RTEMS_EVENT_0 ) != RTEMS_SUCCESSFUL)
        {
            status_code = rtems_event_send( Task_id[TASKID_DUMB], RTEMS_EVENT_3 );
        }
        *nb_sm = 0;
    }
}

unsigned char getSID( rtems_event_set event )
{
    unsigned char sid;

    rtems_event_set eventSetBURST;
    rtems_event_set eventSetSBM;

    //******
    // BURST
    eventSetBURST = RTEMS_EVENT_BURST_BP1_F0
            | RTEMS_EVENT_BURST_BP1_F1
            | RTEMS_EVENT_BURST_BP2_F0
            | RTEMS_EVENT_BURST_BP2_F1;

    //****
    // SBM
    eventSetSBM = RTEMS_EVENT_SBM_BP1_F0
            | RTEMS_EVENT_SBM_BP1_F1
            | RTEMS_EVENT_SBM_BP2_F0
            | RTEMS_EVENT_SBM_BP2_F1;

    if (event & eventSetBURST)
    {
        sid = SID_BURST_BP1_F0;
    }
    else if (event & eventSetSBM)
    {
        sid = SID_SBM1_BP1_F0;
    }
    else
    {
        sid = 0;
    }

    return sid;
}