Discussion :

[AJ_ing]

Bonjour,
je veux calculer le produit des données d'une colonne dans un fichier texte , le format de chaque ligne est le suivant :

P( ..... ) = 0.000802232 logprob = -3.095700 bo_case = 1
P(... ) = 0.25787 logprob = -0.588600 bo_case = 2
P( ...) = 0.000287343 logprob = -3.541600 bo_case = 3-2
P( ... ) = 0.0414477 logprob = -1.382500 bo_case = 3
P(... ) = 0.0880846 logprob = -1.055100 bo_case = 3

je veux calculer le produit de toute les valeurs des lignes se trouvant entre = et logprob
et retourner à chaque fois le nom de fichier suivi de ce produit.

J'ai commencé par extraire cette valeur et la mettre dans un fichier résultat mais ça donne un fichier résultat vide

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
for fic in *.txt 
do 
while read ligne
do       
       val=$(echo "$ligne" | cut -d"=" -f2 | cut -d" " -f1)
       echo "$fich $val" > resultat.txt  
done < $fic
done

merci d'avance

[N_BaH]

tu peux confirmer, s'il te plaît, que tu veux le produit, et non la somme ?

[AJ_ing]

Oui le produit, mais je vois que vu que la première colonne va donner presque des 0,
le produit ou la somme me donne même chose en terme d'information

[frp31]

utilises DC pour faire ça c'est très simple

exemple :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
$ dc
3 4 * p
12
q
$

suffit donc d'afficher un echo du calcul à faire un pipe dans DC pour que ça marche

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
 
$ cd /tmp/d
$ cat > /tmp/toto                                                       
P( ..... ) = 0.000802232 logprob = -3.095700 bo_case = 1
P(... ) = 0.25787 logprob = -0.588600 bo_case = 2
P( ...) = 0.000287343 logprob = -3.541600 bo_case = 3-2
P( ... ) = 0.0414477 logprob = -1.382500 bo_case = 3
P(... ) = 0.0880846 logprob = -1.055100 bo_case = 3
 
$ echo '32k 1 '$(sed 's/.*=\(.*\)logprob.*/\1*/' /tmp/toto ) 'pq' | dc  
.00000000021702104269483793366184
$

le 32k au début est la précision demandée, là j'ai mis 32, mais tu peux utiliser le nombre de décimales que tu veux bien sur.

1 est le premier chiffre pour empêcher d'avoir un stack vide au début du calcul, puisque DC travail en mode de calcul RPN. L'avantage de cette méthode c'est que tous les unix ont DC

la sous commande sed sert à sélectionner la colonne qui nous intéresse...
c'est donc comme si on tapait à la main ceci dans DC

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

 32k 1  0.000802232 * 0.25787 * 0.000287343 * 0.0414477 * 0.0880846 * pq

Quand on a plusieurs fichiers de calculs à traiter il suffit de générer le texte complet de l’opération totale dans un fichier et de faire interpréter ce dernier par DC

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
 
 
         $ cp /tmp/toto /tmp/tutu                                                                                                                                                               
$ cp /tmp/toto /tmp/titi
$ echo '32k 1 ' > /tmp/operation ; for i in /tmp/t?t? ; do echo $(sed 's/.*=\(.*\)logprob.*/\1*/' $i ) >> /tmp/operation ; done ; echo 'pq' >> /tmp/operation ; cat /tmp/operation     
32k 1 
0.000802232 * 0.25787 * 0.000287343 * 0.0414477 * 0.0880846 *
0.000802232 * 0.25787 * 0.000287343 * 0.0414477 * 0.0880846 *
0.000802232 * 0.25787 * 0.000287343 * 0.0414477 * 0.0880846 *
pq

on voit bien que l'opération est bien complête avec les 3 fichiers /tmp/toto /tmp/titi et /tmp/tutu

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
$ echo '32k 1 ' > /tmp/operation ; for i in /tmp/t?t? ; do echo $(sed 's/.*=\(.*\)logprob.*/\1*/' $i ) >> /tmp/operation ; done ; echo 'pq' >> /tmp/operation ; cat /tmp/operation | dc
.00000000000000000000000000001022
$

le calcul est bien bon puisqu'on a bien le résultat attendu.

[AJ_ing]

j'ai aussi ce message d'eureur :

dc: 't' (0164) unimplemented
dc: 't' (0164) unimplemented
.310456
dc: output base must be a number greater than 1
dc: 'b' (0142) unimplemented
dc: value overflows simple integer; punting...
dc: �'g' (0147) unimplemented
dc: stack empty
dc: 't' (0164) unimplemented
dc: 't' (0164) unimplemented
0

j'ai mis le même code :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
for fic in *.txt 
do
  echo '32k 1 '$(sed 's/.*=\(.*\)logprob.*/\1*/' $fic ) 'pq' | dc
  #val=$( { { echo 'scale=20;' ; cat $fic ; } | sed -e 's|.*= *\(.*\) logprob.*|\1 *|' | tr '\n' ' ' ; echo '1' ; } | bc -l)
  #echo "$fich $val" 
done

je pense qu'il vaut mieux calculer la somme, ça donne même information comme valeur
Si je fait la somme c'est comment ?

[frp31]

oui j'ai oublié de préciser qu'il faut lancer ça depuis un repertoire vide à cause de sed qui risque d'interpreter le caractère "*" au lieu de l'utiliser comme caractère.

[frp31]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
$ echo '32k 1 ' > /tmp/operation ; for i in /tmp/t?t? ; do echo $(sed 's/.*=\(.*\)logprob.*/\1+/' $i ) >> /tmp/operation ; done ; echo 'pq' >> /tmp/operation ; cat /tmp/operation     
32k 1 
0.000802232 + 0.25787 + 0.000287343 + 0.0414477 + 0.0880846 +
0.000802232 + 0.25787 + 0.000287343 + 0.0414477 + 0.0880846 +
0.000802232 + 0.25787 + 0.000287343 + 0.0414477 + 0.0880846 +
pq
$ echo '32k 1 ' > /tmp/operation ; for i in /tmp/t?t? ; do echo $(sed 's/.*=\(.*\)logprob.*/\1+/' $i ) >> /tmp/operation ; done ; echo 'pq' >> /tmp/operation ; cat /tmp/operation | dc
2.165475625
$

pour la somme c'est pareil avec un "+"

[jack-ft]

Quelque chose comme ça?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
 
for fic in *.txt ; do
  val=$( { { echo 'scale=20;' ; cat $fic ; } | sed -e 's|.*= *\(.*\) logprob.*|\1 *|' | tr '\n' ' ' ; echo '1' ; } | bc -l)
  echo "$fich $val"
done

Avec le fichier fourni en exemple, on obtient:
.00000000021702104269

)jack(