Those are minor code changes necessary to compile the library with MSVC.
#include <deal.II/base/config.h>
#include <deal.II/base/types.h>
+
+#include <cmath>
#include <complex>
DEAL_II_NAMESPACE_OPEN
typename NumberTraits<number>::real_type
NumberTraits<number>::abs (const number &x)
{
- return std::fabs(x);
+ return std::abs(x);
}
enum EndPoint { left,right };
/// Generate a formula with <tt>n</tt> quadrature points
QGaussRadauChebyshev(const unsigned int n,
- QGaussRadauChebyshev::EndPoint ep=QGaussRadauChebyshev::left);
+ EndPoint ep=QGaussRadauChebyshev::left);
private:
const EndPoint ep;
const unsigned int n_threads = numbers::invalid_unsigned_int,
const types::subdomain_id subdomain_id = numbers::invalid_subdomain_id,
const types::material_id material_id = numbers::invalid_material_id,
- const KellyErrorEstimator<dim,spacedim>::Strategy strategy = cell_diameter_over_24);
+ const Strategy strategy = cell_diameter_over_24);
/**
* Calls the @p estimate function, see above, with
const unsigned int n_threads = numbers::invalid_unsigned_int,
const types::subdomain_id subdomain_id = numbers::invalid_subdomain_id,
const types::material_id material_id = numbers::invalid_material_id,
- const KellyErrorEstimator<dim,spacedim>::Strategy strategy = cell_diameter_over_24);
+ const Strategy strategy = cell_diameter_over_24);
/**
* Same function as above, but accepts more than one solution vector and
const unsigned int n_threads = numbers::invalid_unsigned_int,
const types::subdomain_id subdomain_id = numbers::invalid_subdomain_id,
const types::material_id material_id = numbers::invalid_material_id,
- const KellyErrorEstimator<dim,spacedim>::Strategy strategy = cell_diameter_over_24);
+ const Strategy strategy = cell_diameter_over_24);
/**
* Calls the @p estimate function, see above, with
const unsigned int n_threads = numbers::invalid_unsigned_int,
const types::subdomain_id subdomain_id = numbers::invalid_subdomain_id,
const types::material_id material_id = numbers::invalid_material_id,
- const KellyErrorEstimator<dim,spacedim>::Strategy strategy = cell_diameter_over_24);
+ const Strategy strategy = cell_diameter_over_24);
/**
const unsigned int n_threads = numbers::invalid_unsigned_int,
const types::subdomain_id subdomain_id = numbers::invalid_subdomain_id,
const types::material_id material_id = numbers::invalid_material_id,
- const KellyErrorEstimator<dim,spacedim>::Strategy strategy = cell_diameter_over_24);
+ const Strategy strategy = cell_diameter_over_24);
/**
const unsigned int n_threads = numbers::invalid_unsigned_int,
const types::subdomain_id subdomain_id = numbers::invalid_subdomain_id,
const types::material_id material_id = numbers::invalid_material_id,
- const KellyErrorEstimator<dim,spacedim>::Strategy strategy = cell_diameter_over_24);
+ const Strategy strategy = cell_diameter_over_24);
/**
const unsigned int n_threads = numbers::invalid_unsigned_int,
const types::subdomain_id subdomain_id = numbers::invalid_subdomain_id,
const types::material_id material_id = numbers::invalid_material_id,
- const KellyErrorEstimator<dim,spacedim>::Strategy strategy = cell_diameter_over_24);
+ const Strategy strategy = cell_diameter_over_24);
/**
const unsigned int n_threads = numbers::invalid_unsigned_int,
const types::subdomain_id subdomain_id = numbers::invalid_subdomain_id,
const types::material_id material_id = numbers::invalid_material_id,
- const KellyErrorEstimator<dim,spacedim>::Strategy strategy = cell_diameter_over_24);
+ const Strategy strategy = cell_diameter_over_24);
/**
* Exception
template <>
std::vector<double>
QGaussRadauChebyshev<1>::get_quadrature_points(const unsigned int n,
- QGaussRadauChebyshev::EndPoint ep)
+ EndPoint ep)
{
std::vector<double> points(n);
template <>
std::vector<double>
QGaussRadauChebyshev<1>::get_quadrature_weights(const unsigned int n,
- QGaussRadauChebyshev::EndPoint ep)
+ EndPoint ep)
{
std::vector<double> weights(n);
template <>
QGaussRadauChebyshev<2>::QGaussRadauChebyshev (const unsigned int n,
- QGaussRadauChebyshev::EndPoint ep)
+ EndPoint ep)
:
Quadrature<2> (QGaussRadauChebyshev<1>(n, static_cast<QGaussRadauChebyshev<1>::EndPoint>(ep)),
QGaussRadauChebyshev<1>(n, static_cast<QGaussRadauChebyshev<1>::EndPoint>(ep))),
template <int dim>
QGaussRadauChebyshev<dim>::QGaussRadauChebyshev (const unsigned int n,
- QGaussRadauChebyshev::EndPoint ep)
+ EndPoint ep)
:
Quadrature<dim> (QGaussRadauChebyshev<dim-1>(n,static_cast<typename QGaussRadauChebyshev<dim-1>::EndPoint>(ep)),
QGaussRadauChebyshev<1>(n,static_cast<QGaussRadauChebyshev<1>::EndPoint>(ep))),
const unsigned int n_threads,
const types::subdomain_id subdomain_id,
const types::material_id material_id,
- const typename KellyErrorEstimator<dim,spacedim>::Strategy strategy)
+ const Strategy strategy)
{
// just pass on to the other function
const std::vector<const InputVector *> solutions (1, &solution);
const unsigned int n_threads,
const types::subdomain_id subdomain_id,
const types::material_id material_id,
- const typename KellyErrorEstimator<dim,spacedim>::Strategy strategy)
+ const Strategy strategy)
{
estimate(StaticMappingQ1<dim,spacedim>::mapping, dof_handler, quadrature, neumann_bc, solution,
error, component_mask, coefficients, n_threads,
const unsigned int n_threads,
const types::subdomain_id subdomain_id,
const types::material_id material_id,
- const typename KellyErrorEstimator<dim,spacedim>::Strategy strategy)
+ const Strategy strategy)
{
// just pass on to the other function
const std::vector<const InputVector *> solutions (1, &solution);
const unsigned int n_threads,
const types::subdomain_id subdomain_id,
const types::material_id material_id,
- const typename KellyErrorEstimator<dim,spacedim>::Strategy strategy)
+ const Strategy strategy)
{
estimate(StaticMappingQ1<dim, spacedim>::mapping, dof_handler, quadrature, neumann_bc, solution,
error, component_mask, coefficients, n_threads,
const unsigned int ,
const types::subdomain_id subdomain_id_,
const types::material_id material_id,
- const typename KellyErrorEstimator<dim,spacedim>::Strategy strategy)
+ const Strategy strategy)
{
#ifdef DEAL_II_WITH_P4EST
if (dynamic_cast<const parallel::distributed::Triangulation<dim,spacedim>*>
const unsigned int n_threads,
const types::subdomain_id subdomain_id,
const types::material_id material_id,
- const typename KellyErrorEstimator<dim,spacedim>::Strategy strategy)
+ const Strategy strategy)
{
// forward to the function with the QCollection
estimate (mapping, dof_handler,
const unsigned int n_threads,
const types::subdomain_id subdomain_id,
const types::material_id material_id,
- const typename KellyErrorEstimator<dim,spacedim>::Strategy strategy)
+ const Strategy strategy)
{
estimate(StaticMappingQ1<dim, spacedim>::mapping, dof_handler, quadrature, neumann_bc, solutions,
errors, component_mask, coefficients, n_threads,
const unsigned int n_threads,
const types::subdomain_id subdomain_id,
const types::material_id material_id,
- const typename KellyErrorEstimator<dim,spacedim>::Strategy strategy)
+ const Strategy strategy)
{
estimate(StaticMappingQ1<dim, spacedim>::mapping, dof_handler, quadrature, neumann_bc, solutions,
errors, component_mask, coefficients, n_threads,