AlbPUMI_FMDBMeshStruct.hpp

Go to the documentation of this file.

//*****************************************************************//
//    Albany 2.0:  Copyright 2012 Sandia Corporation               //
//    This Software is released under the BSD license detailed     //
//    in the file "license.txt" in the top-level Albany directory  //
//*****************************************************************//

#ifndef ALBPUMI_FMDBMESHSTRUCT_HPP
#define ALBPUMI_FMDBMESHSTRUCT_HPP

#include "Albany_AbstractMeshStruct.hpp"
#include "AlbPUMI_QPData.hpp"
#include "AlbPUMI_NodeData.hpp"

#include "Teuchos_RCP.hpp"
#include "Teuchos_ParameterList.hpp"
#include "Albany_StateInfoStruct.hpp"
#include "EpetraExt_MultiComm.h"
#include <PHAL_Dimension.hpp>

#include "pumi_mesh.h"
#include "pumi_geom.h"

#ifdef SCOREC_ACIS
#include "pumi_geom_acis.h"
#endif
#ifdef SCOREC_PARASOLID
#include "pumi_geom_parasolid.h"
#endif

#include <apf.h>
#include <apfMesh2.h>
#include <apfPUMI.h>

#define NG_EX_ENTITY_TYPE_MAX 15
#define ENT_DIMS 4

namespace AlbPUMI {

  class FMDBMeshStruct : public Albany::AbstractMeshStruct {

  public:

    FMDBMeshStruct(
                  const Teuchos::RCP<Teuchos::ParameterList>& params,
                  const Teuchos::RCP<const Epetra_Comm>& epetra_comm);

    ~FMDBMeshStruct();

    void setFieldAndBulkData(
                  const Teuchos::RCP<const Epetra_Comm>& comm,
                  const Teuchos::RCP<Teuchos::ParameterList>& params,
                  const unsigned int neq_,
                  const Albany::AbstractFieldContainer::FieldContainerRequirements& req,
                  const Teuchos::RCP<Albany::StateInfoStruct>& sis,
                  const unsigned int worksetSize);

    void splitFields(Teuchos::Array<std::string> fieldLayout);

    Teuchos::ArrayRCP<Teuchos::RCP<Albany::MeshSpecsStruct> >& getMeshSpecs();

    std::vector<Teuchos::RCP<QPData<double, 1> > > scalarValue_states;
    std::vector<Teuchos::RCP<QPData<double, 2> > > qpscalar_states;
    std::vector<Teuchos::RCP<QPData<double, 3> > > qpvector_states;
    std::vector<Teuchos::RCP<QPData<double, 4> > > qptensor_states;

    std::vector<std::string> nsNames;
    std::vector<std::string> ssNames;

    msType meshSpecsType();
    pMeshMdl getMesh() { return mesh; }
    pumi::pGModel getMdl() { return model; }

    // Solution history
    int solutionFieldHistoryDepth;
    void loadSolutionFieldHistory(int step);

    bool useCompositeTet(){ return compositeTet; }

    const Albany::DynamicDataArray<Albany::CellSpecs>::type& getMeshDynamicData() const
        { return meshDynamicData; }

    void setupMeshBlkInfo();

    bool hasRestartSolution;
    double restartDataTime;
    int neq;
    int numDim;
    int cubatureDegree;
    bool interleavedOrdering;
    apf::Mesh2* apfMesh;
    bool solutionInitialized;
    bool residualInitialized;

    Teuchos::Array<std::string> solVectorLayout;
    Teuchos::Array<std::string> resVectorLayout;

    double time;

    // Info to map element block to physics set
    bool allElementBlocksHaveSamePhysics;
    std::map<std::string, int> ebNameToIndex;

    int worksetSize;

    std::string outputFileName;
    int outputInterval;
    int useDistributedMesh;

private:

    Teuchos::RCP<const Teuchos::ParameterList>
      getValidDiscretizationParameters() const;

    const CellTopologyData *getCellTopologyData(const FMDB_EntTopo topo);

    int computeWorksetSize(const int worksetSizeMax, const int ebSizeMax) const;

    Teuchos::RCP<Teuchos::FancyOStream> out;

    Teuchos::ArrayRCP<Teuchos::RCP<Albany::MeshSpecsStruct> > meshSpecs;

    // Information that changes when the mesh adapts
    Albany::DynamicDataArray<Albany::CellSpecs>::type meshDynamicData;

    pumi::pGModel model;
    pMeshMdl mesh;

    bool compositeTet;

#ifndef SCOREC_ACIS
    int PUMI_Geom_RegisterAcis() {
      fprintf(stderr,"ERROR: FMDB Discretization -> Cannot find Acis\n");
      exit(1);
    }
#endif
#ifndef SCOREC_PARASOLID
    int PUMI_Geom_RegisterParasolid() {
      fprintf(stderr,"ERROR: FMDB Discretization -> Cannot find Parasolid\n");
      exit(1);
    }
#endif

  };

}
#endif