src/vtkexport.h

Go to the documentation of this file.

//
// C++ Interface: vtkexport
//
// Description:
//
//
// Author: Roman Geus <roman.geus@psi.ch>, (C) 2004
//
// Copyright: See COPYING file that comes with this distribution
//
//
#ifndef POSTPROCESSVTKEXPORT_H
#define POSTPROCESSVTKEXPORT_H

#include <math.h>
#include <vector>
#include "vtkexport.h"
#include "tetmesh/point.h"
#include "tetmesh/tetmesh.h"
#include "tetmesh/tet.h"
#include "tetmesh/vector3.h"
#include "Epetra_Operator.h"
#include "Epetra_MultiVector.h"
#include "Epetra_SerialDenseVector.h"
#include "Epetra_Vector.h"
#include "nedelecmesh.h"


namespace postprocess {

class VtkExport {

public:
    VtkExport(const NedelecMesh& nedelec_mesh,
              const Epetra_SerialDenseVector& lambda,
              const Epetra_MultiVector& Q);

    void export_mesh(string file);

    void export_eigenfields(std::string file);
    void export_eigenfields2(std::string file, int nof_cells);

    void export_a_to_b(std::string file, mesh::Vector3 a, mesh::Vector3 b, int nof_intervals);

    void integrate_a_to_b(int mode_id, mesh::Vector3 a, mesh::Vector3 b, double& int_e);

    int get_nof_vectors() const
        { return Q_.NumVectors(); }

    int get_nof_dofs() const
        { return Q_.GlobalLength(); }

    void get_eigenpair(int k, double* lambda, colarray::Vector<double>* q);

private: /* private members */
    
    void getdimensions(int N, mesh::Vector3& box_len, int dim[3]);

    const NedelecMesh& nedelec_mesh_;
    const Epetra_SerialDenseVector lambda_;
    const Epetra_MultiVector Q_;
    const Epetra_Comm& comm_;
    int precision_;
}; //end class vtkExport

}; //namespace postprocess

#endif

---