global.h

Go to the documentation of this file.

/*
  Copyright (C) 2011 - 2024 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 <aspect/config.h>
#include <aspect/citation_info.h>
 
#include <deal.II/base/mpi.h>
#include <deal.II/base/exceptions.h>
 
DEAL_II_DISABLE_EXTRA_DIAGNOSTICS
 
#include <boost/archive/binary_oarchive.hpp>
#include <boost/archive/binary_iarchive.hpp>
 
DEAL_II_ENABLE_EXTRA_DIAGNOSTICS
 
#include <boost/container/small_vector.hpp>
 
#ifdef ASPECT_USE_TPETRA
#include <deal.II/lac/trilinos_tpetra_vector.h>
#include <deal.II/lac/trilinos_tpetra_block_vector.h>
#include <deal.II/lac/trilinos_tpetra_sparse_matrix.h>
#include <deal.II/lac/trilinos_tpetra_block_sparse_matrix.h>
#else
#include <deal.II/lac/trilinos_sparse_matrix.h>
#include <deal.II/lac/trilinos_block_sparse_matrix.h>
#endif
 
 
#include <aspect/compat.h>
 
namespace aspect
{
  using namespace dealii;
 
  namespace constants
  {
    constexpr double year_in_seconds = 60*60*24*365.2425;
 
    constexpr double celsius_to_kelvin = 273.15;
 
    constexpr double degree_to_radians = ::numbers::PI / 180.;
 
    constexpr double radians_to_degree = 180. / ::numbers::PI;
 
    constexpr double gas_constant = 8.3144621;
    constexpr double avogadro = 6.02214129e23;
    constexpr double big_g = 6.67430e-11;
 
    namespace earth
    {
 
      namespace masses
      {
        constexpr double planet = 5.9736e24;
        constexpr double core = 1.932e24;
        constexpr double mantle = 4.043e24;
      }
 
      namespace iasp91_radii
      {
        constexpr double inner_core = 1.2171e6;
        constexpr double core = 3.482e6;
        constexpr double lower_mantle = 5.711e6;
        constexpr double planet = 6.371e6;
      }
 
      namespace prem_gravity
      {
        constexpr double icb = 4.4002;
        constexpr double cmb = 10.6823;
        constexpr double ulmb = 10.0143;
        constexpr double surface = 9.8156;
      }
 
      constexpr double surface_gravity = 9.80665;
    }
 
    namespace mars
    {
      namespace radii
      {
        constexpr double planet = 3.3895e6;
        constexpr double core = 1.794e6;
      }
      constexpr double surface_gravity = 3.711;
    }
  }
 
  using constants::year_in_seconds;
 
 
  using iarchive = boost::archive::binary_iarchive;
 
  using oarchive = boost::archive::binary_oarchive;
 
  class QuietException {};
 
  template <class T, unsigned int N = 100>
  using small_vector = boost::container::small_vector<T, N>;
 
  namespace LinearAlgebra
  {
#ifdef ASPECT_USE_TPETRA
    using Vector = ::LinearAlgebra::TpetraWrappers::Vector<double>;
 
    using BlockVector = ::LinearAlgebra::TpetraWrappers::BlockVector<double>;
 
    using SparseMatrix = ::LinearAlgebra::TpetraWrappers::SparseMatrix<double>;
 
    using BlockSparseMatrix = ::LinearAlgebra::TpetraWrappers::BlockSparseMatrix<double>;
#else
    using Vector = ::TrilinosWrappers::MPI::Vector;
 
    using BlockVector = ::TrilinosWrappers::MPI::BlockVector;
 
    using SparseMatrix = ::TrilinosWrappers::SparseMatrix;
 
    using BlockSparseMatrix = ::TrilinosWrappers::BlockSparseMatrix;
#endif
  }
 
  template <class Stream>
  void print_aspect_header(Stream &stream);
}
 
 
 
#endif