embedded.h

/*************************************************************************************
 * MechSys - A C++ library to simulate (Continuum) Mechanical Systems                *
 * Copyright (C) 2005 Dorival de Moraes Pedroso <dorival.pedroso at gmail.com>       *
 * Copyright (C) 2005 Raul Dario Durand Farfan  <raul.durand at gmail.com>           *
 *                                                                                   *
 * This file is part of MechSys.                                                     *
 *                                                                                   *
 * MechSys 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.                                                                          *
 *                                                                                   *
 * MechSys 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      *
 * MechSys; if not, write to the Free Software Foundation, Inc., 51 Franklin Street, *
 * Fifth Floor, Boston, MA 02110-1301, USA                                           *
 *************************************************************************************/

#ifndef MECHSYS_FEM_EMBEDDED_H
#define MECHSYS_FEM_EMBEDDED_H

#ifndef REAL
#define REAL double
#endif

#include "fem/ele/element.h"

namespace FEM
{

class Embedded: public virtual Element
{
public:
    void Jacobian(LinAlg::Matrix<REAL> const & derivs, LinAlg::Matrix<REAL> & J) const;
    void Initialize(Array<Node*> Emb_connects, Element * Tresspassed)
    { 
        _n_emb_nodes  = _connects.size();
        assert(Emb_connects.size()==_n_emb_nodes);
        for (size_t i=0; i<_n_emb_nodes; i++) SetNode(i, Emb_connects[i]);
        _tresspassed  = Tresspassed;
        _n_nodes      = Tresspassed->nNodes();
        _connects     = Tresspassed->Connects();
        _emb_connects = Emb_connects;
    };
    int  nNodes()    const { return _tresspassed->nNodes(); }
    int  nEmbNodes() const { return _n_emb_nodes; }
protected:
    Element      * _tresspassed;
    Array<Node*>   _emb_connects;
    size_t         _n_emb_nodes;
}; // class Emb3Equilib




inline void Embedded::Jacobian(LinAlg::Matrix<REAL> const & derivs, LinAlg::Matrix<REAL> & J) const // {{{
{
    // Calculate a matrix with nodal coordinates
    LinAlg::Matrix<REAL> cmatrix; // size = _n_nodes x 3
    cmatrix.Resize(_n_emb_nodes,3);   // 3 => X,Y,Z (all nodes have X,Y and Z)
    
    // Loop along all nodes of this element
    for (size_t i=0; i<_n_emb_nodes; i++)
    {
        cmatrix(i,0) = _emb_connects[i]->X();
        cmatrix(i,1) = _emb_connects[i]->Y();
        cmatrix(i,2) = _emb_connects[i]->Z();
    }

    // Calculate the Jacobian
    J = derivs*cmatrix;
} // }}}

}; // namespace FEM

#endif // MECHSYS_FEM_EMBEDDED_H

// vim:fdm=marker

---