void mpq_manager<SYNCH>::display_decimal(std::ostream & out, mpq const & a, unsigned prec) { mpz n1, d1, v1; get_numerator(a, n1); get_denominator(a, d1); if (is_neg(a)) { out << "-"; neg(n1); } mpz ten(10); div(n1, d1, v1); display(out, v1); rem(n1, d1, n1); if (is_zero(n1)) goto end; // number is an integer out << "."; for (unsigned i = 0; i < prec; i++) { mul(n1, ten, n1); div(n1, d1, v1); SASSERT(lt(v1, ten)); display(out, v1); rem(n1, d1, n1); if (is_zero(n1)) goto end; // number is precise } out << "?"; end: del(ten); del(n1); del(d1); del(v1); }
//add the two numbers char* add_nums(const char* num1,const char* base1,const char* num2,const char* base2) { const char* base_str[2]; const char* num_str[2]; char *p,*q; int base_gcd; int i,test_res; int base[2]; ln numerator[3],denominator[3]; vector digits; //test bases base_str[0]=base1; base_str[1]=base2; for(i=0;i<2;i++) { base[i]=test_base(base_str[i]); if(base[i]==-1) { printf("Error: invalid base [%s]!\n",base_str[i]); return NULL; } else if(base[i]==-2) { printf("Error: base [%s] is out of range "MACRO2STR(MIN_BASE)"-"MACRO2STR(MAX_BASE)"!\n",base_str[i]); return NULL; } } //test numbers num_str[0]=num1; num_str[1]=num2; for(i=0;i<2;i++) { test_res=test_num(num_str[i],base[i]); if(test_res==0) { printf("Error: invalid number [%s],the number should match the regular expression [0-9]+(.[0-9]+)? !\n",num_str[i]); return NULL; } else if(test_res==2) { printf("Error: invalid number [%s] in base %d !\n",num_str[i],base[i]); return NULL; } } //convert the two numbers into fraction form for(i=0;i<2;i++) { denominator[i]=get_denominator(num_str[i],base[i]); numerator[i]=get_numerator(num_str[i],base[i]); } for(i=0;i<2;i++) { p=ln2str(numerator[i]); q=ln2str(denominator[i]); printf("operand in fraction:%s/%s\n",p,q); fflush(stdout); free(p); free(q); } //compute the fraction denominator[2]=ln_multiply(denominator[0],denominator[1],newln); numerator[0]=ln_multiply(numerator[0],denominator[1],firstln); numerator[1]=ln_multiply(numerator[1],denominator[0],firstln); numerator[2]=ln_add(numerator[0],numerator[1],newln); //free unused numbers for(i=0;i<2;i++) { ln_free(&(numerator[i])); ln_free(&(denominator[i])); } p=ln2str(numerator[2]); q=ln2str(denominator[2]); printf("result in fraction:%s/%s\n",p,q); //simplifying the result fraction base_gcd=gcd(base[0],base[1]); if(base_gcd !=1) { while(ln_divideable_num(numerator[2],base_gcd)==1 && ln_divideable_num(denominator[2],base_gcd)==1) { ln_divide_int(numerator[2],base_gcd,0,trunc_res,firstln); ln_divide_int(denominator[2],base_gcd,0,trunc_res,firstln); } } p=ln2str(numerator[2]); q=ln2str(denominator[2]); printf("result in fraction:%s/%s\n",p,q); fflush(stdout); if(strcmp(q,"1")==0) { free(q); q=(char*)malloc(strlen(p)+4); if(!q) { ln_free(&(numerator[2])); ln_free(&(denominator[2])); printf("Error: coverting fraction to decimal number failed at line %d",__LINE__); return NULL; } sprintf(q,"%s 10",p); free(p); return q; } free(p); free(q); digits=get_decimal_digit(numerator[2],denominator[2]); if(!digits) { ln_free(&(numerator[2])); ln_free(&(denominator[2])); printf("Error: coverting fraction to decimal number failed at line %d",__LINE__); return NULL; } p=get_decimalstr(digits,denominator[2]); if(!p) { ln_free(&(numerator[2])); ln_free(&(denominator[2])); free_digit(digits); printf("Error: coverting fraction to decimal number failed at line %d",__LINE__); return NULL; } return p; }