n_vector.h

Go to the documentation of this file.

// ---------------------------------------------------------------------
//
// Copyright (C) 2020 - 2022 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.md at
// the top level directory of deal.II.
//
// ---------------------------------------------------------------------
 
 
#ifndef dealii_sundials_n_vector_h
#define dealii_sundials_n_vector_h
 
#include <deal.II/base/config.h>
 
#ifdef DEAL_II_WITH_SUNDIALS
#  include <sundials/sundials_nvector.h>
 
#  include <functional>
#  include <memory>
 
DEAL_II_NAMESPACE_OPEN
 
#  ifndef DOXYGEN
namespace SUNDIALS
{
  namespace internal
  {
    template <typename VectorType>
    class NVectorView;
  }
} // namespace SUNDIALS
#  endif
 
namespace SUNDIALS
{
  namespace internal
  {
#  if DEAL_II_SUNDIALS_VERSION_GTE(6, 0, 0)
    template <typename VectorType>
    NVectorView<VectorType>
    make_nvector_view(VectorType &vector, SUNContext nvector_context);
#  else
    template <typename VectorType>
    NVectorView<VectorType>
    make_nvector_view(VectorType &vector);
#  endif
 
    template <typename VectorType>
    VectorType *
    unwrap_nvector(N_Vector v);
 
    template <typename VectorType>
    const VectorType *
    unwrap_nvector_const(N_Vector v);
 
    template <typename VectorType>
    class NVectorView
    {
    public:
      NVectorView() = default;
 
#  if DEAL_II_SUNDIALS_VERSION_GTE(6, 0, 0)
      NVectorView(VectorType &vector, SUNContext nvector_context);
#  else
      NVectorView(VectorType &vector);
#  endif
 
      NVectorView(NVectorView &&) noexcept = default;
 
      NVectorView &
      operator=(NVectorView &&) noexcept = default;
 
      NVectorView(const NVectorView &) = delete;
 
      NVectorView &
      operator=(const NVectorView &) = delete;
 
      ~NVectorView() = default;
 
      operator N_Vector() const;
 
      N_Vector
      operator->() const;
 
    private:
      std::unique_ptr<_generic_N_Vector, std::function<void(N_Vector)>>
        vector_ptr;
    };
  } // namespace internal
} // namespace SUNDIALS
 
DEAL_II_NAMESPACE_CLOSE
 
#endif
#endif