global.h

Go to the documentation of this file.

/*
  Copyright (C) 2011 - 2018 by the authors of the ASPECT code.

  This file is part of ASPECT.

  ASPECT is free software; you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation; either version 2, or (at your option)
  any later version.

  ASPECT is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with ASPECT; see the file LICENSE.  If not see
  <http://www.gnu.org/licenses/>.
*/


#ifndef _aspect_global_h
#define _aspect_global_h

#include <deal.II/base/mpi.h>

DEAL_II_DISABLE_EXTRA_DIAGNOSTICS

#ifdef ASPECT_USE_PETSC
#  include <deal.II/lac/petsc_parallel_block_vector.h>
#  include <deal.II/lac/petsc_parallel_block_sparse_matrix.h>
#  include <deal.II/lac/petsc_precondition.h>
#else
#  include <deal.II/lac/trilinos_parallel_block_vector.h>
#  include <deal.II/lac/trilinos_block_sparse_matrix.h>
#  include <deal.II/lac/trilinos_precondition.h>
#endif

#include <deal.II/lac/generic_linear_algebra.h>

#include <boost/archive/binary_oarchive.hpp>
#include <boost/archive/binary_iarchive.hpp>

DEAL_II_ENABLE_EXTRA_DIAGNOSTICS


#include <aspect/compat.h>

namespace aspect
{
  namespace constants
  {
    extern const double year_in_seconds;

    extern const double celsius_to_kelvin;

    extern const double gas_constant;
    extern const double avogadro;
    extern const double big_g;

    namespace earth
    {

      namespace masses
      {
        extern const double planet;
        extern const double core;
        extern const double mantle;
      }

      namespace iasp91_radii
      {
        extern const double inner_core;
        extern const double core;
        extern const double lower_mantle;
        extern const double planet;
      }

      namespace prem_gravity
      {
        extern const double icb;
        extern const double cmb;
        extern const double ulmb;
        extern const double surface;
      }

      extern const double surface_gravity;
    }

    namespace mars
    {

      namespace radii
      {
        extern const double planet;
        extern const double core;
      }
      extern const double surface_gravity;
    }
  }

  using constants::year_in_seconds;


  typedef boost::archive::binary_iarchive iarchive;

  typedef boost::archive::binary_oarchive oarchive;

  class QuietException {};


  namespace LinearAlgebra
  {
#ifdef ASPECT_USE_PETSC

    typedef ::PETScWrappers::MPI::Vector Vector;

    typedef ::PETScWrappers::MPI::BlockVector BlockVector;

    typedef ::PETScWrappers::MPI::SparseMatrix SparseMatrix;

    typedef ::PETScWrappers::MPI::BlockSparseMatrix BlockSparseMatrix;

    typedef ::PETScWrappers::PreconditionerBase PreconditionBase;

    typedef ::PETScWrappers::PreconditionBoomerAMG PreconditionAMG;

    typedef ::PETScWrappers::PreconditionICC PreconditionIC;

    typedef ::PETScWrappers::PreconditionBlockJacobi PreconditionILU;

    typedef ::PETScWrappers::PreconditionJacobi PreconditionJacobi;

    typedef ::BlockDynamicSparsityPattern BlockDynamicSparsityPattern;

    typedef ::DynamicSparsityPattern DynamicSparsityPattern;
#else

    typedef ::TrilinosWrappers::MPI::Vector Vector;

    typedef ::TrilinosWrappers::MPI::BlockVector BlockVector;

    typedef ::TrilinosWrappers::SparseMatrix SparseMatrix;

    typedef ::TrilinosWrappers::BlockSparseMatrix BlockSparseMatrix;

    typedef ::TrilinosWrappers::PreconditionBase PreconditionBase;

    typedef ::TrilinosWrappers::PreconditionAMG PreconditionAMG;

    typedef ::TrilinosWrappers::PreconditionIC PreconditionIC;

    typedef ::TrilinosWrappers::PreconditionILU PreconditionILU;

    typedef ::TrilinosWrappers::PreconditionJacobi PreconditionJacobi;

    typedef ::TrilinosWrappers::BlockSparsityPattern BlockDynamicSparsityPattern;

    typedef ::TrilinosWrappers::SparsityPattern DynamicSparsityPattern;
#endif
  }
}


template <class Stream>
void print_aspect_header(Stream &stream);

#define ASPECT_INSTANTIATE(INSTANTIATIONS) \
  INSTANTIATIONS(2) \
  INSTANTIATIONS(3)

#endif