public override CDouble EvalDouble(CDouble a, CDouble b)
{
return(new CDouble(a.R + b.R, a.I + b.I));
}
public override CDouble EvalDouble(CDouble a, CDouble b)
{
return(new CDouble(a.R * Math.Cos(b.R), a.R * Math.Sin(b.R)));
}
public override CDouble EvalDouble(CDouble a)
{
return(new CDouble(Math.Sqrt(a.R)));
}
public virtual CDouble EvalDouble(CDouble a, CDouble b)
{
return EvalDouble(a);
}
public override CDouble EvalDouble(CDouble a)
{
return(new CDouble(Math.Atan2(a.R, a.I)));
}
public void Set(int index, CDouble d)
{
if (index >= 0 && index < real.Length) real[index] = d.R;
if (index >= 0 && index < imag.Length) imag[index] = d.I;
}
public virtual CDouble EvalDouble(CDouble a)
{
return a;
}
public virtual CDouble EvalDouble(CDouble a)
{
return(a);
}
public virtual CDouble EvalDouble(CDouble a, CDouble b)
{
return(EvalDouble(a));
}
public override CDouble EvalDouble(CDouble a, CDouble b)
{
return new CDouble(a.R - b.R, a.I - b.I);
}
public override CDouble EvalDouble(CDouble a)
{
return new CDouble(Math.Sqrt(a.R));
}
public override CDouble EvalDouble(CDouble a, CDouble b)
{
return new CDouble(a.R * Math.Cos(b.R),a.R * Math.Sin(b.R));
}
public override CDouble EvalDouble(CDouble a, CDouble b)
{
if (double.IsNaN(a.I) || double.IsNaN(b.I))
return new CDouble(a.R * b.R);
return new CDouble(a.R * b.R - a.I * b.I, a.R * b.I - a.I * b.R);
}
public override CDouble EvalDouble(CDouble a)
{
return new CDouble(Math.Atan2(a.R, a.I));
}
