#include <deal.II/base/function_parser.h>
#include <deal.II/base/utilities.h>
#include <deal.II/base/conditional_ostream.h>
+#include <deal.II/base/std_cxx11/array.h>
#include <deal.II/lac/vector.h>
#include <deal.II/lac/compressed_sparsity_pattern.h>
template <typename InputVector, typename Number>
static
void compute_flux_matrix (const InputVector &W,
- std::array <std::array <Number, dim>, EulerEquations<dim>::n_components > &flux)
+ std_cxx11::array <std_cxx11::array <Number, dim>, EulerEquations<dim>::n_components > &flux)
{
// First compute the pressure that appears in the flux matrix, and then
// compute the first <code>dim</code> columns of the matrix that
const InputVector &Wplus,
const InputVector &Wminus,
const double alpha,
- std::array < Number, n_components> &normal_flux)
+ std_cxx11::array < Number, n_components> &normal_flux)
{
- std::array <std::array <Number, dim>, EulerEquations<dim>::n_components > iflux, oflux;
+ std_cxx11::array <std_cxx11::array <Number, dim>, EulerEquations<dim>::n_components > iflux, oflux;
compute_flux_matrix (Wplus, iflux);
compute_flux_matrix (Wminus, oflux);
template <typename InputVector, typename Number>
static
void compute_forcing_vector (const InputVector &W,
- std::array < Number, n_components> &forcing)
+ std_cxx11::array < Number, n_components> &forcing)
{
const double gravity = -1.0;
// later easily be computed from it:
std::vector <
- std::array <std::array <Sacado::Fad::DFad<double>, dim>, EulerEquations<dim>::n_components >
+ std_cxx11::array <std_cxx11::array <Sacado::Fad::DFad<double>, dim>, EulerEquations<dim>::n_components >
> flux(n_q_points);
std::vector <
- std::array <std::array <double, dim>, EulerEquations<dim>::n_components >
+ std_cxx11::array <std_cxx11::array <double, dim>, EulerEquations<dim>::n_components >
> flux_old(n_q_points);
- std::vector < std::array< Sacado::Fad::DFad<double>, EulerEquations<dim>::n_components> > forcing(n_q_points);
+ std::vector < std_cxx11::array< Sacado::Fad::DFad<double>, EulerEquations<dim>::n_components> > forcing(n_q_points);
- std::vector < std::array< double, EulerEquations<dim>::n_components> > forcing_old(n_q_points);
+ std::vector < std_cxx11::array< double, EulerEquations<dim>::n_components> > forcing_old(n_q_points);
for (unsigned int q=0; q<n_q_points; ++q)
{
// that does so, we also need to determine the Lax-Friedrich's stability
// parameter:
- std::vector< std::array < Sacado::Fad::DFad<double>, EulerEquations<dim>::n_components> > normal_fluxes(n_q_points);
- std::vector< std::array < double, EulerEquations<dim>::n_components> > normal_fluxes_old(n_q_points);
+ std::vector< std_cxx11::array < Sacado::Fad::DFad<double>, EulerEquations<dim>::n_components> > normal_fluxes(n_q_points);
+ std::vector< std_cxx11::array < double, EulerEquations<dim>::n_components> > normal_fluxes_old(n_q_points);
double alpha;