| [637] | 1 | #include "diabolo.h" | 
|---|
|  | 2 | #include "bolo.h" | 
|---|
|  | 3 | #include "acqui.h" | 
|---|
|  | 4 | #include "fichier.h" | 
|---|
|  | 5 | #include "tm.h" | 
|---|
|  | 6 | #include "tache.h" | 
|---|
|  | 7 | #include "regul.h" | 
|---|
|  | 8 | #include "auto.h" | 
|---|
|  | 9 | #include "config_transputer.h" | 
|---|
| [643] | 10 | #include "bit_block.h" | 
|---|
| [637] | 11 |  | 
|---|
|  | 12 |  | 
|---|
|  | 13 |  | 
|---|
|  | 14 |  | 
|---|
|  | 15 | void init_struct(int PCI_actif) | 
|---|
|  | 16 | { | 
|---|
|  | 17 | int i,j; | 
|---|
|  | 18 | def_long_block | 
|---|
|  | 19 | def_nom_block | 
|---|
|  | 20 | printf("memoire pour stglob: %d K  \nmemoire pour  tmtc : %d K \n", | 
|---|
|  | 21 | (sizeof(stglob)>>10),(sizeof(tmtc)>>10)); | 
|---|
|  | 22 |  | 
|---|
|  | 23 | gg=malloc(sizeof(stglob)); | 
|---|
|  | 24 |  | 
|---|
|  | 25 | random(-1); | 
|---|
|  | 26 |  | 
|---|
| [656] | 27 | if(parametr.numero_version!=version_num) {erreur("erreur numero de version");return;} | 
|---|
|  | 28 | if(parametr.n_max_bolo!=nb_max_bolo) {erreur("erreur nb_max_bolo ");return;} | 
|---|
|  | 29 | if(parametr.n_per_block!=nb_per_block) {erreur("erreur nb_per_block");return;} | 
|---|
|  | 30 | if(parametr.n_max_mes_per!=nb_max_mes_per) {erreur("erreur nb_max_mes_per");return;} | 
|---|
|  | 31 | //if(parametr.nb_bolo!=nb_max_bolo) {erreur("erreur nb_bolo_util\n");return;} | 
|---|
| [637] | 32 |  | 
|---|
| [639] | 33 | gg->reglage=reglage_standard[0]; | 
|---|
|  | 34 |  | 
|---|
| [637] | 35 | for(j=0;j<nb_max_bolo;j++) | 
|---|
|  | 36 | { | 
|---|
|  | 37 | //              parametr.bolo[j].bolo_code_util=0; | 
|---|
| [643] | 38 | //              strcpy(parametr.nom_coef[parametr.bolo[j].numero_nom_coef].bolo_nom,""); | 
|---|
| [637] | 39 | //              parametr.bolo[j].bolo_bebo=0; | 
|---|
|  | 40 | //              parametr.bolo[j].bolo_num=0; | 
|---|
|  | 41 | //              parametr.bolo[j].bolo_gain=1000; | 
|---|
|  | 42 | //              parametr.bolo[j].bolo_capa=0; | 
|---|
|  | 43 | //              parametr.bolo[j].bolo_diviseur=0; | 
|---|
|  | 44 | //              parametr.presence_bebo=0; | 
|---|
|  | 45 | //              parametr.coef[nb_max_bolo][7]  non  initialisé | 
|---|
|  | 46 |  | 
|---|
|  | 47 | gg->reglage.bolo[j].mot1=0; | 
|---|
|  | 48 | gg->reglage.bolo[j].mot2=0; | 
|---|
|  | 49 | //              gg->reglage.gainbrut[j]=0; | 
|---|
|  | 50 | //              gg->reglage.gain[j]=1; | 
|---|
|  | 51 | } | 
|---|
|  | 52 | gg->temps_cntl=0; | 
|---|
|  | 53 | gg->depou_en_cours=0; | 
|---|
|  | 54 | gg->trace_ouvert=0; | 
|---|
|  | 55 | for(j=0;j<8;j++) gg->trace_tout[j]=0; | 
|---|
|  | 56 | gg->scan_en_cours=0; | 
|---|
|  | 57 | gg->flag_corel=0; | 
|---|
|  | 58 | gg->flag_tc_reduite=0; | 
|---|
|  | 59 | gg->fichier_dia2_ouvert=0; | 
|---|
|  | 60 | gg->fichier_arch_ouvert=0; | 
|---|
|  | 61 | gg->flag_fou=0; | 
|---|
|  | 62 | gg->temps_origine=0; | 
|---|
|  | 63 |  | 
|---|
|  | 64 | //----------   initialise les structures de la tache     ------------- | 
|---|
|  | 65 | tt=malloc(sizeof(tmtc)); | 
|---|
|  | 66 |  | 
|---|
|  | 67 | tt->PCI_actif=PCI_actif; | 
|---|
|  | 68 | printf("tt->PCI_actif=%d \n",tt->PCI_actif); | 
|---|
|  | 69 |  | 
|---|
|  | 70 | tt->reglage=gg->reglage; | 
|---|
|  | 71 |  | 
|---|
|  | 72 | tt->vi.temps_cntl=0; | 
|---|
|  | 73 |  | 
|---|
|  | 74 | tt->signal_scan_immediat=0; | 
|---|
|  | 75 | // donc toujours zero pour archeops | 
|---|
|  | 76 |  | 
|---|
|  | 77 | for(j=0;j<nb_max_bolo;j++) | 
|---|
|  | 78 | { | 
|---|
|  | 79 | tt->ds.ii[j]=0; | 
|---|
|  | 80 | } | 
|---|
|  | 81 |  | 
|---|
|  | 82 |  | 
|---|
|  | 83 | tt->te.pos_ecrit=0;tt->te.pos_lit=0; | 
|---|
|  | 84 | tt->tc.pos_ecrit=0;tt->tc.pos_lit=0; | 
|---|
|  | 85 | tt->tm.pos_ecrit=0;tt->tm.pos_lit=0; | 
|---|
|  | 86 |  | 
|---|
|  | 87 |  | 
|---|
|  | 88 |  | 
|---|
|  | 89 | j=0; | 
|---|
|  | 90 | for(i=0;i<nb_type_blocks;i++) | 
|---|
|  | 91 | { | 
|---|
|  | 92 | printf(" - size=%5d   block : %s\n",long_block[i],nom_block[i]); | 
|---|
|  | 93 | if (long_block[i]> j ) j=long_block[i]; | 
|---|
|  | 94 | } | 
|---|
|  | 95 | printf(" taille maxi du block de telemesure = %d  //  on reserve %d octets par block \n", | 
|---|
|  | 96 | j,sizeof(block_type_modele)); | 
|---|
|  | 97 | if(j>sizeof(block_type_modele)) erreur("block_type_modele trops petit"); | 
|---|
|  | 98 |  | 
|---|
|  | 99 |  | 
|---|
|  | 100 |  | 
|---|
|  | 101 |  | 
|---|
|  | 102 | init_tache();   // init dependant des taches | 
|---|
|  | 103 |  | 
|---|
|  | 104 | } | 
|---|
|  | 105 |  | 
|---|
|  | 106 |  | 
|---|
|  | 107 |  | 
|---|
|  | 108 | //_________________________________________________________________________________________________________________ | 
|---|
|  | 109 | int     test_erreur(void) | 
|---|
|  | 110 | { | 
|---|
|  | 111 | int err,i; | 
|---|
|  | 112 | DateTimeRec date; | 
|---|
|  | 113 | err=tt->te.pos_ecrit-tt->te.pos_lit;    if (err<0) err+=long_table_err; | 
|---|
|  | 114 | if(err>0) | 
|---|
|  | 115 | { | 
|---|
|  | 116 | GetTime(&date); | 
|---|
|  | 117 | printf("%dh%d:%d  : %d erreurs table: ",date.hour,date.minute,date.second,err); | 
|---|
|  | 118 | for(i=0;(i<err);i++)    { | 
|---|
|  | 119 | if(i<5) switch(tt->te.err[tt->te.pos_lit]) | 
|---|
|  | 120 | { | 
|---|
|  | 121 | case    900+err_sync_repart     : printf(" sync_repart");break; | 
|---|
|  | 122 | //                      case    900+err_fifo_demi       : printf(" fifo_demi");break; | 
|---|
|  | 123 | case    900+err_fifo_pleine     : printf(" fifo_pleine");break; | 
|---|
|  | 124 | case    900+err_sync_bebo       : printf(" sync-bebo");break; | 
|---|
|  | 125 | case    900+err_sync_ctrl       : printf(" sync-ctrl");break; | 
|---|
|  | 126 | case    900+err_lect_pci        : printf(" lect_pci");break; | 
|---|
|  | 127 | case    900+err_lect_ctrl       : printf(" lect_ctrl");break; | 
|---|
|  | 128 | case    900+err_synchro_block   : printf(" synchro_block");break; | 
|---|
|  | 129 | case    900+err_recut_un_block_ok: printf(" err_recut_un_block_ok");break; | 
|---|
|  | 130 | case    900+err_longueur_block  : printf(" err_longueur_block");break; | 
|---|
|  | 131 | case    900+err_sauve_block     : printf(" err_sauve_block");break; | 
|---|
|  | 132 | case    900+err_pile_pleine     : printf(" err_pile_pleine");son(133);break; | 
|---|
|  | 133 |  | 
|---|
| [643] | 134 | case    920+erreur_bit_erreur   : printf(" bit_erreur");break; | 
|---|
|  | 135 | case    920+erreur_bit_vide     : printf(" bit_vide");break; | 
|---|
|  | 136 | case    920+erreur_type         : printf(" bit_type_erreur");break; | 
|---|
|  | 137 | case    920+erreur_longueur     : printf(" bit_erreur_longueur");break; | 
|---|
|  | 138 | case    920+erreur_somme        : printf(" bit_erreur_somme");break; | 
|---|
|  | 139 | case    920+erreur_fin          : printf(" bit_erreur_fin");break; | 
|---|
|  | 140 | case    920+cherche_debut       : printf(" bit_cherche_debut");break; | 
|---|
|  | 141 |  | 
|---|
|  | 142 |  | 
|---|
|  | 143 |  | 
|---|
| [637] | 144 | case    8       : printf(" +");break; | 
|---|
|  | 145 | case    9       : printf(" -");break; | 
|---|
|  | 146 | default : | 
|---|
|  | 147 | if(tt->te.err[tt->te.pos_lit]>=920) | 
|---|
|  | 148 | printf("bit_%d ",tt->te.err[tt->te.pos_lit]-920); | 
|---|
|  | 149 | else    printf("  %d",tt->te.err[tt->te.pos_lit]); | 
|---|
|  | 150 | break; | 
|---|
|  | 151 | } | 
|---|
|  | 152 | tt->te.pos_lit++;       if(     tt->te.pos_lit>=long_table_err) tt->te.pos_lit=0; | 
|---|
|  | 153 | } | 
|---|
|  | 154 | printf("\n"); | 
|---|
|  | 155 | } | 
|---|
|  | 156 | return(err); | 
|---|
|  | 157 | } | 
|---|
|  | 158 | //_________________________________________________________________________________________________________________ | 
|---|
|  | 159 |  | 
|---|
|  | 160 |  | 
|---|
|  | 161 |  | 
|---|
|  | 162 | // appel des bolo pour tracé des courbes individuelles | 
|---|
|  | 163 | void affi_bolo(void) | 
|---|
|  | 164 | { | 
|---|
|  | 165 | static int affi1=0,affi2=0; | 
|---|
|  | 166 | int j; | 
|---|
|  | 167 | static long t=0; | 
|---|
|  | 168 | if(!t) t=TickCount(); | 
|---|
|  | 169 | if(t>TickCount()) return; | 
|---|
|  | 170 | t=TickCount()+60;                //attente au minimum 1 sec entre 2 appels | 
|---|
|  | 171 |  | 
|---|
|  | 172 | for(j=0;j<nb_max_bolo;j++)      if(parametr.bolo[j].bolo_code_util) exec_bolo(j+1,tache_de_fond,0); | 
|---|
|  | 173 | } | 
|---|
|  | 174 |  | 
|---|
|  | 175 |  | 
|---|
|  | 176 |  | 
|---|
|  | 177 |  | 
|---|
|  | 178 |  | 
|---|
|  | 179 |  | 
|---|
|  | 180 | //***************************************************************************************** | 
|---|
|  | 181 | //                                                                                      *** | 
|---|
|  | 182 | #ifdef _archeops        //---------   pour  Archeops    ------------------------------  *** | 
|---|
|  | 183 | //                                                                                      *** | 
|---|
|  | 184 | //***************************************************************************************** | 
|---|
|  | 185 |  | 
|---|
|  | 186 |  | 
|---|
|  | 187 |  | 
|---|
|  | 188 | //_________________________________________________________________________________________________________________ | 
|---|
|  | 189 | void    reset_epld(void) | 
|---|
|  | 190 | { | 
|---|
| [649] | 191 | char mot_tc[10]; | 
|---|
| [637] | 192 | printf("reset epld"); | 
|---|
| [649] | 193 | emission_tc_reduite(tc2_reset,tc3_reset,tc4_reset); | 
|---|
|  | 194 | emission_telecommande(tc_dir_reset_epld,mot_tc); | 
|---|
|  | 195 |  | 
|---|
| [637] | 196 | } | 
|---|
|  | 197 |  | 
|---|
| [649] | 198 | void    reset_trp(void) | 
|---|
|  | 199 | { | 
|---|
|  | 200 | char mot_tc[10]; | 
|---|
|  | 201 | printf("reset trp"); | 
|---|
|  | 202 | emission_tc_reduite(tc2_reset,tc3_reset,tc4_reset); | 
|---|
|  | 203 | emission_telecommande(tc_dir_reset_trp,mot_tc); | 
|---|
|  | 204 | } | 
|---|
| [637] | 205 |  | 
|---|
| [649] | 206 |  | 
|---|
| [637] | 207 | void    calfreq(int flag_on)    //  fabrique l'etat.bolo    gg->reglage.horloge  et  gg->periode_echantillonage | 
|---|
|  | 208 | // le flag_on dit s'il faut commander la sortie des data | 
|---|
|  | 209 | { | 
|---|
|  | 210 | double p,f1,f2,f3; | 
|---|
|  | 211 | int pp,b,tm,flag; | 
|---|
|  | 212 | char mot_tc[10]; | 
|---|
|  | 213 |  | 
|---|
|  | 214 | p=litD(fenetre_acquisition,e_periode,0L); | 
|---|
|  | 215 | if(p<4) p=4;    if(p>16) p=16;  //   butée sur la periode | 
|---|
|  | 216 | pp=p*5.;  p=(double)pp / 5.; | 
|---|
|  | 217 | ecritD(fenetre_acquisition,e_periode,"%6.1f",p);        // reconstruit p | 
|---|
|  | 218 |  | 
|---|
|  | 219 | b=litD(fenetre_acquisition,e_nb_mes,0L); | 
|---|
|  | 220 | if(b<18) {b=18;ecritD(fenetre_acquisition,e_nb_mes,b);} | 
|---|
|  | 221 | if(b>44) {b=44;ecritD(fenetre_acquisition,e_nb_mes,b);} | 
|---|
|  | 222 | b=b*2; | 
|---|
|  | 223 |  | 
|---|
|  | 224 |  | 
|---|
|  | 225 | tm=litD(fenetre_acquisition,e_temp_mort,0L);    // temps mort = 1 .. nbmes par 1/4 per | 
|---|
|  | 226 | if(tm<2) {tm=2;ecritD(fenetre_acquisition,e_temp_mort,tm);} | 
|---|
|  | 227 | if(tm>13) {tm=13;ecritD(fenetre_acquisition,e_temp_mort,tm);} | 
|---|
|  | 228 |  | 
|---|
|  | 229 |  | 
|---|
|  | 230 | flag=flag_on; | 
|---|
|  | 231 | if(litD(fenetre_acquisition,e_flag1,0L))        flag +=2; | 
|---|
|  | 232 | if(litD(fenetre_acquisition,e_flag2,0L))        flag +=4; | 
|---|
|  | 233 |  | 
|---|
|  | 234 | //printf(" b=%d  tm=%d  a=%d   h=%d \n",b,tm,a,h); | 
|---|
|  | 235 |  | 
|---|
|  | 236 | if(gg->flag_tc_reduite) | 
|---|
|  | 237 | { | 
|---|
|  | 238 |  | 
|---|
|  | 239 | if((pp!=gg->reglage.horloge.periode) ou clavier(touche_alt) )    emission_tc_reduite(tc2_horloge,tc3_periode,pp); | 
|---|
|  | 240 | if((b!=gg->reglage.horloge.nb_mesures) ou clavier(touche_alt) )  emission_tc_reduite(tc2_horloge,tc3_nb_mesures,b); | 
|---|
|  | 241 | if((tm!=gg->reglage.horloge.temp_mort) ou clavier(touche_alt) )  emission_tc_reduite(tc2_horloge,tc3_temp_mort,tm); | 
|---|
| [639] | 242 | if(flag ou (flag!=gg->reglage.horloge.flag) ou clavier(touche_alt) )  emission_tc_reduite(tc2_horloge,tc3_flag,flag); | 
|---|
| [637] | 243 | } | 
|---|
|  | 244 |  | 
|---|
|  | 245 | else | 
|---|
|  | 246 | { | 
|---|
|  | 247 | mot_tc[0]=tc_horloge;           //      commande valeur des horloges | 
|---|
|  | 248 | mot_tc[1]=pp;                   //      periode en multiple de 0.2µsec : de 20 a 80  soit  4 a 16 µsec | 
|---|
|  | 249 | mot_tc[2]=b/2 ;                 //      nb pts/demi periode de 23 a 44 | 
|---|
|  | 250 | mot_tc[3]=tm; | 
|---|
|  | 251 | mot_tc[4]=flag;                 //      c'est le flag qui met en route la sortie des mesures | 
|---|
|  | 252 | emission_telecommande(tc_dir_transputer,mot_tc); | 
|---|
|  | 253 |  | 
|---|
|  | 254 | gg->reglage.horloge.periode=pp; | 
|---|
|  | 255 | gg->reglage.horloge.nb_mesures=b;                       //  nombre de mesures dans une periode complete | 
|---|
|  | 256 | gg->reglage.horloge.temp_mort=tm; | 
|---|
|  | 257 | gg->reglage.horloge.flag=flag; | 
|---|
|  | 258 |  | 
|---|
|  | 259 | } | 
|---|
|  | 260 | // affichage de toutes les fréquences | 
|---|
|  | 261 | f1=1000/p;f2=f1/bitmot;f3=f2*1000./(double)(gg->reglage.horloge.nb_mesures); | 
|---|
|  | 262 | ecritD(fenetre_acquisition,e_freq,"freq: %5.1f kHz  --  %5.2f kHz,    modulation: %5.1f Hz",f1,f2,f3); | 
|---|
|  | 263 | printf("freq --> : %5.1f kHz  --  %5.2f kHz,    modulation: %5.1f Hz\n",f1,f2,f3); | 
|---|
|  | 264 | gg->periode_echantillonage=0.5/f3;                                      //  2 fois la frequence de modulation | 
|---|
|  | 265 | printf(" periode d'échantillonage= %g \n",gg->periode_echantillonage); | 
|---|
|  | 266 |  | 
|---|
|  | 267 | tt->reglage=gg->reglage; | 
|---|
|  | 268 | fich_save_partiel(1); | 
|---|
|  | 269 | } | 
|---|
|  | 270 |  | 
|---|
|  | 271 |  | 
|---|
|  | 272 | #endif | 
|---|
|  | 273 | //***************************************************************************************** | 
|---|
|  | 274 | //                                                                                      *** | 
|---|
|  | 275 | #ifdef _diabolo                 //---------   pour  Diabolo    ------------------------------   *** | 
|---|
|  | 276 | //                                                                                      *** | 
|---|
|  | 277 | //***************************************************************************************** | 
|---|
|  | 278 |  | 
|---|
|  | 279 | void    reset_epld(void) | 
|---|
|  | 280 | { | 
|---|
|  | 281 | } | 
|---|
|  | 282 |  | 
|---|
|  | 283 | void    calfreq(int flag_on)                            //  fabrique l'etat.bolo    gg->reglage.horloge.nb_mesures  et  gg->periode_echantillonage | 
|---|
|  | 284 | { | 
|---|
|  | 285 | double p,f1,f2,f3; | 
|---|
|  | 286 | int a,h,pp,b,tm; | 
|---|
|  | 287 | int j; | 
|---|
|  | 288 | char mot_tc[10]; | 
|---|
|  | 289 |  | 
|---|
|  | 290 | p=litD(fenetre_acquisition,e_periode,0L); | 
|---|
|  | 291 | b=litD(fenetre_acquisition,e_nb_mes,0L); | 
|---|
|  | 292 | tm=litD(fenetre_acquisition,e_temp_mort,0L);    // temps mort = 1 .. nbmes par 1/4 per | 
|---|
|  | 293 | if(tm<1) {tm=1;ecritD(fenetre_acquisition,e_temp_mort,tm);} | 
|---|
|  | 294 | if(tm>b) {tm=b;ecritD(fenetre_acquisition,e_temp_mort,b);} | 
|---|
|  | 295 | //  lecture externe ou local | 
|---|
|  | 296 | h=litD(fenetre_acquisition,e_horloge_ext_,0L); | 
|---|
|  | 297 |  | 
|---|
|  | 298 | //if(p<5) p=5;                          //    butée en periode minimum pour éviter de planter le Mac | 
|---|
|  | 299 |  | 
|---|
|  | 300 |  | 
|---|
|  | 301 | //----------------------------------------                              nombre de mesures | 
|---|
|  | 302 |  | 
|---|
|  | 303 | if(b<4) b=4;    if(b>44) b=44; | 
|---|
|  | 304 | gg->reglage.horloge.nb_mesures=(b/2)*4; | 
|---|
|  | 305 |  | 
|---|
|  | 306 | //----------------------------------------                              periode de l'horloge de base | 
|---|
|  | 307 | if (h&1)        p=p*frequ_horloge_externe; | 
|---|
|  | 308 | else            p=p*frequ_horloge_interne; | 
|---|
|  | 309 |  | 
|---|
|  | 310 | a=(p/frequ_divise)-ofset_periode ; | 
|---|
|  | 311 |  | 
|---|
|  | 312 | if(a<0) a=0; | 
|---|
|  | 313 | if(a>1020) a=1020;                      // 10 bits maxi  attention a l'ofset | 
|---|
|  | 314 |  | 
|---|
|  | 315 | pp=a+ofset_periode; | 
|---|
|  | 316 | p=(double)(pp)*frequ_divise; | 
|---|
|  | 317 | if (h&1)         p=p/frequ_horloge_externe; | 
|---|
|  | 318 | else            p=p/frequ_horloge_interne; | 
|---|
|  | 319 |  | 
|---|
|  | 320 | ecritD(fenetre_acquisition,e_periode,"%6.1f",p);  // reconstruit p et pp | 
|---|
|  | 321 |  | 
|---|
|  | 322 |  | 
|---|
|  | 323 | //printf(" b=%d  tm=%d  a=%d   h=%d \n",b,tm,a,h); | 
|---|
|  | 324 |  | 
|---|
|  | 325 | mot_tc[0]=tc_horloge;                   //79..72 | 
|---|
|  | 326 | mot_tc[1]=(tm%7)+ ((b/2)<<3) ;  //71..64 3 bit temps mort et 5 bit  mb (nombre de coups par 1/4 de periode = b/2) | 
|---|
|  | 327 | mot_tc[2]=(a>>8);                       //63..56 poid fort (2 bits) | 
|---|
|  | 328 | mot_tc[3]=a;                            //55..48 poid faible (8 bits) | 
|---|
|  | 329 | mot_tc[4]=h&3;                  //47..40 horloge externe et bolo-iram | 
|---|
|  | 330 | emission_telecommande(0,mot_tc); | 
|---|
|  | 331 |  | 
|---|
|  | 332 |  | 
|---|
|  | 333 | // affichage de toutes les fréquences | 
|---|
|  | 334 | f1=1000/p;f2=f1/bitmot;f3=f2*1000./(double)(gg->reglage.horloge.nb_mesures); | 
|---|
|  | 335 | gg->periode_echantillonage=0.5/f3;              //  2 fois la frequence de modulation | 
|---|
|  | 336 | // periode * 16 * nbdemesure / 2 | 
|---|
|  | 337 | //printf(" durée= %ld \n",gg->periode_echantillonage); | 
|---|
|  | 338 |  | 
|---|
|  | 339 | ecritD(fenetre_acquisition,e_freq,"freq: %5.1f kHz  --  %5.2f kHz,    modulation: %5.1f Hz",f1,f2,f3); | 
|---|
|  | 340 |  | 
|---|
|  | 341 |  | 
|---|
|  | 342 |  | 
|---|
|  | 343 | {//_______________    fabrique les detections synchrones de la modulation électrique | 
|---|
|  | 344 | int i,v,tm; | 
|---|
|  | 345 | tm=litD(fenetre_acquisition,e_temp_mort,0L); | 
|---|
|  | 346 | gg->reglage.temp_mort=tm; | 
|---|
|  | 347 | // detection plus/moins avec tm  points supprimé | 
|---|
|  | 348 | //              ds[0] ==  plus / moins | 
|---|
|  | 349 | //              ds[1] ==  plus / plus | 
|---|
|  | 350 | //              ds[2] ==  + - 0 /  - + 0 | 
|---|
|  | 351 |  | 
|---|
|  | 352 | tt->ds.dsnorme=0; | 
|---|
|  | 353 | i=0; | 
|---|
|  | 354 | for(v=0 ;v<gg->reglage.horloge.nb_mesures/2     ;v++) | 
|---|
|  | 355 | { | 
|---|
|  | 356 | if(v<tm)        tt->ds.dstable[i]=0; | 
|---|
|  | 357 | else            {tt->ds.dstable[i]=1;   tt->ds.dsnorme++;} | 
|---|
|  | 358 | i++; | 
|---|
|  | 359 | } | 
|---|
|  | 360 | for(v=0;v<gg->reglage.horloge.nb_mesures/2      ;v++) | 
|---|
|  | 361 | { | 
|---|
|  | 362 | if(v<tm)        tt->ds.dstable[i]=0; | 
|---|
|  | 363 | else            {tt->ds.dstable[i]=-1;  tt->ds.dsnorme++;} | 
|---|
|  | 364 | i++; | 
|---|
|  | 365 | } | 
|---|
|  | 366 | } | 
|---|
|  | 367 | //printf(" normes = %d  %d \n",tt->ds.dsnorme[0],tt->ds.dsnorme[1]); | 
|---|
|  | 368 | tt->reglage=gg->reglage; | 
|---|
|  | 369 | fich_save_partiel(1); | 
|---|
|  | 370 | } | 
|---|
|  | 371 |  | 
|---|
|  | 372 |  | 
|---|
|  | 373 | //***************************************************************************************** | 
|---|
|  | 374 | //                                                                                      *** | 
|---|
|  | 375 | #endif          //--------------------------------------------------------------------  *** | 
|---|
|  | 376 | //                                                                                      *** | 
|---|
|  | 377 | //***************************************************************************************** | 
|---|
|  | 378 |  | 
|---|
|  | 379 |  | 
|---|
|  | 380 |  | 
|---|
|  | 381 |  | 
|---|
|  | 382 | void    deglitch(void) | 
|---|
|  | 383 | { | 
|---|
|  | 384 | int j; | 
|---|
|  | 385 | double x,dx; | 
|---|
|  | 386 | static int debut=1; | 
|---|
|  | 387 | static double X[nb_max_bolo]; | 
|---|
|  | 388 | static double DX[nb_max_bolo]; | 
|---|
|  | 389 | if(debut) | 
|---|
|  | 390 | { | 
|---|
|  | 391 | for(j=0;j<nb_max_bolo;j++) | 
|---|
|  | 392 | { | 
|---|
|  | 393 | X[j]=0;DX[j]=1000; | 
|---|
|  | 394 | } | 
|---|
|  | 395 | debut=0; | 
|---|
|  | 396 | } | 
|---|
|  | 397 |  | 
|---|
|  | 398 | for(j=0;j<nb_max_bolo;j++) | 
|---|
|  | 399 | { | 
|---|
|  | 400 | x=gg->don.don_bolo[j]; | 
|---|
|  | 401 | dx=x-X[j]; if(dx<0)  dx=-dx; | 
|---|
|  | 402 | if( dx > DX[j] ) | 
|---|
|  | 403 | { | 
|---|
|  | 404 | printf("glitch bolo %d  ancien X=%g  nouveau x=%g  DX=%g \n",j+1,X[j],x,DX[j]); | 
|---|
|  | 405 | gg->don.don_bolo[j]=X[j]; | 
|---|
|  | 406 | DX[j] *=1.1; | 
|---|
|  | 407 | } | 
|---|
|  | 408 | else    { | 
|---|
|  | 409 | DX[j]=0.9 * DX[j] + 0.5 * dx +1 ; | 
|---|
|  | 410 | X[j]=x; | 
|---|
|  | 411 | }               // en moyenne, 5 fois l'interval moyen | 
|---|
|  | 412 | } | 
|---|
|  | 413 |  | 
|---|
|  | 414 | } | 
|---|
|  | 415 |  | 
|---|
|  | 416 |  | 
|---|
|  | 417 |  | 
|---|
|  | 418 |  | 
|---|
|  | 419 |  | 
|---|
|  | 420 |  | 
|---|