Auto status change to Under Review
@@ -2,6 +2,7 | |||
|
2 | 2 | // version 1.0: 31/07/2013 |
|
3 | 3 | // version 1.1: 02/04/2014 |
|
4 | 4 | // version 1.2: 30/04/2014 |
|
5 | // version 1.3: 02/05/2014 | |
|
5 | 6 | |
|
6 | 7 | #include "basic_parameters.h" |
|
7 | 8 | #include <math.h> |
@@ -92,8 +93,7 void init_k_f0( void ) | |||
|
92 | 93 | float alpha_M = 45 * (3.1415927/180); |
|
93 | 94 | |
|
94 | 95 | void BP1_set( float * compressed_spec_mat, uint8_t nb_bins_compressed_spec_mat, uint8_t * lfr_bp1 ){ |
|
95 |
|
|
|
96 | float PSDB; // 32 bits floating point | |
|
96 | float PSDB; // 32-bit floating point | |
|
97 | 97 | float PSDE; |
|
98 | 98 | float tmp; |
|
99 | 99 | float NVEC_V0; |
@@ -110,13 +110,14 void BP1_set( float * compressed_spec_ma | |||
|
110 | 110 | float bx_bx_star; |
|
111 | 111 | float vphi; |
|
112 | 112 | float significand; |
|
113 |
|
|
|
113 | int exponent; // 32-bit signed integer | |
|
114 | uint8_t nbitexp; // 8-bit unsigned integer | |
|
114 | 115 | uint8_t nbitsig; |
|
115 | 116 | uint8_t tmp_uint8; |
|
116 |
uint8_t *pt_uint8; // pointer on unsigned 8-bit |
|
|
117 |
int8_t expmin; // 8 |
|
|
117 | uint8_t *pt_uint8; // pointer on unsigned 8-bit integer | |
|
118 | int8_t expmin; // 8-bit signed integer | |
|
118 | 119 | int8_t expmax; |
|
119 |
int16_t rangesig; |
|
|
120 | uint16_t rangesig; // 16-bit unsigned integer | |
|
120 | 121 | uint16_t psd; |
|
121 | 122 | uint16_t exp; |
|
122 | 123 | uint16_t tmp_uint16; |
@@ -259,16 +260,16 void BP1_set( float * compressed_spec_ma | |||
|
259 | 260 | NVEC_V1 = -compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+4] / tmp; // S13 Im => n2 |
|
260 | 261 | NVEC_V2 = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+2] / tmp; // S12 Im => n3 |
|
261 | 262 | |
|
262 |
lfr_bp1[i*NB_BYTES_BP1+4] = (uint8_t) (NVEC_V0*127.5 + 128); // |
|
|
263 |
lfr_bp1[i*NB_BYTES_BP1+5] = (uint8_t) (NVEC_V1*127.5 + 128); // |
|
|
264 |
pt_uint8 = (uint8_t*) &NVEC_V2; // |
|
|
263 | lfr_bp1[i*NB_BYTES_BP1+4] = (uint8_t) (NVEC_V0*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
264 | lfr_bp1[i*NB_BYTES_BP1+5] = (uint8_t) (NVEC_V1*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
265 | pt_uint8 = (uint8_t*) &NVEC_V2; // Affect an uint8_t pointer with the adress of NVEC_V2 | |
|
265 | 266 | #ifdef LSB_FIRST_TCH |
|
266 |
lfr_bp1[i*NB_BYTES_BP1+6] = pt_uint8[3] & 0x80; // |
|
|
267 |
|
|
|
267 | lfr_bp1[i*NB_BYTES_BP1+6] = pt_uint8[3] & 0x80; // Extract the sign bit of NVEC_V2 (32-bit float, sign bit in the 4th octet:PC convention) | |
|
268 | // Record it at the 8th bit position (from the right to the left) of lfr_bp1[i*NB_BYTES_BP1+6] | |
|
268 | 269 | #endif |
|
269 | 270 | #ifdef MSB_FIRST_TCH |
|
270 |
lfr_bp1[i*NB_BYTES_BP1+6] = pt_uint8[0] & 0x80; // |
|
|
271 |
|
|
|
271 | lfr_bp1[i*NB_BYTES_BP1+6] = pt_uint8[0] & 0x80; // Extract the sign bit of NVEC_V2 (32-bit float, sign bit in the 0th octet:SPARC convention) | |
|
272 | // Record it at the 8th bit position (from the right to the left) of lfr_bp1[i*NB_BYTES_BP1+6] | |
|
272 | 273 | #endif |
|
273 | 274 | #ifdef DEBUG_TCH |
|
274 | 275 | printf("NVEC_V0 : %16.8e\n",NVEC_V0); |
@@ -280,11 +281,11 void BP1_set( float * compressed_spec_ma | |||
|
280 | 281 | #endif |
|
281 | 282 | //======================================================= |
|
282 | 283 | // BP1 ellipticity == PA_LFR_SC_BP1_ELLIP_F0 == 4 bits |
|
283 |
aux = 2*tmp / PSDB; |
|
|
284 | aux = 2*tmp / PSDB; // Compute the ellipticity | |
|
284 | 285 | |
|
285 |
tmp_uint8 = (uint8_t) (aux*15 + 0.5); // |
|
|
286 | tmp_uint8 = (uint8_t) (aux*15 + 0.5); // Shift and cast into a 8-bit uint8_t with rounding | |
|
286 | 287 |
|
|
287 |
lfr_bp1[i*NB_BYTES_BP1+6] = lfr_bp1[i*NB_BYTES_BP1+6] | (tmp_uint8 << 3); // |
|
|
288 | lfr_bp1[i*NB_BYTES_BP1+6] = lfr_bp1[i*NB_BYTES_BP1+6] | (tmp_uint8 << 3); // Put these 4 bits next to the right place | |
|
288 | 289 | // of the sign bit of NVEC_V2 (recorded |
|
289 | 290 | // previously in lfr_bp1[i*NB_BYTES_BP1+6]) |
|
290 | 291 | #ifdef DEBUG_TCH |
@@ -304,11 +305,11 void BP1_set( float * compressed_spec_ma | |||
|
304 | 305 | + 2 * compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+10]*compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+10] |
|
305 | 306 | + 2 * compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+11]*compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+11]; |
|
306 | 307 | aux = PSDB*PSDB; |
|
307 |
tmp = ( 3*tr_SB_SB - aux ) / ( 2 * aux ); |
|
|
308 | tmp = ( 3*tr_SB_SB - aux ) / ( 2 * aux ); // Compute the degree of polarisation | |
|
308 | 309 | |
|
309 |
tmp_uint8 = (uint8_t) (tmp*7 + 0.5); |
|
|
310 | tmp_uint8 = (uint8_t) (tmp*7 + 0.5); // Shift and cast into a 8-bit uint8_t with rounding | |
|
310 | 311 | // where just the first 3 bits are used (0, ..., 7) |
|
311 |
lfr_bp1[i*NB_BYTES_BP1+6] = lfr_bp1[i*NB_BYTES_BP1+6] | tmp_uint8; // |
|
|
312 | lfr_bp1[i*NB_BYTES_BP1+6] = lfr_bp1[i*NB_BYTES_BP1+6] | tmp_uint8; // Record these 3 bits at the 3 first bit positions | |
|
312 | 313 | // (from the right to the left) of lfr_bp1[i*NB_BYTES_BP1+6] |
|
313 | 314 | #ifdef DEBUG_TCH |
|
314 | 315 | printf("DOP : %16.8e\n",tmp); |
@@ -390,7 +391,7 void BP1_set( float * compressed_spec_ma | |||
|
390 | 391 | printf("lfr_bp1[i*NB_BYTES_BP1+7] for ReaSX significand : %u\n",lfr_bp1[i*NB_BYTES_BP1+7]); |
|
391 | 392 | printf("tmp_uint8 for ReaSX exponent : %d\n",tmp_uint8); |
|
392 | 393 | #endif |
|
393 |
lfr_bp1[i*NB_BYTES_BP1+7] = lfr_bp1[i*NB_BYTES_BP1+7] | (tmp_uint8 << 3); // |
|
|
394 | lfr_bp1[i*NB_BYTES_BP1+7] = lfr_bp1[i*NB_BYTES_BP1+7] | (tmp_uint8 << 3); // Shift these 5 bits to the left before logical addition | |
|
394 | 395 | // with lfr_bp1[i*NB_BYTES_BP1+7] |
|
395 | 396 | #ifdef DEBUG_TCH |
|
396 | 397 | printf("lfr_bp1[i*NB_BYTES_BP1+7] for ReaSX exponent + significand : %u\n",lfr_bp1[i*NB_BYTES_BP1+7]); |
@@ -535,12 +536,21 void BP1_set( float * compressed_spec_ma | |||
|
535 | 536 | |
|
536 | 537 | void BP2_set( float * compressed_spec_mat, uint8_t nb_bins_compressed_spec_mat, uint8_t * lfr_bp2 ) |
|
537 | 538 | { |
|
538 | int i, exponent; | |
|
539 | float aux, significand, cross_re, cross_im; | |
|
540 | int8_t nbitexp, nbitsig, expmin, expmax; // 8 bits | |
|
541 | int16_t rangesig; // 16 bits | |
|
542 | uint16_t autocor, tmp_uint16; // 16 bits | |
|
543 |
uint |
|
|
539 | float cross_re; // 32-bit floating point | |
|
540 | float cross_im; | |
|
541 | float aux; | |
|
542 | float significand; | |
|
543 | int exponent; // 32-bit signed integer | |
|
544 | uint8_t nbitexp; // 8-bit unsigned integer | |
|
545 | uint8_t nbitsig; | |
|
546 | uint8_t *pt_uint8; // pointer on unsigned 8-bit integer | |
|
547 | int8_t expmin; // 8-bit signed integer | |
|
548 | int8_t expmax; | |
|
549 | uint16_t rangesig; // 16-bit unsigned integer | |
|
550 | uint16_t autocor; | |
|
551 | uint16_t exp; | |
|
552 | uint16_t tmp_uint16; | |
|
553 | uint16_t i; | |
|
544 | 554 | |
|
545 | 555 | #ifdef DEBUG_TCH |
|
546 | 556 | printf("BP2 : \n"); |
@@ -586,8 +596,8 void BP2_set( float * compressed_spec_ma | |||
|
586 | 596 | aux = sqrt(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX]*compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+9]); |
|
587 | 597 | cross_re = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+1] / aux; |
|
588 | 598 | cross_im = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+2] / aux; |
|
589 |
lfr_bp2[i*NB_BYTES_BP2+10] = (uint8_t) (cross_re*127.5 + 128); // |
|
|
590 |
lfr_bp2[i*NB_BYTES_BP2+20] = (uint8_t) (cross_im*127.5 + 128); // |
|
|
599 | lfr_bp2[i*NB_BYTES_BP2+10] = (uint8_t) (cross_re*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
600 | lfr_bp2[i*NB_BYTES_BP2+20] = (uint8_t) (cross_im*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
591 | 601 | #ifdef DEBUG_TCH |
|
592 | 602 | printf("\nBin number: %d\n", i); |
|
593 | 603 | printf("lfr_bp2[i*NB_BYTES_BP2+10] for cross12_re (%16.8e) : %.3u\n",cross_re, lfr_bp2[i*NB_BYTES_BP2+10]); |
@@ -597,8 +607,8 void BP2_set( float * compressed_spec_ma | |||
|
597 | 607 | aux = sqrt(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX]*compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+16]); |
|
598 | 608 | cross_re = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+3] / aux; |
|
599 | 609 | cross_im = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+4] / aux; |
|
600 |
lfr_bp2[i*NB_BYTES_BP2+11] = (uint8_t) (cross_re*127.5 + 128); // |
|
|
601 |
lfr_bp2[i*NB_BYTES_BP2+21] = (uint8_t) (cross_im*127.5 + 128); // |
|
|
610 | lfr_bp2[i*NB_BYTES_BP2+11] = (uint8_t) (cross_re*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
611 | lfr_bp2[i*NB_BYTES_BP2+21] = (uint8_t) (cross_im*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
602 | 612 | #ifdef DEBUG_TCH |
|
603 | 613 | printf("lfr_bp2[i*NB_BYTES_BP2+11] for cross13_re (%16.8e) : %.3u\n",cross_re, lfr_bp2[i*NB_BYTES_BP2+11]); |
|
604 | 614 | printf("lfr_bp2[i*NB_BYTES_BP2+21] for cross13_im (%16.8e) : %.3u\n",cross_im, lfr_bp2[i*NB_BYTES_BP2+21]); |
@@ -607,8 +617,8 void BP2_set( float * compressed_spec_ma | |||
|
607 | 617 | aux = sqrt(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX]*compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+21]); |
|
608 | 618 | cross_re = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+5] / aux; |
|
609 | 619 | cross_im = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+6] / aux; |
|
610 |
lfr_bp2[i*NB_BYTES_BP2+12] = (uint8_t) (cross_re*127.5 + 128); // |
|
|
611 |
lfr_bp2[i*NB_BYTES_BP2+22] = (uint8_t) (cross_im*127.5 + 128); // |
|
|
620 | lfr_bp2[i*NB_BYTES_BP2+12] = (uint8_t) (cross_re*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
621 | lfr_bp2[i*NB_BYTES_BP2+22] = (uint8_t) (cross_im*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
612 | 622 | #ifdef DEBUG_TCH |
|
613 | 623 | printf("lfr_bp2[i*NB_BYTES_BP2+12] for cross14_re (%16.8e) : %.3u\n",cross_re, lfr_bp2[i*NB_BYTES_BP2+12]); |
|
614 | 624 | printf("lfr_bp2[i*NB_BYTES_BP2+22] for cross14_im (%16.8e) : %.3u\n",cross_im, lfr_bp2[i*NB_BYTES_BP2+22]); |
@@ -617,8 +627,8 void BP2_set( float * compressed_spec_ma | |||
|
617 | 627 | aux = sqrt(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX]*compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+24]); |
|
618 | 628 | cross_re = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+7] / aux; |
|
619 | 629 | cross_im = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+8] / aux; |
|
620 |
lfr_bp2[i*NB_BYTES_BP2+13] = (uint8_t) (cross_re*127.5 + 128); // |
|
|
621 |
lfr_bp2[i*NB_BYTES_BP2+23] = (uint8_t) (cross_im*127.5 + 128); // |
|
|
630 | lfr_bp2[i*NB_BYTES_BP2+13] = (uint8_t) (cross_re*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
631 | lfr_bp2[i*NB_BYTES_BP2+23] = (uint8_t) (cross_im*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
622 | 632 | #ifdef DEBUG_TCH |
|
623 | 633 | printf("lfr_bp2[i*NB_BYTES_BP2+13] for cross15_re (%16.8e) : %.3u\n",cross_re, lfr_bp2[i*NB_BYTES_BP2+13]); |
|
624 | 634 | printf("lfr_bp2[i*NB_BYTES_BP2+23] for cross15_im (%16.8e) : %.3u\n",cross_im, lfr_bp2[i*NB_BYTES_BP2+23]); |
@@ -627,8 +637,8 void BP2_set( float * compressed_spec_ma | |||
|
627 | 637 | aux = sqrt(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+9]*compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+16]); |
|
628 | 638 | cross_re = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+10] / aux; |
|
629 | 639 | cross_im = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+11] / aux; |
|
630 |
lfr_bp2[i*NB_BYTES_BP2+14] = (uint8_t) (cross_re*127.5 + 128); // |
|
|
631 |
lfr_bp2[i*NB_BYTES_BP2+24] = (uint8_t) (cross_im*127.5 + 128); // |
|
|
640 | lfr_bp2[i*NB_BYTES_BP2+14] = (uint8_t) (cross_re*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
641 | lfr_bp2[i*NB_BYTES_BP2+24] = (uint8_t) (cross_im*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
632 | 642 | #ifdef DEBUG_TCH |
|
633 | 643 | printf("lfr_bp2[i*NB_BYTES_BP2+14] for cross23_re (%16.8e) : %.3u\n",cross_re, lfr_bp2[i*NB_BYTES_BP2+14]); |
|
634 | 644 | printf("lfr_bp2[i*NB_BYTES_BP2+24] for cross23_im (%16.8e) : %.3u\n",cross_im, lfr_bp2[i*NB_BYTES_BP2+24]); |
@@ -637,8 +647,8 void BP2_set( float * compressed_spec_ma | |||
|
637 | 647 | aux = sqrt(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+9]*compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+21]); |
|
638 | 648 | cross_re = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+12] / aux; |
|
639 | 649 | cross_im = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+13] / aux; |
|
640 |
lfr_bp2[i*NB_BYTES_BP2+15] = (uint8_t) (cross_re*127.5 + 128); // |
|
|
641 |
lfr_bp2[i*NB_BYTES_BP2+25] = (uint8_t) (cross_im*127.5 + 128); // |
|
|
650 | lfr_bp2[i*NB_BYTES_BP2+15] = (uint8_t) (cross_re*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
651 | lfr_bp2[i*NB_BYTES_BP2+25] = (uint8_t) (cross_im*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
642 | 652 | #ifdef DEBUG_TCH |
|
643 | 653 | printf("lfr_bp2[i*NB_BYTES_BP2+15] for cross24_re (%16.8e) : %.3u\n",cross_re, lfr_bp2[i*NB_BYTES_BP2+15]); |
|
644 | 654 | printf("lfr_bp2[i*NB_BYTES_BP2+25] for cross24_im (%16.8e) : %.3u\n",cross_im, lfr_bp2[i*NB_BYTES_BP2+25]); |
@@ -647,8 +657,8 void BP2_set( float * compressed_spec_ma | |||
|
647 | 657 | aux = sqrt(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+9]*compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+24]); |
|
648 | 658 | cross_re = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+14] / aux; |
|
649 | 659 | cross_im = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+15] / aux; |
|
650 |
lfr_bp2[i*NB_BYTES_BP2+16] = (uint8_t) (cross_re*127.5 + 128); // |
|
|
651 |
lfr_bp2[i*NB_BYTES_BP2+26] = (uint8_t) (cross_im*127.5 + 128); // |
|
|
660 | lfr_bp2[i*NB_BYTES_BP2+16] = (uint8_t) (cross_re*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
661 | lfr_bp2[i*NB_BYTES_BP2+26] = (uint8_t) (cross_im*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
652 | 662 | #ifdef DEBUG_TCH |
|
653 | 663 | printf("lfr_bp2[i*NB_BYTES_BP2+16] for cross25_re (%16.8e) : %.3u\n",cross_re, lfr_bp2[i*NB_BYTES_BP2+16]); |
|
654 | 664 | printf("lfr_bp2[i*NB_BYTES_BP2+26] for cross25_im (%16.8e) : %.3u\n",cross_im, lfr_bp2[i*NB_BYTES_BP2+26]); |
@@ -657,8 +667,8 void BP2_set( float * compressed_spec_ma | |||
|
657 | 667 | aux = sqrt(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+16]*compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+21]); |
|
658 | 668 | cross_re = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+17] / aux; |
|
659 | 669 | cross_im = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+18] / aux; |
|
660 |
lfr_bp2[i*NB_BYTES_BP2+17] = (uint8_t) (cross_re*127.5 + 128); // |
|
|
661 |
lfr_bp2[i*NB_BYTES_BP2+27] = (uint8_t) (cross_im*127.5 + 128); // |
|
|
670 | lfr_bp2[i*NB_BYTES_BP2+17] = (uint8_t) (cross_re*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
671 | lfr_bp2[i*NB_BYTES_BP2+27] = (uint8_t) (cross_im*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
662 | 672 | #ifdef DEBUG_TCH |
|
663 | 673 | printf("lfr_bp2[i*NB_BYTES_BP2+17] for cross34_re (%16.8e) : %.3u\n",cross_re, lfr_bp2[i*NB_BYTES_BP2+17]); |
|
664 | 674 | printf("lfr_bp2[i*NB_BYTES_BP2+27] for cross34_im (%16.8e) : %.3u\n",cross_im, lfr_bp2[i*NB_BYTES_BP2+27]); |
@@ -667,8 +677,8 void BP2_set( float * compressed_spec_ma | |||
|
667 | 677 | aux = sqrt(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+16]*compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+24]); |
|
668 | 678 | cross_re = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+19] / aux; |
|
669 | 679 | cross_im = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+20] / aux; |
|
670 |
lfr_bp2[i*NB_BYTES_BP2+18] = (uint8_t) (cross_re*127.5 + 128); // |
|
|
671 |
lfr_bp2[i*NB_BYTES_BP2+28] = (uint8_t) (cross_im*127.5 + 128); // |
|
|
680 | lfr_bp2[i*NB_BYTES_BP2+18] = (uint8_t) (cross_re*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
681 | lfr_bp2[i*NB_BYTES_BP2+28] = (uint8_t) (cross_im*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
672 | 682 | #ifdef DEBUG_TCH |
|
673 | 683 | printf("lfr_bp2[i*NB_BYTES_BP2+18] for cross35_re (%16.8e) : %.3u\n",cross_re, lfr_bp2[i*NB_BYTES_BP2+18]); |
|
674 | 684 | printf("lfr_bp2[i*NB_BYTES_BP2+28] for cross35_im (%16.8e) : %.3u\n",cross_im, lfr_bp2[i*NB_BYTES_BP2+28]); |
@@ -677,8 +687,8 void BP2_set( float * compressed_spec_ma | |||
|
677 | 687 | aux = sqrt(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+21]*compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+24]); |
|
678 | 688 | cross_re = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+22] / aux; |
|
679 | 689 | cross_im = compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+23] / aux; |
|
680 |
lfr_bp2[i*NB_BYTES_BP2+19] = (uint8_t) (cross_re*127.5 + 128); // |
|
|
681 |
lfr_bp2[i*NB_BYTES_BP2+29] = (uint8_t) (cross_im*127.5 + 128); // |
|
|
690 | lfr_bp2[i*NB_BYTES_BP2+19] = (uint8_t) (cross_re*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
691 | lfr_bp2[i*NB_BYTES_BP2+29] = (uint8_t) (cross_im*127.5 + 128); // Shift and cast into a 8-bit uint8_t (0, ..., 255) with rounding | |
|
682 | 692 | #ifdef DEBUG_TCH |
|
683 | 693 | printf("lfr_bp2[i*NB_BYTES_BP2+19] for cross45_re (%16.8e) : %.3u\n",cross_re, lfr_bp2[i*NB_BYTES_BP2+19]); |
|
684 | 694 | printf("lfr_bp2[i*NB_BYTES_BP2+29] for cross45_im (%16.8e) : %.3u\n",cross_im, lfr_bp2[i*NB_BYTES_BP2+29]); |
@@ -713,19 +723,27 void BP2_set( float * compressed_spec_ma | |||
|
713 | 723 | |
|
714 | 724 | autocor = (uint16_t) ((significand*2-1)*rangesig + 0.5); // Shift and cast into a 16-bit unsigned int with rounding |
|
715 | 725 |
|
|
716 |
|
|
|
717 |
|
|
|
718 | pt_u_short_int = (uint16_t*) &lfr_bp2[i*NB_BYTES_BP2+0]; // Affect an uint16_t pointer with the | |
|
719 |
|
|
|
720 | *pt_u_short_int = autocor | (tmp_uint16 << nbitsig); // Put the exponent bits (nbitexp) next to the | |
|
721 | // left place of the significand bits (nbitsig), making | |
|
722 | // the 16-bit word to be recorded, and record it using the pointer | |
|
726 | exp = (uint16_t) (exponent-expmin); // Shift and cast into a 16-bit unsigned int where just | |
|
727 | // the first nbitexp bits are used (0, ..., 2^nbitexp-1) | |
|
728 | tmp_uint16 = autocor | (exp << nbitsig); // Put the exponent bits (nbitexp) next to the | |
|
729 | // left place of the significand bits (nbitsig), | |
|
730 | // making the 16-bit word to be recorded | |
|
731 | pt_uint8 = (uint8_t*) &tmp_uint16; // Affect an uint8_t pointer with the adress of tmp_uint16 | |
|
732 | #ifdef LSB_FIRST_TCH | |
|
733 | lfr_bp2[i*NB_BYTES_BP2+0] = pt_uint8[0]; // Record LSB of tmp_uint16 | |
|
734 | lfr_bp2[i*NB_BYTES_BP2+1] = pt_uint8[1]; // Record MSB of tmp_uint16 | |
|
735 | #endif | |
|
736 | #ifdef MSB_FIRST_TCH | |
|
737 | lfr_bp2[i*NB_BYTES_BP2+0] = pt_uint8[1]; // Record LSB of tmp_uint16 | |
|
738 | lfr_bp2[i*NB_BYTES_BP2+1] = pt_uint8[0]; // Record MSB of tmp_uint16 | |
|
739 | #endif | |
|
723 | 740 | #ifdef DEBUG_TCH |
|
724 | 741 |
printf("autocor for S11 significand : %u\n",autocor |
|
725 |
printf(" |
|
|
726 |
printf(" |
|
|
727 | printf("lfr_bp2[i*NB_BYTES_BP2+1] : %u or %x\n",lfr_bp2[i*NB_BYTES_BP2+1], lfr_bp2[i*NB_BYTES_BP2+1]); | |
|
728 |
printf("lfr_bp2[i*NB_BYTES_BP2+ |
|
|
742 | printf("exp for S11 exponent : %u\n",exp); | |
|
743 | printf("pt_uint8[1] for S11 exponent + significand : %.3d or %2x\n",pt_uint8[1], pt_uint8[1]); | |
|
744 | printf("pt_uint8[0] for S11 exponent + significand : %.3d or %2x\n",pt_uint8[0], pt_uint8[0]); | |
|
745 | printf("lfr_bp2[i*NB_BYTES_BP2+1] : %3u or %2x\n",lfr_bp2[i*NB_BYTES_BP2+1], lfr_bp2[i*NB_BYTES_BP2+1]); | |
|
746 | printf("lfr_bp2[i*NB_BYTES_BP2+0] : %3u or %2x\n",lfr_bp2[i*NB_BYTES_BP2+0], lfr_bp2[i*NB_BYTES_BP2+0]); | |
|
729 | 747 | #endif |
|
730 | 748 | // S22 |
|
731 | 749 | significand = frexpf(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+9], &exponent); // 0.5 <= significand < 1 |
@@ -750,19 +768,27 void BP2_set( float * compressed_spec_ma | |||
|
750 | 768 | |
|
751 | 769 | autocor = (uint16_t) ((significand*2-1)*rangesig + 0.5); // Shift and cast into a 16-bit unsigned int with rounding |
|
752 | 770 |
|
|
753 |
|
|
|
754 |
|
|
|
755 | pt_u_short_int = (uint16_t*) &lfr_bp2[i*NB_BYTES_BP2+2]; // Affect an uint16_t pointer with the | |
|
756 |
|
|
|
757 | *pt_u_short_int = autocor | (tmp_uint16 << nbitsig); // Put the exponent bits (nbitexp) next to the | |
|
758 | // left place of the significand bits (nbitsig), making | |
|
759 | // the 16-bit word to be recorded, and record it using the pointer | |
|
771 | exp = (uint16_t) (exponent-expmin); // Shift and cast into a 16-bit unsigned int where just | |
|
772 | // the first nbitexp bits are used (0, ..., 2^nbitexp-1) | |
|
773 | tmp_uint16 = autocor | (exp << nbitsig); // Put the exponent bits (nbitexp) next to the | |
|
774 | // left place of the significand bits (nbitsig), | |
|
775 | // making the 16-bit word to be recorded | |
|
776 | pt_uint8 = (uint8_t*) &tmp_uint16; // Affect an uint8_t pointer with the adress of tmp_uint16 | |
|
777 | #ifdef LSB_FIRST_TCH | |
|
778 | lfr_bp2[i*NB_BYTES_BP2+2] = pt_uint8[0]; // Record LSB of tmp_uint16 | |
|
779 | lfr_bp2[i*NB_BYTES_BP2+3] = pt_uint8[1]; // Record MSB of tmp_uint16 | |
|
780 | #endif | |
|
781 | #ifdef MSB_FIRST_TCH | |
|
782 | lfr_bp2[i*NB_BYTES_BP2+2] = pt_uint8[1]; // Record LSB of tmp_uint16 | |
|
783 | lfr_bp2[i*NB_BYTES_BP2+3] = pt_uint8[0]; // Record MSB of tmp_uint16 | |
|
784 | #endif | |
|
760 | 785 | #ifdef DEBUG_TCH |
|
761 |
printf("autocor for S22 significand : % |
|
|
762 |
printf(" |
|
|
763 |
printf(" |
|
|
764 | printf("lfr_bp2[i*NB_BYTES_BP2+3] : %.3d or %x\n",lfr_bp2[i*NB_BYTES_BP2+3], lfr_bp2[i*NB_BYTES_BP2+3]); | |
|
765 |
printf("lfr_bp2[i*NB_BYTES_BP2+ |
|
|
786 | printf("autocor for S22 significand : %u\n",autocor); | |
|
787 | printf("exp for S11 exponent : %u\n",exp); | |
|
788 | printf("pt_uint8[1] for S22 exponent + significand : %.3d or %2x\n",pt_uint8[1], pt_uint8[1]); | |
|
789 | printf("pt_uint8[0] for S22 exponent + significand : %.3d or %2x\n",pt_uint8[0], pt_uint8[0]); | |
|
790 | printf("lfr_bp2[i*NB_BYTES_BP2+3] : %3u or %2x\n",lfr_bp2[i*NB_BYTES_BP2+3], lfr_bp2[i*NB_BYTES_BP2+3]); | |
|
791 | printf("lfr_bp2[i*NB_BYTES_BP2+2] : %3u or %2x\n",lfr_bp2[i*NB_BYTES_BP2+2], lfr_bp2[i*NB_BYTES_BP2+2]); | |
|
766 | 792 | #endif |
|
767 | 793 | // S33 |
|
768 | 794 | significand = frexpf(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+16], &exponent); // 0.5 <= significand < 1 |
@@ -787,19 +813,27 void BP2_set( float * compressed_spec_ma | |||
|
787 | 813 | |
|
788 | 814 | autocor = (uint16_t) ((significand*2-1)*rangesig + 0.5); // Shift and cast into a 16-bit unsigned int with rounding |
|
789 | 815 |
|
|
790 |
|
|
|
791 |
|
|
|
792 | pt_u_short_int = (uint16_t*) &lfr_bp2[i*NB_BYTES_BP2+4]; // Affect an uint16_t pointer with the | |
|
793 |
|
|
|
794 | *pt_u_short_int = autocor | (tmp_uint16 << nbitsig); // Put the exponent bits (nbitexp) next to the | |
|
795 | // left place of the significand bits (nbitsig), making | |
|
796 | // the 16-bit word to be recorded, and record it using the pointer | |
|
816 | exp = (uint16_t) (exponent-expmin); // Shift and cast into a 16-bit unsigned int where just | |
|
817 | // the first nbitexp bits are used (0, ..., 2^nbitexp-1) | |
|
818 | tmp_uint16 = autocor | (exp << nbitsig); // Put the exponent bits (nbitexp) next to the | |
|
819 | // left place of the significand bits (nbitsig), | |
|
820 | // making the 16-bit word to be recorded | |
|
821 | pt_uint8 = (uint8_t*) &tmp_uint16; // Affect an uint8_t pointer with the adress of tmp_uint16 | |
|
822 | #ifdef LSB_FIRST_TCH | |
|
823 | lfr_bp2[i*NB_BYTES_BP2+4] = pt_uint8[0]; // Record LSB of tmp_uint16 | |
|
824 | lfr_bp2[i*NB_BYTES_BP2+5] = pt_uint8[1]; // Record MSB of tmp_uint16 | |
|
825 | #endif | |
|
826 | #ifdef MSB_FIRST_TCH | |
|
827 | lfr_bp2[i*NB_BYTES_BP2+4] = pt_uint8[1]; // Record LSB of tmp_uint16 | |
|
828 | lfr_bp2[i*NB_BYTES_BP2+5] = pt_uint8[0]; // Record MSB of tmp_uint16 | |
|
829 | #endif | |
|
797 | 830 | #ifdef DEBUG_TCH |
|
798 |
printf("autocor for S33 significand : % |
|
|
799 |
printf(" |
|
|
800 |
printf(" |
|
|
801 | printf("lfr_bp2[i*NB_BYTES_BP2+5] : %.3d or %x\n",lfr_bp2[i*NB_BYTES_BP2+5], lfr_bp2[i*NB_BYTES_BP2+5]); | |
|
802 |
printf("lfr_bp2[i*NB_BYTES_BP2+ |
|
|
831 | printf("autocor for S33 significand : %u\n",autocor); | |
|
832 | printf("exp for S33 exponent : %u\n",exp); | |
|
833 | printf("pt_uint8[1] for S33 exponent + significand : %.3d or %2x\n",pt_uint8[1], pt_uint8[1]); | |
|
834 | printf("pt_uint8[0] for S33 exponent + significand : %.3d or %2x\n",pt_uint8[0], pt_uint8[0]); | |
|
835 | printf("lfr_bp2[i*NB_BYTES_BP2+5] : %3u or %2x\n",lfr_bp2[i*NB_BYTES_BP2+5], lfr_bp2[i*NB_BYTES_BP2+5]); | |
|
836 | printf("lfr_bp2[i*NB_BYTES_BP2+4] : %3u or %2x\n",lfr_bp2[i*NB_BYTES_BP2+4], lfr_bp2[i*NB_BYTES_BP2+4]); | |
|
803 | 837 | #endif |
|
804 | 838 | // S44 |
|
805 | 839 | significand = frexpf(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+21], &exponent); // 0.5 <= significand < 1 |
@@ -825,19 +859,27 void BP2_set( float * compressed_spec_ma | |||
|
825 | 859 | |
|
826 | 860 | autocor = (uint16_t) ((significand*2-1)*rangesig + 0.5); // Shift and cast into a 16-bit unsigned int with rounding |
|
827 | 861 |
|
|
828 |
|
|
|
829 |
|
|
|
830 | pt_u_short_int = (uint16_t*) &lfr_bp2[i*NB_BYTES_BP2+6]; // Affect an uint16_t pointer with the | |
|
831 |
|
|
|
832 | *pt_u_short_int = autocor | (tmp_uint16 << nbitsig); // Put the exponent bits (nbitexp) next to the | |
|
833 | // left place of the significand bits (nbitsig), making | |
|
834 | // the 16-bit word to be recorded, and record it using the pointer | |
|
862 | exp = (uint16_t) (exponent-expmin); // Shift and cast into a 16-bit unsigned int where just | |
|
863 | // the first nbitexp bits are used (0, ..., 2^nbitexp-1) | |
|
864 | tmp_uint16 = autocor | (exp << nbitsig); // Put the exponent bits (nbitexp) next to the | |
|
865 | // left place of the significand bits (nbitsig), | |
|
866 | // making the 16-bit word to be recorded | |
|
867 | pt_uint8 = (uint8_t*) &tmp_uint16; // Affect an uint8_t pointer with the adress of tmp_uint16 | |
|
868 | #ifdef LSB_FIRST_TCH | |
|
869 | lfr_bp2[i*NB_BYTES_BP2+6] = pt_uint8[0]; // Record LSB of tmp_uint16 | |
|
870 | lfr_bp2[i*NB_BYTES_BP2+7] = pt_uint8[1]; // Record MSB of tmp_uint16 | |
|
871 | #endif | |
|
872 | #ifdef MSB_FIRST_TCH | |
|
873 | lfr_bp2[i*NB_BYTES_BP2+6] = pt_uint8[1]; // Record LSB of tmp_uint16 | |
|
874 | lfr_bp2[i*NB_BYTES_BP2+7] = pt_uint8[0]; // Record MSB of tmp_uint16 | |
|
875 | #endif | |
|
835 | 876 | #ifdef DEBUG_TCH |
|
836 |
printf("autocor for S44 significand : % |
|
|
837 |
printf(" |
|
|
838 |
printf(" |
|
|
839 | printf("lfr_bp2[i*NB_BYTES_BP2+7] : %.3d or %x\n",lfr_bp2[i*NB_BYTES_BP2+7], lfr_bp2[i*NB_BYTES_BP2+7]); | |
|
840 |
printf("lfr_bp2[i*NB_BYTES_BP2+ |
|
|
877 | printf("autocor for S44 significand : %u\n",autocor); | |
|
878 | printf("exp for S44 exponent : %u\n",exp); | |
|
879 | printf("pt_uint8[1] for S44 exponent + significand : %.3d or %2x\n",pt_uint8[1], pt_uint8[1]); | |
|
880 | printf("pt_uint8[0] for S44 exponent + significand : %.3d or %2x\n",pt_uint8[0], pt_uint8[0]); | |
|
881 | printf("lfr_bp2[i*NB_BYTES_BP2+7] : %3u or %2x\n",lfr_bp2[i*NB_BYTES_BP2+7], lfr_bp2[i*NB_BYTES_BP2+7]); | |
|
882 | printf("lfr_bp2[i*NB_BYTES_BP2+6] : %3u or %2x\n",lfr_bp2[i*NB_BYTES_BP2+6], lfr_bp2[i*NB_BYTES_BP2+6]); | |
|
841 | 883 | #endif |
|
842 | 884 | // S55 |
|
843 | 885 | significand = frexpf(compressed_spec_mat[i*NB_VALUES_PER_SPECTRAL_MATRIX+24], &exponent); // 0.5 <= significand < 1 |
@@ -862,19 +904,27 void BP2_set( float * compressed_spec_ma | |||
|
862 | 904 | |
|
863 | 905 | autocor = (uint16_t) ((significand*2-1)*rangesig + 0.5); // Shift and cast into a 16-bit unsigned int with rounding |
|
864 | 906 |
|
|
865 |
|
|
|
866 |
|
|
|
867 | pt_u_short_int = (uint16_t*) &lfr_bp2[i*NB_BYTES_BP2+8]; // Affect an uint16_t pointer with the | |
|
868 |
|
|
|
869 | *pt_u_short_int = autocor | (tmp_uint16 << nbitsig); // Put the exponent bits (nbitexp) next to the | |
|
870 | // left place of the significand bits (nbitsig), making | |
|
871 | // the 16-bit word to be recorded, and record it using the pointer | |
|
907 | exp = (uint16_t) (exponent-expmin); // Shift and cast into a 16-bit unsigned int where just | |
|
908 | // the first nbitexp bits are used (0, ..., 2^nbitexp-1) | |
|
909 | tmp_uint16 = autocor | (exp << nbitsig); // Put the exponent bits (nbitexp) next to the | |
|
910 | // left place of the significand bits (nbitsig), | |
|
911 | // making the 16-bit word to be recorded | |
|
912 | pt_uint8 = (uint8_t*) &tmp_uint16; // Affect an uint8_t pointer with the adress of tmp_uint16 | |
|
913 | #ifdef LSB_FIRST_TCH | |
|
914 | lfr_bp2[i*NB_BYTES_BP2+8] = pt_uint8[0]; // Record LSB of tmp_uint16 | |
|
915 | lfr_bp2[i*NB_BYTES_BP2+9] = pt_uint8[1]; // Record MSB of tmp_uint16 | |
|
916 | #endif | |
|
917 | #ifdef MSB_FIRST_TCH | |
|
918 | lfr_bp2[i*NB_BYTES_BP2+8] = pt_uint8[1]; // Record LSB of tmp_uint16 | |
|
919 | lfr_bp2[i*NB_BYTES_BP2+9] = pt_uint8[0]; // Record MSB of tmp_uint16 | |
|
920 | #endif | |
|
872 | 921 | #ifdef DEBUG_TCH |
|
873 |
printf("autocor for S55 significand : % |
|
|
874 |
printf(" |
|
|
875 |
printf(" |
|
|
876 | printf("lfr_bp2[i*NB_BYTES_BP2+9] : %.3d or %x\n",lfr_bp2[i*NB_BYTES_BP2+9], lfr_bp2[i*NB_BYTES_BP2+9]); | |
|
877 |
printf("lfr_bp2[i*NB_BYTES_BP2+ |
|
|
922 | printf("autocor for S55 significand : %u\n",autocor); | |
|
923 | printf("exp for S55 exponent : %u\n",exp); | |
|
924 | printf("pt_uint8[1] for S55 exponent + significand : %.3d or %2x\n",pt_uint8[1], pt_uint8[1]); | |
|
925 | printf("pt_uint8[0] for S55 exponent + significand : %.3d or %2x\n",pt_uint8[0], pt_uint8[0]); | |
|
926 | printf("lfr_bp2[i*NB_BYTES_BP2+9] : %3u or %2x\n",lfr_bp2[i*NB_BYTES_BP2+9], lfr_bp2[i*NB_BYTES_BP2+9]); | |
|
927 | printf("lfr_bp2[i*NB_BYTES_BP2+8] : %3u or %2x\n",lfr_bp2[i*NB_BYTES_BP2+8], lfr_bp2[i*NB_BYTES_BP2+8]); | |
|
878 | 928 | #endif |
|
879 | 929 | } |
|
880 | 930 | } |
@@ -2,12 +2,13 | |||
|
2 | 2 | // version 1.0: 31/07/2013 |
|
3 | 3 | // version 1.1: 02/04/2014 |
|
4 | 4 | // version 1.2: 30/04/2014 |
|
5 | // version 1.3: 02/05/2014 | |
|
5 | 6 | |
|
6 | 7 | #ifndef BASIC_PARAMETERS_H_INCLUDED |
|
7 | 8 | #define BASIC_PARAMETERS_H_INCLUDED |
|
8 | 9 | |
|
9 | 10 | #define NB_VALUES_PER_SPECTRAL_MATRIX 25 |
|
10 | #define NB_BINS_COMPRESSED_MATRIX_f0 1 | |
|
11 | #define NB_BINS_COMPRESSED_MATRIX_f0 11 | |
|
11 | 12 | |
|
12 | 13 | #define NB_BYTES_BP1 9 |
|
13 | 14 | #define NB_BYTES_BP2 30 |
@@ -2,6 +2,7 | |||
|
2 | 2 | // version 1.0: 31/07/2013 |
|
3 | 3 | // version 1.1: 02/04/2014 |
|
4 | 4 | // version 1.2: 30/04/2014 |
|
5 | // version 1.3: 02/05/2014 | |
|
5 | 6 | |
|
6 | 7 | #include <file_utilities.h> |
|
7 | 8 | |
@@ -17,10 +18,11 int lecture_file_sm(const char *fileName | |||
|
17 | 18 | } |
|
18 | 19 | (void) fread(compressed_spectral_matrix_f0, sizeof(compressed_spectral_matrix_f0), 1, infile); |
|
19 | 20 | (void) fclose(infile); |
|
20 | printf("size of compressed_spectral_matrix_f0 : %lu\n", sizeof(compressed_spectral_matrix_f0)); | |
|
21 | printf("Number of bins: %d\n\n", NB_BINS_COMPRESSED_MATRIX_f0); | |
|
22 | 21 | |
|
23 |
printf(" |
|
|
22 | printf("Compressed_spectral_matrix_f0 : \n"); | |
|
23 | printf("Number of bins: %d\n", NB_BINS_COMPRESSED_MATRIX_f0); | |
|
24 | printf("Number of values per spectral matrix: %d\n", NB_VALUES_PER_SPECTRAL_MATRIX); | |
|
25 | printf("Size of compressed_spectral_matrix_f0 : %lu\n", sizeof(compressed_spectral_matrix_f0)); | |
|
24 | 26 | |
|
25 | 27 | for(i=0; i<NB_BINS_COMPRESSED_MATRIX_f0; i++){ |
|
26 | 28 |
General Comments 1
You need to be logged in to leave comments.
Login now