for (i=1;i<l;i++) /* the Average Magnitude Difference Function */
    {
      for (j=start;j<=end;j++)
	  {
	    temp[j] = i + j < l ? in[i+j] : 0.0;
	    temp0 = atom_getfloatarg(i, 4096, ctl->otemp);
	    temp0 += i == 0 ? 0.0 : fabs(in[j] - temp[j]);
	  }
      temp0 += ((float)i / (float)l) * ctl->f_sum_abs;
      SETFLOAT(&ctl->otemp[i], temp0);
    }