search

optimad / bitpit

Open source library for scientific HPC
http://optimad.github.io/bitpit/
GNU Lesser General Public License v3.0
117 stars 33 forks source link

discretization: low order cell basis can be evaluated on all type of cells #254

Closed andrea-iob closed 2 years ago

andrea-iob commented 2 years ago

Polynomials of degree 0 and 1 are supported on all type of cells, for pixels and voxels also degree 3 is supported.

Note: coefficients are evaluated as the volume average of the Taylor series expansion coefficients. In order to support higher order degrees on all type of cells, we will need to implement numerical integration of the expansions terms over the volume of the cell. This can be done using Gauss quadrature rules and evaluating the terms on the integration points using the function that evaluate the basis function on a specified point.

andrea-iob commented 2 years ago

For some reasons the review done by @lucask81 is not enough to allow me to merge the pull request. @marcocisternino Can you have a look at these changes?

andrea-iob commented 2 years ago

(We require a review from a user that has write access and @lucask81 doesn't.)