- for(i=0;i<dim;i++){
- double error=work[i];
- if(bs[0]->c->q_sequencep){
- amplitude=a[i]-base;
- base=a[i];
- }else
- amplitude=a[i];
-
- meanamplitude_acc+=fabs(amplitude);
- meanamplitudesq_acc+=amplitude*amplitude;
- meanerror_acc+=fabs(error);
- meanerrorsq_acc+=error*error;
-
- if(amplitude)
- meandev_acc+=fabs(error/amplitude);
- else
- meandev_acc+=fabs(error); /* yeah, yeah */
-
- histogram_errorsq[entry*dim+i]+=error*error;
- histogram_error[entry*dim+i]+=fabs(error);
- if(histogram[entry]==0 || histogram_hi[entry*dim+i]<error)
- histogram_hi[entry*dim+i]=error;
- if(histogram[entry]==0 || histogram_lo[entry*dim+i]>error)
- histogram_lo[entry*dim+i]=error;
- distance+=error*error;
+ if(inter){
+ for(i=0;i<n/dim;i++)
+ base=process_one(b,bi,a+i,dim,n/dim,addmul[bi],base);
+ }else{
+ for(i=0;i<=n-dim;i+=dim)
+ base=process_one(b,bi,a+i,dim,1,addmul[bi],base);
+ }