SundanceSpatiallyHomogeneousDOFMapBase.cpp

Go to the documentation of this file.

/* @HEADER@ */
// ************************************************************************
// 
//                             Sundance
//                 Copyright 2011 Sandia Corporation
// 
// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
// the U.S. Government retains certain rights in this software.
//
// Redistribution and use in source and binary forms, with or without
// modification, are permitted provided that the following conditions are
// met:
//
// 1. Redistributions of source code must retain the above copyright
// notice, this list of conditions and the following disclaimer.
//
// 2. Redistributions in binary form must reproduce the above copyright
// notice, this list of conditions and the following disclaimer in the
// documentation and/or other materials provided with the distribution.
//
// 3. Neither the name of the Corporation nor the names of the
// contributors may be used to endorse or promote products derived from
// this software without specific prior written permission.
//
// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
//
// Questions? Contact Kevin Long (kevin.long@ttu.edu)
// 

/* @HEADER@ */

#include "SundanceSpatiallyHomogeneousDOFMapBase.hpp"
#include "PlayaTabs.hpp"
#include "SundanceOut.hpp"
#include "SundanceMaximalCellFilter.hpp"

using namespace Sundance;
using namespace Teuchos;


SpatiallyHomogeneousDOFMapBase
::SpatiallyHomogeneousDOFMapBase(const Mesh& mesh,
  int nTotalFuncs,
  int setupVerb)
  : DOFMapBase(mesh, setupVerb), allowedFuncs_(), funcDomains_()
{
  RCP<Set<int> > f = rcp(new Set<int>());
  for (int i=0; i<nTotalFuncs; i++) f->put(i);
  allowedFuncs_ = f;
  CellFilter cf = new MaximalCellFilter();
  funcDomains_ = Array<CellFilter>(nTotalFuncs, cf);
}


void SpatiallyHomogeneousDOFMapBase::print(std::ostream& os) const
{
  int myRank = mesh().comm().getRank();

  Tabs tabs;
  int dim = mesh().spatialDim();
  RCP<const MapStructure> s = mapStruct();

  for (int p=0; p<mesh().comm().getNProc(); p++)
  {
    mesh().comm().synchronize();
    mesh().comm().synchronize();
    if (p == myRank)
    {
      os << tabs << 
        "========= DOFMap on proc p=" << p << " =============" << std::endl;
      for (int d=dim; d>=0; d--)
      {
        Tabs tabs1;
        os << tabs1 << "dimension = " << d << std::endl;
        for (int c=0; c<mesh().numCells(d); c++)
        {
          Tabs tabs2;
          os << tabs2 << "Cell d=" << d << " LID=" << c << " GID=" 
             << mesh().mapLIDToGID(d, c);
          if (d==0) 
          {
            os << " x=" << mesh().nodePosition(c) << std::endl;
          }
          else 
          {
            Array<int> facetLIDs;
            Array<int> facetDirs;
            mesh().getFacetArray(d, c, 0, facetLIDs, facetDirs);
            Array<int> facetGIDs(facetLIDs.size());
            for (int v=0; v<facetLIDs.size(); v++)
            {
              facetGIDs[v] = mesh().mapLIDToGID(0, facetLIDs[v]);
            }
            os << " nodes LIDs=" << facetLIDs << " GIDs=" << facetGIDs
               << std::endl;
          }
          for (int b=0; b<s->numBasisChunks(); b++)
          {
            for (int f=0; f<s->funcs(b).size(); f++)
            {
              Tabs tabs3;
              Array<int> dofs;
              getDOFsForCell(d, c, s->funcs(b)[f], dofs);
              os << tabs3 << "f=" << s->funcs(b)[f] << " " 
                 << dofs << std::endl;
              if (false)
              {
                os << tabs3 << "{";
                for (int i=0; i<dofs.size(); i++)
                {
                  if (i != 0) os << ", ";
                  if (isLocalDOF(dofs[i])) os << "L";
                  else os << "R";
                }
                os << "}" << std::endl;
              }
            }
          }
        }
      }
    }
    mesh().comm().synchronize();
    mesh().comm().synchronize();
  }
}