/// <summary>
/// Returns the sampler for a specified rng, data type and distribution spec
/// </summary>
/// <param name="rng">The random stream</param>
/// <param name="spec">The distribution specifier</param>
/// <typeparam name="T">The sample point type</typeparam>
public static ISampler <T> Sampler <T>(this MklRng rng, IDistributionSpec <T> spec)
where T : unmanaged
{
var sampler = default(ISampler <T>);
switch (spec.DistKind)
{
case DistKind.Uniform:
sampler = rng.UniformSampler <T>(UniformSpec.From(spec));
break;
case DistKind.UniformBits:
sampler = rng.UniformBitsSampler <T>(UniformBitsSpec.From(spec));
break;
case DistKind.Bernoulli:
sampler = rng.BernoulliSampler <T>(BernoulliSpec.From(spec));
break;
case DistKind.Gaussian:
sampler = rng.GaussianSampler <T>(GaussianSpec.From(spec));
break;
default:
throw unsupported <T>();
}
return(sampler);
}
public static ISampler <T> UniformBitsSampler <T>(this MklRng rng, UniformBitsSpec <T>?spec = null)
where T : unmanaged
{
var sampler = default(ISampler <T>);
var _spec = spec ?? UniformBitsSpec.Define <T>();
if (typeof(T) == typeof(uint))
{
sampler = samplers.bits(rng, _spec.ToUInt32()) as ISampler <T>;
}
else if (typeof(T) == typeof(ulong))
{
sampler = samplers.bits(rng, _spec.ToUInt64()) as ISampler <T>;
}
else
{
throw unsupported <T>();
}
return(sampler);
}
public static IRngSampler <ulong> bits(MklRng src, UniformBitsSpec <ulong> spec, int?capacity = null) => new UniformBitsSampler <ulong>(src, capacity);
