[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] Non-polynomial basis functions
From: |
Torquil Macdonald Sørensen |
Subject: |
Re: [Getfem-users] Non-polynomial basis functions |
Date: |
Sat, 26 Apr 2014 16:22:40 +0200 |
Hi!
Thanks, I have created a FUNC class with methods eval() and derivative(), and used getfem_fem<FUNC> to create a new mesh_fem. This worked great. I verified that I could solve a PDE using this mesh_fem together with an approximate integration method.
Now I would like to create an integration method that is exact for products between polynomials and my new non-polynomial functions. I've looked at the examples in the source code, e.g. exact_simplex() in getfem_integration.cc. It uses a class "simplex_poly_integration_" that inherits from "poly_integration", and passes it to the constructor of "integration_method".
Looking at "integration_method" in getfem_integration.h, it seems to me that in order to have "im_type = IM_EXACT", my mesh_fem must be polynomial? That is my impression, since the only constructor of "integration_method" that gives im_type = IM_EXACT is the one that takes an argument of type ppoly_integration.
Does this mean that I cannot create an exact integration method for non-polynomial functions?
Best regards
Torquil Sørensen