INNER CODE UNIT · C
avg
aquaskyline/SOAPdenovo2 · fusion/bundle.c:262
long long int avg=sum/weight;
sum = 0;
for(ii=0;ii<weight;ii++){
sum+=((avg-(long long int)tmp->PE[ii])*(avg-(long long int)tmp->PE[ii]));
}
double SD=(sqrt((double)sum/(weight-1)))*3;//just for fast
sum=0;
int num=0;
for(ii=0;ii<weight;ii++){
if(abs(tmp->PE[ii]-avg)<=SD){
sum+=tmp->PE[ii];
num++;
}else{
//ignorePE5++;
counter++;
}
}