Skip to content

A
abc
Commit 9e668f1b by Alan Mishchenko
Options

Synthesis for mesh of LUTs.

parent d66ff2cf
Showing with 525 additions and 0 deletions
abclib.dsp
...@@ -1963,6 +1963,10 @@ SOURCE=.\src\sat\bmc\bmcMaxi.c ...@@ -1963,6 +1963,10 @@ SOURCE=.\src\sat\bmc\bmcMaxi.c
# End Source File # End Source File
# Begin Source File # Begin Source File
SOURCE=.\src\sat\bmc\bmcMesh.c
# End Source File
# Begin Source File
SOURCE=.\src\sat\bmc\bmcMulti.c SOURCE=.\src\sat\bmc\bmcMulti.c
# End Source File # End Source File
# Begin Source File # Begin Source File
......
src/base/abci/abc.c
...@@ -457,6 +457,7 @@ static int Abc_CommandAbc9ReachN ( Abc_Frame_t * pAbc, int argc, cha ...@@ -457,6 +457,7 @@ static int Abc_CommandAbc9ReachN ( Abc_Frame_t * pAbc, int argc, cha
static int Abc_CommandAbc9ReachY ( Abc_Frame_t * pAbc, int argc, char ** argv ); static int Abc_CommandAbc9ReachY ( Abc_Frame_t * pAbc, int argc, char ** argv );
#endif #endif
static int Abc_CommandAbc9Undo ( Abc_Frame_t * pAbc, int argc, char ** argv ); static int Abc_CommandAbc9Undo ( Abc_Frame_t * pAbc, int argc, char ** argv );
static int Abc_CommandAbc9Mesh ( Abc_Frame_t * pAbc, int argc, char ** argv );
static int Abc_CommandAbc9Iso ( Abc_Frame_t * pAbc, int argc, char ** argv ); static int Abc_CommandAbc9Iso ( Abc_Frame_t * pAbc, int argc, char ** argv );
static int Abc_CommandAbc9IsoNpn ( Abc_Frame_t * pAbc, int argc, char ** argv ); static int Abc_CommandAbc9IsoNpn ( Abc_Frame_t * pAbc, int argc, char ** argv );
static int Abc_CommandAbc9IsoSt ( Abc_Frame_t * pAbc, int argc, char ** argv ); static int Abc_CommandAbc9IsoSt ( Abc_Frame_t * pAbc, int argc, char ** argv );
...@@ -1109,6 +1110,7 @@ void Abc_Init( Abc_Frame_t * pAbc ) ...@@ -1109,6 +1110,7 @@ void Abc_Init( Abc_Frame_t * pAbc )
Cmd_CommandAdd( pAbc, "ABC9", "&reachy", Abc_CommandAbc9ReachY, 0 ); Cmd_CommandAdd( pAbc, "ABC9", "&reachy", Abc_CommandAbc9ReachY, 0 );
#endif #endif
Cmd_CommandAdd( pAbc, "ABC9", "&undo", Abc_CommandAbc9Undo, 0 ); Cmd_CommandAdd( pAbc, "ABC9", "&undo", Abc_CommandAbc9Undo, 0 );
Cmd_CommandAdd( pAbc, "ABC9", "&mesh", Abc_CommandAbc9Mesh, 0 );
Cmd_CommandAdd( pAbc, "ABC9", "&iso", Abc_CommandAbc9Iso, 0 ); Cmd_CommandAdd( pAbc, "ABC9", "&iso", Abc_CommandAbc9Iso, 0 );
Cmd_CommandAdd( pAbc, "ABC9", "&isonpn", Abc_CommandAbc9IsoNpn, 0 ); Cmd_CommandAdd( pAbc, "ABC9", "&isonpn", Abc_CommandAbc9IsoNpn, 0 );
Cmd_CommandAdd( pAbc, "ABC9", "&isost", Abc_CommandAbc9IsoSt, 0 ); Cmd_CommandAdd( pAbc, "ABC9", "&isost", Abc_CommandAbc9IsoSt, 0 );
...@@ -38031,6 +38033,106 @@ usage: ...@@ -38031,6 +38033,106 @@ usage:
SeeAlso [] SeeAlso []
***********************************************************************/ ***********************************************************************/
int Abc_CommandAbc9Mesh( Abc_Frame_t * pAbc, int argc, char ** argv )
{
extern void Bmc_MeshTest( Gia_Man_t * p, int X, int Y, int T, int fVerbose );
int X = 4;
int Y = 4;
int T = 3;
int c, fVerbose = 1;
// set defaults
Extra_UtilGetoptReset();
while ( ( c = Extra_UtilGetopt( argc, argv, "XYTh" ) ) != EOF )
{
switch ( c )
{
case 'X':
if ( globalUtilOptind >= argc )
{
Abc_Print( -1, "Command line switch \"-X\" should be followed by an integer.\n" );
goto usage;
}
X = atoi(argv[globalUtilOptind]);
globalUtilOptind++;
if ( X < 3 )
goto usage;
break;
case 'Y':
if ( globalUtilOptind >= argc )
{
Abc_Print( -1, "Command line switch \"-Y\" should be followed by an integer.\n" );
goto usage;
}
Y = atoi(argv[globalUtilOptind]);
globalUtilOptind++;
if ( Y < 3 )
goto usage;
break;
case 'T':
if ( globalUtilOptind >= argc )
{
Abc_Print( -1, "Command line switch \"-T\" should be followed by an integer.\n" );
goto usage;
}
T = atoi(argv[globalUtilOptind]);
globalUtilOptind++;
if ( T < 2 )
goto usage;
break;
case 'v':
fVerbose ^= 1;
break;
case 'h':
goto usage;
default:
goto usage;
}
}
if ( pAbc->pGia == NULL )
{
Abc_Print( -1, "Abc_CommandAbc9Mesh(): There is no design.\n" );
return 1;
}
if ( Gia_ManCoNum(pAbc->pGia) != 1 )
{
Abc_Print( -1, "Currently this command expects AIG with one output.\n" );
return 1;
}
if ( Gia_ManCiNum(pAbc->pGia) > 20 )
{
Abc_Print( -1, "Currently this command expects AIG with no more than 20 nodes.\n" );
return 1;
}
if ( Gia_ManLevelNum(pAbc->pGia) > T )
{
Abc_Print( -1, "The depth of the AIG (%d) cannot be larger than the latency (%d).\n", Gia_ManLevelNum(pAbc->pGia), T );
return 1;
}
Bmc_MeshTest( pAbc->pGia, X, Y, T, fVerbose );
return 0;
usage:
Abc_Print( -2, "usage: &mesh [-XYT num] [-h]\n" );
Abc_Print( -2, "\t creates a mesh architecture for the given AIG\n" );
Abc_Print( -2, "\t-X num : horizontal size of the mesh (X >= 3) [default = %d]\n", X );
Abc_Print( -2, "\t-Y num : vertical size of the mesh (Y >= 3) [default = %d]\n", Y );
Abc_Print( -2, "\t-T num : the latency of the mesh (T >= 2) [default = %d]\n", T );
Abc_Print( -2, "\t-v : toggle printing verbose information [default = %s]\n", fVerbose? "yes": "no" );
Abc_Print( -2, "\t-h : print the command usage\n");
return 1;
}
/**Function*************************************************************
Synopsis []
Description []
SideEffects []
SeeAlso []
***********************************************************************/
int Abc_CommandAbc9Iso( Abc_Frame_t * pAbc, int argc, char ** argv ) int Abc_CommandAbc9Iso( Abc_Frame_t * pAbc, int argc, char ** argv )
{ {
Gia_Man_t * pAig; Gia_Man_t * pAig;
src/sat/bmc/bmcMesh.c 0 → 100644
/**CFile****************************************************************
FileName [bmcMesh.c]
SystemName [ABC: Logic synthesis and verification system.]
PackageName [SAT-based bounded model checking.]
Synopsis [Synthesis for mesh of LUTs.]
Author [Alan Mishchenko]
Affiliation [UC Berkeley]
Date [Ver. 1.0. Started - June 20, 2005.]
Revision [$Id: bmcMesh.c,v 1.00 2005/06/20 00:00:00 alanmi Exp $]
***********************************************************************/
#include "bmc.h"
#include "sat/satoko/satoko.h"
#include "sat/satoko/solver.h"
ABC_NAMESPACE_IMPL_START
////////////////////////////////////////////////////////////////////////
/// DECLARATIONS ///
////////////////////////////////////////////////////////////////////////
#define NCPARS 16
inline int Bmc_MeshTVar( int Me[102][102], int x, int y ) { return Me[x][y]; }
inline int Bmc_MeshGVar( int Me[102][102], int x, int y ) { return Me[x][y] + Me[101][100]; }
inline int Bmc_MeshCVar( int Me[102][102], int x, int y ) { return Me[x][y] + Me[101][100] + Me[101][101]; }
inline int Bmc_MeshUVar( int Me[102][102], int x, int y ) { return Me[x][y] + Me[101][100] + Me[101][101] + NCPARS; }
////////////////////////////////////////////////////////////////////////
/// FUNCTION DEFINITIONS ///
////////////////////////////////////////////////////////////////////////
/**Function*************************************************************
Synopsis []
Description []
SideEffects []
SeeAlso []
***********************************************************************/
static inline Bmc_MeshVarValue( satoko_t * p, int v )
{
int value = var_value(p, v) != VAR_UNASSING ? var_value(p, v) : var_polarity(p, v);
return value == LIT_TRUE;
}
/**Function*************************************************************
Synopsis []
Description []
SideEffects []
SeeAlso []
***********************************************************************/
int Bmc_MeshAddOneHotness( satoko_t * pSat, int iFirst, int iLast )
{
int i, j, v, pVars[100], nVars = 0, nCount = 0;
assert( iFirst < iLast && iFirst + 110 > iLast );
for ( v = iFirst; v < iLast; v++ )
if ( Bmc_MeshVarValue(pSat, v) ) // v = 1
{
assert( nVars < 100 );
pVars[ nVars++ ] = v;
}
if ( nVars <= 1 )
return 0;
for ( i = 0; i < nVars; i++ )
for ( j = i+1; j < nVars; j++ )
{
int pLits[2], RetValue;
pLits[0] = Abc_Var2Lit( pVars[i], 1 );
pLits[1] = Abc_Var2Lit( pVars[j], 1 );
RetValue = satoko_add_clause( pSat, pLits, 2 ); assert( RetValue );
nCount++;
}
return nCount;
}
/**Function*************************************************************
Synopsis []
Description []
SideEffects []
SeeAlso []
***********************************************************************/
void Bmc_MeshTest( Gia_Man_t * p, int X, int Y, int T, int fVerbose )
{
abctime clk = Abc_Clock();
satoko_t * pSat = satoko_create();
Gia_Obj_t * pObj;
int Me[102][102] = {{0}};
int pN[102][2] = {{0}};
int I = Gia_ManPiNum(p);
int G = I + Gia_ManAndNum(p);
int i, x, y, t, g, c, status, RetValue, Lit, iVar, nClauses = 0;
assert( X <= 100 && Y <= 100 && T <= 100 && G <= 100 );
// init the graph
for ( i = 0; i < I; i++ )
pN[i][0] = pN[i][1] = -1;
Gia_ManForEachAnd( p, pObj, i )
{
pN[i-1][0] = Gia_ObjFaninId0(pObj, i)-1;
pN[i-1][1] = Gia_ObjFaninId1(pObj, i)-1;
}
if ( fVerbose )
{
printf( "The graph has %d inputs: ", Gia_ManPiNum(p) );
for ( i = 0; i < I; i++ )
printf( "%c ", 'a' + i );
printf( " and %d nodes: ", Gia_ManAndNum(p) );
for ( i = I; i < G; i++ )
printf( "%c=%c%c ", 'a' + i, 'a' + pN[i][0] , 'a' + pN[i][1] );
printf( "\n" );
}
// init SAT variables (time vars + graph vars + config vars)
// config variables: 16 = 4 buff vars + 12 node vars
iVar = 0;
for ( y = 0; y < Y; y++ )
for ( x = 0; x < X; x++ )
{
printf( "%3d %3d %3d %s", iVar, iVar+T, iVar+T+G, x == X-1 ? "\n":"" );
Me[x][y] = iVar;
iVar += T + G + NCPARS + 1;
}
Me[101][100] = T;
Me[101][101] = G;
if ( fVerbose )
printf( "SAT variable count is %d (%d time vars + %d graph vars + %d config vars + %d aux vars)\n", iVar, X*Y*T, X*Y*G, X*Y*NCPARS, X*Y );
// add constraints
// time 0 and primary inputs only on the boundary
for ( x = 0; x < X; x++ )
for ( y = 0; y < Y; y++ )
{
int iTVar = Bmc_MeshTVar( Me, x, y );
int iGVar = Bmc_MeshGVar( Me, x, y );
if ( x == 0 || x == X-1 || y == 0 || y == Y-1 ) // boundary
{
// time 0 is required
for ( t = 0; t < T; t++ )
{
Lit = Abc_Var2Lit( iTVar+t, (int)(t > 0) );
RetValue = satoko_add_clause( pSat, &Lit, 1 ); assert( RetValue );
}
// internal nodes are not allowed
for ( g = I; g < G; g++ )
{
Lit = Abc_Var2Lit( iGVar+g, 1 );
RetValue = satoko_add_clause( pSat, &Lit, 1 ); assert( RetValue );
}
}
else // not a boundary
{
Lit = Abc_Var2Lit( iTVar, 1 ); // cannot have time 0
RetValue = satoko_add_clause( pSat, &Lit, 1 ); assert( RetValue );
}
}
for ( x = 1; x < X-1; x++ )
for ( y = 1; y < Y-1; y++ )
{
int pLits[100], nLits;
int iTVar = Bmc_MeshTVar( Me, x, y );
int iGVar = Bmc_MeshGVar( Me, x, y );
int iCVar = Bmc_MeshCVar( Me, x, y );
int iUVar = Bmc_MeshUVar( Me, x, y );
// 0=left 1=up 2=right 3=down
int iTVars[4];
int iGVars[4];
iTVars[0] = Bmc_MeshTVar( Me, x-1, y );
iGVars[0] = Bmc_MeshGVar( Me, x-1, y );
iTVars[1] = Bmc_MeshTVar( Me, x, y-1 );
iGVars[1] = Bmc_MeshGVar( Me, x, y-1 );
iTVars[2] = Bmc_MeshTVar( Me, x+1, y );
iGVars[2] = Bmc_MeshGVar( Me, x+1, y );
iTVars[3] = Bmc_MeshTVar( Me, x, y+1 );
iGVars[3] = Bmc_MeshGVar( Me, x, y+1 );
// condition when cell is used
for ( g = 0; g < G; g++ )
{
pLits[0] = Abc_Var2Lit( iGVar+g, 1 );
pLits[1] = Abc_Var2Lit( iUVar, 0 );
RetValue = satoko_add_clause( pSat, pLits, 2 ); assert( RetValue );
}
// at least one time is used
pLits[0] = Abc_Var2Lit( iUVar, 1 );
for ( t = 1; t < T; t++ )
pLits[t] = Abc_Var2Lit( iTVar+t, 0 );
RetValue = satoko_add_clause( pSat, pLits, T ); assert( RetValue );
// at least one config is used
pLits[0] = Abc_Var2Lit( iUVar, 1 );
for ( c = 0; c < NCPARS; c++ )
pLits[c+1] = Abc_Var2Lit( iCVar+c, 0 );
RetValue = satoko_add_clause( pSat, pLits, NCPARS+1 ); assert( RetValue );
// constraints for each time
for ( t = 1; t < T; t++ )
{
int Conf[12][2] = {{0, 1}, {0, 2}, {0, 3}, {1, 2}, {1, 3}, {2, 3}, {1, 0}, {2, 0}, {3, 0}, {2, 1}, {3, 1}, {3, 2}};
// buffer
for ( g = 0; g < G; g++ )
for ( c = 0; c < 4; c++ )
{
nLits = 0;
pLits[ nLits++ ] = Abc_Var2Lit( iTVar+t, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iGVar+g, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iCVar+c, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iTVars[c]+t-1, 0 );
RetValue = satoko_add_clause( pSat, pLits, nLits ); assert( RetValue );
nLits = 0;
pLits[ nLits++ ] = Abc_Var2Lit( iTVar+t, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iGVar+g, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iCVar+c, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iGVars[c]+g, 0 );
RetValue = satoko_add_clause( pSat, pLits, nLits ); assert( RetValue );
nClauses += 2;
}
// node
for ( g = I; g < G; g++ )
for ( c = 0; c < 12; c++ )
{
assert( pN[g][0] >= 0 && pN[g][1] >= 0 );
nLits = 0;
pLits[ nLits++ ] = Abc_Var2Lit( iTVar+t, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iGVar+g, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iCVar+c+4, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iTVars[Conf[c][0]]+t-1, 0 );
RetValue = satoko_add_clause( pSat, pLits, nLits ); assert( RetValue );
nLits = 0;
pLits[ nLits++ ] = Abc_Var2Lit( iTVar+t, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iGVar+g, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iCVar+c+4, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iTVars[Conf[c][1]]+t-1, 0 );
RetValue = satoko_add_clause( pSat, pLits, nLits ); assert( RetValue );
nLits = 0;
pLits[ nLits++ ] = Abc_Var2Lit( iTVar+t, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iGVar+g, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iCVar+c+4, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iGVars[Conf[c][0]]+pN[g][0], 0 );
RetValue = satoko_add_clause( pSat, pLits, nLits ); assert( RetValue );
nLits = 0;
pLits[ nLits++ ] = Abc_Var2Lit( iTVar+t, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iGVar+g, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iCVar+c+4, 1 );
pLits[ nLits++ ] = Abc_Var2Lit( iGVars[Conf[c][1]]+pN[g][1], 0 );
RetValue = satoko_add_clause( pSat, pLits, nLits ); assert( RetValue );
nClauses += 4;
}
}
}
// final condition
{
int iGVar = Bmc_MeshGVar( Me, 1, 1 ) + G-1;
Lit = Abc_Var2Lit( iGVar, 0 );
RetValue = satoko_add_clause( pSat, &Lit, 1 );
if ( RetValue == 0 )
{
printf( "Problem has no solution. " );
Abc_PrintTime( 1, "Time", Abc_Clock() - clk );
satoko_destroy( pSat );
return;
}
}
if ( fVerbose )
printf( "Finished adding %d clauses. Started solving...\n", nClauses );
while ( 1 )
{
int nAddClauses = 0;
status = satoko_solve( pSat );
if ( status == SATOKO_UNSAT )
{
printf( "Problem has no solution. " );
break;
}
if ( status == SATOKO_UNDEC )
{
printf( "Computation timed out. " );
break;
}
assert( status == SATOKO_SAT );
// check if the solution is valid and add constraints
for ( x = 0; x < X; x++ )
for ( y = 0; y < Y; y++ )
{
if ( x == 0 || x == X-1 || y == 0 || y == Y-1 ) // boundary
{
int iGVar = Bmc_MeshGVar( Me, x, y );
nAddClauses += Bmc_MeshAddOneHotness( pSat, iGVar, iGVar + G );
}
else
{
int iTVar = Bmc_MeshTVar( Me, x, y );
int iGVar = Bmc_MeshGVar( Me, x, y );
int iCVar = Bmc_MeshCVar( Me, x, y );
nAddClauses += Bmc_MeshAddOneHotness( pSat, iTVar, iTVar + T );
nAddClauses += Bmc_MeshAddOneHotness( pSat, iGVar, iGVar + G );
nAddClauses += Bmc_MeshAddOneHotness( pSat, iCVar, iCVar + NCPARS );
}
}
if ( nAddClauses > 0 )
{
printf( "Adding %d one-hotness clauses.\n", nAddClauses );
continue;
}
printf( "Satisfying solution found. " );
iVar = solver_varnum(pSat);
for ( i = 0; i < iVar; i++ )
if ( Bmc_MeshVarValue(pSat, i) )
printf( "%d ", i );
printf( "\n" );
break;
}
Abc_PrintTime( 1, "Time", Abc_Clock() - clk );
if ( status == SATOKO_SAT )
{
// count the number of nodes and buffers
int nBuffs = 0, nNodes = 0;
for ( x = 1; x < X-1; x++ )
for ( y = 1; y < Y-1; y++ )
{
int iCVar = Bmc_MeshCVar( Me, x, y );
for ( c = 0; c < 4; c++ )
if ( Bmc_MeshVarValue(pSat, iCVar+c) )
{
printf( "Buffer x=%d y=%d (var = %d; config = %d)\n", x, y, iCVar+c, c );
nBuffs++;
}
for ( c = 4; c < NCPARS; c++ )
if ( Bmc_MeshVarValue(pSat, iCVar+c) )
{
printf( "Node x=%d y=%d (var = %d; config = %d)\n", x, y, iCVar+c, c );
nNodes++;
}
}
printf( "The %d x %d mesh with latency %d with %d active cells (%d nodes and %d buffers):\n", X, Y, T, nNodes+nBuffs, nNodes, nBuffs );
// print mesh
printf( " Y\\X " );
for ( x = 0; x < X; x++ )
printf( " %-2d ", x );
printf( "\n" );
for ( y = 0; y < Y; y++ )
{
printf( " %-2d ", y );
for ( x = 0; x < X; x++ )
{
int iTVar = Bmc_MeshTVar( Me, x, y );
int iGVar = Bmc_MeshGVar( Me, x, y );
int fFound = 0; ;
for ( t = 0; t < T; t++ )
for ( g = 0; g < G; g++ )
if ( Bmc_MeshVarValue(pSat, iTVar+t) && Bmc_MeshVarValue(pSat, iGVar+g) )
{
printf( " %c%-2d ", 'a' + g, t );
fFound = 1;
}
if ( !fFound )
printf( " " );
}
printf( "\n" );
}
}
satoko_destroy( pSat );
}
////////////////////////////////////////////////////////////////////////
/// END OF FILE ///
////////////////////////////////////////////////////////////////////////
ABC_NAMESPACE_IMPL_END
src/sat/bmc/module.make
...@@ -21,5 +21,6 @@ SRC += src/sat/bmc/bmcBCore.c \ ...@@ -21,5 +21,6 @@ SRC += src/sat/bmc/bmcBCore.c \
src/sat/bmc/bmcInse.c \ src/sat/bmc/bmcInse.c \
src/sat/bmc/bmcLoad.c \ src/sat/bmc/bmcLoad.c \
src/sat/bmc/bmcMaxi.c \ src/sat/bmc/bmcMaxi.c \
src/sat/bmc/bmcMesh.c \
src/sat/bmc/bmcMulti.c \ src/sat/bmc/bmcMulti.c \
src/sat/bmc/bmcUnroll.c src/sat/bmc/bmcUnroll.c
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment