#include "petscvec.h"
#include <complex>
#include <vector>

template<typename T>
class VectorP
{
  public:
  void dot_product(VectorP<T>& vecA, VectorP<T>& vecB, T& prod)
  {
    VecDot(vecA.x, vecB.x, &prod);
  }

  private:
  Vec x;
};

int main()
{
  VectorP<std::complex<double> > a, b;
  std::vector<PetscScalar> p(2, 0.0);


  a.dot_product(a, b, p[1]);
  return 0;
}