public double ChooseAplpha_LipMethodForConditionalGradient(Matrix KSI, Matrix f, double alpha, double tau, double h)
{
double rho = ChooseAlpha_LipMethod();
int N = f.Length.n;
int M = f.Length.m;
double NORM = 0;
Matrix KSI2 = new Matrix(N, M);
Matrix f2 = new Matrix(N, M);
for (int i = 0; i < N; i++)
{
for (int j = 0; j < M; j++)
{
KSI2[i, j] = Math.Pow(Math.Abs(KSI[i, j]), 2);
}
}
NORM = Math.Sqrt(MyMath.IntegrateSurface(KSI2, tau, h));
double f_prime = 0, temp = 0;
for (int i = 0; i < N; i++)
{
for (int j = 0; j < M; j++)
{
//R = 10;
f_prime = -R * KSI[i, j] / NORM;
temp = f_prime - f[i, j];
f2[i, j] = temp * temp;
}
}
NORM = MyMath.IntegrateSurface(f2, tau, h);
double dJ = DifferentialEquation.dJ_f(KSI, f, h, tau);
//MessageBox.Show("dJ = " + dJ + " NORM = " + NORM);
return(min(1d, rho * dJ / (NORM)));
}
