double *add_3d_vectors1(double a[3], double b[3])
{
  double sum[3];
  sum[0] = a[0] + b[0];
  sum[1] = a[1] + b[1];
  sum[2] = a[2] + b[2];
  return sum;   /* DON'T DO THIS!!! */
}