matrix_shop_cell.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//    $Id: matrix_shop_cell.h 1348 2013-08-16 00:45:07Z secanell $
//
//    Copyright (C) 2006, 2007, 2008, 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_cell_h
#define __appshop__matrix_shop_cell_h

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

using namespace AppFrame;

namespace AppShop
{
  namespace Matrix
  {
    namespace Cell
    {
      template <int dim>
      void mass (FullMatrix<double>& M,
                 const FEValuesBase<dim>& fe,
                 const double a = 1.);
      template <int dim>
      void mass (FullMatrix<double>& M,
                 const FEValuesBase<dim>& fe,
                 const FEValuesBase<dim>& fetest,
                 const unsigned int component,
                 const double a = 1.);
      template <int dim>
      void mass_scaled (FullMatrix<double>& M,
                        const FEValuesBase<dim>& fe,
                        const std::vector<double>& a,
                        bool divide);

      template <int dim>
      void mass_scaled (FullMatrix<double>& M,
                        const FEValuesBase<dim>& fe,
                        const FEValuesBase<dim>& fetest,
                        const unsigned int select,
                        const std::vector<double>& factor,
                        bool divide = false);
      
      template <int dim>
      void advection(
        FullMatrix<double>& M,
        const FEValuesBase<dim>& fe,
        const std::vector<Point<dim> >& velocity);

      template <int dim>
      void advection(
        FullMatrix<double>& M,
        const FEValuesBase<dim>& fe,
        const VectorSlice<const std::vector<std::vector<double> > >& velocity,
        double factor = 1.);

      template <int dim>
      void div_primal(
        std::vector<MatrixBlock<FullMatrix<double> > >& M,
        const std::vector<unsigned int>& indices,
        const FEValuesBase<dim>& fe,
        const FEValuesBase<dim>& fetest,
        double factor = 1.);

      template <int dim>
      void div(
        std::vector<MatrixBlock<FullMatrix<double> > >& M,
        const std::vector<unsigned int>& indices,
        const FEValuesBase<dim>& fe,
        const FEValuesBase<dim>& fetest,
        double factor = 1.);

      template <int dim>
      void div_scaled (
        FullMatrix<double>& M,
        const FEValuesBase<dim>& fe,
        const FEValuesBase<dim>& fetest, 
        const std::vector<double>& factor);

      template <int dim>
      void di(
        FullMatrix<double>& M,
        const FEValuesBase<dim>& fe,
        const FEValuesBase<dim>& fetest,
        unsigned int i);

      template <int dim>
      void di_scaled (
        FullMatrix<double>& M,
        const FEValuesBase<dim>& fe,
        const FEValuesBase<dim>& fetest,
        unsigned int d,
        std::vector<double>& factor);

      template <int dim>
      void eps_div(
        FullMatrix<double>& M,
        const FEValuesBase<dim>& fe,
        const FEValuesBase<dim>& fetest);

      template <int dim>
      void laplacian(
        FullMatrix<double>& M,
        const FEValuesBase<dim>& fe,
        double nu = 1.);
      template <int dim>
      void laplacian(
        FullMatrix<double>& M,
        const FEValuesBase<dim>& fe,
        Tensor<2,dim> nu);

      template <int dim>
      void laplacian_scaled(
        FullMatrix<double>& M,
        const FEValuesBase<dim>& fe,
        const std::vector<double>& nu,
        bool divide = false);

      template <int dim>
      void laplacian_scaled(
        FullMatrix<double>& M,
        const FEValuesBase<dim>& fe,
        const double nu,
        bool divide = false);      
      
      template <int dim>
        void div_div (FullMatrix<double>& M,
                      const FEValuesBase<dim>& fe,
                      const double factor = 1.);

    }
  }
}

#endif