Логарифм по произвольному основанию

double logarithm(double b, double a, int *err)
{   double res = 0;
   *err = 0;
   if (b==1 || b<0) *err = 1;
   else if (a<0) *err = 2;
   else res = log(a) / log(b);
   return res;
}

#include <iostream>
#include <cmath>
using namespace std;
static inline unsigned long long int tick()
{
unsigned long long int d;
__asm__ __volatile__ ("rdtsc" : "=A" (d) );
return d;
}
 
const int AIM = 1000; // кол-во проходов
 
unsigned long long test()
{
    unsigned long long int a = 0;
    unsigned long long int b = 0;
    a=tick();
    log10(99999999999999999);      //тестируемая функция
    b=tick();
    return b-a;
}
 
int main()
{
    unsigned long long res = 0;
    for (int i = 0; i < AIM; ++i)
        res += test();
    cout << res/AIM << endl;
}