История изменений
Исправление cvs-255, (текущая версия) :
Это какие именно, где про них почитать? И чем они лучше расчета по 5 точкам?
Если что, вот мой код расчета:
tensor2 dGi(Lvector p, int n)
{
int L = p.dim();
tensor2 res(L);
double coef[5] = {1/(12*H), -2/(3*H), 0, 2/(3*H), -1/(12*H)};
Lvector h(L);
int i;
for (i = 0; i < L; i++)
{
h[i] = (i==n)*H;
}
for (i = 0; i < 5; i++)
{
if (i == 2)
continue;
Lvector dp = h * (i-2);
tensor2 metr = Metric(p + dp);
res += metr*(coef[i]);
}
return res;
}
Исправление cvs-255, :
Это какие именно, где про них почитать? И чем они лучше расчета по 5 точкам?
Если что, вот мой код расчета:
tensor2 dGi(Lvector p, int n)
{
int L = p.dim();
tensor2 res(L);
double coef[5] = {1/(12*H), -2/(3*H), 0, 2/(3*H), -1/(12*H)};
Lvector h(L);
int i;
for (i = 0; i < L; i++)
{
h[i] = (i==n)*H;
}
for (i = 0; i < 5; i++)
{
Lvector dp = h * (i-2);
tensor2 metr = Metric(p + dp);
res += metr*(coef[i]);
}
return res;
}
Исправление cvs-255, :
Это какие именно, где про них почитать? И чем они лучше расчета по 5 точкам?
Исходная версия cvs-255, :
Это какие именно, где про них почитать?