//C++ TO C# CONVERTER WARNING: 'const' methods are not available in C#:
//ORIGINAL LINE: virtual double FuncDeriv(const Vector & x, const Vector & dir, double & deriv) const
public override double FuncDeriv(Vector x, Vector dir, ref double deriv)
{
// from 2d:
int j;
int k;
int lpi;
int gpi;
Vec <3> n, vgrad;
Point <3> pp1;
Vec2d g1 = new Vec2d();
Vec2d vdir = new Vec2d();
double badness;
double hbad;
double hderiv;
vgrad = 0;
badness = 0;
ld.meshthis.GetNormalVector(ld.surfi, ld.sp1, ld.gi1, n);
// pp1 = sp1;
// pp1.Add2 (x.Get(1), t1, x.Get(2), t2);
pp1 = ld.sp1 + x(0) * ld.t1.functorMethod + x(1) * ld.t2.functorMethod;
//C++ TO C# CONVERTER NOTE: This static local variable declaration (not allowed in C#) has been moved just prior to the method:
// static Array<Point2d> pts2d;
FuncDeriv_pts2d.SetSize(mesh.GetNP());
deriv = 0;
for (j = 1; j <= ld.locelements.Size(); j++)
{
lpi = ld.locrots.Get(j);
Element2d bel = mesh[ld.locelements.Get(j)];
gpi = bel.PNum(lpi);
for (k = 1; k <= bel.GetNP(); k++)
{
PointIndex pi = bel.PNum(k);
FuncDeriv_pts2d.Elem(pi) = new Point2d(ld.t1.functorMethod * (new mesh.Point(pi) - ld.sp1), ld.t2.functorMethod * (new mesh.Point(pi) - ld.sp1));
}
FuncDeriv_pts2d.Elem(gpi) = new Point2d(x(0), x(1));
vdir = new Vec2d(dir(0), dir(1));
hbad = bel.CalcJacobianBadnessDirDeriv(FuncDeriv_pts2d, lpi, vdir, ref hderiv);
deriv += hderiv;
badness += hbad;
}
return(badness);
}
//C++ TO C# CONVERTER WARNING: 'const' methods are not available in C#:
//ORIGINAL LINE: virtual double FuncGrad(const Vector & x, Vector & grad) const
public override double FuncGrad(Vector x, ref Vector grad)
{
// from 2d:
int lpi;
int gpi;
Vec <3> n, vgrad;
Point <3> pp1;
Vec2d g1 = new Vec2d();
Vec2d vdir = new Vec2d();
double badness;
double hbad;
double hderiv;
vgrad = 0;
badness = 0;
ld.meshthis.GetNormalVector(ld.surfi, ld.sp1, ld.gi1, n);
pp1 = ld.sp1 + x(0) * ld.t1.functorMethod + x(1) * ld.t2.functorMethod;
// meshthis -> ProjectPoint (surfi, pp1);
// meshthis -> GetNormalVector (surfi, pp1, n);
//C++ TO C# CONVERTER NOTE: This static local variable declaration (not allowed in C#) has been moved just prior to the method:
// static Array<Point2d> pts2d;
FuncGrad_pts2d.SetSize(mesh.GetNP());
grad = 0;
for (int j = 1; j <= ld.locelements.Size(); j++)
{
lpi = ld.locrots.Get(j);
Element2d bel = mesh[ld.locelements.Get(j)];
gpi = bel.PNum(lpi);
for (int k = 1; k <= bel.GetNP(); k++)
{
PointIndex pi = bel.PNum(k);
FuncGrad_pts2d.Elem(pi) = new Point2d(ld.t1.functorMethod * (new mesh.Point(pi) - ld.sp1), ld.t2.functorMethod * (new mesh.Point(pi) - ld.sp1));
}
FuncGrad_pts2d.Elem(gpi) = new Point2d(x(0), x(1));
for (int k = 1; k <= 2; k++)
{
if (k == 1)
{
vdir = new Vec2d(1, 0);
}
else
{
vdir = new Vec2d(0, 1);
}
hbad = bel.CalcJacobianBadnessDirDeriv(FuncGrad_pts2d, lpi, vdir, ref hderiv);
grad(k - 1) += hderiv;
if (k == 1)
{
badness += hbad;
}
}
}
/*
* vgrad.Add (-(vgrad * n), n);
*
* grad.Elem(1) = vgrad * t1;
* grad.Elem(2) = vgrad * t2;
*/
return(badness);
}
