Skip to content
Snippets Groups Projects
Commit 2624ac39 authored by Porrmann, Maik's avatar Porrmann, Maik
Browse files

remove testmibasis

parent 87fac121
No related branches found
No related tags found
No related merge requests found
add_executable("C1Elements" C1Elements.cc)
target_link_dune_default_libraries("C1Elements")
add_executable("testMIBasis" testMIBasis.cc)
target_link_dune_default_libraries("testMIBasis")
\ No newline at end of file
#include <fstream>
#include <utility>
#include <dune/geometry/type.hh>
#include <dune/localfunctions/utility/polynomialbasis.hh>
#include <dune/localfunctions/hermite.hh>
#include <dune/localfunctions/argyris.hh>
int main(int argc, char *argv[])
{
constexpr unsigned int dim = 2;
typedef Dune::MultiIndex<dim, double> MI;
typedef Dune::MonomialBasis<Dune::GeometryTypes::simplex(dim).toId(), MI> MIBasis;
MIBasis basis(5);
Dune::FieldVector<MI, dim> x;
/*
* Init MultiIndices
* x[0]=(1,0,0) x
* x[1]=(0,1,0) y
* x[2]=(0,0,1) z
*/
for (unsigned int i = 0; i < dim; ++i)
x[i].set(i, 1);
std::vector<MI> val(basis.size());
// val now contains all monomials in $n$ dimensions with degree $\leq order+1$
basis.evaluate(x, val);
std::cout << "Monomial Basis " << std::endl
<< val << std::endl;
Dune::Impl::HermiteLocalBasis<MI, MI, 2> hermite;
val.resize(hermite.size());
hermite.evaluate(x, val);
std::cout << "Hermite Basis 2d " << std::endl << val << std::endl;
Dune::Impl::ArgyrisLocalBasis<MI, MI> argyris;
val.resize(argyris.size());
argyris.evaluate(x, val);
std::cout << "Argyris Basis " << std::endl << val << std::endl;
std::vector<decltype(argyris)::Traits::JacobianType> jacs;
argyris.evaluateJacobian(x, jacs);
std::cout << "Argyris Basis Jacobian " << std::endl << jacs << std::endl;
return 0;
}
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment