public static RpropBufferData BufferTGaInit()
{
var result = new RpropBufferData();
if (!VariablePointOp_RpropGamma.UseMean)
{
result.lowerBound = 0;
}
return(result);
}
public static RpropBufferData Buffer0Init()
{
var result = new RpropBufferData();
if (!UseMean)
{
result.lowerBound = 0;
}
return(result);
}
/// <include file='FactorDocs.xml' path='factor_docs/message_op_class[@name="VariablePointOp_RpropBeta"]/message_doc[@name="MarginalAverageConditional(Beta, Beta, RpropBufferData, Beta)"]/*'/>
public static Beta MarginalAverageConditional([IgnoreDependency] Beta use, [IgnoreDependency] Beta def, [RequiredArgument] RpropBufferData bufferBeta, Beta result)
{
result.Point = bufferBeta.nextPoint;
return(result);
}
public static RpropBufferData BufferBeta([NoInit] Beta use, Beta def, Beta to_marginal, RpropBufferData bufferBeta)
{
var currDist = use * def;
if (currDist.IsPointMass)
{
bufferBeta.nextPoint = currDist.Point;
return(bufferBeta);
}
// cannot use buffer.nextPoint as currPoint since Marginal could be initialized by user
double currPoint;
if (to_marginal.IsPointMass)
{
currPoint = to_marginal.Point;
}
else
{
currPoint = currDist.GetMean();
}
double currDeriv, currDeriv2;
currDist.GetDerivatives(currPoint, out currDeriv, out currDeriv2);
bufferBeta.SetNextPoint(currPoint, currDeriv);
return(bufferBeta);
}
/// <include file='FactorDocs.xml' path='factor_docs/message_op_class[@name="VariablePointOp_RpropGamma"]/message_doc[@name="MarginalAverageConditional(Gamma, Gamma, RpropBufferData, Gamma)"]/*'/>
public static Gamma MarginalAverageConditional([IgnoreDependency] Gamma use, [IgnoreDependency] Gamma def, [RequiredArgument] RpropBufferData buffer0, Gamma result)
{
if (UseMean)
{
result.Point = Math.Exp(buffer0.nextPoint);
}
else
{
result.Point = buffer0.nextPoint;
}
return(result);
}
public static RpropBufferData Buffer0([NoInit] Gamma use, Gamma def, Gamma to_marginal, RpropBufferData buffer0)
{
var currDist = use * def;
if (currDist.IsPointMass)
{
if (UseMean)
{
buffer0.nextPoint = Math.Log(currDist.Point);
}
else
{
buffer0.nextPoint = currDist.Point;
}
return(buffer0);
}
// cannot use buffer.nextPoint as currPoint since Marginal could be initialized by user
double currPoint;
if (to_marginal.IsPointMass)
{
currPoint = to_marginal.Point;
}
else
{
currPoint = currDist.GetMean();
}
double currDeriv, currDeriv2;
if (UseMean)
{
currDeriv = currDist.Shape - currDist.Rate * currPoint;
if (currPoint <= 0)
{
throw new ArgumentException($"currPoint ({currPoint}) <= 0");
}
buffer0.SetNextPoint(Math.Log(currPoint), currDeriv);
}
else
{
currDist.GetDerivatives(currPoint, out currDeriv, out currDeriv2);
buffer0.SetNextPoint(currPoint, currDeriv);
}
return(buffer0);
}
/// <include file='FactorDocs.xml' path='factor_docs/message_op_class[@name="VariablePointOp_RpropTruncatedGamma"]/message_doc[@name="MarginalAverageConditional(TruncatedGamma, TruncatedGamma, RpropBufferData, TruncatedGamma)"]/*'/>
public static TruncatedGamma MarginalAverageConditional([IgnoreDependency] TruncatedGamma use, [IgnoreDependency] TruncatedGamma def, [RequiredArgument] RpropBufferData bufferTGa, TruncatedGamma result)
{
if (VariablePointOp_RpropGamma.UseMean)
{
result.Point = Math.Exp(bufferTGa.nextPoint);
}
else
{
result.Point = bufferTGa.nextPoint;
}
return(result);
}
public static RpropBufferData BufferTGa([NoInit] TruncatedGamma use, TruncatedGamma def, TruncatedGamma to_marginal, RpropBufferData bufferTGa)
{
var currDist = use * def;
if (currDist.IsPointMass)
{
if (double.IsInfinity(currDist.Point) || double.IsNaN(currDist.Point))
{
throw new ArgumentOutOfRangeException();
}
if (VariablePointOp_RpropGamma.UseMean)
{
bufferTGa.nextPoint = Math.Log(currDist.Point);
}
else
{
bufferTGa.nextPoint = currDist.Point;
}
return(bufferTGa);
}
// cannot use buffer.nextPoint as currPoint since Marginal could be initialized by user
double currPoint;
if (to_marginal.IsPointMass)
{
currPoint = to_marginal.Point;
}
else
{
currPoint = currDist.GetMean();
}
double currDeriv, currDeriv2;
if (VariablePointOp_RpropGamma.UseMean)
{
bufferTGa.lowerBound = Math.Log(currDist.LowerBound);
bufferTGa.upperBound = Math.Log(currDist.UpperBound);
currDeriv = currDist.Gamma.Shape - currDist.Gamma.Rate * currPoint;
//Trace.WriteLine($"use deriv = {(use.Gamma.Shape-1) - use.Gamma.Rate*currPoint} def deriv = {def.Gamma.Shape - def.Gamma.Rate*currPoint} total deriv = {currDeriv}");
if (currPoint <= 0)
{
throw new ArgumentException($"currPoint ({currPoint}) <= 0");
}
bufferTGa.SetNextPoint(Math.Log(currPoint), currDeriv);
}
else
{
bufferTGa.lowerBound = currDist.LowerBound;
bufferTGa.upperBound = currDist.UpperBound;
currDist.Gamma.GetDerivatives(currPoint, out currDeriv, out currDeriv2);
bufferTGa.SetNextPoint(currPoint, currDeriv);
}
return(bufferTGa);
}
/// <include file='FactorDocs.xml' path='factor_docs/message_op_class[@name="VariablePointOp_RpropTruncatedGaussian"]/message_doc[@name="MarginalAverageConditional(TruncatedGaussian, TruncatedGaussian, RpropBufferData, TruncatedGaussian)"]/*'/>
public static TruncatedGaussian MarginalAverageConditional([IgnoreDependency] TruncatedGaussian use, [IgnoreDependency] TruncatedGaussian def, [RequiredArgument] RpropBufferData bufferTG, TruncatedGaussian result)
{
result.Point = bufferTG.nextPoint;
return(result);
}
public static RpropBufferData BufferTG([NoInit] TruncatedGaussian use, TruncatedGaussian def, TruncatedGaussian to_marginal, RpropBufferData bufferTG)
{
var currDist = use * def;
if (currDist.IsPointMass)
{
if (double.IsInfinity(currDist.Point))
{
throw new ArgumentOutOfRangeException(nameof(use), "infinite point mass");
}
bufferTG.nextPoint = currDist.Point;
return(bufferTG);
}
// cannot use buffer.nextPoint as currPoint since Marginal could be initialized by user
double currPoint;
if (to_marginal.IsPointMass)
{
currPoint = to_marginal.Point;
}
else
{
currPoint = currDist.GetMean();
}
// deriv of -0.5*prec*x^2+pm*x
// is -prec*x + pm
double currDeriv = currDist.Gaussian.MeanTimesPrecision - currDist.Gaussian.Precision * currPoint;
bufferTG.lowerBound = currDist.LowerBound;
bufferTG.upperBound = currDist.UpperBound;
bufferTG.SetNextPoint(currPoint, currDeriv);
return(bufferTG);
}
public static RpropBufferData Buffer([NoInit] Gaussian use, Gaussian def, Gaussian to_marginal, RpropBufferData buffer)
{
var currDist = use * def;
if (currDist.IsPointMass)
{
buffer.nextPoint = currDist.Point;
return(buffer);
}
// cannot use buffer.nextPoint as currPoint since Marginal could be initialized by user
double currPoint;
if (to_marginal.IsPointMass)
{
currPoint = to_marginal.Point;
}
else
{
currPoint = currDist.GetMean();
}
// deriv of -0.5*prec*x^2+pm*x
// is -prec*x + pm
double currDeriv = currDist.MeanTimesPrecision - currDist.Precision * currPoint;
buffer.SetNextPoint(currPoint, currDeriv);
return(buffer);
}
