void uvadd(int n, double *a, double *b, double *c)
{
   int i;
   for (i=0; i<=n-1; i++)
      c[i]=a[i]+b[i];
}   

void uvsub(int n, double *a, double *b, double *c)
{
   int i;
   for (i=0; i<=n-1; i++)
      c[i]=a[i]-b[i];
}   

void uvsadd(int n, double *a, double s, double *b, double *c)
{
   int i;
   for (i=0; i<=n-1; i++)
      c[i]=s*a[i]+b[i];
}