4 releases
0.2.2 | May 15, 2024 |
---|---|
0.2.1 | Apr 7, 2024 |
0.2.0 | Apr 4, 2024 |
0.1.0 | Mar 19, 2024 |
#127 in Science
28 downloads per month
12KB
137 lines
Rugfield
Rugfield is a Rust library for generating Gaussian Random Fields (GRFs) using the circulant embedding method. It provides an efficient and easy-to-use implementation for simulating GRFs with various kernel functions.
Features
- Generate GRFs using the circulant embedding method
- Support for various kernels.
SquaredExponential
Matern
LocalPeriodic
RationalQuadratic
- Efficient computation using the
rustfft
library - Optional serialization support with the
serde
feature
Installation
Add the following to your Cargo.toml
file:
[dependencies]
rugfield = "0.2.2"
Usage
Here's a simple example of how to use Rugfield to generate a GRF with a squared exponential kernel:
use rugfield::{grf, Kernel};
fn main() {
let n = 100;
let kernel = Kernel::SquaredExponential(0.1);
let grf_data = grf(n, kernel);
// Plot the GRF data
// ...
}
This code generates a GRF with a squared exponential kernel and a length scale of 0.1. The resulting GRF data is stored in the grf_data
vector.
For a complete example, see the examples/squared_exponential.rs
file:
use peroxide::fuga::*;
use rugfield::{grf, Kernel::SquaredExponential};
fn main() -> Result<(), Box<dyn Error>> {
let x_max = 100.0;
let x_min = 0.0;
let sigma = 0.1;
let n = 1000;
let samples = 8;
let kernel = SquaredExponential(sigma);
let x = linspace_with_precision(x_min, x_max, n, 2);
let grfs = (0..samples).map(|_| grf(n, kernel)).collect::<Vec<_>>();
// Plot the GRF data
// ...
Ok(())
}
With a specific random number generator:
use peroxide::fuga::*;
use rugfield::{grf_with_rng, Kernel::SquaredExponential};
fn main() -> Result<(), Box<dyn Error>> {
let n = 100;
let kernel = SquaredExponential(0.1);
let mut rng = stdrng_from_seed(42);
let grf_data = grf_with_rng(&mut rng, n, kernel);
// ...
}
The above code generates multiple GRFs with a squared exponential kernel and plots the resulting data. Here's an example output:
Documentation
For detailed documentation and API reference, please refer to the Rustdoc documentation.
License
Rugfield is licensed under the MIT License.
Contributing
Contributions are welcome! If you find any issues or have suggestions for improvements, please open an issue or submit a pull request on the GitHub repository.
Acknowledgments
Rugfield was inspired by the paper "An Effective Method for Simulating Gaussian Random Fields" by Grace Chan (1999).
We would like to express our gratitude to the authors of the rustfft
library for their excellent work, which has been instrumental in the development of Rugfield.
Dependencies
~7.5MB
~148K SLoC