gsBarrierPatch.h

Go to the documentation of this file.

#pragma once

#include <gsIO/gsOptionList.h>

#ifdef gsHLBFGS_ENABLED
#include <gsHLBFGS/gsHLBFGS.h>
#endif

using namespace gismo;

namespace gismo
{
template<short_t d, typename T=real_t>
class gsBarrierPatch
{
 public:
  explicit gsBarrierPatch(const gsMultiPatch<T> &mp, const gsDofMapper &mapper);

  explicit gsBarrierPatch(const gsMultiPatch<T> &mp, bool patchWise = true);

  void setMapper(const gsDofMapper &mapper) { m_mapper = mapper; };

  void compute();

  const gsMultiPatch<T> &result() const { return m_mp; };

  gsOptionList &options() { return m_options; };

  void defaultOptions();

 private:
  void _makeMapper();

  gsDofMapper _makeMapperOnePatch(const gsGeometry<T>&currPatch) const;

  void _makeMapperGlobalPatches();

  void _makeMapperLocalPatches();

  void logMapperInformation();

  mutable gsExprEvaluator<T> m_evaluator;
  mutable gsExprAssembler<T> m_assembler;
  gsDofMapper m_mapper;

  mutable gsMultiPatch<T> m_mp;
  gsMultiBasis<T> m_mb;
  gsMultiPatch<T> m_bRep;

  std::string m_filename;
  T m_boxsize = 1.0;

  gsVector<T, d> m_boundingBoxLeftBottomCorner;
  gsVector<T, d> m_scalingVec;

  size_t m_freeInterface = 1;
  gsOptionList m_options;
};


#ifdef GISMO_WITH_PYBIND11

  void pybind11_init_gsBarrierPatch2(pybind11::module &m);
  void pybind11_init_gsBarrierPatch3(pybind11::module &m);

#endif // GISMO_WITH_PYBIND11

}// namespace gismo

#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsBarrierPatch.hpp)
#endif