public void Axpy()
 {
   DoubleVector a = new DoubleVector(new double[4]{0,1,2,3});
   Double scal = 3;
   DoubleVector b = new DoubleVector(4);
   
   b.Axpy(scal,a);
   a.Scale(scal);
   
   Assert.AreEqual(a[0],b[0]);
   Assert.AreEqual(a[1],b[1]);
   Assert.AreEqual(a[2],b[2]);
   Assert.AreEqual(a[3],b[3]);
 }