void expand_bvpprod(bvexp_table_t *table, bvarith_buffer_t *buffer, pp_buffer_t *p, uint32_t n, uint32_t *c) { bv_vartable_t *vtbl; bvmlist_t *q; pp_buffer_t *aux; pprod_t *r; uint32_t *a; uint32_t i, m, d, e, k; thvar_t x; assert(buffer->store == &table->store && buffer->ptbl == &table->pprods); bvarith_buffer_prepare(buffer, n); bvarith_buffer_set_one(buffer); aux = p; vtbl = table->vtbl; if (total_degree_test(table, vtbl, p)) { aux = &table->pp; pp_buffer_reset(aux); // make a copy of c in the internal bvconst buffer bvconstant_copy(&table->bvconst, n, c); c = table->bvconst.data; k = (n + 31) >> 5; m = p->len; for (i=0; i<m; i++) { x = p->prod[i].var; d = p->prod[i].exp; if (bvvar_is_const(vtbl, x)) { a = bvvar_val(vtbl, x); bvconst_mulpower(c, k, a, d); } else { q = bvexp_def(table, x); if (q != NULL && mlist_is_short(q, &e) && d * e <= BVEXP_DEGREE_LIMIT) { bvarith_buffer_mul_mlist_power(buffer, q, d, &table->aux); } else { pp_buffer_mul_varexp(aux, x, d); } } } // normalize bvconst_normalize(c, n); pp_buffer_normalize(aux); } /* * The result is c * aux * buffer */ r = pprod_from_buffer(&table->pprods, aux); bvarith_buffer_mul_mono(buffer, c, r); bvarith_buffer_normalize(buffer); }
void expand_bvpoly(bvexp_table_t *table, bvarith_buffer_t *buffer, bvpoly_buffer_t *p) { bv_vartable_t *vtbl; bvmlist_t *q; uint32_t *c; uint32_t i, n; thvar_t x; assert(buffer->store == &table->store && buffer->ptbl == &table->pprods); bvarith_buffer_prepare(buffer, bvpoly_buffer_bitsize(p)); n = bvpoly_buffer_num_terms(p); if (n > 0) { vtbl = table->vtbl; i = 0; // constant term of p if (bvpoly_buffer_var(p, 0) == const_idx) { bvarith_buffer_add_const(buffer, bvpoly_buffer_coeff(p, 0)); i ++; } // non-constant terms while (i < n) { x = bvpoly_buffer_var(p, i); c = bvpoly_buffer_coeff(p, i); i ++; if (bvvar_is_const(vtbl, x)) { bvarith_buffer_add_const_times_const(buffer, c, bvvar_val(vtbl, x)); } else { q = bvexp_def(table, x); if (q != NULL) { bvarith_buffer_add_const_times_mlist(buffer, q, c); } else { bvarith_buffer_add_varmono(buffer, c, x); } } } bvarith_buffer_normalize(buffer); } }
/* * Print the definition of x in vtbl */ static void print_bv_vardef(FILE *f, bv_vartable_t *vtbl, thvar_t x) { uint32_t nbits; assert(valid_bvvar(vtbl, x)); nbits = bvvar_bitsize(vtbl, x); print_bvvar(f, x); fprintf(f, ":bv[%"PRIu32"] = ", nbits); switch (bvvar_tag(vtbl, x)) { case BVTAG_VAR: fputs("var", f); break; case BVTAG_CONST64: bvconst64_print(f, bvvar_val64(vtbl, x), nbits); break; case BVTAG_CONST: bvconst_print(f, bvvar_val(vtbl, x), nbits); break; case BVTAG_POLY64: print_bv_poly64(f, bvvar_poly64_def(vtbl, x)); break; case BVTAG_POLY: print_bv_poly(f, bvvar_poly_def(vtbl, x)); break; case BVTAG_PPROD: print_bv_product(f, bvvar_pprod_def(vtbl, x)); break; case BVTAG_BIT_ARRAY: print_litarray(f, nbits, bvvar_bvarray_def(vtbl, x)); break; case BVTAG_ITE: print_bv_ite(f, bvvar_ite_def(vtbl, x)); break; case BVTAG_UDIV: print_bv_binop(f, "div", bvvar_binop(vtbl, x)); break; case BVTAG_UREM: print_bv_binop(f, "rem", bvvar_binop(vtbl, x)); break; case BVTAG_SDIV: print_bv_binop(f, "sdiv", bvvar_binop(vtbl, x)); break; case BVTAG_SREM: print_bv_binop(f, "srem", bvvar_binop(vtbl, x)); break; case BVTAG_SMOD: print_bv_binop(f, "smod", bvvar_binop(vtbl, x)); break; case BVTAG_SHL: print_bv_binop(f, "shl", bvvar_binop(vtbl, x)); break; case BVTAG_LSHR: print_bv_binop(f, "lshr", bvvar_binop(vtbl, x)); break; case BVTAG_ASHR: print_bv_binop(f, "ashr", bvvar_binop(vtbl, x)); break; case BVTAG_ADD: print_bv_binop(f, "add", bvvar_binop(vtbl, x)); break; case BVTAG_SUB: print_bv_binop(f, "sub", bvvar_binop(vtbl, x)); break; case BVTAG_MUL: print_bv_binop(f, "mul", bvvar_binop(vtbl, x)); break; case BVTAG_NEG: print_bvneg(f, bvvar_binop(vtbl, x)); break; } }