Comments (4)
Thanks, can you PR the fix and include a source? The lineage here seems to be SciML/Surrogates.jl#126 and #20, but there's no source. So it would be good to have that in the docstring.
from quasimontecarlo.jl.
Thanks, can you PR the fix and include a source? The lineage here seems to be SciML/Surrogates.jl#126 and #20, but there's no source. So it would be good to have that in the docstring.
I can definitely add a source for using them in 1 or 2 dimensions, but there seem to be no references establishing the adequacy of these sequences in more than 2 dimensions, apart from a single blog post by a non-mathematician which contains no formal proofs.
There is, however, literature on using them in either 1 or 2 dimensions, where they're conjectured to have optimal properties. It's not at all clear the generalization in the blog post has these properties, and the sequences also appear to have exceedingly poor 2d projections. Here it is struggling with 8 dimensions and 64 points:
While Sobol has very little trouble:
And Faure handles it easily as well.
from quasimontecarlo.jl.
Can you paste the links to these blog posts and such? This library needs to capture more of this kind of information in the docstrings (for all of the sampling methods 😅)
from quasimontecarlo.jl.
In an extremely funny turn of events I have since learned about component-by-component constructions for QMC lattices. The original, buggy code probably produced better samples than the new "correct" code that follows the method laid out here.
from quasimontecarlo.jl.
Related Issues (20)
- Better Sampling Interface HOT 15
- Calculating discrepancies (esp. L2) and searching for optimized nets
- Error when the number of samples is zero
- https://github.com/SciML/QuasiMonteCarlo.jl/pull/50 breaks downstream HOT 2
- Remove `Distributions` from dependencies HOT 7
- Surrogates model breaks downstream HOT 4
- Base `b` Gray code
- Lazy iteration over points HOT 3
- Evenly distributed samples? HOT 6
- skip=n argument for parallelization? HOT 3
- Halton Sample not defined HOT 4
- Sobol' Sequence Missing First 4 Points HOT 2
- Error in sample() usage HOT 3
- Incorrect Halton Sequence HOT 5
- Integrate QMCGenerators.jl Features
- Estimators and Statistics for RQMC Integraiton HOT 4
- Citation HOT 1
- Improve docstrings
- GridSample appears to sample on a line HOT 1
- Randomized Halton sequence? HOT 2
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 quasimontecarlo.jl.