public List <double> GetActualArgs(Lens len, double h1)
{
var nextRefraction = len.Refractiond;
var sinI = (this.L - len.Radius) / len.Radius * Math.Sin(this.U);
if (this.L >= Utility.Infinity)
{
sinI = h1 / len.Radius;
}
var sinINext = this.NowRefraction / nextRefraction * sinI;
var uNext = U + Math.Asin(sinI) - Math.Asin(sinINext);
var lNext = len.Radius + len.Radius * sinINext / Math.Sin(uNext);
return(new List <double> {
Math.Asin(sinI), Math.Asin(sinINext), uNext, lNext
});
}
public Light Paraxial(Lens len, double a, char ty)
{
double nextRefraction = 1;
switch (ty)
{
case 'd':
{
nextRefraction = len.Refractiond;
break;
}
case 'C':
{
nextRefraction = len.RefractionC;
break;
}
case 'F':
{
nextRefraction = len.RefractionF;
break;
}
}
double i = (this.L - len.Radius) / len.Radius * this.U;
// MessageBox.Show(L.ToString() +" "+len.Radius+" "+U+" "+i);
if (this.L >= Utility.Infinity)
{
i = a / len.Radius;
}
// ReSharper disable once InconsistentNaming
double i_ = this.NowRefraction / nextRefraction * i;
// ReSharper disable once InconsistentNaming
double u_ = this.U + i - i_;
// ReSharper disable once InconsistentNaming
double l_ = len.Radius + len.Radius * i_ / u_;
return(new Light(nextRefraction, l_ - len.Thickness, u_));
}
public Light Actual(Lens len, double h1, char ty)
{
double nextRefraction = 1;
switch (ty)
{
case 'd':
{
nextRefraction = len.Refractiond;
break;
}
case 'C':
{
nextRefraction = len.RefractionC;
break;
}
case 'F':
{
nextRefraction = len.RefractionF;
break;
}
}
var sinI = (this.L - len.Radius) / len.Radius * Math.Sin(this.U);
if (L >= Utility.Infinity)
{
sinI = h1 / len.Radius;
}
var sinINext = this.NowRefraction / nextRefraction * sinI;
var uNext = U + Math.Asin(sinI) - Math.Asin(sinINext);
var lNext = len.Radius + len.Radius * sinINext / Math.Sin(uNext);
return(new Light(nextRefraction, lNext - len.Thickness, uNext));
}