/// <summary>
/// default plotting
/// </summary>
protected override void PlotCurrentState(double phystime, TimestepNumber timestepNo, int superSampling = 0)
{
string caseStr = "";
DGField[] Fields = new DGField[0];
SinglePhaseField phiDist = ComputeDistanceField();
Fields = Fields.Cat(this.phi, this.mu, this.Velocity, this.gradPhi0, this.Curvature, this.DCurvature, phiDist);
BoSSS.Solution.Tecplot.Tecplot.PlotFields(Fields, "Phasefield-" + timestepNo + caseStr, phystime, superSampling);
}
//protected override IEnumerable<DGField> InstantiateParameterFields() {
// DGField[] ParameterFields = Velocity.ToArray().Cat(phi0);
// return ParameterFields;
//}
public override IEnumerable <DGField> InstantiateResidualFields()
{
DGField[] ResidualFields = new DGField[] { phi_Resi };
switch (this.Control.ModTyp)
{
case PhasefieldControl.ModelType.modelB:
ResidualFields = ResidualFields.Cat(mu_Resi);
break;
case PhasefieldControl.ModelType.modelA:
case PhasefieldControl.ModelType.modelC:
default:
break;
}
if (this.Control.CurvatureCorrectionType != PhasefieldControl.CurvatureCorrection.None)
{
ResidualFields.Cat(curvature_Resi);
}
return(ResidualFields);
}
protected override IEnumerable <DGField> InstantiateSolutionFields()
{
DGField[] SolutionFields = new DGField[] { phi };
switch (this.Control.ModTyp)
{
case PhasefieldControl.ModelType.modelB:
SolutionFields = SolutionFields.Cat(mu);
break;
case PhasefieldControl.ModelType.modelA:
case PhasefieldControl.ModelType.modelC:
default:
break;
}
if (this.Control.CurvatureCorrectionType != PhasefieldControl.CurvatureCorrection.None)
{
SolutionFields.Cat(Curvature);
}
return(SolutionFields);
}
