/// <summary>
/// Test on a Cartesian grid, with an exact polynomial solution.
/// </summary>
public static ippControl TestCartesian1(int xRes = 32, double xStretch = 1.0, int yRes = 16, double yStretch = 1.01, int pDG = 2)
{
var RR = new ippControl();
RR.ProjectName = "ipPoison/cartesian";
RR.savetodb = false;
RR.FieldOptions.Add("T", new FieldOpts()
{
Degree = pDG, SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
RR.FieldOptions.Add("Tex", new FieldOpts()
{
Degree = pDG
});
RR.InitialValues_Evaluators.Add("RHS", X => 1.0);
RR.InitialValues_Evaluators.Add("Tex", X => (0.5 * X[0].Pow2() - 10 * X[0]));
RR.ExactSolution_provided = true;
RR.GridFunc = delegate() {
double[] xNodes = CreateNodes(xRes, xStretch, 0, 10);
double[] yNodes = CreateNodes(yRes, yStretch, -1, +1);
var grd = Grid2D.Cartesian2DGrid(xNodes, yNodes);
//Console.WriteLine("Achtung: Dreieck-gitter.");
//var grd = Grid2D.UnstructuredTriangleGrid(xNodes, yNodes);
return(grd);
};
RR.IsDirichlet = delegate(CommonParamsBnd inp) {
return(Math.Abs(inp.X[0] - 0.0) <= 1.0e-6);
};
RR.g_Diri = delegate(CommonParamsBnd inp) {
return(0.0);
};
RR.g_Neum = delegate(CommonParamsBnd inp) {
return(0.0);
};
//RR.solver_name = "direct";
RR.solver_name = null;
RR.GridPartType = BoSSS.Foundation.Grid.GridPartType.none;
return(RR);
}
/// <summary>
/// Test on a Cartesian grid, with an exact polynomial solution.
/// </summary>
public static ippControl TestCartesian3D(int xRes = 32, double xStretch = 1.0, int yRes = 16, double yStretch = 1.0, int zRes = 16, double zStretch = 1.0)
{
var R = new ippControl();
R.ProjectName = "ipPoison/cartesian";
R.savetodb = false;
R.FieldOptions.Add("T", new FieldOpts()
{
Degree = 6, SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.FieldOptions.Add("Tex", new FieldOpts()
{
Degree = 6
});
R.InitialValues_Evaluators.Add("RHS", X => 1.0);
R.InitialValues_Evaluators.Add("Tex", X => (0.5 * X[0].Pow2() - 10 * X[0]));
R.ExactSolution_provided = true;
R.GridFunc = delegate() {
double[] xNodes = CreateNodes(xRes, xStretch, 0, 10);
double[] yNodes = CreateNodes(yRes, yStretch, -1, +1);
double[] zNodes = CreateNodes(zRes, zStretch, -1, +1);
var grd = Grid3D.Cartesian3DGrid(xNodes, yNodes, zNodes);
return(grd);
};
R.IsDirichlet = delegate(CommonParamsBnd inp) {
return(Math.Abs(inp.X[0] - 0.0) <= 1.0e-6);
};
R.g_Diri = delegate(CommonParamsBnd inp) {
return(0.0);
};
R.g_Neum = delegate(CommonParamsBnd inp) {
return(0.0);
};
R.solver_name = null;
return(R);
}
/// <summary>
/// Test on a curved grid.
/// </summary>
public static ippControl TestCurved()
{
var R = new ippControl();
R.ProjectName = "ipPoison/curved";
R.savetodb = false;
R.FieldOptions.Add("T", new FieldOpts()
{
Degree = 2, SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.FieldOptions.Add("Tex", new FieldOpts()
{
Degree = 15
});
R.InitialValues_Evaluators.Add("RHS", X => 0.0);
R.InitialValues_Evaluators.Add("Tex", X => (Math.Log(X[0].Pow2() + X[1].Pow2()) / Math.Log(4.0)) + 1.0);
R.ExactSolution_provided = true;
R.GridFunc = delegate() {
var grd = Grid2D.CurvedSquareGrid(GenericBlas.Linspace(1, 2, 3), GenericBlas.Linspace(0, 1, 11), CellType.Square_9, true);
return(grd);
};
R.IsDirichlet = delegate(CommonParamsBnd inp) {
return(true);
};
R.g_Diri = delegate(CommonParamsBnd inp) {
double x = inp.X[0], y = inp.X[1];
return(Math.Sqrt(x * x + y * y));
};
R.solver_name = null;
return(R);
}
public ipFlux(double penalty_const, MultidimensionalArray cj, ippControl __ctrl)
: base(penalty_const, cj, "T") //
{
ctrl = __ctrl;
}
/// <summary>
/// Poisson Equation on a (-1,1)x(-1,1), Dirichlet everywhere
/// </summary>
public static ippControl Square(int xRes = 21, int yRes = 16, int deg = 5)
{
//Func<double[], double> exRhs = X => 2 * X[0] * X[0] + 2 * X[1] * X[1] - 4;
//Func<double[], double> exSol = X => (1.0 - X[0] * X[0]) * (1.0 - X[1] * X[1]);
//Func<double[], double> exSol = X => (1.0 - X[1]);
//Func<double[], double> exRhs = X => 0.0;
Func <double[], double> exSol = X => - Math.Cos(X[0] * Math.PI * 0.5) * Math.Cos(X[1] * Math.PI * 0.5);
Func <double[], double> exRhs = X => (Math.PI * Math.PI * 0.5 * Math.Cos(X[0] * Math.PI * 0.5) * Math.Cos(X[1] * Math.PI * 0.5)); // == - /\ exSol
var R = new ippControl();
R.ProjectName = "ipPoison/square";
R.savetodb = false;
//R.DbPath = "D:\\BoSSS-db";
R.FieldOptions.Add("T", new FieldOpts()
{
Degree = deg, SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.FieldOptions.Add("Tex", new FieldOpts()
{
Degree = 4
});
R.InitialValues_Evaluators.Add("RHS", exRhs);
R.InitialValues_Evaluators.Add("Tex", exSol);
R.ExactSolution_provided = true;
R.GridFunc = delegate() {
double[] xNodes = GenericBlas.Linspace(-1, 1, xRes);
double[] yNodes = GenericBlas.Linspace(-1, 1, yRes);
var grd = Grid2D.Cartesian2DGrid(xNodes, yNodes);
return(grd);
};
R.IsDirichlet = delegate(CommonParamsBnd inp) {
double x = inp.X[0];
double y = inp.X[1];
return(true);
//return !((Math.Abs(y - 1) < 1.0e-8) || (Math.Abs(x - 1) < 1.0e-8));
};
R.g_Diri = delegate(CommonParamsBnd inp) {
return(exSol(inp.X));
};
R.g_Neum = delegate(CommonParamsBnd inp) {
double x = inp.X[0];
double y = inp.X[1];
// y == 1
if (Math.Abs(y - 1) < 1.0e-8)
{
return(0.5 * Math.PI * Math.Cos(0.5 * x * Math.PI));
}
// x == 1
if (Math.Abs(x - 1) < 1.0e-8)
{
return(0.5 * Math.PI * Math.Cos(0.5 * y * Math.PI));
}
return(double.NaN);
};
R.solver_name = null;
R.NoOfSolverRuns = 1;
return(R);
}
/// <summary>
/// Test on a Cartesian grid, with a sinusodial solution.
/// </summary>
/// <param name="Res">
/// Grid resolution, 2 entries for 2D test, 3 entries for 3D test.
/// </param>
/// <param name="Stretch">
/// Grid stretching factors.
/// </param>
/// <param name="solver_name">
/// Name of solver to use.
/// </param>
public static ippControl TestCartesian2(int[] Res, double[] Stretch = null, string solver_name = "direct", int deg = 3)
{
if (Res.Length != 2 && Res.Length != 3)
{
throw new ArgumentOutOfRangeException();
}
if (Stretch == null)
{
Stretch = new double[Res.Length];
Stretch.SetAll(1.0);
}
else
{
if (Res.Length != Stretch.Length)
{
throw new ArgumentException();
}
}
var R = new ippControl();
R.ProjectName = "ipPoison/cartesian";
R.savetodb = false;
R.FieldOptions.Add("T", new FieldOpts()
{
Degree = deg, SaveToDB = FieldOpts.SaveToDBOpt.TRUE
});
R.FieldOptions.Add("Tex", new FieldOpts()
{
Degree = deg
});
R.InitialValues_Evaluators.Add("RHS", X => - Math.Sin(X[0]));
R.InitialValues_Evaluators.Add("Tex", X => Math.Sin(X[0]));
R.ExactSolution_provided = true;
R.GridFunc = delegate() {
if (Res.Length == 2)
{
double[] xNodes = CreateNodes(Res[0], Stretch[0], 0, 10);
double[] yNodes = CreateNodes(Res[1], Stretch[1], -1, +1);
var grd = Grid2D.Cartesian2DGrid(xNodes, yNodes);
return(grd);
}
else if (Res.Length == 3)
{
double[] xNodes = CreateNodes(Res[0], Stretch[0], 0, 10);
double[] yNodes = CreateNodes(Res[1], Stretch[1], -1, +1);
double[] zNodes = CreateNodes(Res[2], Stretch[2], -1, +1);
var grd = Grid3D.Cartesian3DGrid(xNodes, yNodes, zNodes);
return(grd);
}
else
{
throw new NotSupportedException();
}
};
R.IsDirichlet = delegate(CommonParamsBnd inp) {
return(Math.Abs(inp.X[0] - 0.0) <= 1.0e-6);
};
R.g_Diri = delegate(CommonParamsBnd inp) {
return(0.0);
};
R.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);
}
if (inp.X.Length > 2 && (Math.Abs(inp.X[2] - 1.0) < 1.0e-8 || Math.Abs(inp.X[2] + 1.0) < 1.0e-8))
{
return(0);
}
return(Math.Cos(10.0));
};
R.solver_name = solver_name;
return(R);
}