Method LoadNodeType

comp/gridfunction.cpp:653–788

  template <class SCAL>  template <int N, NODE_TYPE NTYPE>
  void  S_GridFunction<SCAL> :: LoadNodeType (istream & ist, int mdcomp) 
  {
    static Timer t ("GridFunctionLoad node type"); RegionTimer r (t);    
#ifdef PARALLEL
    const auto& ma = fespace->GetMeshAccess();
    auto comm = ma->GetCommunicator();
    int id = comm.Rank();
    int ntasks = comm.Size();
    
    const FESpace & fes = *GetFESpace();
    shared_ptr<ParallelDofs> par = fes.GetParallelDofs ();
    
    if(id > 0)
      {
	int nnodes = ma->GetNNodes (NTYPE);
	
	Array<Vec<N+1, int> > nodekeys;
	Array<int> master_nodes;
	Array<DofId> dnums, pnums;
	
	for(int i = 0; i < nnodes; i++)
	  {
	    fes.GetDofNrs (NodeId(NTYPE, i),  dnums);
	    
	    if (dnums.Size() == 0) continue;
	    if (!par->IsMasterDof (dnums[0])) continue;
	    
	    master_nodes.Append(i);

	    switch (NTYPE)
	      {
	      case NT_VERTEX: pnums.SetSize(1); pnums[0] = i; break;
                // case NT_EDGE: ma->GetEdgePNums (i, pnums); break;
              case NT_EDGE: pnums = ma->GetEdgePNums (i); break;
                // case NT_FACE: ma->GetFacePNums (i, pnums); break;
              case NT_FACE: pnums = ma->GetFacePNums (i); break;
                // case NT_CELL: ma->GetElVertices (i, pnums); break;
              case NT_CELL: pnums = ma->GetElVertices (ElementId(VOL,i)); break;
	      }

	    Vec<N+1, int> key;
	    key = -1;
	    for (int j = 0; j < pnums.Size(); j++)
	      key[j] = ma->GetGlobalVertexNum (pnums[j]);
	    key[N] = dnums.Size();
	    
	    nodekeys.Append (key);	
	  }
	
	comm.Send (nodekeys, 0, 12);
	
	Array<SCAL> loc_data;
	comm.Recv (loc_data, 0, 13);
	
	for (int i = 0, cnt = 0; i < master_nodes.Size(); i++)
	  {
	    // fes.GetNodeDofNrs (NTYPE, master_nodes[i], dnums);
            fes.GetDofNrs (NodeId(NTYPE, master_nodes[i]), dnums);