Method BlockVector

linalg/basevector.cpp:799–823

  BlockVector :: BlockVector (const Array<shared_ptr<BaseVector>> & avecs)
    : vecs(avecs), ispar(vecs.Size())
  {
    size = 0;
    for (auto & vec:vecs)
      size += vec->Size();
    // #ifdef PARALLEL
    ispar.Clear();
    for (size_t k = 0; k < vecs.Size(); k++)
      if (vecs[k]->GetParallelStatus() != NOT_PARALLEL)
        {
          ispar.SetBit(k);
          comm = dynamic_pointer_cast<ParallelBaseVector> (vecs[k])->GetParallelDofs()->GetCommunicator();
        }
    /*
    for (size_t k = 0; k<vecs.Size(); k++) {
      auto stat = vecs[k]->GetParallelStatus();
      if ( stat==NOT_PARALLEL ) continue;
      ispar.SetBit(k);
      auto * pv = dynamic_cast_ParallelBaseVector(vecs[k].get());
      comm = pv->GetParallelDofs()->GetCommunicator();
    }
    */
    // #endif
  }

  void * BlockVector :: Memory () const
  { throw Exception("BlockVector::Memory is not useful"); }