vecellipt.cc

#include <config.h>

#include <iostream>

#include <amdis/AMDiS.hpp>
#include <amdis/LocalOperators.hpp>
#include <amdis/ProblemStat.hpp>

using namespace AMDiS;

using ElliptParam   = YaspGridBasis<GRIDDIM, 2,2>;
using ElliptProblem = ProblemStat<ElliptParam>;

int main(int argc, char** argv)
{
  Environment env(argc, argv);

  ElliptProblem prob("ellipt");
  prob.initialize(INIT_ALL);

  AdaptInfo adaptInfo("adapt");

  // ⟨∇u₀,∇v₀⟩
  auto opL = makeOperator(tag::gradtest_gradtrial{}, 1.0);
  prob.addMatrixOperator(opL, 0, 0);

  // ⟨u₁,v₀⟩
  auto opIminus = makeOperator(tag::test_trial{}, -1.0);
  prob.addMatrixOperator(opIminus, 1, 0);

  // ⟨u₁,v₁⟩
  auto opI = makeOperator(tag::test_trial{}, 1.0);
  prob.addMatrixOperator(opI, 1, 1);

  // (-1.0,v₀)
  auto opForce = makeOperator(tag::test{}, [](auto const& x) { return -1.0; }, 0);
  prob.addVectorOperator(opForce, 0);


  // 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, 0, 0, dbcValues);

  prob.buildAfterAdapt(adaptInfo, Flag(0));
  prob.solve(adaptInfo);
  prob.writeFiles(adaptInfo);

  return 0;
}