mapping_p1.hh

Go to the documentation of this file.

/*!
 *
 * Copyright (C) 2015 Technical University of Liberec.  All rights reserved.
 * 
 * This program is free software; you can redistribute it and/or modify it under
 * the terms of the GNU General Public License version 3 as published by the
 * Free Software Foundation. (http://www.gnu.org/licenses/gpl-3.0.en.html)
 * 
 * This program 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.
 *
 * 
 * @file    mapping_p1.hh
 * @brief   Class MappingP1 implements the affine transformation of
 *          the unit cell onto the actual cell.
 * @author  Jan Stebel
 */

#ifndef MAPPING_P1_HH_
#define MAPPING_P1_HH_

#include <armadillo>
#include "fem/mapping.hh"                      // for MappingInternalData (p...
#include "fem/update_flags.hh"                 // for operator&, operator|
#include "mesh/accessors.hh"                     // for ElementAccessor

class Quadrature;


/**
 * Auxiliary templates to resolve nonzero matrix sizes for small dimensions. 
 * ( some compilers do not accept ternary operator in template parameters.)
 */ 
template<unsigned int dim>
class MatrixSizes {
public:  
  static const unsigned int dim_minus_one = dim-1;
};  

template<>
class MatrixSizes<0> {
public:  
  static const unsigned int dim_minus_one = 0;
};  


/**
 * @brief Affine mapping between reference and actual cell.
 *
 * Class MappingP1 implements the affine transformation of
 * the reference cell onto the actual cell.
 *
 * @param dim Dimension of the cells.
 * @param spacedim Dimension of the Euclidean space.
 */
template<unsigned int dim, unsigned int spacedim = 3>
class MappingP1
{
public:

    typedef arma::vec::fixed<dim+1> BaryPoint;
    typedef arma::vec::fixed<spacedim> RealPoint;
    typedef arma::mat::fixed<spacedim, dim+1> ElementMap;
    
    /**
     * @brief Determines which additional quantities have to be computed.
     *
     * @param flags Update flags for required quantities.
     * @return All necessary flags.
     */
    static UpdateFlags update_each(UpdateFlags flags);
    
    /**
     * Map from reference element (barycentric coords) to global coord system.
     * Matrix(3, dim+1) M: x_real = M * x_bary;
     * M columns are real coordinates of nodes.
     */
    static ElementMap element_map(ElementAccessor<3> elm);

    /**
     * Compute jacobian matrix for an element given by the @p coords element map.
     */
    static arma::mat::fixed<spacedim,dim> jacobian(const ElementMap &coords);

    /**
     * Project given point in real coordinates to reference element (barycentic coordinates).
     * Result vector have dimension dim()+1.
     * Use RefElement<dim>::bary_to_local() to get local coordinates.
     */
    static BaryPoint project_real_to_unit(const RealPoint &point, const ElementMap &map);
    
    /**
     * Project given point from reference element (barycentic coordinates) to real coordinates.
     * Use RefElement<dim>::local_to_bary() to get barycentric coordinates in input.
     */
    static RealPoint project_unit_to_real(const BaryPoint &point, const ElementMap &map);

    /**
     * Clip a point given by barycentric cocordinates to the element.
     * If the point is out of the element the closest point
     * projection to the element surface is used.
     */
    static BaryPoint clip_to_element(BaryPoint &barycentric);

    /// Test if element contains given point.
    static bool contains_point(arma::vec point, ElementAccessor<3> elm);


};








#endif /* MAPPING_P1_HH_ */