fsw_processing.h
332 lines
| 12.7 KiB
| text/x-c
|
CLexer
paul
|
r126 | #ifndef FSW_PROCESSING_H_INCLUDED | ||
#define FSW_PROCESSING_H_INCLUDED | ||||
#include <rtems.h> | ||||
#include <grspw.h> | ||||
#include <math.h> | ||||
#include <stdlib.h> // abs() is in the stdlib | ||||
paul
|
r227 | #include <stdio.h> | ||
paul
|
r126 | #include <math.h> | ||
paul
|
r175 | #include <grlib_regs.h> | ||
paul
|
r126 | |||
#include "fsw_params.h" | ||||
typedef struct ring_node_asm | ||||
{ | ||||
struct ring_node_asm *next; | ||||
float matrix[ TOTAL_SIZE_SM ]; | ||||
unsigned int status; | ||||
} ring_node_asm; | ||||
paul
|
r133 | typedef struct | ||
paul
|
r126 | { | ||
paul
|
r179 | unsigned char targetLogicalAddress; | ||
unsigned char protocolIdentifier; | ||||
unsigned char reserved; | ||||
unsigned char userApplication; | ||||
unsigned char packetID[2]; | ||||
unsigned char packetSequenceControl[2]; | ||||
unsigned char packetLength[2]; | ||||
// DATA FIELD HEADER | ||||
unsigned char spare1_pusVersion_spare2; | ||||
unsigned char serviceType; | ||||
unsigned char serviceSubType; | ||||
unsigned char destinationID; | ||||
unsigned char time[6]; | ||||
// AUXILIARY HEADER | ||||
unsigned char sid; | ||||
unsigned char biaStatusInfo; | ||||
paul
|
r192 | unsigned char sy_lfr_common_parameters_spare; | ||
unsigned char sy_lfr_common_parameters; | ||||
paul
|
r179 | unsigned char acquisitionTime[6]; | ||
unsigned char pa_lfr_bp_blk_nr[2]; | ||||
// SOURCE DATA | ||||
paul
|
r181 | unsigned char data[ 780 ]; // MAX size is 26 bins * 30 Bytes [TM_LFR_SCIENCE_BURST_BP2_F1] | ||
paul
|
r126 | } bp_packet; | ||
paul
|
r133 | typedef struct | ||
paul
|
r126 | { | ||
paul
|
r181 | unsigned char targetLogicalAddress; | ||
unsigned char protocolIdentifier; | ||||
unsigned char reserved; | ||||
unsigned char userApplication; | ||||
unsigned char packetID[2]; | ||||
unsigned char packetSequenceControl[2]; | ||||
unsigned char packetLength[2]; | ||||
// DATA FIELD HEADER | ||||
unsigned char spare1_pusVersion_spare2; | ||||
unsigned char serviceType; | ||||
unsigned char serviceSubType; | ||||
unsigned char destinationID; | ||||
unsigned char time[6]; | ||||
// AUXILIARY HEADER | ||||
unsigned char sid; | ||||
unsigned char biaStatusInfo; | ||||
paul
|
r192 | unsigned char sy_lfr_common_parameters_spare; | ||
unsigned char sy_lfr_common_parameters; | ||||
paul
|
r181 | unsigned char acquisitionTime[6]; | ||
unsigned char source_data_spare; | ||||
unsigned char pa_lfr_bp_blk_nr[2]; | ||||
// SOURCE DATA | ||||
paul
|
r216 | unsigned char data[ 143 ]; // 13 bins * 11 Bytes | ||
} bp_packet_with_spare; // only for TM_LFR_SCIENCE_NORMAL_BP1_F0 and F1 | ||||
paul
|
r126 | |||
paul
|
r197 | typedef struct asm_msg | ||
paul
|
r126 | { | ||
ring_node_asm *norm; | ||||
ring_node_asm *burst_sbm; | ||||
rtems_event_set event; | ||||
paul
|
r179 | unsigned int coarseTimeNORM; | ||
unsigned int fineTimeNORM; | ||||
unsigned int coarseTimeSBM; | ||||
unsigned int fineTimeSBM; | ||||
paul
|
r126 | } asm_msg; | ||
paul
|
r259 | extern unsigned char thisIsAnASMRestart; | ||
paul
|
r126 | extern volatile int sm_f0[ ]; | ||
extern volatile int sm_f1[ ]; | ||||
extern volatile int sm_f2[ ]; | ||||
// parameters | ||||
extern struct param_local_str param_local; | ||||
paul
|
r195 | extern Packet_TM_LFR_PARAMETER_DUMP_t parameter_dump_packet; | ||
paul
|
r126 | |||
// registers | ||||
extern time_management_regs_t *time_management_regs; | ||||
paul
|
r179 | extern volatile spectral_matrix_regs_t *spectral_matrix_regs; | ||
paul
|
r126 | |||
extern rtems_name misc_name[5]; | ||||
extern rtems_id Task_id[20]; /* array of task ids */ | ||||
paul
|
r179 | ring_node * getRingNodeForAveraging( unsigned char frequencyChannel); | ||
paul
|
r126 | // ISR | ||
rtems_isr spectral_matrices_isr( rtems_vector_number vector ); | ||||
//****************** | ||||
// Spectral Matrices | ||||
void reset_nb_sm( void ); | ||||
// SM | ||||
void SM_init_rings( void ); | ||||
void SM_reset_current_ring_nodes( void ); | ||||
// ASM | ||||
void ASM_generic_init_ring(ring_node_asm *ring, unsigned char nbNodes ); | ||||
//***************** | ||||
// Basic Parameters | ||||
void BP_reset_current_ring_nodes( void ); | ||||
paul
|
r181 | void BP_init_header(bp_packet *packet, | ||
paul
|
r243 | unsigned int apid, unsigned char sid, | ||
unsigned int packetLength , unsigned char blkNr); | ||||
paul
|
r181 | void BP_init_header_with_spare(bp_packet_with_spare *packet, | ||
paul
|
r243 | unsigned int apid, unsigned char sid, | ||
unsigned int packetLength, unsigned char blkNr ); | ||||
paul
|
r133 | void BP_send( char *data, | ||
paul
|
r243 | rtems_id queue_id, | ||
paul
|
r133 | unsigned int nbBytesToSend , unsigned int sid ); | ||
paul
|
r243 | void BP_send_s1_s2(char *data, | ||
rtems_id queue_id, | ||||
unsigned int nbBytesToSend, unsigned int sid ); | ||||
paul
|
r126 | |||
//****************** | ||||
// general functions | ||||
paul
|
r171 | void reset_sm_status( void ); | ||
paul
|
r126 | void reset_spectral_matrix_regs( void ); | ||
void set_time(unsigned char *time, unsigned char *timeInBuffer ); | ||||
paul
|
r139 | unsigned long long int get_acquisition_time( unsigned char *timePtr ); | ||
paul
|
r149 | unsigned char getSID( rtems_event_set event ); | ||
paul
|
r126 | |||
extern rtems_status_code get_message_queue_id_prc1( rtems_id *queue_id ); | ||||
extern rtems_status_code get_message_queue_id_prc2( rtems_id *queue_id ); | ||||
//*************************************** | ||||
// DEFINITIONS OF STATIC INLINE FUNCTIONS | ||||
paul
|
r172 | static inline void SM_average(float *averaged_spec_mat_NORM, float *averaged_spec_mat_SBM, | ||
paul
|
r179 | ring_node *ring_node_tab[], | ||
unsigned int nbAverageNORM, unsigned int nbAverageSBM, | ||||
asm_msg *msgForMATR ); | ||||
paul
|
r195 | |||
paul
|
r185 | static inline void SM_average_debug(float *averaged_spec_mat_NORM, float *averaged_spec_mat_SBM, | ||
paul
|
r243 | ring_node *ring_node_tab[], | ||
unsigned int nbAverageNORM, unsigned int nbAverageSBM, | ||||
asm_msg *msgForMATR ); | ||||
paul
|
r182 | |||
void ASM_patch( float *inputASM, float *outputASM ); | ||||
paul
|
r195 | |||
paul
|
r182 | void extractReImVectors(float *inputASM, float *outputASM, unsigned int asmComponent ); | ||
paul
|
r126 | static inline void ASM_reorganize_and_divide(float *averaged_spec_mat, float *averaged_spec_mat_reorganized, | ||
paul
|
r243 | float divider ); | ||
paul
|
r195 | |||
paul
|
r126 | static inline void ASM_compress_reorganize_and_divide(float *averaged_spec_mat, float *compressed_spec_mat, | ||
paul
|
r243 | float divider, | ||
unsigned char nbBinsCompressedMatrix, unsigned char nbBinsToAverage , unsigned char ASMIndexStart); | ||||
paul
|
r195 | |||
paul
|
r126 | static inline void ASM_convert(volatile float *input_matrix, char *output_matrix); | ||
paul
|
r173 | void SM_average( float *averaged_spec_mat_NORM, float *averaged_spec_mat_SBM, | ||
paul
|
r179 | ring_node *ring_node_tab[], | ||
unsigned int nbAverageNORM, unsigned int nbAverageSBM, | ||||
asm_msg *msgForMATR ) | ||||
paul
|
r126 | { | ||
float sum; | ||||
unsigned int i; | ||||
for(i=0; i<TOTAL_SIZE_SM; i++) | ||||
{ | ||||
sum = ( (int *) (ring_node_tab[0]->buffer_address) ) [ i ] | ||||
+ ( (int *) (ring_node_tab[1]->buffer_address) ) [ i ] | ||||
+ ( (int *) (ring_node_tab[2]->buffer_address) ) [ i ] | ||||
+ ( (int *) (ring_node_tab[3]->buffer_address) ) [ i ] | ||||
+ ( (int *) (ring_node_tab[4]->buffer_address) ) [ i ] | ||||
+ ( (int *) (ring_node_tab[5]->buffer_address) ) [ i ] | ||||
+ ( (int *) (ring_node_tab[6]->buffer_address) ) [ i ] | ||||
+ ( (int *) (ring_node_tab[7]->buffer_address) ) [ i ]; | ||||
paul
|
r150 | if ( (nbAverageNORM == 0) && (nbAverageSBM == 0) ) | ||
paul
|
r126 | { | ||
paul
|
r150 | averaged_spec_mat_NORM[ i ] = sum; | ||
averaged_spec_mat_SBM[ i ] = sum; | ||||
paul
|
r179 | msgForMATR->coarseTimeNORM = ring_node_tab[0]->coarseTime; | ||
msgForMATR->fineTimeNORM = ring_node_tab[0]->fineTime; | ||||
msgForMATR->coarseTimeSBM = ring_node_tab[0]->coarseTime; | ||||
msgForMATR->fineTimeSBM = ring_node_tab[0]->fineTime; | ||||
paul
|
r126 | } | ||
paul
|
r150 | else if ( (nbAverageNORM != 0) && (nbAverageSBM != 0) ) | ||
paul
|
r126 | { | ||
paul
|
r150 | averaged_spec_mat_NORM[ i ] = ( averaged_spec_mat_NORM[ i ] + sum ); | ||
averaged_spec_mat_SBM[ i ] = ( averaged_spec_mat_SBM[ i ] + sum ); | ||||
paul
|
r126 | } | ||
paul
|
r150 | else if ( (nbAverageNORM != 0) && (nbAverageSBM == 0) ) | ||
paul
|
r126 | { | ||
paul
|
r150 | averaged_spec_mat_NORM[ i ] = ( averaged_spec_mat_NORM[ i ] + sum ); | ||
averaged_spec_mat_SBM[ i ] = sum; | ||||
paul
|
r179 | msgForMATR->coarseTimeSBM = ring_node_tab[0]->coarseTime; | ||
msgForMATR->fineTimeSBM = ring_node_tab[0]->fineTime; | ||||
paul
|
r126 | } | ||
else | ||||
{ | ||||
paul
|
r223 | averaged_spec_mat_NORM[ i ] = sum; | ||
averaged_spec_mat_SBM[ i ] = ( averaged_spec_mat_SBM[ i ] + sum ); | ||||
paul
|
r231 | msgForMATR->coarseTimeNORM = ring_node_tab[0]->coarseTime; | ||
msgForMATR->fineTimeNORM = ring_node_tab[0]->fineTime; | ||||
paul
|
r243 | // PRINTF2("ERR *** in SM_average *** unexpected parameters %d %d\n", nbAverageNORM, nbAverageSBM) | ||
paul
|
r126 | } | ||
} | ||||
} | ||||
paul
|
r173 | void SM_average_debug( float *averaged_spec_mat_NORM, float *averaged_spec_mat_SBM, | ||
paul
|
r243 | ring_node *ring_node_tab[], | ||
unsigned int nbAverageNORM, unsigned int nbAverageSBM, | ||||
asm_msg *msgForMATR ) | ||||
paul
|
r172 | { | ||
float sum; | ||||
unsigned int i; | ||||
for(i=0; i<TOTAL_SIZE_SM; i++) | ||||
{ | ||||
sum = ( (int *) (ring_node_tab[0]->buffer_address) ) [ i ]; | ||||
paul
|
r185 | averaged_spec_mat_NORM[ i ] = sum; | ||
averaged_spec_mat_SBM[ i ] = sum; | ||||
msgForMATR->coarseTimeNORM = ring_node_tab[0]->coarseTime; | ||||
msgForMATR->fineTimeNORM = ring_node_tab[0]->fineTime; | ||||
msgForMATR->coarseTimeSBM = ring_node_tab[0]->coarseTime; | ||||
msgForMATR->fineTimeSBM = ring_node_tab[0]->fineTime; | ||||
paul
|
r172 | } | ||
} | ||||
paul
|
r126 | void ASM_reorganize_and_divide( float *averaged_spec_mat, float *averaged_spec_mat_reorganized, float divider ) | ||
{ | ||||
int frequencyBin; | ||||
int asmComponent; | ||||
paul
|
r179 | unsigned int offsetASM; | ||
unsigned int offsetASMReorganized; | ||||
paul
|
r126 | |||
paul
|
r179 | // BUILD DATA | ||
paul
|
r126 | for (asmComponent = 0; asmComponent < NB_VALUES_PER_SM; asmComponent++) | ||
{ | ||||
for( frequencyBin = 0; frequencyBin < NB_BINS_PER_SM; frequencyBin++ ) | ||||
{ | ||||
paul
|
r179 | offsetASMReorganized = | ||
paul
|
r126 | frequencyBin * NB_VALUES_PER_SM | ||
+ asmComponent; | ||||
paul
|
r179 | offsetASM = | ||
paul
|
r126 | asmComponent * NB_BINS_PER_SM | ||
+ frequencyBin; | ||||
paul
|
r179 | averaged_spec_mat_reorganized[offsetASMReorganized ] = | ||
averaged_spec_mat[ offsetASM ] / divider; | ||||
paul
|
r126 | } | ||
} | ||||
} | ||||
void ASM_compress_reorganize_and_divide(float *averaged_spec_mat, float *compressed_spec_mat , float divider, | ||||
paul
|
r243 | unsigned char nbBinsCompressedMatrix, unsigned char nbBinsToAverage, unsigned char ASMIndexStart ) | ||
paul
|
r126 | { | ||
int frequencyBin; | ||||
int asmComponent; | ||||
int offsetASM; | ||||
int offsetCompressed; | ||||
int k; | ||||
paul
|
r179 | // BUILD DATA | ||
paul
|
r126 | for (asmComponent = 0; asmComponent < NB_VALUES_PER_SM; asmComponent++) | ||
{ | ||||
for( frequencyBin = 0; frequencyBin < nbBinsCompressedMatrix; frequencyBin++ ) | ||||
{ | ||||
offsetCompressed = // NO TIME OFFSET | ||||
frequencyBin * NB_VALUES_PER_SM | ||||
+ asmComponent; | ||||
offsetASM = // NO TIME OFFSET | ||||
asmComponent * NB_BINS_PER_SM | ||||
+ ASMIndexStart | ||||
+ frequencyBin * nbBinsToAverage; | ||||
compressed_spec_mat[ offsetCompressed ] = 0; | ||||
for ( k = 0; k < nbBinsToAverage; k++ ) | ||||
{ | ||||
compressed_spec_mat[offsetCompressed ] = | ||||
( compressed_spec_mat[ offsetCompressed ] | ||||
paul
|
r243 | + averaged_spec_mat[ offsetASM + k ] ); | ||
paul
|
r126 | } | ||
paul
|
r179 | compressed_spec_mat[ offsetCompressed ] = | ||
compressed_spec_mat[ offsetCompressed ] / (divider * nbBinsToAverage); | ||||
paul
|
r126 | } | ||
} | ||||
} | ||||
void ASM_convert( volatile float *input_matrix, char *output_matrix) | ||||
{ | ||||
unsigned int frequencyBin; | ||||
unsigned int asmComponent; | ||||
char * pt_char_input; | ||||
char * pt_char_output; | ||||
unsigned int offsetInput; | ||||
unsigned int offsetOutput; | ||||
pt_char_input = (char*) &input_matrix; | ||||
pt_char_output = (char*) &output_matrix; | ||||
// convert all other data | ||||
for( frequencyBin=0; frequencyBin<NB_BINS_PER_SM; frequencyBin++) | ||||
{ | ||||
for ( asmComponent=0; asmComponent<NB_VALUES_PER_SM; asmComponent++) | ||||
{ | ||||
offsetInput = (frequencyBin*NB_VALUES_PER_SM) + asmComponent ; | ||||
offsetOutput = 2 * ( (frequencyBin*NB_VALUES_PER_SM) + asmComponent ) ; | ||||
pt_char_input = (char*) &input_matrix [ offsetInput ]; | ||||
pt_char_output = (char*) &output_matrix[ offsetOutput ]; | ||||
pt_char_output[0] = pt_char_input[0]; // bits 31 downto 24 of the float | ||||
pt_char_output[1] = pt_char_input[1]; // bits 23 downto 16 of the float | ||||
} | ||||
} | ||||
} | ||||
paul
|
r195 | void ASM_compress_reorganize_and_divide_mask(float *averaged_spec_mat, float *compressed_spec_mat, | ||
paul
|
r243 | float divider, | ||
unsigned char nbBinsCompressedMatrix, unsigned char nbBinsToAverage , unsigned char ASMIndexStart, unsigned char channel); | ||||
paul
|
r195 | |||
paul
|
r236 | int getFBinMask(int k, unsigned char channel); | ||
paul
|
r195 | |||
paul
|
r214 | void init_kcoeff_sbm_from_kcoeff_norm( float *input_kcoeff, float *output_kcoeff, unsigned char nb_bins_norm); | ||
paul
|
r126 | #endif // FSW_PROCESSING_H_INCLUDED | ||