vecellipt.cc

 
#include <config.h>
 
#include <iostream>
 
#include <amdis/AMDiS.hpp>
#include <amdis/LocalOperators.hpp>
#include <amdis/ProblemStat.hpp>
 
using namespace AMDiS;
 
int main(int argc, char** argv)
{
  Environment env(argc, argv);
 
  using Grid = Dune::YaspGrid<GRIDDIM>;
  auto grid = MeshCreator<Grid>{"elliptMesh"}.create();
 
  using namespace Dune::Functions::BasisFactory;
  ProblemStat prob("ellipt", *grid, power<2>(lagrange<2>()));
  prob.initialize(INIT_ALL);
 
  AdaptInfo adaptInfo("adapt");
 
  // ⟨∇u₁,∇v₁⟩
  auto opL = makeOperator(tag::gradtest_gradtrial{}, 1.0);
  prob.addMatrixOperator(opL, 1, 1);
 
  // ⟨u₀,v₀⟩ and ⟨u₁,v₀⟩
  auto opM = makeOperator(tag::test_trial{}, 1.0);
  prob.addMatrixOperator(opM, 0, 0);
  auto opM2 = makeOperator(tag::test_trial{}, -1.0);
  prob.addMatrixOperator(opM2, 0, 1);
 
  // (f,v₁) with f ≡ -1
  auto opForce = makeOperator(tag::test{}, [](auto const& x) { return -1.0; }, 0);
  prob.addVectorOperator(opForce, 1);
 
 
  // set boundary condition
  auto predicate = [](auto const& x){ return x[0] < 1.e-8 || x[1] < 1.e-8; }; // define boundary
  auto dbcValues = [](auto const& x){ return 0.0; }; // set value
  prob.addDirichletBC(predicate, 1, 1, dbcValues);
 
  prob.buildAfterAdapt(adaptInfo, Flag(0));
  prob.solve(adaptInfo);
  prob.writeFiles(adaptInfo);
 
  return 0;
}