InteractionToStressConverter.h

Go to the documentation of this file.

//                                                         //
// Copyright (c) 2003-2017 by The University of Queensland //
// Centre for Geoscience Computing                         //
// http://earth.uq.edu.au/centre-geoscience-computing      //
//                                                         //
// Primary Business: Brisbane, Queensland, Australia       //
// Licensed under the Open Software License version 3.0    //
// http://www.apache.org/licenses/LICENSE-2.0              //
//                                                         //
 
 
#ifndef ESYS_LSMINTERACTIONTOSTRESSCONVERTOR_H
#define ESYS_LSMINTERACTIONTOSTRESSCONVERTOR_H
 
#include "Tools/StressCalculator/StressTensorCalculator.h"
#include "Tools/StressCalculator/StressTensorCollection.h"
#include "Tools/StressCalculator/VtkUnstructuredGrid.h"
#include "Tools/StressCalculator/CartesianGrid.h"
#include "Tools/StressCalculator/Vec3Comparer.h"
 
#include "Foundation/BoundingBox.h"
 
#include <iostream>
#include <iomanip>
 
namespace std
{
  template <>
  struct less<Vec3>
  {
    bool operator()(const Vec3 &v1, const Vec3 &v2) const
    {
      return esys::lsm::Vec3XyzComparer()(v1, v2);
    }
  };
}
 
namespace esys
{
  namespace lsm
  {
    typedef CartesianGrid<double>         DoubleGrid;
    typedef CartesianGrid<StressTensor *> StressTensorPtrGrid;
    typedef CartesianGrid<Tensor>         TensorGrid;
 
    class InteractionToStressConverter
    {
    public:
      typedef StressTensorCollection<ContactPtTensorCalculator> StressTensCollection;
      typedef StressTensCollection::StressCalculator          StressTensorCalculator;
 
      InteractionToStressConverter(const BoundingBox &box, double gridSpacing);
 
      void addRaw2Interactions(std::istream &iStream);
 
      void writeVtkUnstructuredXml(std::ostream &oStream);
 
      void writeVtkUnstructuredXmlGridInformation(std::ostream &oStream);
 
      void writeVtkStructuredXml(std::ostream &oStream);
 
      StressTensorPtrGrid &getTensorIrregularGrid();
      void calcTensorIrregularGrid();
 
      TensorGrid &getTensorRegularGrid();
 
      void calcTensorRegularGrid();
 
      DoubleGrid &getDevRegularGrid();
 
      void calcDevRegularGrid();
 
      void writeFlatStructured(std::ostream &oStream);
 
      void writeFlatUnstructured
      (
        std::ostream &oStream
      );
 
      void writeUnstructuredDx(std::ostream &oStream);
 
      double getMaxRadius();
 
    protected:
      double getRealDevStress(const Tensor &stressTensor) const;
      double getNormDevStress(const Tensor &stressTensor) const;
 
      bool is3d() const;
      
    private:
      double                 m_gridSpacing;
      BoundingBox            m_bBox;
      StressTensorCalculator m_stressCalculator;
      StressTensCollection m_stressTensorCollection;
      TensorGrid             m_regTensorGrid;
      DoubleGrid             m_regDevStressGrid;
      StressTensorPtrGrid    m_irrStressTensorGrid;
    };
  }
}
#endif