Comments (2)
I think your question is a little too broad to spark any new ideas. A solver library can be created in a lot of different ways. SpectralDNS is one library that initially only contained Navier-Stokes, and then grew from that. I'm sure it's not the ideal library, but it seems to be working:-) One thing is that I have considered rewriting it with a nicer interface to integrators, as I don't really like the way it is now. The integrators that are in Shenfun I find nicer. Did you look at those yet? I also want to add some adaptive integrators here eventually. The good thing with these is that they separate the right hand side of the equations into linear and nonlinear parts completely generically.
When it comes to gracefully switching between bases I'm not so sure this should be possible. If you have a solution in a Chebyshev basis and want it in Legendre, then you need to make a projection. There's really no other switch allowed. One thing that actually is allowed is to switch directly from a Dirichlet basis to a regular Chebyshev basis, because these bases can use the same quadrature points in physical space. There is no simple mechanism implemented to do this switch yet though. You need to create both bases and then copy the data from one physical array in one basis, to an array in the other one.
One option is to implement solvers in spectralDNS in the current framework. Have you thought about that possibility?
from shenfun.
yeah, I agree, this issue is a bit vague. It did steered some useful replies from your part though ;)
I'll look at integrators more closely, this looks interesting.
Regarding the bases switching issue, I am afraid I was not accurate enough in my initial post.
In a most simple situation, imagine we want to write a solver for an elliptic problem and we would like to let the basis choice as an option to the user. My intuition is that this may have an important impact on the architecture of the library but it's tough for me to estimate what this would be precisely.
I like your last suggestion about trying spectralDNS architecture.
Would it be difficult to build a minimalist version of spectralDNS?
from shenfun.
Related Issues (20)
- Dx not working with 1D mixed function space HOT 4
- Issue on Stokes equations with doubly periodic boundary conditions. HOT 2
- Error when evaluating derivatives for nonstandard domain HOT 7
- Installation problem: error in "cimport fastgl_wrap" HOT 2
- Old shenfun script for 2D Poisson equation dosn't work with new version HOT 2
- 1d poisson HOT 1
- Bug in FunctionSpace with bc as dictionary HOT 4
- Solve DGL with non-constant (discrete) variable? HOT 5
- AttributeError: 'tuple' object has no attribute 'copy' HOT 1
- Install problem HOT 1
- problems in solving 2D Allen Cahn equation HOT 2
- Two questions in the demo MixedPoisson.py HOT 7
- How to compute the numerical solution on the general points HOT 4
- 1. Solver for multidimensional problem 2. Product with variable coefficients HOT 4
- Two questions about the Legendre-Galerkin method for the Cahn-Hilliard equation HOT 2
- How to compute the integral about the nonlinear term in the Navier-Stokes equation? HOT 2
- How to assemble the right-hand side vector? HOT 3
- Understanding the chebyshev solver in poisson2ND.py HOT 3
- Time dependent boundary conditions in segments for Rayleigh-Benard 2D HOT 2
- Implementing the one-dimensional Kuramoto-Sivashinsky system HOT 3
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from shenfun.