/// <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> UniformSampler <T>(this MklRng rng, UniformSpec <T>?spec = null)
where T : unmanaged
{
var _spec = spec ?? UniformSpec.Define <T>(PrimalInfo.minval <T>(), PrimalInfo.maxval <T>());
var sampler = default(ISampler <T>);
if (typeof(T) == typeof(int))
{
sampler = samplers.uniform(rng, _spec.ToInt32()) as ISampler <T>;
}
else if (typeof(T) == typeof(float))
{
sampler = samplers.uniform(rng, _spec.ToFloat32()) as ISampler <T>;
}
else if (typeof(T) == typeof(double))
{
sampler = samplers.uniform(rng, _spec.ToFloat64()) as ISampler <T>;
}
else
{
throw unsupported <T>();
}
return(sampler);
}
public static IRngSampler <T> UniformSampler <T>(this MklRng rng, UniformSpec <T>?spec = null)
where T : unmanaged
{
var _spec = spec ?? Distributions.uniform <T>(Numeric.minval <T>(), Numeric.maxval <T>());
var sampler = default(IRngSampler <T>);
if (typeof(T) == typeof(int))
{
sampler = samplers.uniform(rng, _spec.ToInt32()) as IRngSampler <T>;
}
else if (typeof(T) == typeof(float))
{
sampler = samplers.uniform(rng, _spec.ToFloat32()) as IRngSampler <T>;
}
else if (typeof(T) == typeof(double))
{
sampler = samplers.uniform(rng, _spec.ToFloat64()) as IRngSampler <T>;
}
else
{
throw no <T>();
}
return(sampler);
}
public UniformSampler(MklRng src, UniformSpec <T> distspec, int?buferLen = null)
: base(src, distspec, buferLen)
{
}
