// Copyright 2010 The Code Cavern

void	bda_adddiag(bda_ptr z,bda_srcptr x,bda_len n)
{bda_len i;bda_digit h,l;bda_bit c=0;

for(i=0;i<n;i++)
   {bda_digit_mul(h,l,x[i],x[i]);
    bda_digit_adc(c,z[2*i],z[2*i],l,c);
    bda_digit_adc(c,z[2*i+1],z[2*i+1],h,c);
   }
return;}