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

Orthogonalise polynomial sets #495

Merged
merged 9 commits into from
May 17, 2022
Merged

Orthogonalise polynomial sets #495

merged 9 commits into from
May 17, 2022

Conversation

mscroggs
Copy link
Member

@mscroggs mscroggs commented May 13, 2022

Orthogonalise the polynomial basis during element construction.

This reduces the condition number of the dual matrices. eg for serendipity elements (built-in = current elements on main, orthogonalised = after this PR):

image

This has a knock on effect of reducing the error when solving problems, eg:

image

@mscroggs mscroggs marked this pull request as draft May 13, 2022 09:14
@mscroggs mscroggs changed the title Add empty function Orthogonalise polynomial sets May 13, 2022
@mscroggs mscroggs marked this pull request as ready for review May 13, 2022 09:29
@mscroggs mscroggs merged commit b257d86 into main May 17, 2022
@mscroggs mscroggs deleted the mscroggs/orthogonalise branch May 17, 2022 07:03
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

Successfully merging this pull request may close these issues.

2 participants