public double Rdiff(IModelParam param)
{
CPHParam v = param as CPHParam;
double o = NMath.Abs(omega - v.omega) / v.omega;
double a = Blas.Rdiff(ndim, alpha, v.alpha);
double b = Blas.Rdiff(ndim, rate, v.rate);
return(NMath.Max(o, a, b));
}
public double Adiff(IModelParam param)
{
HErlangParam v = param as HErlangParam;
double o = NMath.Abs(omega - v.omega);
double a = Blas.Adiff(m, alpha, v.alpha);
double b = Blas.Adiff(m, rate, v.rate);
return(NMath.Max(o, a, b));
}
public void CopyTo(IModelParam other)
{
OriginalParam v = other as OriginalParam;
if (n == v.n)
{
Blas.Dcopy(n, param, v.param);
}
else
{
throw new InvalidCastException();
}
}
public void CopyTo(IModelParam p)
{
HErlangParam v = p as HErlangParam;
if (m == v.m)
{
v.omega = omega;
Blas.Dcopy(m, alpha, v.alpha);
Array.Copy(shape, v.shape, m);
Blas.Dcopy(m, rate, v.rate);
}
else
{
throw new InvalidCastException();
}
}
public void CopyTo(IModelParam p)
{
CPHParam v = p as CPHParam;
if (ndim == v.ndim)
{
v.omega = omega;
v.lambda = lambda;
Blas.Dcopy(ndim, alpha, v.alpha);
Blas.Dcopy(ndim, rate, v.rate);
Blas.Dcopy(ndim, scaledRate, v.scaledRate);
}
else
{
throw new InvalidCastException();
}
}
public double Rdiff(IModelParam other)
{
OriginalParam v = other as OriginalParam;
return(Blas.Rdiff(n, param, v.param));
}
private void Init()
{
prev_param = model.GetParam().Create();
}
