hermitetest.cc

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#include <config.h>

#include <iostream>

#include <dune/common/exceptions.hh>
#include <dune/common/parallel/mpihelper.hh>

#include <dune/grid/yaspgrid.hh>
#include <dune/grid/uggrid.hh>
#include <dune/grid/onedgrid.hh>
#include <dune/grid/io/file/vtk/subsamplingvtkwriter.hh>
#include <dune/grid/utility/structuredgridfactory.hh>
#include <dune/grid/uggrid.hh>
#include <dune/functions/functionspacebases/lagrangebasis.hh>
#include <dune/functions/hermitebasis.hh>

#include <dune/functions/functionspacebases/test/basistest.hh>
#include <dune/functions/gridfunctions/discreteglobalbasisfunction.hh>

using namespace Dune;
using namespace Dune::Functions;



int main (int argc, char* argv[])
{
  Dune::MPIHelper::instance(argc, argv);

  Dune::TestSuite test;



  using namespace Dune::Functions::BasisFactory;

  { // 1d
    std::unique_ptr<OneDGrid> grid
      = StructuredGridFactory<OneDGrid>::createSimplexGrid({0.}, {1.}, {10});

    auto gridView = grid->levelGridView(0);

    {
      auto preBasisFactory = hermite();
      auto preBasis = preBasisFactory.template makePreBasis<FlatMultiIndex<std::size_t>>(gridView);
      using PB=decltype(preBasis);
      PB::Node node(preBasis.makeNode());
      auto node2 = preBasis.makeNode();
      // node.complain();
      // node2.complain();

    }
    {
      auto basis = makeBasis(gridView, lagrange<2>());
      test.subTest(checkBasis(basis, EnableContinuityCheck()));

      
    }


  }
  /*
  {   //2d
    auto grid = StructuredGridFactory<UGGrid<2> >::createSimplexGrid({0.,0.}, {1.,1.}, {10,10});

    auto gridView = grid->leafGridView();
    using GridView = decltype(gridView);
    {
      using namespace Dune::Functions::BasisFactory;
      // auto basis = makeBasis(gridView, hermite());
      // test.subTest(checkBasis(basis, EnableContinuityCheck()));
    }
  }
  */

  return test.exit();
}