Method RestrictionOperator

comp/fespace.cpp:3701–3710

  
  shared_ptr<BaseMatrix> CompoundFESpace :: RestrictionOperator (int spacenr) const
  {
    shared_ptr<BaseMatrix> emb = make_shared<EmbeddingTranspose> (GetNDof(), GetRange(spacenr), IsComplex());
    if (IsParallel())
      emb = make_shared<ParallelMatrix> (emb,
                                         GetParallelDofs(),
                                         spaces[spacenr] -> GetParallelDofs(),
                                         C2C);
    return emb;
  }