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

This commit adds two jupyter notebooks exhibiting the Ultraspherical polynomials. #476

Merged
merged 4 commits into from
Jun 20, 2023

Conversation

kyle-augustson
Copy link
Contributor

The file Scalar_Diffusion.ipnb has some mathematical details about the Ultraspherical polynomials regarding how to build sparse spectral representations of operators such as derivatives. It also details how to maintain sparsity of those operators even if they are multiplied by non-constant coefficients. The example problem is the 1D diffusion equation both with and without a spatially variable diffusion coefficient for scalar quantity q. This uses the second order Crank-Nicholson method for time stepping. The file Reaction_Diffusion.ipynb shows how to construct a non-linear PDE using the Newell–Whitehead-Segel equation as an example.

…nb has some

mathematical details about the Ultraspherical polynomials regarding how to build
sparse spectral representations of operators such as derivatives.  It also details
how to maintain sparsity of those operators even if they are multiplied by
non-constant coefficients. The example problem is the 1D diffusion equation both
with and without a spatially variable diffusion coefficient for scalar quantity q.
This uses the second order Crank-Nicholson method for time stepping. The file
Reaction_Diffusion.ipynb shows how to construct a non-linear PDE using the
Newell–Whitehead-Segel equation as an example.
… pointer to ffmpeg).

Scalar_Diffusion.ipynb has quite a few more comments to help guide the reader as to what is going on locally.
Copy link
Contributor

@feathern feathern left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Awesome. This should serve as the basis for some solid documentation once we get this formalism implemented into the main code.

@feathern feathern merged commit 232d067 into geodynamics:main Jun 20, 2023
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