/usr/include/deal.II/base/function_bessel.h is in libdeal.ii-dev 8.1.0-6ubuntu1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 | // ---------------------------------------------------------------------
// $Id: function_bessel.h 30036 2013-07-18 16:55:32Z maier $
//
// Copyright (C) 2010 - 2013 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, 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.
// The full text of the license can be found in the file LICENSE at
// the top level of the deal.II distribution.
//
// ---------------------------------------------------------------------
#ifndef __deal2__function_bessel_h
#define __deal2__function_bessel_h
#include <deal.II/base/config.h>
#include <deal.II/base/function.h>
#include <deal.II/base/point.h>
DEAL_II_NAMESPACE_OPEN
namespace Functions
{
/**
* The Bessel functions of first kind or positive integer order.
*
* @author Guido Kanschat
* @date 2010
*/
template <int dim>
class Bessel1 : public Function<dim>
{
public:
Bessel1(const unsigned int order,
const double wave_number,
const Point<dim> center = Point<dim>());
virtual double value (const Point<dim> &points, const unsigned int component) const;
virtual void value_list (const std::vector<Point<dim> > &points,
std::vector<double> &values,
const unsigned int component = 0) const;
virtual Tensor<1,dim> gradient (const Point<dim> &p,
const unsigned int component = 0) const;
virtual void gradient_list (const std::vector<Point<dim> > &points,
std::vector<Tensor<1,dim> > &gradients,
const unsigned int component = 0) const;
private:
unsigned int order;
double wave_number;
Point<dim> center;
};
}
DEAL_II_NAMESPACE_CLOSE
#endif
|