HG_REPOSITORIES/LPP/INSTRUMENTATION/USERS/JEANDET/LFR_FSW Files · src/basic_parameters.c

commit before release

paul - - Load All Authors

File last commit:

r0:452a68f7426f default


                r25:68625add1e50

default

Download file

             basic_parameters.c
        
                    194 lines
            
             | 11.4 KiB
            
                | text/x-c
            
             |
                CLexer
            
             / src / basic_parameters.c
          
                    History
                
                 |
                  Source
                 | Raw
                 |Copy content
                 |Copy permalink

        paul
    
Initial commit.

              r0
            
      #include <..\header\basic_parameters.h>

      #include<math.h>

      float k14_re = 1;

      float k14_im = 1;

      float k14_bis_re = 1;

      float k14_bis_im = 1;

      float k14_tris_re = 1;

      float k14_tris_im = 1;

      float k15_re = 1;

      float k15_im = 1;

      float k15_bis_re = 1;

      float k15_bis_im = 1;

      float k24_re = 1;

      float k24_im = 1;

      float k24_bis_re = 1;

      float k24_bis_im = 1;

      float k24_tris_re = 1;

      float k24_tris_im = 1;

      float k25_re = 1;

      float k25_im = 1;

      float k25_bis_re = 1;

      float k25_bis_im = 1;

      float k34_re = 1;

      float k34_im = 1;

      float k44 = 1;

      float k55 = 1;

      float k45_re = 1;

      float k45_im = 1;

      float alpha_M = M_PI/4;

      void BP1_set(){

          int i, j;

          unsigned char tmp_u_char;

          unsigned char * pt_char;

          float PSDB, PSDE;

          float NVEC_V0, NVEC_V1, NVEC_V2;

          float significand;

          int exponent;

          float aux, tr_SB_SB, tmp;

          float e_cross_b_re, e_cross_b_im;

          float n_cross_e_scal_b_re = 0, n_cross_e_scal_b_im = 0;

          float nx = 0, ny = 0;

          float bz_bz_star = 0;

          for(i=0; i<NB_BINS_COMPRESSED_MATRIX_f0; i++){

              //==============================================

              // BP1 PSD == B PAR_LFR_SC_BP1_PE_FL0 == 16 bits

              PSDB = compressed_spectral_matrix_f0[i*30]      // S11

                  + compressed_spectral_matrix_f0[i*30+10]    // S22

                  + compressed_spectral_matrix_f0[i*30+18];   // S33

              significand = frexp(PSDB, &exponent);

              pt_char = (unsigned char*) &PSDB;

              LFR_BP1_F0[i*9+8] = pt_char[0]; // bits 31 downto 24 of the float

              LFR_BP1_F0[i*9+7] = pt_char[1];  // bits 23 downto 16 of the float

              //==============================================

              // BP1 PSD == E PAR_LFR_SC_BP1_PB_FL0 == 16 bits

              PSDE = compressed_spectral_matrix_f0[i*30+24] * k44     // S44

                  + compressed_spectral_matrix_f0[i*30+28] * k55      // S55

                  + compressed_spectral_matrix_f0[i*30+26] * k45_re   // S45

                  - compressed_spectral_matrix_f0[i*30+27] * k45_im;  // S45

              pt_char = (unsigned char*) &PSDE;

              LFR_BP1_F0[i*9+6] = pt_char[0]; // bits 31 downto 24 of the float

              LFR_BP1_F0[i*9+5] = pt_char[1]; // bits 23 downto 16 of the float

              //==============================================================================

              // BP1 normal wave vector == PAR_LFR_SC_BP1_NVEC_V0_F0 == 8 bits

                                     // == PAR_LFR_SC_BP1_NVEC_V1_F0 == 8 bits

                                     // == PAR_LFR_SC_BP1_NVEC_V2_F0 == 1 bits

              tmp = sqrt(

                          compressed_spectral_matrix_f0[i*30+3]*compressed_spectral_matrix_f0[i*30+3]     //Im S12

                          +compressed_spectral_matrix_f0[i*30+5]*compressed_spectral_matrix_f0[i*30+5]    //Im S13

                          +compressed_spectral_matrix_f0[i*30+5]*compressed_spectral_matrix_f0[i*30+13]   //Im S23

                          );

              NVEC_V0 = compressed_spectral_matrix_f0[i*30+13] / tmp;  // Im S23

              NVEC_V1 = -compressed_spectral_matrix_f0[i*30+5] / tmp;  // Im S13

              NVEC_V2 = compressed_spectral_matrix_f0[i*30+1] / tmp;   // Im S12

              LFR_BP1_F0[i*9+4] = (char) (NVEC_V0*256);

              LFR_BP1_F0[i*9+3] = (char) (NVEC_V1*256);

              pt_char = (unsigned char*) &NVEC_V2;

              LFR_BP1_F0[i*9+2] = pt_char[0] & 0x80;  // extract the sign of NVEC_V2

              //=======================================================

              // BP1 ellipticity == PAR_LFR_SC_BP1_ELLIP_F0   == 4 bits

              aux = 2*tmp / PSDB;                                             // compute the ellipticity

              tmp_u_char = (unsigned char) (aux*(16-1));                      // convert the ellipticity

              LFR_BP1_F0[i*9+2] = LFR_BP1_F0[i*9+2] | ((tmp_u_char&0x0f)<<3); // keeps 4 bits of the resulting unsigned char

              //==============================================================

              // BP1 degree of polarization == PAR_LFR_SC_BP1_DOP_F0 == 3 bits

              for(j = 0; j<NB_VALUES_PER_SPECTRAL_MATRIX;j++){

                  tr_SB_SB = compressed_spectral_matrix_f0[i*30]*compressed_spectral_matrix_f0[i*30]

                          + compressed_spectral_matrix_f0[i*30+10]*compressed_spectral_matrix_f0[i*30+10]

                          + compressed_spectral_matrix_f0[i*30+18]*compressed_spectral_matrix_f0[i*30+18]

                          + 2 * compressed_spectral_matrix_f0[i*30+2]*compressed_spectral_matrix_f0[i*30+2]

                          + 2 * compressed_spectral_matrix_f0[i*30+3]*compressed_spectral_matrix_f0[i*30+3]

                          + 2 * compressed_spectral_matrix_f0[i*30+4]*compressed_spectral_matrix_f0[i*30+4]

                          + 2 * compressed_spectral_matrix_f0[i*30+5]*compressed_spectral_matrix_f0[i*30+5]

                          + 2 * compressed_spectral_matrix_f0[i*30+12]*compressed_spectral_matrix_f0[i*30+12]

                          + 2 * compressed_spectral_matrix_f0[i*30+13]*compressed_spectral_matrix_f0[i*30+13];

              }

              aux = PSDB*PSDB;

              tmp = ( 3*tr_SB_SB - aux ) / ( 2 * aux );

              tmp_u_char = (unsigned char) (NVEC_V0*(8-1));

              LFR_BP1_F0[i*9+2] = LFR_BP1_F0[i*9+2] | ((tmp_u_char&0x07)); // keeps 3 bits of the resulting unsigned char

              //=======================================================================================

              // BP1 z-component of the normalized Poynting flux == PAR_LFR_SC_BP1_SZ_F0 == 8 bits (7+1)

              e_cross_b_re = compressed_spectral_matrix_f0[i*30+20]*k34_re

                              + compressed_spectral_matrix_f0[i*30+6]*k14_re

                              + compressed_spectral_matrix_f0[i*30+8]*k15_re

                              + compressed_spectral_matrix_f0[i*30+14]*k24_re

                              + compressed_spectral_matrix_f0[i*30+16]*k25_re;

              e_cross_b_im = compressed_spectral_matrix_f0[i*30+21]*k34_im

                              + compressed_spectral_matrix_f0[i*30+7]*k14_im

                              + compressed_spectral_matrix_f0[i*30+9]*k15_im

                              + compressed_spectral_matrix_f0[i*30+15]*k24_im

                              + compressed_spectral_matrix_f0[i*30+17]*k25_im;

              tmp = e_cross_b_re / PSDE; // compute ReaSz

              LFR_BP1_F0[i*9+1] = ((unsigned char) (tmp * 128)) & 0x7f; // is it always positive?

              tmp = e_cross_b_re * e_cross_b_im;

              pt_char = (unsigned char*) &tmp;

              LFR_BP1_F0[i*9+1] = LFR_BP1_F0[i*9+1] | (pt_char[0] & 0x80);  // extract the sign of ArgSz

              //======================================================================

              // BP1 phase velocity estimator == PAR_LFR_SC_BP1_VPHI_F0 == 8 bits (7+1)

              nx = -sin(alpha_M)*NVEC_V0 - cos(alpha_M)*NVEC_V1;

              ny = NVEC_V2;

              bz_bz_star = cos(alpha_M) * cos(alpha_M) * compressed_spectral_matrix_f0[i*30]              // re S11

                              + sin(alpha_M) * sin(alpha_M) * compressed_spectral_matrix_f0[i*30+10]      // re S22

                              - 2 * sin(alpha_M) * cos(alpha_M) * compressed_spectral_matrix_f0[i*30+2];  // re S12

              n_cross_e_scal_b_re = nx * (compressed_spectral_matrix_f0[i*30+8]*k15_bis_re

                                              +compressed_spectral_matrix_f0[i*30+6]*k14_bis_re

                                              +compressed_spectral_matrix_f0[i*30+16]*k25_bis_re

                                              +compressed_spectral_matrix_f0[i*30+14]*k24_bis_re)

                                      + ny * (compressed_spectral_matrix_f0[i*30+6]*k14_tris_re

                                              +compressed_spectral_matrix_f0[i*30+14]*k24_tris_re);

              n_cross_e_scal_b_im = nx * (compressed_spectral_matrix_f0[i*30+8]*k15_bis_im

                                              +compressed_spectral_matrix_f0[i*30+6]*k14_bis_im

                                              +compressed_spectral_matrix_f0[i*30+16]*k25_bis_im

                                              +compressed_spectral_matrix_f0[i*30+14]*k24_bis_im)

                                      + ny * (compressed_spectral_matrix_f0[i*30+6]*k14_tris_im

                                              +compressed_spectral_matrix_f0[i*30+14]*k24_tris_im);

              tmp = n_cross_e_scal_b_re / bz_bz_star;

              LFR_BP1_F0[i*9+0] = ((unsigned char) (tmp * 128)) & 0x7f; // is it always positive?

              tmp = n_cross_e_scal_b_re * n_cross_e_scal_b_im;

              pt_char = (unsigned char*) &tmp;

              LFR_BP1_F0[i*9+1] = LFR_BP1_F0[i*9+0] | (pt_char[0] & 0x80);  // extract the sign of ArgV

          }

      }

      void BP2_set(){

          // BP2 autocorrelation

          int i, aux = 0;

          for(i = 0; i<NB_BINS_COMPRESSED_MATRIX_f0; i++){

              // S12

              aux = sqrt(compressed_spectral_matrix_f0[i*30]*compressed_spectral_matrix_f0[i*30+10]);

              compressed_spectral_matrix_f0[i*30+2] = compressed_spectral_matrix_f0[i*30+2] / aux;

              compressed_spectral_matrix_f0[i*30+3] = compressed_spectral_matrix_f0[i*30+3] / aux;

              // S13

              aux = sqrt(compressed_spectral_matrix_f0[i*30]*compressed_spectral_matrix_f0[i*30+18]);

              compressed_spectral_matrix_f0[i*30+4] = compressed_spectral_matrix_f0[i*30+4] / aux;

              compressed_spectral_matrix_f0[i*30+5] = compressed_spectral_matrix_f0[i*30+5] / aux;

              // S23

              aux = sqrt(compressed_spectral_matrix_f0[i*30+12]*compressed_spectral_matrix_f0[i*30+18]);

              compressed_spectral_matrix_f0[i*30+12] = compressed_spectral_matrix_f0[i*30+12] / aux;

              compressed_spectral_matrix_f0[i*30+13] = compressed_spectral_matrix_f0[i*30+13] / aux;

              // S45

              aux = sqrt(compressed_spectral_matrix_f0[i*30+24]*compressed_spectral_matrix_f0[i*30+28]);

              compressed_spectral_matrix_f0[i*30+26] = compressed_spectral_matrix_f0[i*30+26] / aux;

              compressed_spectral_matrix_f0[i*30+27] = compressed_spectral_matrix_f0[i*30+27] / aux;

              // S14

              aux = sqrt(compressed_spectral_matrix_f0[i*30]*compressed_spectral_matrix_f0[i*30+24]);

              compressed_spectral_matrix_f0[i*30+6] = compressed_spectral_matrix_f0[i*30+6] / aux;

              compressed_spectral_matrix_f0[i*30+7] = compressed_spectral_matrix_f0[i*30+7] / aux;

              // S15

              aux = sqrt(compressed_spectral_matrix_f0[i*30]*compressed_spectral_matrix_f0[i*30+28]);

              compressed_spectral_matrix_f0[i*30+8] = compressed_spectral_matrix_f0[i*30+8] / aux;

              compressed_spectral_matrix_f0[i*30+9] = compressed_spectral_matrix_f0[i*30+9] / aux;

              // S24

              aux = sqrt(compressed_spectral_matrix_f0[i*10]*compressed_spectral_matrix_f0[i*30+24]);

              compressed_spectral_matrix_f0[i*30+14] = compressed_spectral_matrix_f0[i*30+14] / aux;

              compressed_spectral_matrix_f0[i*30+15] = compressed_spectral_matrix_f0[i*30+15] / aux;

              // S25

              aux = sqrt(compressed_spectral_matrix_f0[i*10]*compressed_spectral_matrix_f0[i*30+28]);

              compressed_spectral_matrix_f0[i*30+16] = compressed_spectral_matrix_f0[i*30+16] / aux;

              compressed_spectral_matrix_f0[i*30+17] = compressed_spectral_matrix_f0[i*30+17] / aux;

              // S34

              aux = sqrt(compressed_spectral_matrix_f0[i*18]*compressed_spectral_matrix_f0[i*30+24]);

              compressed_spectral_matrix_f0[i*30+20] = compressed_spectral_matrix_f0[i*30+20] / aux;

              compressed_spectral_matrix_f0[i*30+21] = compressed_spectral_matrix_f0[i*30+21] / aux;

              // S35

              aux = sqrt(compressed_spectral_matrix_f0[i*18]*compressed_spectral_matrix_f0[i*30+28]);

              compressed_spectral_matrix_f0[i*30+22] = compressed_spectral_matrix_f0[i*30+22] / aux;

              compressed_spectral_matrix_f0[i*30+23] = compressed_spectral_matrix_f0[i*30+23] / aux;

          }

      }

	Site-wide shortcuts
/	Use quick search box
g h	Goto home page
g g	Goto my private gists page
g G	Goto my public gists page
g 0-9	Goto bookmarked items from 0-9
n r	New repository page
n g	New gist page

	Repositories
g s	Goto summary page
g c	Goto changelog page
g f	Goto files page
g F	Goto files page with file search activated
g p	Goto pull requests page
g o	Goto repository settings
g O	Goto repository access permissions settings
t s	Toggle sidebar on some pages

paul Initial commit.	r0	#include <..\header\basic_parameters.h>
		#include<math.h>

		float k14_re = 1;
		float k14_im = 1;
		float k14_bis_re = 1;
		float k14_bis_im = 1;
		float k14_tris_re = 1;
		float k14_tris_im = 1;
		float k15_re = 1;
		float k15_im = 1;
		float k15_bis_re = 1;
		float k15_bis_im = 1;
		float k24_re = 1;
		float k24_im = 1;
		float k24_bis_re = 1;
		float k24_bis_im = 1;
		float k24_tris_re = 1;
		float k24_tris_im = 1;
		float k25_re = 1;
		float k25_im = 1;
		float k25_bis_re = 1;
		float k25_bis_im = 1;
		float k34_re = 1;
		float k34_im = 1;
		float k44 = 1;
		float k55 = 1;
		float k45_re = 1;
		float k45_im = 1;
		float alpha_M = M_PI/4;

		void BP1_set(){
		int i, j;
		unsigned char tmp_u_char;
		unsigned char * pt_char;
		float PSDB, PSDE;
		float NVEC_V0, NVEC_V1, NVEC_V2;
		float significand;
		int exponent;
		float aux, tr_SB_SB, tmp;
		float e_cross_b_re, e_cross_b_im;
		float n_cross_e_scal_b_re = 0, n_cross_e_scal_b_im = 0;
		float nx = 0, ny = 0;
		float bz_bz_star = 0;
		for(i=0; i<NB_BINS_COMPRESSED_MATRIX_f0; i++){
		//==============================================
		// BP1 PSD == B PAR_LFR_SC_BP1_PE_FL0 == 16 bits
		PSDB = compressed_spectral_matrix_f0[i*30] // S11
		+ compressed_spectral_matrix_f0[i*30+10] // S22
		+ compressed_spectral_matrix_f0[i*30+18]; // S33
		significand = frexp(PSDB, &exponent);
		pt_char = (unsigned char*) &PSDB;
		LFR_BP1_F0[i*9+8] = pt_char[0]; // bits 31 downto 24 of the float
		LFR_BP1_F0[i*9+7] = pt_char[1]; // bits 23 downto 16 of the float
		//==============================================
		// BP1 PSD == E PAR_LFR_SC_BP1_PB_FL0 == 16 bits
		PSDE = compressed_spectral_matrix_f0[i30+24] k44 // S44
		+ compressed_spectral_matrix_f0[i30+28] k55 // S55
		+ compressed_spectral_matrix_f0[i30+26] k45_re // S45
		- compressed_spectral_matrix_f0[i30+27] k45_im; // S45
		pt_char = (unsigned char*) &PSDE;
		LFR_BP1_F0[i*9+6] = pt_char[0]; // bits 31 downto 24 of the float
		LFR_BP1_F0[i*9+5] = pt_char[1]; // bits 23 downto 16 of the float
		//==============================================================================
		// BP1 normal wave vector == PAR_LFR_SC_BP1_NVEC_V0_F0 == 8 bits
		// == PAR_LFR_SC_BP1_NVEC_V1_F0 == 8 bits
		// == PAR_LFR_SC_BP1_NVEC_V2_F0 == 1 bits
		tmp = sqrt(
		compressed_spectral_matrix_f0[i30+3]compressed_spectral_matrix_f0[i*30+3] //Im S12
		+compressed_spectral_matrix_f0[i30+5]compressed_spectral_matrix_f0[i*30+5] //Im S13
		+compressed_spectral_matrix_f0[i30+5]compressed_spectral_matrix_f0[i*30+13] //Im S23
		);
		NVEC_V0 = compressed_spectral_matrix_f0[i*30+13] / tmp; // Im S23
		NVEC_V1 = -compressed_spectral_matrix_f0[i*30+5] / tmp; // Im S13
		NVEC_V2 = compressed_spectral_matrix_f0[i*30+1] / tmp; // Im S12
		LFR_BP1_F0[i9+4] = (char) (NVEC_V0256);
		LFR_BP1_F0[i9+3] = (char) (NVEC_V1256);
		pt_char = (unsigned char*) &NVEC_V2;
		LFR_BP1_F0[i*9+2] = pt_char[0] & 0x80; // extract the sign of NVEC_V2
		//=======================================================
		// BP1 ellipticity == PAR_LFR_SC_BP1_ELLIP_F0 == 4 bits
		aux = 2*tmp / PSDB; // compute the ellipticity
		tmp_u_char = (unsigned char) (aux*(16-1)); // convert the ellipticity
		LFR_BP1_F0[i9+2] = LFR_BP1_F0[i9+2] \| ((tmp_u_char&0x0f)<<3); // keeps 4 bits of the resulting unsigned char
		//==============================================================
		// BP1 degree of polarization == PAR_LFR_SC_BP1_DOP_F0 == 3 bits
		for(j = 0; j<NB_VALUES_PER_SPECTRAL_MATRIX;j++){
		tr_SB_SB = compressed_spectral_matrix_f0[i30]compressed_spectral_matrix_f0[i*30]
		+ compressed_spectral_matrix_f0[i30+10]compressed_spectral_matrix_f0[i*30+10]
		+ compressed_spectral_matrix_f0[i30+18]compressed_spectral_matrix_f0[i*30+18]
		+ 2 * compressed_spectral_matrix_f0[i30+2]compressed_spectral_matrix_f0[i*30+2]
		+ 2 * compressed_spectral_matrix_f0[i30+3]compressed_spectral_matrix_f0[i*30+3]
		+ 2 * compressed_spectral_matrix_f0[i30+4]compressed_spectral_matrix_f0[i*30+4]
		+ 2 * compressed_spectral_matrix_f0[i30+5]compressed_spectral_matrix_f0[i*30+5]
		+ 2 * compressed_spectral_matrix_f0[i30+12]compressed_spectral_matrix_f0[i*30+12]
		+ 2 * compressed_spectral_matrix_f0[i30+13]compressed_spectral_matrix_f0[i*30+13];
		}
		aux = PSDB*PSDB;
		tmp = ( 3tr_SB_SB - aux ) / ( 2 aux );
		tmp_u_char = (unsigned char) (NVEC_V0*(8-1));
		LFR_BP1_F0[i9+2] = LFR_BP1_F0[i9+2] \| ((tmp_u_char&0x07)); // keeps 3 bits of the resulting unsigned char
		//=======================================================================================
		// BP1 z-component of the normalized Poynting flux == PAR_LFR_SC_BP1_SZ_F0 == 8 bits (7+1)
		e_cross_b_re = compressed_spectral_matrix_f0[i30+20]k34_re
		+ compressed_spectral_matrix_f0[i30+6]k14_re
		+ compressed_spectral_matrix_f0[i30+8]k15_re
		+ compressed_spectral_matrix_f0[i30+14]k24_re
		+ compressed_spectral_matrix_f0[i30+16]k25_re;
		e_cross_b_im = compressed_spectral_matrix_f0[i30+21]k34_im
		+ compressed_spectral_matrix_f0[i30+7]k14_im
		+ compressed_spectral_matrix_f0[i30+9]k15_im
		+ compressed_spectral_matrix_f0[i30+15]k24_im
		+ compressed_spectral_matrix_f0[i30+17]k25_im;
		tmp = e_cross_b_re / PSDE; // compute ReaSz
		LFR_BP1_F0[i9+1] = ((unsigned char) (tmp 128)) & 0x7f; // is it always positive?
		tmp = e_cross_b_re * e_cross_b_im;
		pt_char = (unsigned char*) &tmp;
		LFR_BP1_F0[i9+1] = LFR_BP1_F0[i9+1] \| (pt_char[0] & 0x80); // extract the sign of ArgSz
		//======================================================================
		// BP1 phase velocity estimator == PAR_LFR_SC_BP1_VPHI_F0 == 8 bits (7+1)
		nx = -sin(alpha_M)NVEC_V0 - cos(alpha_M)NVEC_V1;
		ny = NVEC_V2;
		bz_bz_star = cos(alpha_M) * cos(alpha_M) * compressed_spectral_matrix_f0[i*30] // re S11
		+ sin(alpha_M) * sin(alpha_M) * compressed_spectral_matrix_f0[i*30+10] // re S22
		- 2 * sin(alpha_M) * cos(alpha_M) * compressed_spectral_matrix_f0[i*30+2]; // re S12
		n_cross_e_scal_b_re = nx * (compressed_spectral_matrix_f0[i30+8]k15_bis_re
		+compressed_spectral_matrix_f0[i30+6]k14_bis_re
		+compressed_spectral_matrix_f0[i30+16]k25_bis_re
		+compressed_spectral_matrix_f0[i30+14]k24_bis_re)
		+ ny * (compressed_spectral_matrix_f0[i30+6]k14_tris_re
		+compressed_spectral_matrix_f0[i30+14]k24_tris_re);
		n_cross_e_scal_b_im = nx * (compressed_spectral_matrix_f0[i30+8]k15_bis_im
		+compressed_spectral_matrix_f0[i30+6]k14_bis_im
		+compressed_spectral_matrix_f0[i30+16]k25_bis_im
		+compressed_spectral_matrix_f0[i30+14]k24_bis_im)
		+ ny * (compressed_spectral_matrix_f0[i30+6]k14_tris_im
		+compressed_spectral_matrix_f0[i30+14]k24_tris_im);
		tmp = n_cross_e_scal_b_re / bz_bz_star;
		LFR_BP1_F0[i9+0] = ((unsigned char) (tmp 128)) & 0x7f; // is it always positive?
		tmp = n_cross_e_scal_b_re * n_cross_e_scal_b_im;
		pt_char = (unsigned char*) &tmp;
		LFR_BP1_F0[i9+1] = LFR_BP1_F0[i9+0] \| (pt_char[0] & 0x80); // extract the sign of ArgV
		}

		}

		void BP2_set(){
		// BP2 autocorrelation
		int i, aux = 0;
		for(i = 0; i<NB_BINS_COMPRESSED_MATRIX_f0; i++){
		// S12
		aux = sqrt(compressed_spectral_matrix_f0[i30]compressed_spectral_matrix_f0[i*30+10]);
		compressed_spectral_matrix_f0[i30+2] = compressed_spectral_matrix_f0[i30+2] / aux;
		compressed_spectral_matrix_f0[i30+3] = compressed_spectral_matrix_f0[i30+3] / aux;
		// S13
		aux = sqrt(compressed_spectral_matrix_f0[i30]compressed_spectral_matrix_f0[i*30+18]);
		compressed_spectral_matrix_f0[i30+4] = compressed_spectral_matrix_f0[i30+4] / aux;
		compressed_spectral_matrix_f0[i30+5] = compressed_spectral_matrix_f0[i30+5] / aux;
		// S23
		aux = sqrt(compressed_spectral_matrix_f0[i30+12]compressed_spectral_matrix_f0[i*30+18]);
		compressed_spectral_matrix_f0[i30+12] = compressed_spectral_matrix_f0[i30+12] / aux;
		compressed_spectral_matrix_f0[i30+13] = compressed_spectral_matrix_f0[i30+13] / aux;
		// S45
		aux = sqrt(compressed_spectral_matrix_f0[i30+24]compressed_spectral_matrix_f0[i*30+28]);
		compressed_spectral_matrix_f0[i30+26] = compressed_spectral_matrix_f0[i30+26] / aux;
		compressed_spectral_matrix_f0[i30+27] = compressed_spectral_matrix_f0[i30+27] / aux;
		// S14
		aux = sqrt(compressed_spectral_matrix_f0[i30]compressed_spectral_matrix_f0[i*30+24]);
		compressed_spectral_matrix_f0[i30+6] = compressed_spectral_matrix_f0[i30+6] / aux;
		compressed_spectral_matrix_f0[i30+7] = compressed_spectral_matrix_f0[i30+7] / aux;
		// S15
		aux = sqrt(compressed_spectral_matrix_f0[i30]compressed_spectral_matrix_f0[i*30+28]);
		compressed_spectral_matrix_f0[i30+8] = compressed_spectral_matrix_f0[i30+8] / aux;
		compressed_spectral_matrix_f0[i30+9] = compressed_spectral_matrix_f0[i30+9] / aux;
		// S24
		aux = sqrt(compressed_spectral_matrix_f0[i10]compressed_spectral_matrix_f0[i*30+24]);
		compressed_spectral_matrix_f0[i30+14] = compressed_spectral_matrix_f0[i30+14] / aux;
		compressed_spectral_matrix_f0[i30+15] = compressed_spectral_matrix_f0[i30+15] / aux;
		// S25
		aux = sqrt(compressed_spectral_matrix_f0[i10]compressed_spectral_matrix_f0[i*30+28]);
		compressed_spectral_matrix_f0[i30+16] = compressed_spectral_matrix_f0[i30+16] / aux;
		compressed_spectral_matrix_f0[i30+17] = compressed_spectral_matrix_f0[i30+17] / aux;
		// S34
		aux = sqrt(compressed_spectral_matrix_f0[i18]compressed_spectral_matrix_f0[i*30+24]);
		compressed_spectral_matrix_f0[i30+20] = compressed_spectral_matrix_f0[i30+20] / aux;
		compressed_spectral_matrix_f0[i30+21] = compressed_spectral_matrix_f0[i30+21] / aux;
		// S35
		aux = sqrt(compressed_spectral_matrix_f0[i18]compressed_spectral_matrix_f0[i*30+28]);
		compressed_spectral_matrix_f0[i30+22] = compressed_spectral_matrix_f0[i30+22] / aux;
		compressed_spectral_matrix_f0[i30+23] = compressed_spectral_matrix_f0[i30+23] / aux;
		}
		}