Changeset b39d4d in git for Singular/mpr_complex.cc
- Timestamp:
- Dec 18, 2000, 4:44:42 PM (23 years ago)
- Branches:
- (u'spielwiese', 'a7324b6e0b44a1a8ed3fa4d9ca3e2ff210ddd52c')
- Children:
- e3dc1cc75eb8813142d38ad4d5b8638f5748228a
- Parents:
- 6a9b5a892f6ba5ac58be777a540985149140afae
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/mpr_complex.cc
r6a9b5a rb39d4d 2 2 * Computer Algebra System SINGULAR * 3 3 ****************************************/ 4 /* $Id: mpr_complex.cc,v 1.2 8 2000-11-21 15:35:44 SingularExp $ */4 /* $Id: mpr_complex.cc,v 1.29 2000-12-18 15:44:41 obachman Exp $ */ 5 5 6 6 /* … … 63 63 gmp_output_digits= digits; 64 64 mpf_set_default_prec( db ); 65 mpf_set_prec( (mpf_ptr)*(diff.mpfp()),32);66 mpf_set_prec( (mpf_ptr)*gmpRel.mpfp(),32);67 mpf_set_d( (mpf_ptr)*gmpRel.mpfp(),0.1);68 mpf_pow_ui( (mpf_ptr)*gmpRel.mpfp(),(mpf_ptr)*gmpRel.mpfp(),digits);65 mpf_set_prec(*diff._mpfp(),32); 66 mpf_set_prec(*gmpRel._mpfp(),32); 67 mpf_set_d(*gmpRel._mpfp(),0.1); 68 mpf_pow_ui(*gmpRel._mpfp(),*gmpRel._mpfp(),digits); 69 69 } 70 70 … … 706 706 gmp_float ai=this->imag(); 707 707 if (ar.isZero() || ai.isZero()) return; 708 mpf_abs( (mpf_ptr)*ar.mpfp(), (mpf_ptr)*ar.mpfp());709 mpf_abs( (mpf_ptr)*ai.mpfp(), (mpf_ptr)*ai.mpfp());710 mpf_set_prec( (mpf_ptr)*ar.mpfp(), 32);711 mpf_set_prec( (mpf_ptr)*ai.mpfp(), 32);708 mpf_abs(*ar._mpfp(), *ar._mpfp()); 709 mpf_abs(*ai._mpfp(), *ai._mpfp()); 710 mpf_set_prec(*ar._mpfp(), 32); 711 mpf_set_prec(*ai._mpfp(), 32); 712 712 if (ar > ai) 713 713 { 714 mpf_div( (mpf_ptr)*ai.mpfp(), (mpf_ptr)*ai.mpfp(), (mpf_ptr)*ar.mpfp());714 mpf_div(*ai._mpfp(), *ai._mpfp(), *ar._mpfp()); 715 715 if (ai < gmpRel) this->imag(0.0); 716 716 } 717 717 else 718 718 { 719 mpf_div( (mpf_ptr)*ar.mpfp(), (mpf_ptr)*ar.mpfp(), (mpf_ptr)*ai.mpfp());719 mpf_div(*ar._mpfp(), *ar._mpfp(), *ai._mpfp()); 720 720 if (ar < gmpRel) this->real(0.0); 721 721 }
Note: See TracChangeset
for help on using the changeset viewer.