public void OneDemExampleTest()
{
var s = new OneDemExample();
double RoTst = s.AllParticles.Sum(p =>
p.M * KernelF.W(s.AllParticles[50].X - p.X, s.h));
Assert.AreEqual(1.0, RoTst, 0.00001);
}
public void Nevermind()
{
var s = new OneDemExample();
var lstfail = new List <Tuple <double, double> >(s.Particles.Count);
var lstRo = new List <Tuple <double, double> >(s.Particles.Count);
s.h /= 2;
foreach (var particle in s.Particles)
{
var ro = s.AllParticles.Sum(p => p.M * KernelF.W(particle.X - p.X, s.h));
lstRo.Add(new Tuple <double, double>(particle.X, ro));
if (Math.Abs(particle.Ro - ro) > 0.0001)
{
lstfail.Add(new Tuple <double, double>(particle.X, ro));
}
}
Assert.AreEqual(1, 1, 0.1);
}