matrix_shop_bdry.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//    $Id: matrix_shop_bdry.h 1348 2013-08-16 00:45:07Z secanell $
//
//    Copyright (C) 2006, 2007, 2009 by Guido Kanschat
//
//    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 __appshop__matrix_shop_bdry_h
#define __appshop__matrix_shop_bdry_h

#include <appframe/app_shop.h>
#include <appframe/matrix_base.h>

using namespace AppFrame;


namespace AppShop
{
  namespace Matrix
  {
    namespace Boundary
    {
      template <int dim>
      void advection(
        FullMatrix<double>& M,
        const FEFaceValuesBase<dim>& fe,
        const std::vector<Point<dim> >& velocity,
        double factor = 1.);
      
      template <int dim>
      void advection(
        FullMatrix<double>& M,
        const FEFaceValuesBase<dim>& fe,
        const VectorSlice<const std::vector<std::vector<double> > >& velocity,
        double factor = 1.);
      
      template <int dim>
      void u_dot_n (
        std::vector<MatrixBlock<FullMatrix<double> > >& M,
        const std::vector<unsigned int>& indices,
        const FEFaceValuesBase<dim>& fe,
        const FEFaceValuesBase<dim>& fetest,
        double factor = 1.);
      
      template <int dim>
      void u_dot_n_scaled (
        std::vector<MatrixBlock<FullMatrix<double> > >& M,
        const std::vector<unsigned int>& indices,
        const FEFaceValuesBase<dim>& fe,
        const FEFaceValuesBase<dim>& fetest,
        const std::vector<double>& factor,
        bool divide = false);
      
      template <int dim>
      void u_ni (
        FullMatrix<double>& M,
        const FEFaceValuesBase<dim>& fe,
        const FEFaceValuesBase<dim>& fetest,
        unsigned int i);
      
      template <int dim>
      void jump (
        FullMatrix<double>& M,
        const FEFaceValuesBase<dim>& fe,
        const double a = 1.);
      template <int dim>
      void jump_scaled (
        FullMatrix<double>& M,
        const FEFaceValuesBase<dim>& fe,
        const std::vector<double>& factor,
        const double factor2 = 1.,
        bool divide = false);
      template <int dim>
      void nitsche (
        FullMatrix<double>& M,
        const FEFaceValuesBase<dim>& fe,
        double penalty,
        double factor = 1.);
      
      template <int dim>
      void nitsche_scaled (
        FullMatrix<double>& M,
        const FEFaceValuesBase<dim>& fe,
        const std::vector<double>& factor,
        double penalty);
      
      template <int dim>
      void nitsche_div (
        FullMatrix<double>& M,
        const FEFaceValuesBase<dim>& fe,
        double penalty);
    }
  }
}


#endif