projects
/
physik
/
nlsop.git
/ commitdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
| commitdiff |
tree
raw
|
patch
|
inline
| side by side (parent:
f8eab70
)
calc average
author
hackbard
<hackbard>
Thu, 30 Sep 2004 14:38:05 +0000
(14:38 +0000)
committer
hackbard
<hackbard>
Thu, 30 Sep 2004 14:38:05 +0000
(14:38 +0000)
parse_trim_collision.c
patch
|
blob
|
history
diff --git
a/parse_trim_collision.c
b/parse_trim_collision.c
index
5469bfa
..
eed12c5
100644
(file)
--- a/
parse_trim_collision.c
+++ b/
parse_trim_collision.c
@@
-30,6
+30,7
@@
int main(int argc,char **argv) {
int fd,ion,i,j;
int z,skipped;
int fd,ion,i,j;
int z,skipped;
+ int avg1,avg2;
double en;
double max_e;
int max_chfi,max_th;
double en;
double max_e;
int max_chfi,max_th;
@@
-49,6
+50,8
@@
int main(int argc,char **argv) {
ion=1;
skipped=0;
ion=1;
skipped=0;
+ avg1=0;
+ avg2=0;
for(i=0;i<MAXZ;i++) {
Z[i].cell_hits_from_ions=0;
Z[i].total_hits=0;
for(i=0;i<MAXZ;i++) {
Z[i].cell_hits_from_ions=0;
Z[i].total_hits=0;
@@
-73,12
+76,15
@@
int main(int argc,char **argv) {
if(Z[i].energy>max_e) max_e=Z[i].energy;
if(Z[i].total_hits>max_th) max_th=Z[i].total_hits;
if(Z[i].cell_hits_from_ions>max_chfi) max_chfi=Z[i].cell_hits_from_ions;
if(Z[i].energy>max_e) max_e=Z[i].energy;
if(Z[i].total_hits>max_th) max_th=Z[i].total_hits;
if(Z[i].cell_hits_from_ions>max_chfi) max_chfi=Z[i].cell_hits_from_ions;
+ avg1+=Z[i].total_hits;
+ avg2+=Z[i].cell_hits_from_ions;
}
for(i=0;i<MAXZ;i++)
printf("%d %f %f %f\n",i*3,1.0*Z[i].total_hits/max_th,
Z[i].energy/max_e,
1.0*Z[i].cell_hits_from_ions/max_chfi);
printf("skipped = %d\n",skipped);
}
for(i=0;i<MAXZ;i++)
printf("%d %f %f %f\n",i*3,1.0*Z[i].total_hits/max_th,
Z[i].energy/max_e,
1.0*Z[i].cell_hits_from_ions/max_chfi);
printf("skipped = %d\n",skipped);
+ printf("average hits per ion %d / %d\n",avg1/ion,avg2/ion);
return 1;
}
if(buf[i]=='\n') break;
return 1;
}
if(buf[i]=='\n') break;