Oskar,
I'm sure these are all related. I did not mean to imply that your
compiler was broken! Only that as I myself only have one, the only
testing I have done is with that one; and since Sage has started
shipping with eclib (which is about 6 months) I know that lots of
other compilers out the
Hi John,
*JUST FOR TESTING OUT MY PREVIOUS SUSPICION*, I changed to
bigcomplex cubic::hess_root() const
{
bigfloat discr = I2bigfloat(disc());
// if(!is_positive(disc()))
//{
// cout<<"Error: hess_root called with negative dicriminant!\n";
// return to_bigfloat(0);
//}
bi
Oskar,
I think you are probably right there. Of course that code is never
executed in practice since the way the code is organised this function
is only called when disc is negative, but of course that is not going
to keep the compiler happy. So the line
return to_bigfloat(0);
should be repla
Hi John,
thanks for your answer. Maybe it's of no relevance for the described
problem but I noticed the following anyway:
In procs/interface.h, I have
...
#include
#define bigfloat RR
...
#include
typedef complex CC;
#define bigcomplex CC
...
inline RR to_bigfloat(const int& n) {return to_RR(
Well, as the author of eclib I hope I can help. It is a C++ library.
It uses the NTL library, whcih includes a multiprecision real number
type called RR. MTL does not have multiprecision complexes, so
implement those using RR for the real and imaginary parts, and a
standard C++ library template
