/// <summary>
/// A circular interface in the 2D domain \f$ (3/2,3/2)^2 \f$.
/// </summary>
public static XdgPoisson3Control Circle(int Resolution = 16, int p = 1, string DBPath = null, LinearSolverCode solver = LinearSolverCode.classic_pardiso)
{
//BoSSS.Application.XdgPoisson3.HardCodedControl.Circle(Resolution: 8);
XdgPoisson3Control R = new XdgPoisson3Control();
R.ProjectName = "XdgPoisson3/Circle";
if (DBPath == null)
{
R.savetodb = false;
}
else
{
R.savetodb = false;
R.DbPath = DBPath;
}
R.SetDGdegree(p);
R.GridFunc = delegate() {
return(Grid2D.Cartesian2DGrid(GenericBlas.Linspace(-1.5, 1.5, Resolution + 1), GenericBlas.Linspace(-1.5, 1.5, Resolution + 1)));
};
double RADIUS = 0.8;
R.InitialValues_Evaluators.Add("Phi", X => - (X[0] - 0.0).Pow2() - (X[1] - 0.0).Pow2() + (RADIUS).Pow2());
R.ExcactSolSupported = false;
R.InitialValues_Evaluators.Add("uEx#A", X => 0.0);
R.InitialValues_Evaluators.Add("uEx#B", X => 0.0);
R.InitialValues_Evaluators.Add("u#A", X => 0.0);
R.InitialValues_Evaluators.Add("u#B", X => 0.0);
R.InitialValues_Evaluators.Add("rhs#A", X => + 1.0);
R.InitialValues_Evaluators.Add("rhs#B", X => + 1.0);
R.MU_A = -1.0;
R.MU_B = -1000.0;
R.xLaplaceBCs.g_Diri = (X => 0.0);
R.xLaplaceBCs.IsDirichlet = (inp => true);
R.LinearSolver.SolverCode = solver;
R.LinearSolver.SolverCode = LinearSolverCode.classic_pardiso;
R.TimesteppingMode = AppControl._TimesteppingMode.Steady;
R.AgglomerationThreshold = 0.1;
R.PrePreCond = MultigridOperator.Mode.SymPart_DiagBlockEquilib;
R.LinearSolver.NoOfMultigridLevels = 5;
return(R);
}
/// <summary>
/// A 45-degree interface in the 2D domain \f$ (-2,2)^2 \f$.
/// </summary>
public static XdgPoisson3Control Schraeg()
{
XdgPoisson3Control R = new XdgPoisson3Control();
R.ProjectName = "XdgPoisson3/Schräg";
R.savetodb = false;
R.FieldOptions.Add("Phi", new FieldOpts()
{
Degree = 3,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.FieldOptions.Add("u", new FieldOpts()
{
Degree = 3,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.GridFunc = delegate() {
return(Grid2D.Cartesian2DGrid(GenericBlas.Linspace(-2, 2, 8), GenericBlas.Linspace(-2, 2, 8)));
};
// set the level-set
R.InitialValues_Evaluators.Add("Phi", X => (X[0] + X[1]) + 0.075);
//R.InitialValues.Add("Phi", X => X[0] + +X[1]*0.001);
R.ExcactSolSupported = true;
R.InitialValues_Evaluators.Add("uEx#A", X => X[1]);
R.InitialValues_Evaluators.Add("uEx#B", X => X[1]);
R.InitialValues_Evaluators.Add("rhs#A", X => 0.0);
R.InitialValues_Evaluators.Add("rhs#B", X => 0.0);
R.MU_A = -1.0;
R.MU_B = -1.0;
R.xLaplaceBCs.g_Diri = delegate(CommonParamsBnd inp) {
double y = inp.X[1];
return(y);
};
R.xLaplaceBCs.IsDirichlet = (inp => true);
R.LinearSolver.SolverCode = LinearSolverConfig.Code.classic_pardiso;//R.solverName = "direct";
R.AgglomerationThreshold = 0.0;
R.PrePreCond = MultigridOperator.Mode.IdMass;
R.penalty_multiplyer = 1.1;
R.ViscosityMode = XLaplace_Interface.Mode.SIP;
return(R);
}
/// <summary>
/// A spherical interface in the 3D domain \f$ (-2, 2)^3 \f$.
/// </summary>
public static XdgPoisson3Control Ball3D()
{
XdgPoisson3Control R = new XdgPoisson3Control();
R.ProjectName = "XdgPoisson3/Ball3D";
R.savetodb = false;
R.FieldOptions.Add("Phi", new FieldOpts()
{
Degree = 1,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.FieldOptions.Add("u", new FieldOpts()
{
Degree = 2,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.GridFunc = delegate() {
return(Grid3D.Cartesian3DGrid(GenericBlas.Linspace(-2, 2, 6), GenericBlas.Linspace(-2, 2, 6), GenericBlas.Linspace(-2, 2, 6)));
};
// set the level-set
R.InitialValues_Evaluators.Add("Phi", X => (X[0].Pow2() + X[1].Pow2() + X[2].Pow2()) - 1.0.Pow2());
//R.InitialValues.Add("Phi", X => X[0] + 0.1);
R.ExcactSolSupported = false;
//R.InitialValues.Add("uEx#A", X => X[1]);
//R.InitialValues.Add("uEx#B", X => X[1]);
R.InitialValues_Evaluators.Add("rhs#A", X => 1.0);
R.InitialValues_Evaluators.Add("rhs#B", X => 1.0);
R.MU_A = -1.0;
R.MU_B = -1.0;
R.xLaplaceBCs.g_Diri = ((CommonParamsBnd inp) => 0.0);
R.xLaplaceBCs.IsDirichlet = (inp => true);
R.LinearSolver.SolverCode = LinearSolverConfig.Code.classic_pardiso;//R.solverName = "direct";
R.AgglomerationThreshold = 0.1;
R.PrePreCond = MultigridOperator.Mode.DiagBlockEquilib;
R.penalty_multiplyer = 1.1;
R.ViscosityMode = XLaplace_Interface.Mode.SIP;
return(R);
}
/// <summary>
/// A circular interface in the 2D domain \f$ (3/2,3/2)^2 \f$.
/// </summary>
public static XdgPoisson3Control Circle(int Resolution = 16, int p = 1, string DBPath = null, LinearSolverConfig.Code solver = LinearSolverConfig.Code.classic_pardiso)
{
XdgPoisson3Control R = new XdgPoisson3Control();
R.ProjectName = "XdgPoisson3/Circle";
if (DBPath == null)
{
R.savetodb = false;
}
else
{
R.savetodb = false;
R.DbPath = DBPath;
}
R.FieldOptions.Add("Phi", new FieldOpts()
{
Degree = p,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.FieldOptions.Add("u", new FieldOpts()
{
Degree = p,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.GridFunc = delegate() {
return(Grid2D.Cartesian2DGrid(GenericBlas.Linspace(-1.5, 1.5, Resolution), GenericBlas.Linspace(-1.5, 1.5, Resolution)));
};
double RADIUS = 0.8;
R.InitialValues_Evaluators.Add("Phi", X => - X[0].Pow2() - X[1].Pow2() + (RADIUS).Pow2());
R.ExcactSolSupported = false;
R.InitialValues_Evaluators.Add("uEx#A", X => 0.0);
R.InitialValues_Evaluators.Add("uEx#B", X => 0.0);
R.InitialValues_Evaluators.Add("u#A", X => 0.0);
R.InitialValues_Evaluators.Add("u#B", X => 0.0);
R.InitialValues_Evaluators.Add("rhs#A", X => + 1.0);
R.InitialValues_Evaluators.Add("rhs#B", X => + 1.0);
R.MU_A = -1.0;
R.MU_B = -1000.0;
R.xLaplaceBCs.g_Diri = (X => 0.0);
R.xLaplaceBCs.IsDirichlet = (inp => true);
R.LinearSolver.SolverCode = solver;
R.LinearSolver.SolverCode = LinearSolverConfig.Code.exp_softpcg_mg;
R.dtMax = 0.1;
R.dtMin = 0.1;
R.NoOfTimesteps = 10;
R.Endtime = 100000.0;
R.AgglomerationThreshold = 0.1;
R.PrePreCond = MultigridOperator.Mode.SymPart_DiagBlockEquilib_DropIndefinite;
R.LinearSolver.NoOfMultigridLevels = 5;
return(R);
}
/// <summary>
/// A circular interface within the 2D domain \f$ (-1,1)^2 \f$, with Dirichlet boundary conditions at \f$ x = -1 \f$ and Neumann boundary conditions elsewhere.
/// </summary>
public static XdgPoisson3Control CircleNeum(int Resolution = 13)
{
XdgPoisson3Control C = new XdgPoisson3Control();
C.DbPath = @"D:\\BoSSS-db";
C.savetodb = false;
C.ProjectName = "XDGPoisson/Circle";
C.ProjectDescription = "XDG Poisson Circle";
//Diffusion Parameter
C.MU_A = -1.0;
C.MU_B = -1.0;
//Funktion on the RHS of the Poisson-Equation
C.InitialValues_Evaluators.Add("rhs#A", X => Math.Sin((X[0] + 1.0) * 3));
C.InitialValues_Evaluators.Add("rhs#B", X => Math.Sin((X[0] + 1.0) * 3));
// Problem Definition
C.FieldOptions.Add("Phi", new FieldOpts()
{
Degree = 3,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
double RADIUS = 0.8;
C.InitialValues_Evaluators.Add("Phi", X => X[0].Pow2() + X[1].Pow2() - (RADIUS).Pow2());
C.FieldOptions.Add("u", new FieldOpts()
{
Degree = 3,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
//Definition of Grid - Only Dirichlet Boundaries
C.GridFunc = delegate() {
return(Grid2D.Cartesian2DGrid(GenericBlas.Linspace(-1, 1, Resolution), GenericBlas.Linspace(-1, 1, Resolution)));
};
C.xLaplaceBCs.g_Diri = (inp => 0.0);
C.xLaplaceBCs.g_Neum = delegate(CommonParamsBnd inp) {
if (Math.Abs(inp.X[1] - 1.0) < 1.0e-8 || Math.Abs(inp.X[1] + 1.0) < 1.0e-8)
{
return(0);
}
return(Math.Cos(2.0 * 3.0) * (1.0 / 3.0));
};
C.xLaplaceBCs.IsDirichlet = (inp => Math.Abs(inp.X[0] + 1.0) <= 1.0e-8);
/// Exact Solution
C.ExcactSolSupported = true;
C.FieldOptions.Add("uEx", new FieldOpts()
{
Degree = 3,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
C.InitialValues_Evaluators.Add("uEx#A", X => Math.Sin((X[0] + 1.0) * 3) * (1.0 / 9.0));
C.InitialValues_Evaluators.Add("uEx#B", X => Math.Sin((X[0] + 1.0) * 3) * (1.0 / 9.0));
/// Solver Parameters
//C.solverName = "pcg+mg+schwarz";
C.LinearSolver.SolverCode = LinearSolverConfig.Code.classic_pardiso;//R.solverName = "direct";
/// Discretization Parameters
C.ViscosityMode = XLaplace_Interface.Mode.SIP;
return(C);
}
/// <summary>
/// This is a testrun similar to the one in \public\doc\handbook\apdx-NodeSolverPerformance\XDGPoisson\Part1-Calculations.bws
/// phase A: 3D sphere, pahse B: [-1,1]^3\sphere
/// </summary>
/// <param name="myDB"></param>
/// <returns></returns>
public static XdgPoisson3Control TestOrTreat(int solver = 1, int blocksize = 10000, string myDB = null)
{
XdgPoisson3Control C = new XdgPoisson3Control();
switch (solver)
{
case 1:
C.LinearSolver.SolverCode = LinearSolverConfig.Code.exp_softpcg_mg;
break;
case 2:
C.LinearSolver.SolverCode = LinearSolverConfig.Code.exp_softpcg_schwarz_directcoarse;
break;
default:
throw new NotImplementedException("guess again");
}
C.savetodb = false;
//C.DbPath = @"E:\\XdgPerformance";
int Res = 8;
C.GridFunc = delegate() {
double[] xNodes = GenericBlas.Linspace(-1, +1, Res + 1);
double[] yNodes = GenericBlas.Linspace(-1, +1, Res + 1);
double[] zNodes = GenericBlas.Linspace(-1, +1, Res + 1);
int J = (xNodes.Length - 1) * (yNodes.Length - 1) * (zNodes.Length - 1);
var grid = Grid3D.Cartesian3DGrid(xNodes, yNodes, zNodes);
grid.Name = "thisisatestgrid";
grid.EdgeTagNames.Add(1, "Dirichlet");
grid.DefineEdgeTags(delegate(double[] X) {
return(1);
});
return(grid);
};
//these are parameters for batchprocessing. They are here for testing ...
//C.PerformanceModeON = true;
//C.SuppressExceptionPrompt = true;
C.LinearSolver.TargetBlockSize = blocksize;
C.SetDGdegree(2);
//C.LinearSolver.SolverCode = LinearSolverConfig.Code.exp_softpcg_mg;
//C.LinearSolver.SolverCode = LinearSolverConfig.Code.exp_softpcg_schwarz_directcoarse;
//C.LinearSolver.SolverCode = LinearSolverConfig.Code.exp_softpcg_jacobi_mg;
C.LinearSolver.NoOfMultigridLevels = 10;
C.LinearSolver.ConvergenceCriterion = 1e-6;
C.ExcactSolSupported = false;
double radius = 0.7;
C.InitialValues_Evaluators.Add("Phi", X => X[0].Pow2() + X[1].Pow2() + X[2].Pow2() - radius.Pow2());
C.MU_A = -1;
C.MU_B = -1000;
C.InitialValues_Evaluators.Add("rhs#A", X => 1.0);
C.InitialValues_Evaluators.Add("rhs#B", X => 1.0);
C.InitialValues_Evaluators.Add("u#A", X => 0);
C.InitialValues_Evaluators.Add("u#B", X => 0);
C.CutCellQuadratureType = XQuadFactoryHelper.MomentFittingVariants.Classic;
C.SetDefaultDiriBndCnd = true;
//C.xLaplaceBCs.g_Diri = ((CommonParamsBnd inp) => 0.0);
//C.xLaplaceBCs.IsDirichlet = (inp => true);
C.ViscosityMode = XLaplace_Interface.Mode.SIP;
C.AgglomerationThreshold = 0.1;
return(C);
}
/// <summary>
/// a piecewise parabolic solution.
/// </summary>
public static XdgPoisson3Control PiecewiseParabola(double delta = 0.0, double muA = -20, double muB = -1)
{
XdgPoisson3Control R = new XdgPoisson3Control();
if (delta < -1.5 || delta > 1.5)
{
throw new ArgumentOutOfRangeException();
}
R.ProjectName = "XdgPoisson3/PiecewiseParabola";
R.savetodb = false;
R.FieldOptions.Add("Phi", new FieldOpts()
{
Degree = 3,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.FieldOptions.Add("u", new FieldOpts()
{
Degree = 3,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.GridFunc = delegate() {
double[] xNodes = new double[] { -6, -5, -4, -3, -2, -1, 1, 2, 3, 4, 5, 6 };
double[] yNodes = GenericBlas.Linspace(-3, 3, 7);
//double[] xNodes = new double[] { -6, -4, -2, -1, 1, 6 };
//double[] yNodes = GenericBlas.Linspace(-3, 3, 7);
return(Grid2D.Cartesian2DGrid(xNodes, yNodes, periodicY: true));
};
R.MU_A = muA;
R.MU_B = muB;
double A0 = 36 / (muB + muA);
double A1 = (3 * (-muB + muA)) / (muA * (muB + muA));
double A2 = -1 / (2 * muA);
double B0 = 36 / (muB + muA);
double B1 = (3 * (-muB + muA)) / ((muB + muA) * muB);
double B2 = -1 / (2 * muB);
Func <double[], double> uEx_A = (X => - (A2 * (X[0] - delta).Pow2() + A1 * (X[0] - delta) + A0));
Func <double[], double> uEx_B = (X => - (B2 * (X[0] - delta).Pow2() + B1 * (X[0] - delta) + B0));
// set the level-set
R.InitialValues_Evaluators.Add("Phi", X => (X[0] - delta));
R.ExcactSolSupported = true;
R.InitialValues_Evaluators.Add("uEx#A", uEx_A);
R.InitialValues_Evaluators.Add("uEx#B", uEx_B);
R.InitialValues_Evaluators.Add("rhs#A", X => + 1.0);
R.InitialValues_Evaluators.Add("rhs#B", X => + 1.0);
double Diri_left = uEx_A(new double[] { -6, 0 });
double Diri_rigt = uEx_B(new double[] { +6, 0 });
R.xLaplaceBCs.g_Diri = (inp => ((inp.X[0] < 0) ? Diri_left : Diri_rigt));
R.xLaplaceBCs.g_Neum = (inp => 0.0);
R.xLaplaceBCs.IsDirichlet = delegate(CommonParamsBnd inp) {
double x = inp.X[0];
if (Math.Abs(x + 6) < 1.0e-6 || Math.Abs(x - 6) < 1.0e-6)
{
return(true);
}
else
{
return(false);
}
};
R.LinearSolver.SolverCode = LinearSolverConfig.Code.classic_pardiso;//R.solverName = "direct";
R.AgglomerationThreshold = 0.0;
R.PrePreCond = MultigridOperator.Mode.SymPart_DiagBlockEquilib;
return(R);
}
/// <summary>
/// A piecewise linear solution.
/// </summary>
public static XdgPoisson3Control MultigridVsAggregation(double delta = 0.1)
{
XdgPoisson3Control R = new XdgPoisson3Control();
R.ProjectName = "XdgPoisson3/PiecewiseLinear";
R.savetodb = false;
//R.DbPath = "C:\\BoSSS-db";
R.FieldOptions.Add("Phi", new FieldOpts()
{
Degree = 1,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.FieldOptions.Add("u", new FieldOpts()
{
Degree = 1,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.FieldOptions.Add("uEx", new FieldOpts()
{
Degree = 1,
SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
//var cutout = new double[2, 2] { { 0, 0 }, { 3, 3 } };
//var tt = new BoundingBox[] { new BoundingBox(2) };
R.GridFunc = delegate() {
double[] xNodes = GenericBlas.Linspace(-3, 3, 10);
double[] yNodes = GenericBlas.Linspace(-3, 3, 10);
return(Grid2D.Cartesian2DGrid(xNodes, yNodes, CutOuts: new BoundingBox[] { new BoundingBox(new double[2, 2] {
{ 0, 0 }, { 3, 3 }
}) }));
};
R.LinearSolver.NoOfMultigridLevels = 3;
const double kA = -1;
const double kB = -2;
// set the level-set
R.InitialValues_Evaluators.Add("Phi", X => (X[0] - delta));
R.ExcactSolSupported = true;
R.InitialValues_Evaluators.Add("uEx#A", X => kA * (X[0] - delta));
R.InitialValues_Evaluators.Add("uEx#B", X => kB * (X[0] - delta));
R.InitialValues_Evaluators.Add("rhs#A", X => 0.0);
R.InitialValues_Evaluators.Add("rhs#B", X => 0.0);
R.ExcactSolSupported = true;
R.MU_A = -1.0 / kA;
R.MU_B = -1.0 / kB;
R.xLaplaceBCs.g_Diri = delegate(CommonParamsBnd inp) {
double x = inp.X[0];
if (x < 0)
{
return((x - delta) * kA);
}
else
{
return((x - delta) * kB);
}
};
R.xLaplaceBCs.g_Neum = (inp => 0.0);
R.xLaplaceBCs.IsDirichlet = delegate(CommonParamsBnd inp) {
double x = inp.X[0];
if (Math.Abs(x + 3) < 1.0e-6 || Math.Abs(x - 3) < 1.0e-6)
{
return(true);
}
else
{
return(false);
}
};
R.LinearSolver.SolverCode = LinearSolverConfig.Code.classic_pardiso;//R.solverName = "direct";
R.AgglomerationThreshold = 0.2;
return(R);
}
/// <summary>
/// A piecewise linear solution.
/// </summary>
public static XdgPoisson3Control PiecewiseLinear(double delta = 0.5)
{
XdgPoisson3Control R = new XdgPoisson3Control();
R.ProjectName = "XdgPoisson3/PiecewiseLinear";
R.savetodb = false;
//R.DbPath = "C:\\BoSSS-db";
R.SetDGdegree(1);
R.GridFunc = delegate() {
double[] xNodes = new double[] { -6, -3, -2, -1, 1, 2, 3, 6 };
double[] yNodes = GenericBlas.Linspace(-3, 3, 2);
return(Grid2D.Cartesian2DGrid(xNodes, yNodes));
};
const double kA = -1;
const double kB = -1;
// set the level-set
R.InitialValues_Evaluators.Add("Phi", X => (X[0] - delta));
R.ExcactSolSupported = true;
R.InitialValues_Evaluators.Add("uEx#A", X => kA * (X[0] - delta));
R.InitialValues_Evaluators.Add("uEx#B", X => kB * (X[0] - delta));
R.InitialValues_Evaluators.Add("rhs#A", X => 0.0);
R.InitialValues_Evaluators.Add("rhs#B", X => 0.0);
R.MU_A = -1.0 / kA;
R.MU_B = -1.0 / kB;
R.xLaplaceBCs.g_Diri = delegate(CommonParamsBnd inp) {
double x = inp.X[0];
if (x < 0)
{
return((x - delta) * kA);
}
else
{
return((x - delta) * kB);
}
};
R.xLaplaceBCs.g_Neum = (inp => 0.0);
R.xLaplaceBCs.IsDirichlet = delegate(CommonParamsBnd inp) {
double x = inp.X[0];
if (Math.Abs(x + 6) < 1.0e-6 || Math.Abs(x - 6) < 1.0e-6)
{
return(true);
}
else
{
return(false);
}
};
R.LinearSolver.SolverCode = LinearSolverConfig.Code.exp_softpcg_schwarz;//R.solverName = "pcg+schwarz";
R.LinearSolver.NoOfMultigridLevels = 2;
R.AgglomerationThreshold = 0.0;
return(R);
}
/// <summary>
/// This is a testrun similar to the one in \public\doc\handbook\apdx-NodeSolverPerformance\XDGPoisson\Part1-Calculations.bws
/// phase A: 3D sphere, pahse B: [-1,1]^3\sphere
/// </summary>
/// <param name="myDB"></param>
/// <returns></returns>
public static XdgPoisson3Control TestOrTreat(int solver = 1, int blocksize = 10000, string myDB = null)
{
XdgPoisson3Control C = new XdgPoisson3Control();
switch (solver)
{
case 0:
C.LinearSolver.SolverCode = LinearSolverCode.classic_pardiso;
break;
case 1:
C.LinearSolver.SolverCode = LinearSolverCode.exp_Kcycle_schwarz;
break;
case 3:
C.LinearSolver.SolverCode = LinearSolverCode.exp_gmres_levelpmg;
break;
case 4:
C.LinearSolver.SolverCode = LinearSolverCode.exp_OrthoS_pMG;
break;
default:
throw new NotImplementedException("guess again");
}
C.savetodb = false;
//C.DbPath = @"D:\trash_db";
int Res = 8;
C.GridFunc = delegate() {
double[] xNodes = GenericBlas.Linspace(-1, +1, Res + 1);
double[] yNodes = GenericBlas.Linspace(-1, +1, Res + 1);
double[] zNodes = GenericBlas.Linspace(-1, +1, Res + 1);
int J = (xNodes.Length - 1) * (yNodes.Length - 1) * (zNodes.Length - 1);
var grid = Grid3D.Cartesian3DGrid(xNodes, yNodes, zNodes);
grid.Name = "thisisatestgrid";
grid.EdgeTagNames.Add(1, "Dirichlet");
grid.DefineEdgeTags(delegate(double[] X) {
return(1);
});
return(grid);
};
C.SessionName = String.Format("XDGPoison_solver{0}_blsz{1}_Xdg2lowB", solver, blocksize);
C.ProjectName = "PoisonTest";
C.GridPartType = GridPartType.Hilbert;
C.LinearSolver.TargetBlockSize = blocksize / 2;
C.SetDGdegree(2);
C.LinearSolver.NoOfMultigridLevels = 2;
C.LinearSolver.ConvergenceCriterion = 1e-8;
C.LinearSolver.MaxSolverIterations = 1000;
C.LinearSolver.MaxKrylovDim = 50;
//C.LinearSolver.TargetBlockSize = 79;
C.ExcactSolSupported = false;
double radius = 0.7;
C.InitialValues_Evaluators.Add("Phi", X => X[0].Pow2() + X[1].Pow2() + X[2].Pow2() - radius.Pow2());
C.MU_A = -1;
C.MU_B = -1000;
C.InitialValues_Evaluators.Add("rhs#A", X => 1.0);
C.InitialValues_Evaluators.Add("rhs#B", X => 1.0);
C.InitialValues_Evaluators.Add("u#A", X => 0);
C.InitialValues_Evaluators.Add("u#B", X => 0);
//C.CutCellQuadratureType = XQuadFactoryHelper.MomentFittingVariants.Classic;
C.CutCellQuadratureType = XQuadFactoryHelper.MomentFittingVariants.Saye;
C.SetDefaultDiriBndCnd = true; //which means ...
//C.xLaplaceBCs.g_Diri = ((CommonParamsBnd inp) => 0.0);
//C.xLaplaceBCs.IsDirichlet = (inp => true);
// ... but stuff is not serializable, therefore this workaround.
C.ViscosityMode = XLaplace_Interface.Mode.SIP;
C.AgglomerationThreshold = 0.1;
return(C);
}