public void DerivativeTest()
{
double[] expected =
{
0.5, 0.497502, 0.490033, 0.477668, 0.46053, 0.438791,
0.412668, 0.382421, 0.348353, 0.310805, 0.270151
};
SinLinkFunction target = new SinLinkFunction();
for (int i = 0; i < 11; i++)
{
double x = i / 10.0;
double y = target.Inverse(x);
double d1 = target.Derivative(x);
double d2 = target.Derivative2(y);
Assert.AreEqual(expected[i], d1, 1e-6);
Assert.AreEqual(expected[i], d2, 1e-6);
Assert.IsFalse(Double.IsNaN(d1));
Assert.IsFalse(Double.IsNaN(d2));
}
}
public static ILinkFunction Get(LinkFunctionType link)
{
ILinkFunction result = null;
switch (link)
{
case LinkFunctionType.Absolute: result = new AbsoluteLinkFunction(); break;
case LinkFunctionType.Cauchit: result = new CauchitLinkFunction(); break;
case LinkFunctionType.Identity: result = new IdentityLinkFunction(); break;
case LinkFunctionType.Inverse: result = new InverseLinkFunction(); break;
case LinkFunctionType.InverseSquared: result = new InverseSquaredLinkFunction(); break;
case LinkFunctionType.Logit: result = new LogitLinkFunction(); break;
case LinkFunctionType.Log: result = new LogLinkFunction(); break;
case LinkFunctionType.LogLog: result = new LogLogLinkFunction(); break;
case LinkFunctionType.Probit: result = new ProbitLinkFunction(); break;
case LinkFunctionType.Sin: result = new SinLinkFunction(); break;
case LinkFunctionType.Threshold: result = new ThresholdLinkFunction(); break;
}
return(result);
}
public void SinLinkFunctionConstructorTest()
{
SinLinkFunction target = new SinLinkFunction();
Assert.AreEqual(-1, target.A);
Assert.AreEqual(2, target.B);
for (int i = 0; i < 11; i++)
{
double x = i / 10.0;
double y = (Math.Sin(x) + 1) / 2.0;
Assert.AreEqual(y, target.Function(x), 1e-10);
Assert.AreEqual(x, target.Inverse(y), 1e-10);
}
}
