optimize calculation speed

solaoi · solaoi · commit 1a41807cc5ed · 2022-03-28T19:25:24.000+09:00
diff --git a/colc.ts b/colc.ts
@@ -118,7 +118,7 @@ if (binSize === null) {
     }
     bash.push("| sort -n | awk");
     bash.push(
-      `'BEGIN{OFMT="%.6f"}NR==1{min=$1}{sum+=$1;d[NR]=$1}END{avg=sum/NR;for(i in d)s+=(d[i]-avg)^2;stddev=sqrt(s/(NR-1));q1=(3*d[int((NR-1)/4)+1]+d[int((NR-1)/4)+2])/4;q3=(d[int(3*(NR-1)/4)+1]+3*d[int(3*(NR-1)/4)+2])/4;iqr=q3-q1;stur=1+log(NR)/log(2);sturi=int(stur);sturges=stur>sturi?sturi+1:sturi;max=d[NR];print stddev,avg,sum,NR,max,min,sqrt(s/(NR-1))/sqrt(NR),s/(NR-1),(NR%2)?d[(NR+1)/2]:(d[NR/2]+d[NR/2+1])/2,avg+stddev,avg-stddev,avg+2*stddev,avg-2*stddev,avg+3*stddev,avg-3*stddev,(max-min)/sturges,(3.5*stddev)/(exp(log(NR)/3)),q1,q3,iqr,q1-1.5*iqr,q3+1.5*iqr,(max-min)/sqrt(NR),2*iqr/(exp(log(NR)/3))}'`,
+      `'BEGIN{OFMT="%.6f"}NR==1{min=$1}{sum+=$1;d[NR]=$1}END{avg=sum/NR;for(i in d)s+=(d[i]-avg)^2;stddev=sqrt(s/(NR-1));q1=(3*d[int((NR-1)/4)+1]+d[int((NR-1)/4)+2])/4;q3=(d[int(3*(NR-1)/4)+1]+3*d[int(3*(NR-1)/4)+2])/4;iqr=q3-q1;stur=1+log(NR)/log(2);sturi=int(stur);sturges=stur>sturi?sturi+1:sturi;max=d[NR];range=max-min;sqrtnr=sqrt(NR);threerootnr=exp(log(NR)/3);print stddev,avg,sum,NR,max,min,sqrt(s/(NR-1))/sqrtnr,s/(NR-1),(NR%2)?d[(NR+1)/2]:(d[NR/2]+d[NR/2+1])/2,avg+stddev,avg-stddev,avg+2*stddev,avg-2*stddev,avg+3*stddev,avg-3*stddev,range/sturges,(3.5*stddev)/threerootnr,q1,q3,iqr,q1-1.5*iqr,q3+1.5*iqr,range/sqrtnr,2*iqr/threerootnr}'`,
     );
     return bash.join(" ");
   })();

Original file line number	Diff line number	Diff line change
`@@ -118,7 +118,7 @@ if (binSize === null) {`
`118`	`118`	`}`
`119`	`119`	`bash.push("\| sort -n \| awk");`
`120`	`120`	`bash.push(`
`121`		- `'BEGIN{OFMT="%.6f"}NR==1{min=$1}{sum+=$1;d[NR]=$1}END{avg=sum/NR;for(i in d)s+=(d[i]-avg)^2;stddev=sqrt(s/(NR-1));q1=(3d[int((NR-1)/4)+1]+d[int((NR-1)/4)+2])/4;q3=(d[int(3(NR-1)/4)+1]+3d[int(3(NR-1)/4)+2])/4;iqr=q3-q1;stur=1+log(NR)/log(2);sturi=int(stur);sturges=stur>sturi?sturi+1:sturi;max=d[NR];print stddev,avg,sum,NR,max,min,sqrt(s/(NR-1))/sqrt(NR),s/(NR-1),(NR%2)?d[(NR+1)/2]:(d[NR/2]+d[NR/2+1])/2,avg+stddev,avg-stddev,avg+2stddev,avg-2stddev,avg+3stddev,avg-3stddev,(max-min)/sturges,(3.5stddev)/(exp(log(NR)/3)),q1,q3,iqr,q1-1.5iqr,q3+1.5iqr,(max-min)/sqrt(NR),2iqr/(exp(log(NR)/3))}'`,
	`121`	+ `'BEGIN{OFMT="%.6f"}NR==1{min=$1}{sum+=$1;d[NR]=$1}END{avg=sum/NR;for(i in d)s+=(d[i]-avg)^2;stddev=sqrt(s/(NR-1));q1=(3d[int((NR-1)/4)+1]+d[int((NR-1)/4)+2])/4;q3=(d[int(3(NR-1)/4)+1]+3d[int(3(NR-1)/4)+2])/4;iqr=q3-q1;stur=1+log(NR)/log(2);sturi=int(stur);sturges=stur>sturi?sturi+1:sturi;max=d[NR];range=max-min;sqrtnr=sqrt(NR);threerootnr=exp(log(NR)/3);print stddev,avg,sum,NR,max,min,sqrt(s/(NR-1))/sqrtnr,s/(NR-1),(NR%2)?d[(NR+1)/2]:(d[NR/2]+d[NR/2+1])/2,avg+stddev,avg-stddev,avg+2stddev,avg-2stddev,avg+3stddev,avg-3stddev,range/sturges,(3.5stddev)/threerootnr,q1,q3,iqr,q1-1.5iqr,q3+1.5iqr,range/sqrtnr,2iqr/threerootnr}'`,
`122`	`122`	`);`
`123`	`123`	`return bash.join(" ");`
`124`	`124`	`})();`