Comments (4)
Guys SciML/DifferentialEquations.jl#260 (comment)
Chris has done most of the work for you. So writing this up cleanly and checking that this makes sense is the goal. Don't worry about writing up the connection to the spectral and finite element methods.
from diffeqoperators.jl.
@stevenzhangdx I have merged our comments on the two expressions that we both do not fully understand (and your original text is still in the file, just in case you donβt like the merged version. It is on the top of the merged comment but between \iffalse \fi opperators in order to not show up).
from diffeqoperators.jl.
Some comments from editing the document:
- Need a notation for the extension of the domain and the domain itself. This will become especailly important with jumps
- Use
$i=1,... I$ as the grid in the interior. So the solution to the discretized function is$u \in R^I$ if stationary - The total size if
$\hat{I} > I$ . - In general, use a hat when looking at the extension. Then
$\hat{u} \in R^{\hat{I}}$ is the extension of$u \in R^I$ . Replace the current$v$ with$u$ - For the steady state, just drop the
$t$ for now. - When dealing with a univariate function with 2 sides, can denote the extra points as
$I_L$ and$I_H$ so that$\hat{I} = I + I_L + I_H$ . For example, with a univariate diffusion operator,$I_L = I_H = 1$ -
$R$ is a$I \times \hat{I}$ matrix, etc. - Change the basic notation to stack the BC1 and BC2 on the bottom instead of on the top and bottom.
- Rename BC1 and BC2.
- Use
- When we start writing up the examples, we can make the main description more generic
- As we go through the first example, we can convert the main description to be fully generic in terms of
$\hat{I}, I,$ etc. - Move the Dirichlet, etc. to the examples sections to map into the B, Q, etc.
- I think (6) and (7) will be clearer with the generic notation.
- Try to get rid of the
$v_t$ , etc. notation for the discretization?. Too easy to confuse with a partial... Can abuse the$u$ vector and the$u$ function notation? - Lets get rid of the
$f$ and just use$u$ ...
- As we go through the first example, we can convert the main description to be fully generic in terms of
from diffeqoperators.jl.
Steps:
- Go through the current description and update the obvious notation from above, be very clear on the sizes of all vectors/matrices
- Start with the examples 3.1, 3.2 and write down the
$B$ ,$R$ , etc. and move anything not generic out of the main description. - Add in a new example, which is a drift only process, no diffusion. There,
$\hat{I} = I + 1$ which side you add the point on depends on the drift direction if upwind. What do the R and B look like? - For (new) 3.3 and 3.4, write the operators for the drift and diffusion as separate linear operators. So a Q_drift and a Q_diffusion. But a single R and B! The R and B will depend on the upwind direction.
from diffeqoperators.jl.
Related Issues (20)
- Performance issues with nonlinear_diffusion! HOT 4
- Fix symbolic arrays test to allow MTK v5.21.0 HOT 7
- UndefVarError: dereference not defined when precompiling DiffEqOperators (Win10, julia 1.6.1-2) HOT 7
- Version Compatibility? HOT 1
- Error defining second-order differential HOT 8
- The latest versions of LoopVectorization (0.12.80-81) are not compatible and cause errors HOT 1
- Noob Help: Simple, Linear, Inhomogenous PDE HOT 3
- concretization of High Dimension PDE HOT 2
- LoopVectorization break HOT 5
- Unexpected concretization of Laplace operator HOT 3
- v4.35.0 failed to release HOT 2
- Feature request: Support Summation-By-Parts operators HOT 4
- Support Unitful HOT 6
- Incorrect boundary padded vector with composed PeriodicBC on 2d data HOT 6
- Implementing Drift-Diffusion model in higher dimensions HOT 3
- Uninitialized field `opnorm` in JacVecOperator HOT 4
- SplitODEProblem broken for MatrixFreeOperator HOT 3
- ERROR: LoadError: Some tests did not pass: 30 passed, 0 failed, 4 errored, 4 broken. HOT 18
- Warning related to "LoopVectorization.check_args"
- Lots of precompilation errors and noise HOT 7
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 diffeqoperators.jl.