src/fem/nedelecelement1.cpp

Go to the documentation of this file.

/***************************************************************************
                          nedelecelement1.cpp  -  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.                                   *
 *                                                                         *
 ***************************************************************************/

#include "nedelecelement1.h"

using namespace colarray;
using namespace mesh;

NedelecElement1::NedelecElement1()
{ }

NedelecElement1::~NedelecElement1()
{ }

const int* NedelecElement1::get_dof_ids(mesh::TetMesh* m, mesh::Tet* tet) const {
    static int dof_ids[6];
    for (int i = 0; i < 6; i ++) {
        dof_ids[i] = tet->get_edge_id(i);
    }
    return dof_ids;
}

const bool* NedelecElement1::get_eliminated_dofs(mesh::TetMesh* m,
                                                 mesh::Tet* tet,
                                                 unsigned int sym_plane_config) const {
    static bool eliminated_dofs[6];
    for (int i = 0; i < 6; i ++) {
        if (tet->get_edge(i)->on_boundary()) {
            eliminated_dofs[i] = true;      
        } else if ((sym_plane_config ^ 0x7) & tet->get_edge(i)->on_symmetry_planes()) {
            eliminated_dofs[i] = true;
        } else {
            eliminated_dofs[i] = false;
        }
    }
    return eliminated_dofs;
}
      
void NedelecElement1::get_element_matrix(mesh::Tet* tet, Matrix<double>& Ae, Matrix<double>& Me) {

    if (Ae.get_m() != 6 || Ae.get_n() != 6) {
        // inappropriate Matrix shape: reassign
        Ae = Matrix<double>(6, 6);
    }
    if (Me.get_m() != 6 || Me.get_n() != 6) {
        // inappropriate Matrix shape: reassign
        Me = Matrix<double>(6, 6);
    }

#ifndef GAUSSIANQUADRATURE  
    Tet::GradientMatrix J(*tet);
    get_Ae(tet, 1, J, Ae);
    get_Me(tet, 1, J, Me);
#else
    get_AeMe(tet, 1, Ae, Me);
#endif
}

mesh::Vector3 NedelecElement1::eval(Tet* tet, const Vector3& x, const Vector<double>& dof) {
    mesh::Tet::GradientMatrix J(*tet);
    return eval_cartesian(tet, 1, dof, J, x);  
}
mesh::Vector3 NedelecElement1::eval_curl(Tet* tet, const Vector3& x, const Vector<double>& dof) {
    mesh::Tet::GradientMatrix J(*tet);
    return eval_curl_cartesian(tet, 1, dof, J, x);  
}

---