Add an assertation

Note: This assumes that any vector used to construct a BasisGridFunction
already has the proper size. If you want to be able to write

  Vector vec; // unitialised
  BasisGridFunction func(vec);
  initialize(vec);

for some reason, then please revert this change and let me know.