Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Fix initialization of MatrixFreeFunctions::ShapeInfo for many quadrature points #7690

Closed
kronbichler opened this issue Feb 4, 2019 · 1 comment

Comments

@kronbichler
Copy link
Member

As pointed out in #7687, our Lagrange polynomials underflow in case we have more than a few hundreds of points. We should take one of the following three steps:

  • Equip polynomials with variable-length evaluation, i.e., break the Lagrange polynomial sequence (x_i-x_0)(x_i-x_1)...(x_i-x_n) into several parts if we use more than 100 points or so. It would add quite some extensive logic to what is supposed to be a simple problem. If nobody sees a need for polynomials of such a high degree, I would not like to do this.
  • Disable constructing the collocation polynomial if the number of quadrature points is very high, and make sure that we do not use this path anywhere in FEEvaluation.
@masterleinad
Copy link
Member

I was discussing this problem with @stefanmeggendorfer and would also prefer the second solution.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants