template <typename T>
  inline
  MSpline<T>::MSpline(const DVector<Vector<T,3>> &c, int d) {
    _d = d;
    _makeKnotVector(c.getDim());
    _C = c;
    for (int i=0;i<_C.getDim();i++){
        Selector<T,3>* s = new Selector<T,3>(_C[i],i,this);
        this->insert(s);

    }
    auto sk = new SelectorGridVisualizer<T>;
    sk->setSelectors(_C,0,isClosed());
    this->insertVisualizer(sk);

  }