level_set_intersected_elem.h

Go to the documentation of this file.

/*
 * MAST: Multidisciplinary-design Adaptation and Sensitivity Toolkit
 * Copyright (C) 2013-2020  Manav Bhatia and MAST authors
 *
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * This library 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
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with this library; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 */


#ifndef __mast_level_set_intersected_elem_h__
#define __mast_level_set_intersected_elem_h__

// MAST includes
#include "mesh/geom_elem.h"


namespace MAST {
    
    // Forward declerations
    class LevelSetIntersection;
    
    class LevelSetIntersectedElem: public MAST::GeomElem  {
        
    public:
        LevelSetIntersectedElem();
        
        virtual ~LevelSetIntersectedElem();
        
        virtual const libMesh::Elem& get_quadrature_elem() const;
        
        virtual const libMesh::Elem& get_quadrature_local_elem() const;

        bool if_elem_has_level_set_boundary() const;

        
        bool if_subelem_has_side_on_level_set_boundary() const;

        
        int get_subelem_side_on_level_set_boundary() const;

        virtual void init(const libMesh::Elem& elem,
                          const MAST::SystemInitialization& sys_init) {
            
            libmesh_error(); // should not get called
        }
        
        virtual void init(const libMesh::Elem& elem,
                          const MAST::SystemInitialization& sys_init,
                          MAST::LevelSetIntersection& intersection);

        
        virtual std::unique_ptr<MAST::FEBase>
        init_fe(bool init_grads,
                bool init_second_order_derivative,
                int extra_quadrature_order = 0) const;
        
        
        virtual std::unique_ptr<MAST::FEBase>
        init_side_fe(unsigned int s,
                     bool init_grads,
                     bool init_second_order_derivative,
                     int extra_quadrature_order = 0) const;

        
    protected:
        
        const libMesh::Elem                   *_sub_elem;
        libMesh::Elem                         *_local_sub_elem;

        MAST::LevelSetIntersection            *_intersection;

        std::vector<libMesh::Node*>           _local_subelem_nodes;
    };
}

#endif  // __mast_level_set_intersected_elem_h__