// ==========================
// Check boundary conditions
// ==========================
// ========================
// Check level set movement
// ========================
//[Test] Deactivated, because failing & to much variations
public static void LineMovementTest(
[Values(LevelSetEvolution.FastMarching, LevelSetEvolution.ExtensionVelocity, LevelSetEvolution.ScalarConvection, LevelSetEvolution.Fourier)] LevelSetEvolution lsEvo,
[Values(LevelSetHandling.Coupled_Once, LevelSetHandling.LieSplitting, LevelSetHandling.Coupled_Iterative)] LevelSetHandling lsHandl,
[Values(TimeSteppingScheme.ImplicitEuler, TimeSteppingScheme.CrankNicolson, TimeSteppingScheme.BDF2)] TimeSteppingScheme tsScheme,
[Values(0.5)] double cLength)
{
var C = PhysicalBasedTestcases.ChannelFlow.CF_LevelSetMovementTest(1, cLength, lsEvo, lsHandl, tsScheme);
using (var solver = new XNSE_SolverMain())
{
solver.Init(C);
solver.RunSolverMode();
double[] BmQ_LL = solver.ComputeBenchmarkQuantities_LineInterface();
double err_thrsld = 1e-4;
// length of contact-line
double err = Math.Abs(2 - BmQ_LL[0]);
Assert.Less(err, err_thrsld, "error interface length");
Console.WriteLine("error in interface length = {0}", err);
// area of species
err = Math.Abs(cLength * 1 - BmQ_LL[1]);
Assert.Less(err, err_thrsld, "error in area");
Console.WriteLine("error in area = {0}", err);
}
}
