-void print_real_number(ostream & os, const cl_R & num)
-{
- cl_print_flags ourflags;
- if (::instanceof(num, ::cl_RA_ring)) {
- // case 1: integer or rational, nothing special to do:
- ::print_real(os, ourflags, num);
- } else {
- // case 2: float
- // make CLN believe this number has default_float_format, so it prints
- // 'E' as exponent marker instead of 'L':
- ourflags.default_float_format = ::cl_float_format(The(cl_F)(num));
- ::print_real(os, ourflags, num);
- }
- return;
-}
-
-/** This method adds to the output so it blends more consistently together
- * with the other routines and produces something compatible to ginsh input.
- *
- * @see print_real_number() */
-void numeric::print(ostream & os, unsigned upper_precedence) const
-{
- debugmsg("numeric print", LOGLEVEL_PRINT);
- if (this->is_real()) {
- // case 1, real: x or -x
- if ((precedence<=upper_precedence) && (!this->is_nonneg_integer())) {
- os << "(";
- print_real_number(os, The(cl_R)(*value));
- os << ")";
- } else {
- print_real_number(os, The(cl_R)(*value));
- }
- } else {
- // case 2, imaginary: y*I or -y*I
- if (::realpart(*value) == 0) {
- if ((precedence<=upper_precedence) && (::imagpart(*value) < 0)) {
- if (::imagpart(*value) == -1) {
- os << "(-I)";
- } else {
- os << "(";
- print_real_number(os, The(cl_R)(::imagpart(*value)));
- os << "*I)";
- }
- } else {
- if (::imagpart(*value) == 1) {
- os << "I";
- } else {
- if (::imagpart (*value) == -1) {
- os << "-I";
- } else {
- print_real_number(os, The(cl_R)(::imagpart(*value)));
- os << "*I";
- }
- }
- }
- } else {
- // case 3, complex: x+y*I or x-y*I or -x+y*I or -x-y*I
- if (precedence <= upper_precedence)
- os << "(";
- print_real_number(os, The(cl_R)(::realpart(*value)));
- if (::imagpart(*value) < 0) {
- if (::imagpart(*value) == -1) {
- os << "-I";
- } else {
- print_real_number(os, The(cl_R)(::imagpart(*value)));
- os << "*I";
- }
- } else {
- if (::imagpart(*value) == 1) {
- os << "+I";
- } else {
- os << "+";
- print_real_number(os, The(cl_R)(::imagpart(*value)));
- os << "*I";
- }
- }
- if (precedence <= upper_precedence)
- os << ")";
- }
- }
-}
-
-
-void numeric::printraw(ostream & os) const
-{
- // The method printraw doesn't do much, it simply uses CLN's operator<<()
- // for output, which is ugly but reliable. e.g: 2+2i
- debugmsg("numeric printraw", LOGLEVEL_PRINT);
- os << "numeric(" << *value << ")";
-}
-
-
-void numeric::printtree(ostream & os, unsigned indent) const
-{
- debugmsg("numeric printtree", LOGLEVEL_PRINT);
- os << string(indent,' ') << *value
- << " (numeric): "
- << "hash=" << hashvalue << " (0x" << hex << hashvalue << dec << ")"
- << ", flags=" << flags << endl;
-}
-
-
-void numeric::printcsrc(ostream & os, unsigned type, unsigned upper_precedence) const
-{
- debugmsg("numeric print csrc", LOGLEVEL_PRINT);
- ios::fmtflags oldflags = os.flags();
- os.setf(ios::scientific);
- if (this->is_rational() && !this->is_integer()) {
- if (compare(_num0()) > 0) {
- os << "(";
- if (type == csrc_types::ctype_cl_N)
- os << "cl_F(\"" << numer().evalf() << "\")";
- else
- os << numer().to_double();
- } else {
- os << "-(";
- if (type == csrc_types::ctype_cl_N)
- os << "cl_F(\"" << -numer().evalf() << "\")";
- else
- os << -numer().to_double();
- }
- os << "/";
- if (type == csrc_types::ctype_cl_N)
- os << "cl_F(\"" << denom().evalf() << "\")";
- else
- os << denom().to_double();
- os << ")";
- } else {
- if (type == csrc_types::ctype_cl_N)
- os << "cl_F(\"" << evalf() << "\")";
- else
- os << to_double();
- }
- os.flags(oldflags);
+static void print_real_cl_N(const print_context & c, const cln::cl_R & x)
+{
+ if (cln::instanceof(x, cln::cl_I_ring)) {
+
+ // Integer number
+ c.s << "cln::cl_I(\"";
+ print_real_number(c, x);
+ c.s << "\")";
+
+ } else if (cln::instanceof(x, cln::cl_RA_ring)) {
+
+ // Rational number
+ cln::cl_print_flags ourflags;
+ c.s << "cln::cl_RA(\"";
+ cln::print_rational(c.s, ourflags, cln::the<cln::cl_RA>(x));
+ c.s << "\")";
+
+ } else {
+
+ // Anything else
+ c.s << "cln::cl_F(\"";
+ print_real_number(c, cln::cl_float(1.0, cln::default_float_format) * x);
+ c.s << "_" << Digits << "\")";
+ }
+}
+
+void numeric::print_numeric(const print_context & c, const char *par_open, const char *par_close, const char *imag_sym, const char *mul_sym, unsigned level) const
+{
+ const cln::cl_R r = cln::realpart(value);
+ const cln::cl_R i = cln::imagpart(value);
+
+ if (cln::zerop(i)) {
+
+ // case 1, real: x or -x
+ if ((precedence() <= level) && (!this->is_nonneg_integer())) {
+ c.s << par_open;
+ print_real_number(c, r);
+ c.s << par_close;
+ } else {
+ print_real_number(c, r);
+ }
+
+ } else {
+ if (cln::zerop(r)) {
+
+ // case 2, imaginary: y*I or -y*I
+ if (i == 1)
+ c.s << imag_sym;
+ else {
+ if (precedence()<=level)
+ c.s << par_open;
+ if (i == -1)
+ c.s << "-" << imag_sym;
+ else {
+ print_real_number(c, i);
+ c.s << mul_sym << imag_sym;
+ }
+ if (precedence()<=level)
+ c.s << par_close;
+ }
+
+ } else {
+
+ // case 3, complex: x+y*I or x-y*I or -x+y*I or -x-y*I
+ if (precedence() <= level)
+ c.s << par_open;
+ print_real_number(c, r);
+ if (i < 0) {
+ if (i == -1) {
+ c.s << "-" << imag_sym;
+ } else {
+ print_real_number(c, i);
+ c.s << mul_sym << imag_sym;
+ }
+ } else {
+ if (i == 1) {
+ c.s << "+" << imag_sym;
+ } else {
+ c.s << "+";
+ print_real_number(c, i);
+ c.s << mul_sym << imag_sym;
+ }
+ }
+ if (precedence() <= level)
+ c.s << par_close;
+ }
+ }