7u83/stc_rpncalc
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

NUM_A_ARRAY cannot be less than num_lsu

Browse Source
This commit is contained in:
Jeff Wang 2019-10-02 22:45:02 -04:00
parent ac1aaac096
commit 23e49a0405
1 changed files with 7 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/decn/proto/ln_mfp.cpp
View File

@ -14,7 +14,7 @@ static const int BSIZE = 200;
char buf[BSIZE]; char buf[BSIZE];
char pbuf[BSIZE]; char pbuf[BSIZE];
static const int NUM_A_ARR = 10; static const int NUM_A_ARR = 9;
mpfr_float a_arr[NUM_A_ARR]; mpfr_float a_arr[NUM_A_ARR];
mpfr_float ln_a_arr[NUM_A_ARR]; mpfr_float ln_a_arr[NUM_A_ARR];
@ -23,6 +23,9 @@ static const unsigned int CALC_PRECISION = 18;
int main(void){ int main(void){
cout << std::scientific << std::setprecision(CALC_PRECISION); cout << std::scientific << std::setprecision(CALC_PRECISION);
mpfr_float a(10, CALC_PRECISION);
a = log(a);
cout << "ln(10):" << a << endl;
//initiallize a_arr and ln_a_arr //initiallize a_arr and ln_a_arr
for (int i = 0; i < NUM_A_ARR; i++){ for (int i = 0; i < NUM_A_ARR; i++){
// a[i] = (1 + 10^-i) // a[i] = (1 + 10^-i)
@ -47,6 +50,7 @@ int main(void){
actual = log(A); actual = log(A);
int num_times[NUM_A_ARR]; int num_times[NUM_A_ARR];
//track number of times multiplied by a_arr[j]
for (int j = 0; j < NUM_A_ARR; j++){ for (int j = 0; j < NUM_A_ARR; j++){
if (j != 0){ if (j != 0){
b_j *= 10.0; b_j *= 10.0;
@ -57,6 +61,8 @@ int main(void){
b_j = accum; b_j = accum;
accum *= a_arr[j]; accum *= a_arr[j];
accum -= 1.0; accum -= 1.0;
// printf(" %d: ", k_j);
// cout << b_j << endl;
k_j++; k_j++;
} }
num_times[j] = k_j - 1; num_times[j] = k_j - 1;