src/fem/nedelecmesh.h

Go to the documentation of this file.

/***************************************************************************
                          nedelecmesh.h  -  description
                             -------------------
    begin                : Fri Dec 19 2003
    copyright            : (C) 2003 by Roman Geus
    email                : roman.geus@psi.ch
***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program 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 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/

#ifndef NEDELECMESH_H
#define NEDELECMESH_H

#include "Epetra_CrsMatrix.h"
#include "Epetra_FECrsMatrix.h"
#include "tetmesh/paralleltetmesh.h"
#include <vector>

#include "tetmesh/tetmesh.h"
#include "tetmesh/vector3.h"
#include "colarray/vector.h"
#include "nedelecelement.h"

#include "tetmesh/entity.h"

class NedelecMesh {
public:
  NedelecMesh(mesh::TetMesh* mesh, int order, unsigned int sym_plane_config);
    ~NedelecMesh();
    
    NedelecElement* get_element(mesh::id_t i) const
    { return _element; }
    
    void assembleAM(Epetra_FECrsMatrix* A, Epetra_FECrsMatrix* M);
    
    mesh::id_t get_num_my_dofs(int order) const;
    
    // int get_nof_mapped_dofs(int order) const;
    
    mesh::id_t get_num_global_dofs(int order) const;
    
    mesh::id_t get_num_global_mdofs(int order) const;
    
    mesh::id_t gdof(mesh::id_t ldof) const;
    
    bool is_owned_gdof(mesh::id_t gdof) const;
    
    mesh::id_t map_dof(mesh::id_t orig_dof) const;
    
    int map_dof_rev(int mapped_dof) const;
    
    void set_map(unsigned int sym_plane_config);
    
    mesh::Vector3 eval(const mesh::Vector3& x, const colarray::Vector<double>& q) const;
    
    mesh::Vector3 eval_curl(const mesh::Vector3& x, const colarray::Vector<double>& q) const;

        /* quality factor calculation for specific eigenpar where lambda is the
         * eigenvalue and q ist the eigenvector; offers seamless parallel and serial usage */
        double NedelecMesh::quality(double lambda, const double* q) const;


    double q_factor(double lambda, const double* q) const;

    //return pointer to the tetmesh of this nedelec_mesh
    mesh::TetMesh* getTetMesh()const
    {return _mesh;}
    
    void generate_sorted_AM(const Epetra_FECrsMatrix* A, const Epetra_FECrsMatrix* M, Epetra_FECrsMatrix* A_s, Epetra_FECrsMatrix* M_s);
  
private:

    typedef std::pair<mesh::id_t, mesh::id_t> nodepair;   
    typedef utility::triple<nodepair, mesh::id_t, mesh::id_t> dof_with_ids;
    void sort_matrix(const Epetra_FECrsMatrix* M, Epetra_FECrsMatrix* M_s, const vector<dof_with_ids>& all_dofs);
    
    mesh::TetMesh* _mesh;
    mesh::ParallelTetMesh* _pmesh;

    NedelecElement* _element;
    int _order;
    std::map<mesh::id_t, mesh::id_t> _mgid_map;
    std::vector<mesh::id_t> _map_rev;
    std::vector<mesh::id_t> _nof_orig_dofs;
    std::vector<mesh::id_t> _nof_mapped_dofs;
    std::vector<mesh::id_t> _num_global_dofs;
};

#endif

---