mesh_worker_workers.h

Go to the documentation of this file.

// ------------------------------------------------------------------
//
// $Id: mesh_worker_workers.h 1348 2013-08-16 00:45:07Z secanell $
//
// Copyright (C) 2006, 2007, 2008, 2009 by Guido Kanschat.
// Copyright (C) 2012 by Valentin N. Zingan, University of Alberta.
//
// This file is subject to QPL and may not be  distributed
// without copyright and license information. Please refer
// to the file deal.II/doc/license.html for the  text  and
// further information on this license.
//
// ------------------------------------------------------------------

#ifndef _DEALII_APPFRAME_MESH_WORKER_WORKERS_H_
#define _DEALII_APPFRAME_MESH_WORKER_WORKERS_H_

#include <appframe/mesh_worker_info.h>

// ################################################################################################

// ##### CATEGORY-2 ################################################# WORKER_OBJECTS ##############
//
// ##### MeshWorker::WorkerObjects::LocalWorker       #####
// ##### MeshWorker::WorkerObjects::IntegrationWorker #####

namespace MeshWorker // # 2
{

namespace WorkerObjects
{

// --- 1 ---

template<int dim>
class LocalWorker
{
public:

  LocalWorker()
  {
    boundary_fluxes = true;
    interior_fluxes = true;
  }

  void cell(const InfoObjects::DoFInfo<dim>& cell);

  void bdry(const InfoObjects::DoFInfo<dim>& face);

  void face(const InfoObjects::DoFInfo<dim>& face1,
            const InfoObjects::DoFInfo<dim>& face2);

  bool boundary_fluxes;

  bool interior_fluxes;
};

// --- 2 ---

template<int dim>
class IntegrationWorker : public LocalWorker<dim>
{
public:

  IntegrationWorker();

  void initialize_selectors(const VectorSelector& cs,
                            const VectorSelector& bs,
                            const VectorSelector& fs);

  void add_selector(const std::string& name,
                    bool               values,
                    bool               gradients,
                    bool               hessians,
                    bool               cell,
                    bool               bdry,
                    bool               face);

  void add_update_flags(const UpdateFlags flags,
                        bool              cell = true,
                        bool              bdry = true,
                        bool              face = true,
                        bool              ngbr = true)
  {
    if(cell) cell_flags |= flags;
    if(bdry) bdry_flags |= flags;
    if(face) face_flags |= flags;
    if(ngbr) ngbr_flags |= flags;
  }

  void initialize_gauss_quadrature(unsigned int n_cell_points,
                                   unsigned int n_bdry_points,
                                   unsigned int n_face_points)
  {
    cell_quadrature = QGauss<dim>(n_cell_points);
    bdry_quadrature = QGauss<dim-1>(n_bdry_points);
    face_quadrature = QGauss<dim-1>(n_face_points);
  }

  MeshWorker::VectorSelector cell_selector;

  MeshWorker::VectorSelector bdry_selector;

  MeshWorker::VectorSelector face_selector;

  Quadrature<dim> cell_quadrature;

  Quadrature<dim-1> bdry_quadrature;

  Quadrature<dim-1> face_quadrature;

  UpdateFlags cell_flags;

  UpdateFlags bdry_flags;

  UpdateFlags face_flags;

  UpdateFlags ngbr_flags;
};

// --- implementation ---

template<int dim>
inline
IntegrationWorker<dim>::IntegrationWorker()
:
LocalWorker<dim>()
{
  cell_flags = update_JxW_values;
  bdry_flags = UpdateFlags(update_normal_vectors | update_JxW_values);
  face_flags = bdry_flags;
  ngbr_flags = update_default;
}

template<int dim>
inline
void
IntegrationWorker<dim>::initialize_selectors(const VectorSelector& cs,
                                             const VectorSelector& bs,
                                             const VectorSelector& fs)
{
  cell_selector = cs;
  bdry_selector = bs;
  face_selector = fs;

  if(cell_selector.has_values()    != 0) cell_flags |= update_values;
  if(cell_selector.has_gradients() != 0) cell_flags |= update_gradients;
  if(cell_selector.has_hessians()  != 0) cell_flags |= update_hessians;

  if(bdry_selector.has_values()    != 0) bdry_flags |= update_values;
  if(bdry_selector.has_gradients() != 0) bdry_flags |= update_gradients;
  if(bdry_selector.has_hessians()  != 0) bdry_flags |= update_hessians;

  if(face_selector.has_values()    != 0) face_flags |= update_values;
  if(face_selector.has_gradients() != 0) face_flags |= update_gradients;
  if(face_selector.has_hessians()  != 0) face_flags |= update_hessians;

  if(face_selector.has_values()    != 0) ngbr_flags |= update_values;
  if(face_selector.has_gradients() != 0) ngbr_flags |= update_gradients;
  if(face_selector.has_hessians()  != 0) ngbr_flags |= update_hessians;
}

template<int dim>
inline
void
IntegrationWorker<dim>::add_selector(const std::string& name,
                                     bool               values,
                                     bool               gradients,
                                     bool               hessians,
                                     bool               cell,
                                     bool               bdry,
                                     bool               face)
{
  if(cell) cell_selector.add_vector(name, values, gradients, hessians);
  if(bdry) bdry_selector.add_vector(name, values, gradients, hessians);
  if(face) face_selector.add_vector(name, values, gradients, hessians);

  if(cell_selector.has_values()    != 0) cell_flags |= update_values;
  if(cell_selector.has_gradients() != 0) cell_flags |= update_gradients;
  if(cell_selector.has_hessians()  != 0) cell_flags |= update_hessians;

  if(bdry_selector.has_values()    != 0) bdry_flags |= update_values;
  if(bdry_selector.has_gradients() != 0) bdry_flags |= update_gradients;
  if(bdry_selector.has_hessians()  != 0) bdry_flags |= update_hessians;

  if(face_selector.has_values()    != 0) face_flags |= update_values;
  if(face_selector.has_gradients() != 0) face_flags |= update_gradients;
  if(face_selector.has_hessians()  != 0) face_flags |= update_hessians;

  if(face_selector.has_values()    != 0) ngbr_flags |= update_values;
  if(face_selector.has_gradients() != 0) ngbr_flags |= update_gradients;
  if(face_selector.has_hessians()  != 0) ngbr_flags |= update_hessians;
}

} // namespace WorkerObjects

} // namespace MeshWorker # 2

#endif