     __complex__ float g_b;
     g_a = s_a.real() + s_a.imag()*__I__;
     g_b = s_b.real() + s_b.imag()*__I__;
-    // Ok here we screwed big time
     b = cusp::acos(a);
     g_b = cacosf(g_a);
     s_b = std::complex<ValueType>(creal(g_b),cimag(g_b));


    g_a = s_a.real() + s_a.imag()*__I__;
    g_b = s_b.real() + s_b.imag()*__I__;
    // Ok here we screwed big time
    b = cusp::atan(a);
    g_b = catanf(g_a);
    s_b = std::complex<ValueType>(creal(g_b),cimag(g_b));

  ASSERT_COMPLEX_ALMOST_EQUAL(b,s_b);
  b = cusp::log10(a);
  s_b = std::log10(s_a);
  ASSERT_COMPLEX_ALMOST_EQUAL(b,s_b);

  b = cusp::pow(a,b);
  s_b = std::pow(s_a,s_b);
  ASSERT_COMPLEX_ALMOST_EQUAL(b,s_b);
  b = cusp::pow(a,ValueType(1.3));
  s_b = std::pow(s_a,ValueType(1.3));
