/* recursively convert an AST expression into a BDD */
DdNode* bdd_expr(DdManager* m, symbol** symtab, ast_node* expr) {
int i;
DdNode* ret, *tmp1, *tmp2;
switch (expr->tag) {
case _id_expr: /* simply return that variable's BDD */
if ((i = symtab_lookup(symtab, expr->name)) < 0) {
printf("Error: variable '%s' not found in symbol table\n", expr->name);
exit(1);
}
ret = Cudd_bddIthVar(m, symtab_lookup(symtab, expr->name));
break;
case _lit_expr: /* return a BDD 1 or 0 */
ret = expr->val ? Cudd_ReadOne(m) : Cudd_ReadLogicZero(m);
break;
case _not_expr: /* return complement BDD */
tmp1 = bdd_expr(m, symtab, expr->children[EXPR1]);
ret = Cudd_Not(tmp1);
Cudd_Ref(ret);
Cudd_RecursiveDeref(m, tmp1);
break;
case _and_expr: /* return logical AND */
tmp1 = bdd_expr(m, symtab, expr->children[EXPR1]);
tmp2 = bdd_expr(m, symtab, expr->children[EXPR2]);
ret = Cudd_bddAnd(m, tmp1, tmp2);
Cudd_Ref(ret);
Cudd_RecursiveDeref(m, tmp1);
Cudd_RecursiveDeref(m, tmp2);
break;
case _or_expr: /* return logical OR */
tmp1 = bdd_expr(m, symtab, expr->children[EXPR1]);
tmp2 = bdd_expr(m, symtab, expr->children[EXPR2]);
ret = Cudd_bddOr(m, tmp1, tmp2);
Cudd_Ref(ret);
Cudd_RecursiveDeref(m, tmp1);
Cudd_RecursiveDeref(m, tmp2);
break;
case _eq_expr: /* check for equality, return 1 or 0 */
tmp1 = bdd_expr(m, symtab, expr->children[EXPR1]);
tmp2 = bdd_expr(m, symtab, expr->children[EXPR2]);
if (Cudd_EquivDC(m, tmp1, tmp2, Cudd_ReadLogicZero(m)))
ret = Cudd_ReadOne(m);
else
ret = Cudd_ReadLogicZero(m);
Cudd_RecursiveDeref(m, tmp1);
Cudd_RecursiveDeref(m, tmp2);
break;
case _impl_expr: /* return logical implication (!a|b) */
tmp1 = Cudd_Not(bdd_expr(m, symtab, expr->children[EXPR1]));
tmp2 = bdd_expr(m, symtab, expr->children[EXPR2]);
ret = Cudd_bddOr(m, tmp1, tmp2);
Cudd_Ref(ret);
Cudd_RecursiveDeref(m, tmp1);
Cudd_RecursiveDeref(m, tmp2);
break;
default:
printf("Error: invalid expression\n");
exit(1);
}
return ret;
}
static double Expectation(DdNode **nodes_ex, int lenNodes) {
int i;
double rootProb, CLL = 0;
for (i = 0; i < lenNodes; i++) {
if (!Cudd_IsConstant(nodes_ex[i])) {
nodesB = init_table(boolVars_ex[i]);
nodesF = init_table(boolVars_ex[i]);
Forward(nodes_ex[i], i);
rootProb = GetOutsideExpe(nodes_ex[i], example_prob[i], i);
if (rootProb <= 0.0)
CLL = CLL + LOGZERO * example_prob[i];
else
CLL = CLL + log(rootProb) * example_prob[i];
nodes_probs_ex[i] = rootProb;
destroy_table(nodesB, boolVars_ex[i]);
destroy_table(nodesF, boolVars_ex[i]);
} else if (nodes_ex[i] == Cudd_ReadLogicZero(mgr_ex[i])) {
CLL = CLL + LOGZERO * example_prob[i];
nodes_probs_ex[i] = 0.0;
} else
nodes_probs_ex[i] = 1.0;
}
return CLL;
}
BddBuilder::BddBuilder(BddBuilder * pA){
__pddWireHead = __pddWireTail = NULL;
__pddOutputWireHead = __pddOutputWireTail = NULL;
__pddOutputNodes = NULL;
__pddGateHead = __pddGateTail = NULL;
__pddInputNodes = NULL;
__pddManager = pA->__pddManager;
__Vcc = Cudd_ReadOne(__pddManager);
__GND = Cudd_ReadLogicZero(__pddManager);
char * name_vcc = new char[4];
char * name_gnd = new char[4];
sprintf(name_vcc, "Vcc");
sprintf(name_gnd, "GND");
__VccWire = new DdWire(__pddManager, name_vcc, 0, 0);
__GNDWire = new DdWire(__pddManager, name_gnd, 0, 0);
__VccWire->setDdNode(__Vcc);
__GNDWire->setDdNode(__GND);
if(__pddManager == NULL){
perror("DdManager initializing error.");
}
//__inputWireCnt = 2; // Vcc & GND
__inputWireCnt = -1;
__outputWireCnt = -1;
__pddInputWireHead = __VccWire;
__VccWire->setInputListNext(__GNDWire);
__pddInputWireTail = __GNDWire;
}
int zeroc(void)
{
DdNode * node;
node=Cudd_ReadLogicZero(mgr);
Cudd_Ref(node);
return (int) node;
}
ref_t shadow_zero(shadow_mgr mgr) {
if (do_ref(mgr))
return REF_ZERO;
else {
DdNode * n = Cudd_ReadLogicZero(mgr->bdd_manager);
reference_dd(mgr, n);
return dd2ref(n, IS_BDD);
}
}
/* General conversion from ref to dd */
static DdNode *get_ddnode(shadow_mgr mgr, ref_t r) {
DdNode *n;
if (!do_ref(mgr)) {
n = ref2dd(mgr, r);
} else if (REF_IS_INVALID(r)) {
err(fatal, "No node associated with invalid ref");
n = Cudd_ReadLogicZero(mgr->bdd_manager);
reference_dd(mgr, n);
} else if (!keyvalue_find(mgr->r2c_table, (word_t ) r, (word_t *) &n)) {
char buf[24];
shadow_show(mgr, r, buf);
err(fatal, "No node associated with ref %s (0x%llx)", buf, r);
n = mgr->do_cudd ? Cudd_ReadLogicZero(mgr->bdd_manager) :
(DdNode *) REF_ZERO;
reference_dd(mgr, n);
}
return n;
}
static YAP_Bool zero(void) {
YAP_Term arg, out;
DdNode *node;
arg = YAP_ARG1;
node = Cudd_ReadLogicZero(mgr_ex[ex]);
Cudd_Ref(node);
out = YAP_MkIntTerm((YAP_Int)node);
return (YAP_Unify(out, arg));
}
DdNode* compute_EF(DdManager* m, pos* postab, DdNode* R, DdNode* p) {
DdNode* last = Cudd_ReadOne(m);
DdNode* current = p;
while (!Cudd_EquivDC(m, last, current, Cudd_ReadLogicZero(m))) {
// printf("iterating\n");
last = Cudd_bddAnd(m, current, Cudd_ReadOne(m));
current = Cudd_bddOr(m, p, compute_EX(m, postab, R, current));
}
return current;
}
/* wrap the CFG to BDD conversion */
DdNode* encode_prog(DdManager* m, pos* postab, cfg_node* proc_list_head) {
DdNode* tmp, *proc, *init;
DdNode* R = Cudd_ReadLogicZero(m);
Cudd_Ref(R);
cfg_node* next = proc_list_head;
int proc_no = 0; /* all process begin with an 'init' */
while (next) { /* skip node that has id = 0 */
// init = Cudd_ReadOne(m);
// init = bdd_encode_pc(m, init, postab->pc, postab->pc_size, proc_no, 0);
init = bdd_mk_skip(m, postab, 0, proc_no);
proc = bdd_convert_cfg(m, postab, proc_list_head, proc_no, init, R, 0);
tmp = Cudd_bddOr(m, R, proc); /* add process to transition relation */
Cudd_Ref(tmp);
Cudd_RecursiveDeref(m, R);
Cudd_RecursiveDeref(m, proc);
R = tmp;
next = next->next_proc; /* next process */
++proc_no;
}
return R; /* return complete transition relation */
}
BddBuilder::BddBuilder(){
__pddWireHead = __pddWireTail = NULL;
__pddOutputWireHead = __pddOutputWireTail = NULL;
__pddOutputNodes = NULL;
__pddGateHead = __pddGateTail = NULL;
__pddInputNodes = NULL;
__pddManager = Cudd_Init(0, 0,CUDD_UNIQUE_SLOTS,CUDD_CACHE_SLOTS , 0);
Cudd_ReduceHeap(__pddManager,CUDD_REORDER_RANDOM_PIVOT,0);
Cudd_AutodynEnable(__pddManager,CUDD_REORDER_RANDOM_PIVOT);
__Vcc = Cudd_ReadOne(__pddManager);
__GND = Cudd_ReadLogicZero(__pddManager);
char * name_vcc = new char[4];
char * name_gnd = new char[4];
sprintf(name_vcc, "Vcc");
sprintf(name_gnd, "GND");
__VccWire = new DdWire(__pddManager, name_vcc, 0, 0);
__GNDWire = new DdWire(__pddManager, name_gnd, 0, 0);
__VccWire->setDdNode(__Vcc);
__GNDWire->setDdNode(__GND);
Cudd_Ref(__Vcc);
Cudd_Ref(__GND);
if(__pddManager == NULL){
perror("DdManager initializing error.");
}
//__inputWireCnt = 2; // Vcc & GND
__inputWireCnt = -1;
__outputWireCnt = -1;
__pddInputWireHead = __VccWire;
__VccWire->setInputListNext(__GNDWire);
__pddInputWireTail = __GNDWire;
}
/**Function*************************************************************
Synopsis [Attempts to solve the miter using a number of tricks.]
Description [Returns -1 if timed out; 0 if SAT; 1 if UNSAT. Returns
a simplified version of the original network (or a constant 0 network).
In case the network is not a constant zero and a SAT assignment is found,
pNtk->pModel contains a satisfying assignment.]
SideEffects []
SeeAlso []
***********************************************************************/
int Abc_NtkMiterProve( Abc_Ntk_t ** ppNtk, void * pPars )
{
Prove_Params_t * pParams = pPars;
Abc_Ntk_t * pNtk, * pNtkTemp;
int RetValue, nIter, nSatFails, Counter, clk, timeStart = clock();
sint64 nSatConfs, nSatInspects, nInspectLimit;
// get the starting network
pNtk = *ppNtk;
assert( Abc_NtkIsStrash(pNtk) );
assert( Abc_NtkPoNum(pNtk) == 1 );
if ( pParams->fVerbose )
{
printf( "RESOURCE LIMITS: Iterations = %d. Rewriting = %s. Fraiging = %s.\n",
pParams->nItersMax, pParams->fUseRewriting? "yes":"no", pParams->fUseFraiging? "yes":"no" );
printf( "Mitering = %d (%3.1f). Rewriting = %d (%3.1f). Fraiging = %d (%3.1f).\n",
pParams->nMiteringLimitStart, pParams->nMiteringLimitMulti,
pParams->nRewritingLimitStart, pParams->nRewritingLimitMulti,
pParams->nFraigingLimitStart, pParams->nFraigingLimitMulti );
printf( "Mitering last = %d.\n",
pParams->nMiteringLimitLast );
}
// if SAT only, solve without iteration
if ( !pParams->fUseRewriting && !pParams->fUseFraiging )
{
clk = clock();
RetValue = Abc_NtkMiterSat( pNtk, (sint64)pParams->nMiteringLimitLast, (sint64)0, 0, NULL, NULL );
Abc_NtkMiterPrint( pNtk, "SAT solving", clk, pParams->fVerbose );
*ppNtk = pNtk;
return RetValue;
}
// check the current resource limits
for ( nIter = 0; nIter < pParams->nItersMax; nIter++ )
{
if ( pParams->fVerbose )
{
printf( "ITERATION %2d : Confs = %6d. FraigBTL = %3d. \n", nIter+1,
(int)(pParams->nMiteringLimitStart * pow(pParams->nMiteringLimitMulti,nIter)),
(int)(pParams->nFraigingLimitStart * pow(pParams->nFraigingLimitMulti,nIter)) );
fflush( stdout );
}
// try brute-force SAT
clk = clock();
nInspectLimit = pParams->nTotalInspectLimit? pParams->nTotalInspectLimit - pParams->nTotalInspectsMade : 0;
RetValue = Abc_NtkMiterSat( pNtk, (sint64)(pParams->nMiteringLimitStart * pow(pParams->nMiteringLimitMulti,nIter)), (sint64)nInspectLimit, 0, &nSatConfs, &nSatInspects );
Abc_NtkMiterPrint( pNtk, "SAT solving", clk, pParams->fVerbose );
if ( RetValue >= 0 )
break;
// add to the number of backtracks and inspects
pParams->nTotalBacktracksMade += nSatConfs;
pParams->nTotalInspectsMade += nSatInspects;
// check if global resource limit is reached
if ( (pParams->nTotalBacktrackLimit && pParams->nTotalBacktracksMade >= pParams->nTotalBacktrackLimit) ||
(pParams->nTotalInspectLimit && pParams->nTotalInspectsMade >= pParams->nTotalInspectLimit) )
{
printf( "Reached global limit on conflicts/inspects. Quitting.\n" );
*ppNtk = pNtk;
return -1;
}
// try rewriting
if ( pParams->fUseRewriting )
{
clk = clock();
Counter = (int)(pParams->nRewritingLimitStart * pow(pParams->nRewritingLimitMulti,nIter));
// Counter = 1;
while ( 1 )
{
/*
extern Abc_Ntk_t * Abc_NtkIvyResyn( Abc_Ntk_t * pNtk, int fUpdateLevel, int fVerbose );
pNtk = Abc_NtkIvyResyn( pNtkTemp = pNtk, 0, 0 ); Abc_NtkDelete( pNtkTemp );
if ( (RetValue = Abc_NtkMiterIsConstant(pNtk)) >= 0 )
break;
if ( --Counter == 0 )
break;
*/
/*
Abc_NtkRewrite( pNtk, 0, 0, 0, 0, 0 );
if ( (RetValue = Abc_NtkMiterIsConstant(pNtk)) >= 0 )
break;
if ( --Counter == 0 )
break;
*/
Abc_NtkRewrite( pNtk, 0, 0, 0, 0, 0 );
if ( (RetValue = Abc_NtkMiterIsConstant(pNtk)) >= 0 )
break;
if ( --Counter == 0 )
break;
Abc_NtkRefactor( pNtk, 10, 16, 0, 0, 0, 0 );
if ( (RetValue = Abc_NtkMiterIsConstant(pNtk)) >= 0 )
break;
if ( --Counter == 0 )
break;
pNtk = Abc_NtkBalance( pNtkTemp = pNtk, 0, 0, 0 ); Abc_NtkDelete( pNtkTemp );
if ( (RetValue = Abc_NtkMiterIsConstant(pNtk)) >= 0 )
break;
if ( --Counter == 0 )
break;
}
Abc_NtkMiterPrint( pNtk, "Rewriting ", clk, pParams->fVerbose );
}
if ( pParams->fUseFraiging )
{
// try FRAIGing
clk = clock();
nInspectLimit = pParams->nTotalInspectLimit? pParams->nTotalInspectLimit - pParams->nTotalInspectsMade : 0;
pNtk = Abc_NtkMiterFraig( pNtkTemp = pNtk, (int)(pParams->nFraigingLimitStart * pow(pParams->nFraigingLimitMulti,nIter)), nInspectLimit, &RetValue, &nSatFails, &nSatConfs, &nSatInspects ); Abc_NtkDelete( pNtkTemp );
Abc_NtkMiterPrint( pNtk, "FRAIGing ", clk, pParams->fVerbose );
// printf( "NumFails = %d\n", nSatFails );
if ( RetValue >= 0 )
break;
// add to the number of backtracks and inspects
pParams->nTotalBacktracksMade += nSatConfs;
pParams->nTotalInspectsMade += nSatInspects;
// check if global resource limit is reached
if ( (pParams->nTotalBacktrackLimit && pParams->nTotalBacktracksMade >= pParams->nTotalBacktrackLimit) ||
(pParams->nTotalInspectLimit && pParams->nTotalInspectsMade >= pParams->nTotalInspectLimit) )
{
printf( "Reached global limit on conflicts/inspects. Quitting.\n" );
*ppNtk = pNtk;
return -1;
}
}
}
// try to prove it using brute force SAT
if ( RetValue < 0 && pParams->fUseBdds )
{
if ( pParams->fVerbose )
{
printf( "Attempting BDDs with node limit %d ...\n", pParams->nBddSizeLimit );
fflush( stdout );
}
clk = clock();
pNtk = Abc_NtkCollapse( pNtkTemp = pNtk, pParams->nBddSizeLimit, 0, pParams->fBddReorder, 0 );
if ( pNtk )
{
Abc_NtkDelete( pNtkTemp );
RetValue = ( (Abc_NtkNodeNum(pNtk) == 1) && (Abc_ObjFanin0(Abc_NtkPo(pNtk,0))->pData == Cudd_ReadLogicZero(pNtk->pManFunc)) );
}
else
pNtk = pNtkTemp;
Abc_NtkMiterPrint( pNtk, "BDD building", clk, pParams->fVerbose );
}
if ( RetValue < 0 )
{
if ( pParams->fVerbose )
{
printf( "Attempting SAT with conflict limit %d ...\n", pParams->nMiteringLimitLast );
fflush( stdout );
}
clk = clock();
nInspectLimit = pParams->nTotalInspectLimit? pParams->nTotalInspectLimit - pParams->nTotalInspectsMade : 0;
RetValue = Abc_NtkMiterSat( pNtk, (sint64)pParams->nMiteringLimitLast, (sint64)nInspectLimit, 0, NULL, NULL );
Abc_NtkMiterPrint( pNtk, "SAT solving", clk, pParams->fVerbose );
}
// assign the model if it was proved by rewriting (const 1 miter)
if ( RetValue == 0 && pNtk->pModel == NULL )
{
pNtk->pModel = ALLOC( int, Abc_NtkCiNum(pNtk) );
memset( pNtk->pModel, 0, sizeof(int) * Abc_NtkCiNum(pNtk) );
}
/**Function*************************************************************
Synopsis []
Description []
SideEffects []
SeeAlso []
***********************************************************************/
void Abc_MvRead( Mv_Man_t * p )
{
FILE * pFile;
char Buffer[1000], * pLine;
DdNode * bCube, * bTemp, * bProd, * bVar0, * bVar1, * bCubeSum;
int i, v;
// start the cube
bCubeSum = Cudd_ReadLogicZero(p->dd); Cudd_Ref( bCubeSum );
// start the values
for ( i = 0; i < 15; i++ )
for ( v = 0; v < 4; v++ )
{
p->bValues[i][v] = Cudd_ReadLogicZero(p->dd); Cudd_Ref( p->bValues[i][v] );
p->bValueDcs[i][v] = Cudd_ReadLogicZero(p->dd); Cudd_Ref( p->bValueDcs[i][v] );
}
// read the file
pFile = fopen( "input.pla", "r" );
while ( fgets( Buffer, 1000, pFile ) )
{
if ( Buffer[0] == '#' )
continue;
if ( Buffer[0] == '.' )
{
if ( Buffer[1] == 'e' )
break;
continue;
}
// get the cube
bCube = Abc_MvReadCube( p->dd, Buffer, 18 ); Cudd_Ref( bCube );
// add it to the values of the output functions
pLine = Buffer + 19;
for ( i = 0; i < 15; i++ )
{
if ( pLine[2*i] == '-' && pLine[2*i+1] == '-' )
{
for ( v = 0; v < 4; v++ )
{
p->bValueDcs[i][v] = Cudd_bddOr( p->dd, bTemp = p->bValueDcs[i][v], bCube ); Cudd_Ref( p->bValueDcs[i][v] );
Cudd_RecursiveDeref( p->dd, bTemp );
}
continue;
}
else if ( pLine[2*i] == '0' && pLine[2*i+1] == '0' ) // 0
v = 0;
else if ( pLine[2*i] == '1' && pLine[2*i+1] == '0' ) // 1
v = 1;
else if ( pLine[2*i] == '0' && pLine[2*i+1] == '1' ) // 2
v = 2;
else if ( pLine[2*i] == '1' && pLine[2*i+1] == '1' ) // 3
v = 3;
else assert( 0 );
// add the value
p->bValues[i][v] = Cudd_bddOr( p->dd, bTemp = p->bValues[i][v], bCube ); Cudd_Ref( p->bValues[i][v] );
Cudd_RecursiveDeref( p->dd, bTemp );
}
// add the cube
bCubeSum = Cudd_bddOr( p->dd, bTemp = bCubeSum, bCube ); Cudd_Ref( bCubeSum );
Cudd_RecursiveDeref( p->dd, bTemp );
Cudd_RecursiveDeref( p->dd, bCube );
}
// add the complement of the domain to all values
for ( i = 0; i < 15; i++ )
for ( v = 0; v < 4; v++ )
{
if ( p->bValues[i][v] == Cudd_Not(Cudd_ReadOne(p->dd)) )
continue;
p->bValues[i][v] = Cudd_bddOr( p->dd, bTemp = p->bValues[i][v], p->bValueDcs[i][v] ); Cudd_Ref( p->bValues[i][v] );
Cudd_RecursiveDeref( p->dd, bTemp );
p->bValues[i][v] = Cudd_bddOr( p->dd, bTemp = p->bValues[i][v], Cudd_Not(bCubeSum) ); Cudd_Ref( p->bValues[i][v] );
Cudd_RecursiveDeref( p->dd, bTemp );
}
printf( "Domain = %5.2f %%.\n", 100.0*Cudd_CountMinterm(p->dd, bCubeSum, 32)/Cudd_CountMinterm(p->dd, Cudd_ReadOne(p->dd), 32) );
Cudd_RecursiveDeref( p->dd, bCubeSum );
// create each output function
for ( i = 0; i < 15; i++ )
{
p->bFuncs[i] = Cudd_ReadLogicZero(p->dd); Cudd_Ref( p->bFuncs[i] );
for ( v = 0; v < 4; v++ )
{
bVar0 = Cudd_NotCond( Cudd_bddIthVar(p->dd, 30), ((v & 1) == 0) );
bVar1 = Cudd_NotCond( Cudd_bddIthVar(p->dd, 31), ((v & 2) == 0) );
bCube = Cudd_bddAnd( p->dd, bVar0, bVar1 ); Cudd_Ref( bCube );
bProd = Cudd_bddAnd( p->dd, p->bValues[i][v], bCube ); Cudd_Ref( bProd );
Cudd_RecursiveDeref( p->dd, bCube );
// add the value
p->bFuncs[i] = Cudd_bddOr( p->dd, bTemp = p->bFuncs[i], bProd ); Cudd_Ref( p->bFuncs[i] );
Cudd_RecursiveDeref( p->dd, bTemp );
Cudd_RecursiveDeref( p->dd, bProd );
}
}
}
void
I_BDD_new_mutex_domain (DdManager *dd, DdNode *node[], unsigned int size)
{
unsigned int shbuf, vars, ctr;
DdNode **ddVar, *one, *newnode, *oldnode, *bitnode;
int i,j, extra;
/* base cases */
if (size == 0)
{
return;
}
one = Cudd_ReadOne(dd);
if (size == 1)
{
node[0] = one;
Cudd_Ref (node[0]);
return;
}
for(shbuf = size - 1, vars = 0; shbuf != 0; shbuf >>= 1, vars++);
extra = (1 << vars) - size;
ddVar = malloc(vars * sizeof(DdNode *));
for (j = 0; j < vars; j++)
{
ddVar[j] = bitnode = Cudd_bddNewVar (dd);
Cudd_Ref (bitnode);
}
ctr = 0;
for (i = 0; i < size; i++, ctr++)
{
newnode = one;
Cudd_Ref (newnode);
#ifdef I_BDD_MUTEX_DEBUG
printf("i%d:", i);
#endif
for (j = (extra > 0); j < vars; j++)
{
bitnode = (ctr & (1 << j)) ? ddVar[j] : Cudd_Not(ddVar[j]);
newnode = Cudd_bddAnd (dd, bitnode, oldnode = newnode);
Cudd_Ref (newnode);
Cudd_RecursiveDeref (dd, oldnode);
#ifdef I_BDD_MUTEX_DEBUG
if (ctr & (1 << j))
printf("+ v%d ", j);
else
printf("+ !v%d ", j);
#endif
}
#ifdef I_BDD_MUTEX_DEBUG
printf ("\n");
#endif
node[i] = newnode;
if (extra > 0)
{
extra--;
ctr++;
}
}
#ifdef I_BDD_MUTEX_TEST
{
DdNode *test, *old_test;
int i;
test = Cudd_ReadLogicZero(dd);
Cudd_Ref (test);
for (i = 0; i < size; i++)
{
if (Cudd_bddIte (dd, test, node[i], zero) != zero)
I_punt("I_BDD_new_mutex_complex: not mutex!");
test = Cudd_bddIte(dd, node[i], one, old_test = test);
Cudd_Ref (test);
Cudd_RecursiveDeref (dd, old_test);
}
if (test != one)
I_punt("I_BDD_new_mutex_complex: not exhaustive!");
Cudd_RecursiveDeref (dd, test);
}
#endif
for (j = 0; j < vars; j++)
Cudd_RecursiveDeref (dd, ddVar[j]);
free (ddVar);
return;
}
shadow_mgr new_shadow_mgr(bool do_cudd, bool do_local, bool do_dist, chaining_t chaining) {
if (!(do_cudd || do_local || do_dist)) {
err(true, "Must have at least one active evaluation mode");
}
shadow_mgr mgr = (shadow_mgr) malloc_or_fail(sizeof(shadow_ele),
"new_shadow_mgr");
mgr->do_cudd = do_cudd;
mgr->do_local = do_local;
mgr->do_dist = do_dist;
mgr->c2r_table = word_keyvalue_new();
mgr->r2c_table = word_keyvalue_new();
mgr->nvars = 0;
mgr->nzvars = 0;
ref_t r = REF_ZERO;
DdNode *n = NULL;
if (do_cudd) {
/* Modified CUDD Parameters */
unsigned int numVars = 0; /* Default 0 */
unsigned int numVarsZ = 0; /* Default 0 */
unsigned int numSlots = 1u<<18; /* Default 256 */
unsigned int cacheSize = 1u<<22; /* Default 262144 */
/* Default 67,108,864 */
unsigned long int maxMemory = (1u<<31) + 32UL * 1024 * 1024 * 1024;
#if 0
// Use defaults
numSlots = 256;
cacheSize = 262144;
maxMemory = 67108864;
#endif
mgr->bdd_manager = Cudd_Init(numVars, numVarsZ, numSlots, cacheSize, maxMemory);
Cudd_AutodynDisable(mgr->bdd_manager);
Cudd_AutodynDisableZdd(mgr->bdd_manager);
#ifndef NO_CHAINING
Cudd_ChainingType ct = CUDD_CHAIN_NONE;
switch (chaining) {
case CHAIN_NONE:
ct = CUDD_CHAIN_NONE;
report(0, "No chaining enabled");
break;
case CHAIN_CONSTANT:
ct = CUDD_CHAIN_CONSTANT;
report(0, "Constant chaining enabled");
break;
case CHAIN_ALL:
ct = CUDD_CHAIN_ALL;
report(0, "Or chaining enabled");
break;
default:
err(true, "Invalid chaining mode %d\n", chaining);
}
Cudd_SetChaining(mgr->bdd_manager, ct);
#endif
n = Cudd_ReadLogicZero(mgr->bdd_manager);
report(0, "Using CUDD Version %s", CUDD_VERSION);
}
if (do_ref(mgr)) {
mgr->ref_mgr = new_ref_mgr();
if (!do_cudd) {
n = ref2dd(mgr, r);
}
} else {
r = dd2ref(n, IS_BDD);
}
reference_dd(mgr, n);
add_ref(mgr, r, n);
return mgr;
}
static DdNode *tree2BDD(DdManager *mgr, ACCExpr *expr, VarMap &varMap)
{
std::string op = expr->value;
DdNode *ret = nullptr;
if (op == "&&")
op = "&";
else if (op == "||")
op = "|";
if (checkInteger(expr, "1"))
ret = Cudd_ReadOne(mgr);
else if (checkInteger(expr, "0"))
ret = Cudd_ReadLogicZero(mgr);
else if (op == "!")
return Cudd_Not(tree2BDD(mgr, expr->operands.front(), varMap)); // Not passes through ref count
else if (op != "&" && op != "|" && op != "^") {
if ((op == "!=" || op == "==")) {
ACCExpr *lhs = getRHS(expr, 0);
if (boolPossible(lhs) && boolPossible(getRHS(expr,1)))
goto next; // we can analyze relops on booleans
if (trace_bool)
printf("[%s:%d] boolnot %d %d = %s\n", __FUNCTION__, __LINE__, boolPossible(getRHS(expr,0)), boolPossible(getRHS(expr,1)), tree2str(expr).c_str());
if (isIdChar(lhs->value[0])) {
if (trace_bool)
printf("[%s:%d] name %s type %s\n", __FUNCTION__, __LINE__, lhs->value.c_str(), refList[lhs->value].type.c_str());
}
}
if (op == "!=") // normalize comparison strings
expr->value = "==";
std::string name = "( " + tree2str(expr) + " )";
if (!varMap[name]) {
varMap[name] = new MapItem;
varMap[name]->index = varMap.size();
varMap[name]->node = Cudd_bddIthVar(mgr, varMap[name]->index);
}
ret = varMap[name]->node;
if (op == "!=") { // normalize comparison strings
expr->value = op; // restore
ret = Cudd_Not(ret);
}
}
if (ret) {
Cudd_Ref(ret);
return ret;
}
next:;
for (auto item: expr->operands) {
DdNode *operand = tree2BDD(mgr, item, varMap), *next;
if (!ret)
ret = operand;
else {
if (op == "&")
next = Cudd_bddAnd(mgr, ret, operand);
else if (op == "|")
next = Cudd_bddOr(mgr, ret, operand);
else if (op == "^" || op == "!=")
next = Cudd_bddXor(mgr, ret, operand);
else if (op == "==")
next = Cudd_bddXnor(mgr, ret, operand);
else {
printf("[%s:%d] unknown operator\n", __FUNCTION__, __LINE__);
exit(-1);
}
Cudd_Ref(next);
Cudd_RecursiveDeref(mgr, operand);
Cudd_RecursiveDeref(mgr, ret);
ret = next;
}
}
return ret;
}
