MFEMVectorUtils.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#ifdef MOOSE_MFEM_ENABLED

#pragma once
#include <vector>

#include "libmesh/point.h"
#include "libmesh/ignore_warnings.h"
#include "mfem.hpp"
#include "libmesh/restore_warnings.h"

namespace Moose::MFEM
{
inline std::size_t
MFEMIndex(const std::size_t i_dim,
          const std::size_t i_point,
          const std::size_t num_dims,
          const std::size_t num_points,
          const mfem::Ordering::Type ordering)
{
  if (ordering == mfem::Ordering::byNODES)
    return mfem::Ordering::Map<mfem::Ordering::byNODES>(num_points, num_dims, i_point, i_dim);
  else
    return mfem::Ordering::Map<mfem::Ordering::byVDIM>(num_points, num_dims, i_point, i_dim);
}

libMesh::Point libMeshPointFromMFEMVector(const mfem::Vector & vec);

mfem::Vector libMeshPointsToMFEMVector(const std::vector<libMesh::Point> & points,
                                       const unsigned int num_dims,
                                       const mfem::Ordering::Type ordering);

}

#endif