🚀 go-pugleaf

W wolnej chwili zmajstrowałem sobie prostą "biblioteczkę" z podstawowymi funkcjami matematycznymi. Posłużyłem się przy tym "gotowcami" zamieszczonymi na stronie: https://gcc.gnu.org/onlinedocs/gcc/Extended-Asm.html

Sqrt liczy się tylko symbolicznie szybciej niż sqrt z biblioteki math.h, pozostałe funkcje są szybsze od 2x wzwyż.

Za Chiny nie mogę rozkminić tych literek między dwukropkami, ale wprowadziłem pewne korekty, per analogiam - i (chyba) działa... ;-)

W przypadku "jednoargumentowych" wyrażeń sugerowana jest składnia:

asm ("foo" : "=&t" (a) : "f" (b));

Przy okazji testów szybkości zauważyłem jednak, że jeśli w ten sposób zakoduję funkcje Sin() i Cos(), to w przypadku wywołania w kodzie programu:

sin = Sin(a);
cos = Cos(a);

poprawnie liczy się kosinus, ale w przypadku sinusa wychodzą jakieś głupoty.
Wywołanie typu:

tan = Sin(a) / Cos(a)

w ogóle się wywala, a przecież mam tego używać...

Poprzez analogię do przykładu z fsincos i fyl2xp1 do zmieniłem "f" na "0" i wywaliłem ampersanda sprzed t:

asm ("foo" : "=t" (a) : "0" (b));


No i teraz wszystko gra!

A dlaczego? Nie mam pojęcia... ;-)





//------------------------------------------------------------------------------
inline double Sqrt(double x)
{
double s;

   asm ("fsqrt" : "=t" (s) : "0" (x));

   return s;
}
//------------------------------------------------------------------------------
inline double Sin(double a)
{
double s;

   asm ("fsin" : "=t" (s) : "0" (a));

   return s;
}
//------------------------------------------------------------------------------
inline double Cos(double a)
{
double c;

   asm ("fcos" : "=t" (c) : "0" (a));

   return c;
}
//------------------------------------------------------------------------------
inline void SinCos(double a, double *sin, double *cos)
{
double s, c;

   asm ("fsincos" : "=t" (c), "=u" (s) : "0" (a));

   *sin = s;
   *cos = c;
}
//------------------------------------------------------------------------------
inline double Tan(double a)
{
double t, o;

   asm ("fptan" : "=t" (o), "=u" (t) : "0" (a));

   return t;
}
//------------------------------------------------------------------------------
inline double ATan(double y, double x)
{
double a;

   asm ("fpatan" : "=t" (a) : "0" (x), "u" (y) : "st(1)");

   return a;
}
//------------------------------------------------------------------------------