EllipsoidVtkOutput.cpp 3.77 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35
//======================================================================================================================
//
//  This file is part of waLBerla. waLBerla is free software: you can
//  redistribute it and/or modify it under the terms of the GNU General Public
//  License as published by the Free Software Foundation, either version 3 of
//  the License, or (at your option) any later version.
//
//  waLBerla is distributed in the hope that it will be useful, but WITHOUT
//  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
//  FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
//  for more details.
//
//  You should have received a copy of the GNU General Public License along
//  with waLBerla (see COPYING.txt). If not, see <http://www.gnu.org/licenses/>.
//
//! \file EllipsoidVtkOutput.cpp
//! \author Christoph Rettinger <christoph.rettinger@fau.de>
//
//======================================================================================================================

#include "core/DataTypes.h"
#include "EllipsoidVtkOutput.h"
#include "pe/rigidbody/BodyStorage.h"

#include "core/selectable/IsSetSelected.h"
#include "core/uid/GlobalState.h"


namespace walberla {
namespace pe {

std::vector< EllipsoidVtkOutput::Attributes > EllipsoidVtkOutput::getAttributes() const
{
   std::vector< Attributes > attributes;
   attributes.push_back( Attributes( vtk::typeToString< float >(), "mass", uint_c(1) ) );
36
   attributes.push_back( Attributes( vtk::typeToString< float >(), "tensorGlyph", uint_c(6) ) );
37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71
   attributes.push_back( Attributes( vtk::typeToString< float >(), "velocity", uint_c(3) ) );
   attributes.push_back( Attributes( vtk::typeToString< int >(),   "rank", uint_c(1) ) );
   attributes.push_back( Attributes( vtk::typeToString< id_t >(),  "id", uint_c(1) ) );
   attributes.push_back( Attributes( vtk::typeToString< id_t >(),  "uid", uint_c(1) ) );

   return attributes;
}

void EllipsoidVtkOutput::configure()
{
   bodies_.clear();
   tensorGlyphs_.clear();

   for( auto blockIt = blockStorage_.begin(); blockIt != blockStorage_.end(); ++blockIt )
   {

      const Storage& bs = *(blockIt->getData<const Storage>( storageID_ ));

      for( auto it = bs[0].begin(); it != bs[0].end(); ++it )
      {
         if (it->getTypeID() == Ellipsoid::getStaticTypeID())
         {
            auto ellipsoid = static_cast<ConstEllipsoidID> (*it);
            bodies_.push_back(ellipsoid);

            // compute tensor glyph for visualization with ParaView (tensorGlyph)
            Mat3 rotMat = ellipsoid->getRotation();
            Vector3<real_t> directionVectorX(rotMat[0], rotMat[3], rotMat[6]);
            Vector3<real_t> directionVectorY(rotMat[1], rotMat[4], rotMat[7]);
            Vector3<real_t> directionVectorZ(rotMat[2], rotMat[5], rotMat[8]);
            Vector3<real_t> semiAxes = ellipsoid->getSemiAxes();
            Mat3 axa = math::dyadicProduct(directionVectorX, directionVectorX);
            Mat3 bxb = math::dyadicProduct(directionVectorY, directionVectorY);
            Mat3 cxc = math::dyadicProduct(directionVectorZ, directionVectorZ);
            Mat3 tensor = axa * semiAxes[0] + bxb * semiAxes[1] + cxc * semiAxes[2];
72
            // use symmetry to only write 6 of the 9 elements: XX YY ZZ XY YZ XZ
73 74
            tensorGlyphs_.push_back({{tensor(0,0), tensor(1,1), tensor(2,2),
                                      tensor(0,1), tensor(1,2), tensor(0,2)}});
75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93
         }
      }
   }
}

std::vector< math::Vector3< real_t > > EllipsoidVtkOutput::getPoints()
{
   std::vector< math::Vector3< real_t > > result( bodies_.size() );
   auto resultIt = result.begin();
   for( auto it = bodies_.begin(); it != bodies_.end(); ++it, ++resultIt )
   {
      *resultIt = (*it)->getPosition();
   }
   return result;
}

} // namespace pe
} // namespace walberla