   ASSERT_COMPLEX_ALMOST_EQUAL(b,s_b);
   b = cusp::tan(a);  
   s_b = std::tan(s_a);
-  // There's a difference here and it's a close call 
-  // but I think the cusp value came closer
-  // to the value calculated with arbitrary precision
   ASSERT_COMPLEX_ALMOST_EQUAL(b,s_b);
   b = cusp::tanh(a);
   s_b = std::tanh(s_a);
-  // Here we're both clearly the wrong
   ASSERT_COMPLEX_ALMOST_EQUAL(b,s_b);
 
   // Inverse trigonometic functions not part of standard 

    __complex__ double 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 = cacos(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__;
    b = cusp::asin(a);
    g_b = casinf(g_a);
    s_b = std::complex<ValueType>(creal(g_b),cimag(g_b));
