namespace PKMA_
{
 [HandlerCategory("Kaufman's")]
 public class PKMA : IDouble2DoubleHandler
 {
 [HandlerParameter(true, "5", Min = "2", Max = "100", Step = "1")]
 public int Period { get; set; }

 public IList<double> Execute(IList<double> source)
 {
 var P = source;
 var CUM = new double[P.Count];
 var PKMA = new double[P.Count];
 var SMA_ = new double[P.Count];

 for(int i = 0; i < P.Count; i++)
 {
 for(int j = 0; j <=i; j++)
 CUM[i]+=1;
 }
 for(int i = 0; i < P.Count; i++)
 {
 if(i<Period) SMA_[i]=0;
 else
 {
 for(int j = i-Period+1; j <=i; j++)
 SMA_[i]+=P[j];
 }
 SMA_[i]=SMA_[i]/(Period+0.0);
 }

 for(int i = 0; i < P.Count; i++)
 { 
 if(i<Period) PKMA[i]=0;
 else
 {
 var Day = CUM[i]+1;
 var MV = 1.0/(Period+0.0); 
 if(Day<(Period+2)) PKMA[i] = P[i];
 else
 {
 if(Day==(Period+2)) PKMA[i] = SMA_[i];
 else PKMA[i] = PKMA[i-1]+(MV*(P[i]-PKMA[i-1]));
 }
 } 
 }
 return PKMA;
 }
 }
}


namespace PKAMAII
{ 
 [HandlerCategory("Kaufman's Adaptive Moving Average II")]
 public class PKAMAII : IDouble2DoubleHandler
 {
 [HandlerParameter(true, "10", Min = "1", Max = "32767", Step = "1")]
 public int Period { get; set; }

 [HandlerParameter(true, "0.666", Min = "0", Max = "1", Step = "0.001")]
 public double FEndF { get; set; }

 [HandlerParameter(true, "0.0645", Min = "0", Max = "1", Step = "0.0001")]
 public double SEndF { get; set; }

 public IList<double> Execute(IList<double> source)
 {
 var P = source;
 var Dnoise = new double[P.Count];
 var Noise = new double[P.Count];
 var PKAMAII = new double[P.Count];

 for(int i = 0; i < P.Count; i++)
 {
 if(i<1) Dnoise[i]=0;
 else Dnoise[i] = Math.Abs(P[i]-P[i-1]);
 }

 for(int i = 0; i < P.Count; i++)
 {
 if(i<Period) Noise[i]=0;
 else
 {
 for(int j = i-Period+1; j <= i; j++)
 Noise[i]+=Dnoise[j];
 }
 }

 for(int i = 0; i < P.Count; i++)
 {
 if(i<Period) PKAMAII[i]=0;
 else
 {
 var Signal = Math.Abs(P[i]-P[i-Period]);
 var EffRatio = Noise[i]>0?(Signal/Noise[i]):0;
 var FERatio = FEndF*EffRatio+SEndF*(1-EffRatio);
 var SmoothF = Math.Pow(FERatio,2);
 PKAMAII[i]=P[i]*SmoothF+(1-SmoothF)*PKAMAII[i-1];
 }
 }
 return PKAMAII;
 }
 }