ParaViewModule.h

Go to the documentation of this file.

 
#pragma once
 
#include <jem/base/Array.h>
#include <jem/base/CString.h>
#include <jem/base/IllegalArgumentException.h>
#include <jem/base/Slice.h>
#include <jem/base/System.h>
#include <jem/io/FileName.h>
#include <jem/io/FileOpenException.h>
#include <jem/io/FileWriter.h>
#include <jem/io/PrintWriter.h>
#include <jem/util/ArrayBuffer.h>
#include <jem/util/Properties.h>
 
#include <jive/app/Module.h>
#include <jive/app/ModuleFactory.h>
#include <jive/app/Names.h>
#include <jive/fem/ElementGroup.h>
#include <jive/fem/ElementSet.h>
#include <jive/fem/NodeSet.h>
#include <jive/model/Actions.h>
#include <jive/model/Model.h>
#include <jive/model/StateVector.h>
#include <jive/util/Assignable.h>
#include <jive/util/DofSpace.h>
#include <jive/util/FuncUtils.h>
#include <jive/util/Globdat.h>
#include <jive/util/ItemGroup.h>
#include <jive/util/ItemSet.h>
#include <jive/util/SparseTable.h>
#include <jive/util/XTable.h>
 
#include <filesystem>
 
using jem::ALL;
using jem::Array;
using jem::idx_t;
using jem::IllegalArgumentException;
using jem::Limits;
using jem::newInstance;
using jem::Ref;
using jem::SliceFrom;
using jem::SliceFromTo;
using jem::SliceTo;
using jem::String;
using jem::sum;
using jem::where;
using jem::io::endl;
using jem::io::FileFlags;
using jem::io::FileName;
using jem::io::FileOpenException;
using jem::io::FileWriter;
using jem::io::PrintWriter;
using jem::io::Writer;
using jem::util::ArrayBuffer;
using jem::util::Properties;
 
using jive::IdxVector;
using jive::Matrix;
using jive::StringVector;
using jive::Vector;
using jive::app::Module;
using jive::app::PropNames;
using jive::fem::ElementGroup;
using jive::fem::ElementSet;
using jive::fem::NodeSet;
using jive::model::Actions;
using jive::model::Model;
using jive::model::StateVector;
using jive::util::Assignable;
using jive::util::DofSpace;
using jive::util::Function;
using jive::util::Globdat;
using jive::util::ItemGroup;
using jive::util::ItemSet;
using jive::util::SparseTable;
using jive::util::XTable;
 
class ParaViewModule : public Module
{
protected:
  struct ElementInfo
  {
    String name;           
    String shape;          
    StringVector elemData; 
    StringVector nodeData; 
    StringVector dispData; 
    StringVector dofData;  
  };
 
public:
  static const char *TYPE_NAME; 
  static const char *SPACING;   
 
  explicit ParaViewModule(const String &name = "paraView");
 
  virtual Status init(const Properties &conf, const Properties &props,
                      const Properties &globdat) override;
 
  virtual Status run(const Properties &globdat) override;
 
  virtual void shutdown(const Properties &globdat) override;
 
  static Ref<Module> makeNew(const String &name, const Properties &conf,
                             const Properties &props, const Properties &globdat);
 
  static void declare();
 
  static idx_t nameToVTKNum(const String &name);
 
  static IdxVector gmsh2ParaNodeOrder(const IdxVector elNodes, const String &name);
 
private:
 
  void writeFile_(const String &fileName, const Properties &globdat);
 
  void writePiece_(const Ref<PrintWriter> &file, const Assignable<NodeSet> &points,
                   const Assignable<ElementSet> &cells,
                   const Assignable<ElementGroup> &group, const Vector &disp,
                   const Vector &velo, const Vector &acce, const Ref<DofSpace> &dofs,
                   const Ref<Model> &model, const Properties &globdat,
                   const ElementInfo &info);
 
  void writeDataArray_(const Ref<PrintWriter> &file, const Matrix &data,
                       const String &type, const String &name);
 
  void writeDataArray_(const Ref<PrintWriter> &file, const Ref<XTable> &data,
                       const IdxVector &rows, const String &type, const String &name);
 
  void writeDataArray_(const Ref<PrintWriter> &file, const Vector &data,
                       const String &type, const String &name);
 
private:
  String nameFormat_;          
  String fileType_;            
  StringVector elemSets_;      
  Array<ElementInfo> setInfo_; 
  Ref<Function> sampleCond_;   
  Ref<Function> sampleInfo_;   
 
  idx_t outNum_;                      
  bool pvdPrint_;                     
  String pvdName_;                    
  ArrayBuffer<double> pvdTimeBuffer_; 
  ArrayBuffer<String> pvdNameBuffer_; 
};