[an error occurred while processing this directive]
|
за какое время выполняются приведенные ф-и (тактов скажем).
Нужно для осознанного выбора процессора (AD или Ti), а софта для Ti под рукой нету (с AD работал).
Math_eq3(1.0/6.0, 0.5, 50, -100, 3); // ret-1.9381661282967
EQtype1(0, 0.7777777, 9, 1.0/9.0, 20, 40); // ret-0.22616019410722
#define TYPE double
#define MAX_ERR 0.000001
TYPE Math_eq3(TYPE a, TYPE b, TYPE c, TYPE d, TYPE t1)
{
// int i;
TYPE t0, t, y0, y1, de;
t0=0.0;
y0=a*t0*t0*t0 + b*t0*t0 + c*t0 + d;
do{
y1=a*t1*t1*t1 + b*t1*t1 + c*t1 + d;
t=t1-(y1*(t0-t1))/(y0-y1);
de=fabs(t1-t);
t0=t1;
t1=t;
y0=y1;
}while(de>=MAX_ERR);
return t;
}
TYPE EQtype1(TYPE t0, TYPE t1, TYPE p, TYPE j, TYPE K1, TYPE K2)
{
TYPE y0, y1, t, de, g;
g=sqrt(K1*j+t0*t0); g=g*g*g;
y0=p*t0*t0*t0 + p*g + 2*K1*t0 - K2;
do{
g=sqrt(K1*j+t1*t1); g=g*g*g;
y1=p*t1*t1*t1 + p*g + 2*K1*t1 - K2;
t=t1-(y1*(t1-t0))/(y1-y0);
de=fabs(t1-t);
t0=t1;
t1=t;
y0=y1;
}while(de>=MAX_ERR);
return t;
}
E-mail: info@telesys.ru