InteractionToStressConverter.h

Go to the documentation of this file.

//                                                         //
// Copyright (c) 2003-2014 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.opensource.org/licenses/osl-3.0.php          //
//                                                         //


#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