Method SaveNodeType

comp/gridfunction.cpp:878–998

  template <class SCAL>  template <int N, NODE_TYPE NTYPE>
  void S_GridFunction<SCAL> :: SaveNodeType (ostream & ost, int mdcomp) const
  {
    static Timer t ("GridFunctionSave notetype"); 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> > nodenums;
	Array<SCAL> data;
    
	Array<DofId> dnums;
        Array<int> pnums;
    
	for (size_t i = 0; i < nnodes; i++)
	  {
            fes.GetDofNrs (NodeId(NTYPE,i),  dnums);
	    
	    if (dnums.Size() == 0) continue;
	    if (!par->IsMasterDof (dnums[0])) continue;      

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

	    Vec<N+1, int> points;
	    points = -1;
	    for (int j = 0; j < pnums.Size(); j++)
	      points[j] = ma->GetGlobalVertexNum (pnums[j]);
	    points[N] = dnums.Size();
	    
	    nodenums.Append(points);
	    
	    Vector<SCAL> elvec(dnums.Size()*fes.GetDimension());
	    GetElementVector (mdcomp, dnums, elvec);
	    
	    for (int j = 0; j < elvec.Size(); j++)
	      data.Append(elvec(j));
	  }    

	comm.Gather (nodenums.Size());
	comm.Gather (data.Size());
        
	comm.Send(nodenums,0,22);
	comm.Send(data,0,23);
      }