namespace clarify

  aWeights = new Scalar [N];

  for (size_t i = 0; i < N; i++)

    aWeights[i] = 1.0 / N;

  Alloc2D(3, 3, &R);

  Alloc2D(N, 3, &aaXf_shifted);

  Alloc2D(N, 3, &aaXm_shifted);

  MathEigen::Alloc2D(3, 3, &R);

  MathEigen::Alloc2D(N, 3, &aaXf_shifted);

  MathEigen::Alloc2D(N, 3, &aaXm_shifted);

}

template<typename Scalar, typename ConstArrayOfCoords, typename ConstArray>

void Superpose3D<Scalar, ConstArrayOfCoords, ConstArray>::

Dealloc() {

  if (R)

    Dealloc2D(&R);

    MathEigen::Dealloc2D(&R);

  if (aWeights)

    delete [] aWeights;

  if (aaXf_shifted)

    Dealloc2D(&aaXf_shifted);

    MathEigen::Dealloc2D(&aaXf_shifted);

  if (aaXm_shifted)

    Dealloc2D(&aaXm_shifted);

    MathEigen::Dealloc2D(&aaXm_shifted);

}

// memory management: copy and move constructor, swap, and assignment operator: