FyeldGenerator repository
This package provides with a quick way of generating random field having a specified power spectrum.
Example
from FyeldGenerator import generate_field
import matplotlib.pyplot as plt
import numpy as np
# Helper that generates power-law power spectrum
def Pkgen(n):
def Pk(k):
return np.power(k, -n)
return Pk
# Draw samples from a normal distribution
def distrib(shape):
a = np.random.normal(loc=0, scale=1, size=shape)
b = np.random.normal(loc=0, scale=1, size=shape)
return a + 1j * b
shape = (512, 512)
field = generate_field(distrib, Pkgen(2), shape)
plt.imshow(field, cmap="seismic")
Install
It is now on pypi! For the "official" release, use:
pip install FyeldGenerator
For the latest release:
pip install -e git+https://github.com/cphyc/FyeldGenerator.git#egg=FyeldGenerator
License
This work is licensed under the CC-BY-SA license. You are allowed to copy, modify and distribute it as long as you keed the license. See more in the LICENSE file.