Comments (3)
@chenwilliam77 A less focused project than EconPDEs working on these things is https://github.com/JuliaDiffEq/DiffEqOperators.jl/ which implements fornberg, different upwinding, etc. and the intention is to get very solid lazy operator composition as well. It is implemented largely for matrix-free methods using convolutions (which is a little different than the stencils that @matthieugomez is using) but there are concretizations to structured matrices where appropriate
If you are interested in doing any julia code on these topics, we would love to have you work on using these for finance applications.
from econpdes.jl.
Hi Matthieu,
Perhaps another alternative to the way you're currently implementing the finite differences is to use the Fornberg algorithm. It generates weights for one-dimensional finite differences of arbitrary order and accuracy, for both uniform and non-uniform grids. Then to approximate a mixed derivative, you'd just apply the one-dimensional finite difference in both directions. In the case of a uniform grid, this will yield the standard formula for a mixed derivative approximation. Note that I'm fairly certain this is the way Mathematica implements its finite differences (or so it seemed when I implemented a replication of Di Tella (2017) in Matlab recently).
A Matlab implementation of the Fornberg algorithm is posted here, and if you're interested, then I can try coding it up in Julia for you because I would like to use your tools. My personal research uses continuous-time macrofinance methods, and I'd prefer not to re-invent the wheel!
Best,
William
from econpdes.jl.
Thanks for the example of mixed derivatives. Before I thought the monotonicity condition of Barles-Souganidis' framework is violated by the mixed derivative, so it would've been a problem for the finite difference method.
So it seems with a fully implicit algorithm we can still handle the cross derivative term (at least when it's not too large)?
from econpdes.jl.
Related Issues (20)
- Vector vs Elementwise
- Aiyagari Continuous Time eample using EconPDEs HOT 2
- Simulate results from EconPDEs HOT 1
- Finite horizon boundary problem HOT 3
- CambpellCochrane.jl HOT 3
- Finite horizon consumption-saving problem HOT 2
- Bug in AchdouHanLasryLionsMoll one asset example HOT 2
- TagBot trigger issue HOT 4
- Does this work on the latest version of Julia? HOT 1
- Simple Firm Investment Problem HOT 2
- Neoclassical Growth Model Example
- Syntax error? HOT 2
- Cite this package
- More than two states HOT 1
- Order vs named
- DiTella Model HOT 3
- Allow a way to specify different boundary conditions
- Solve Kolmogorov Forward Equation together with an HJB HOT 6
- Register HOT 1
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 econpdes.jl.