/**
* Constructs a post through zero tracks for backwards procedure which is
* a little bit different to handle.
*
* @param aut: base automaton
* @param state: initial state we are getting zero post for
* @param level: level of macro inception
* @param prefix: list of variables for projection
* @return zero post of initial @p state
*/
MacroStateSet* GetZeroMacroPost(Automaton & aut, TStateSet*& state, unsigned level, PrefixListType & prefix) {
if(level == 0) {
return GetZeroPost(aut, state, level, prefix);
} else {
if(level == 1 && (((MacroStateSet*)state)->getMacroStates()).size() == 0) {
return new MacroStateSet();
} else {
const MacroTransMTBDD & transPost = GetMTBDDForPost(aut, state, level, prefix);
int projecting = getProjectionVariable(level, prefix);
#if USE_PRUNED_UNION_FUNCTOR
MacroPrunedUnionFunctor muf(level-1);
#else
MacroUnionFunctor muf;
#endif
MacroStateDeterminizatorFunctor msdf;
MacroTransMTBDD projectedMtbdd = (msdf(transPost)).Project(
[&transPost, projecting](size_t var) { return var < projecting;}, muf);
#if (DEBUG_BDDS == true)
std::cout << "[?] BDD: \n";
std::cout << MacroTransMTBDD::DumpToDot({&projectedMtbdd}) << "\n\n";
#endif
MacroStateSet *postStates = projectedMtbdd.GetValue(constructUniversalTrack());
return postStates;
}
}
}
int changepars(void)
{
extern int use_rxvt;
extern int cluster;
extern int shortqsonr;
extern int searchflg;
extern int demode;
extern int contest;
extern int announcefilter;
extern int showscore_flag;
extern int zonedisplay;
extern int trxmode;
extern char hiscall[];
extern int rit;
extern int trx_control;
extern int editor;
extern int packetinterface;
extern int nopacket;
extern int cqdelay;
extern int ctcomp;
extern SCREEN *mainscreen;
extern char *config_file;
extern int miniterm;
extern char buffer[];
#ifdef HAVE_LIBHAMLIB
extern freq_t outfreq;
#else
extern int outfreq;
#endif
extern int simulator;
extern int keyerport;
extern char synclogfile[];
extern char sc_volume[];
extern int cwstart;
char parameterstring[20];
char parameters[51][19];
char cmdstring[80];
int i, k, x, nopar = 0;
int maxpar = 50;
int volumebuffer;
int currentmode = 0;
int rc;
strcpy(parameters[0], "SPOT");
strcpy(parameters[1], "MAP");
strcpy(parameters[2], "CLOFF");
strcpy(parameters[3], "CLUSTER");
strcpy(parameters[4], "SHORT");
strcpy(parameters[5], "LONG");
strcpy(parameters[6], "MESSAGE");
strcpy(parameters[7], "LIST");
strcpy(parameters[8], "CHECK");
strcpy(parameters[9], "NOCHECK");
strcpy(parameters[10], "TONE");
strcpy(parameters[11], "EDIT");
strcpy(parameters[12], "VIEW");
strcpy(parameters[13], "HELP");
strcpy(parameters[14], "DEMODE");
strcpy(parameters[15], "CONTEST");
strcpy(parameters[16], "FILTER");
strcpy(parameters[17], "SCORE");
strcpy(parameters[18], "WRITE");
strcpy(parameters[19], "EXIT");
strcpy(parameters[20], "TXFILE");
strcpy(parameters[21], "ZONES");
strcpy(parameters[22], "CTY");
strcpy(parameters[23], "MODE");
strcpy(parameters[24], "SET");
strcpy(parameters[25], "MULTI");
strcpy(parameters[26], "PROP");
strcpy(parameters[27], "RITCLEAR");
strcpy(parameters[28], "TRXCONTROL");
strcpy(parameters[29], "CFG");
//strcpy(parameters[30], "CWMODE");
strcpy(parameters[30], "CW");
strcpy(parameters[31], "SSBMODE");
strcpy(parameters[32], "DIGIMODE");
strcpy(parameters[33], "PACKET");
strcpy(parameters[34], "SIMULATOR");
strcpy(parameters[35], "INFO");
strcpy(parameters[36], "FREQ");
strcpy(parameters[37], "RECONNECT");
strcpy(parameters[38], "QUIT");
strcpy(parameters[39], "CQDELAY");
strcpy(parameters[40], "ADIF");
strcpy(parameters[41], "SYNC");
strcpy(parameters[42], "RESCORE");
strcpy(parameters[43], "SCVOLUME");
strcpy(parameters[44], "SCAN");
strcpy(parameters[45], "DEBUG");
strcpy(parameters[46], "MINITERM");
strcpy(parameters[47], "RTTY");
strcpy(parameters[48], "SOUND");
strcpy(parameters[49], "CWMODE");
strcpy(parameters[50], "CHARS");
nopar = 0;
attroff(A_STANDOUT);
attron(COLOR_PAIR(C_HEADER));
mvprintw(12, 29, "PARAMETER? ");
refreshp();
mvprintw(12, 29, " ");
mvprintw(12, 29, "");
refreshp();
echo();
getstr(parameterstring);
noecho();
for (k = 0; parameterstring[k]; k++)
parameterstring[k] = toupper(parameterstring[k]);
for (i = 0; i <= maxpar; i++) {
if (strncmp(parameterstring, parameters[i], 3) == 0) {
break;
}
}
switch (i) {
case 0: /* SPOTS) */
{
/* SPOTS not supported anymore
* - default to MAP*/
cluster = MAP;
break;
}
case 1: /* BANDMAP */
{
cluster = MAP;
break;
}
case 2: /* CLOFF */
{
cluster = NOCLUSTER;
break;
}
case 3: /* CLUSTER */
{
cluster = CLUSTER;
announcefilter = FILTER_ALL;
break;
}
case 4: /* SHORTNR */
{
shortqsonr = SHORTCW;
break;
}
case 5: /* LONGNR */
{
shortqsonr = LONGCW;
break;
}
case 6: /* MESSAGE */
{
message_change(i);
break;
}
case 7: /* LIST */
{
listmessages();
break;
}
case 8: /* CHECK */
{
searchflg = SEARCHWINDOW;
break;
}
case 9: /* NOCHECK */
{
searchflg = 0;
break;
}
case 10: /* TONE */
{
set_tone();
break;
}
case 11: /* EDIT */
{
logedit();
break;
}
case 12: /* VIEW */
{
logview();
break;
}
case 13: /* HELP */
{
show_help();
break;
}
case 14: /* DEMODE */
{
if (demode == SEND_DE)
demode = 0;
else
demode = SEND_DE;
mvprintw(13, 29, "DE-mode is %d", demode);
refreshp();
sleep(1);
break;
}
case 15: /* CONTEST */
{
if (contest == CONTEST)
contest = 0;
else {
contest = CONTEST;
searchflg = SEARCHWINDOW;
}
mvprintw(13, 29, "CONTEST-mode is %d", contest);
refreshp();
sleep(1);
break;
}
case 16: /* FILTER */
{
announcefilter++;
if (announcefilter > 3)
announcefilter = 0;
mvprintw(13, 29, "FILTER-mode is %d", announcefilter);
refreshp();
sleep(1);
break;
}
case 17: /* SCORE */
{
if (showscore_flag == 0)
showscore_flag = 1;
else {
showscore_flag = 0;
}
mvprintw(13, 29, "Show score-mode is %d", showscore_flag);
refreshp();
sleep(1);
break;
}
case 18: /* WRITE CABRILLO FILE */
{
int old_cluster = cluster;
cluster = NOCLUSTER;
write_cabrillo();
cluster = old_cluster;
break;
}
case 19: /* EXIT */
{
writeparas();
clear();
cleanup_telnet();
endwin();
puts("\n\nThanks for using TLF.. 73\n");
exit(0);
break;
}
case 20: /* TXFILE */
{
break;
}
case 21: /* ZONES */
{
if (zonedisplay == 0)
zonedisplay = 1;
else {
zonedisplay = 0;
}
break;
}
case 22: /* COUNTRIES */
{
show_mults();
refreshp();
sleep(1);
break;
}
case 23: /* MODE */
{
if (trxmode == CWMODE)
trxmode = SSBMODE;
else if (trxmode == SSBMODE)
trxmode = DIGIMODE;
else
trxmode = CWMODE;
if (trxmode == CWMODE) {
mvprintw(13, 29, "TRXMODE = CW");
} else if (trxmode == SSBMODE)
mvprintw(13, 29, "TRXMODE = SSB");
else
mvprintw(13, 29, "TRXMODE = DIG");
refreshp();
sleep(1);
break;
}
case 24: /* SET PARAMETERS */
case 29: /* CFG PARAMETERS */
{
clear();
if (editor == EDITOR_JOE) {
strcpy(cmdstring, "joe ");
} else if (editor == EDITOR_VI) {
strcpy(cmdstring, "vi ");
} else if (editor == EDITOR_MC) {
strcpy(cmdstring, "mcedit ");
} else {
strcpy(cmdstring, "e3 ");
}
strcat(cmdstring, config_file);
rc = system(cmdstring);
read_logcfg();
read_rules(); /* also reread rules file */
writeparas();
mvprintw(24, 0, "Logcfg.dat loaded, parameters written..");
refreshp();
clear_display();
break;
}
case 25: /* MULTI */
{
multiplierinfo();
break;
}
case 26: /* PROPAGATION */
{
muf();
clear_display();
break;
}
case 27: /* RITCLEAR */
{
if (rit == RITCLEAR)
rit = 0;
else {
rit = RITCLEAR;
}
if (rit == RITCLEAR) {
mvprintw(13, 29, "RIT clear on");
} else {
mvprintw(13, 29, "RIT clear off");
}
refreshp();
sleep(1);
break;
}
case 28: /* trx ctl */
{
if (trx_control == 1)
trx_control = 0;
else {
trx_control = 1;
}
if (trx_control == 1) {
mvprintw(13, 29, "TRX control on");
} else {
mvprintw(13, 29, "TRX control off");
}
refreshp();
sleep(1);
break;
}
case 30: /* CW */
case 49:
{
if (keyerport == MFJ1278_KEYER) {
strcpy(buffer, "MODE CW,30");
buffer[7] = '\015';
buffer[8] = 'K';
buffer[9] = '\015';
buffer[10] = '\0';
sendbuf();
}
trxmode = CWMODE;
if (trx_control == 1)
outfreq = SETCWMODE;
break;
}
case 31: /* SSBMODE */
{
trxmode = SSBMODE;
outfreq = SETSSBMODE;
break;
}
case 32: /* DIGIMODE */
{
trxmode = DIGIMODE;
break;
}
case 33: /* PACKET */
{
if ((nopacket == 0) && (packetinterface > 0))
packet();
break;
}
case 34: /* SIMULATOR */
{
if (simulator == 0) {
simulator = 1;
if (ctcomp == 1) {
mvprintw(13, 19,
"The simulator only works in TRmode. Switching to TRmode");
ctcomp = 0;
} else
mvprintw(13, 29, "Simulator on");
refreshp();
if (keyerport == NET_KEYER) {
if (netkeyer(K_WORDMODE, NULL) < 0) {
mvprintw(24, 0,
"keyer not active; switching to SSB");
trxmode = SSBMODE;
clear_display();
}
}
} else {
simulator = 0;
mvprintw(13, 29, "Simulator off");
refreshp();
if (keyerport == NET_KEYER) {
if (netkeyer(K_RESET, NULL) < 0) {
mvprintw(24, 0,
"keyer not active; switching to SSB");
trxmode = SSBMODE;
clear_display();
}
}
}
break;
}
case 35: /* INFO */
{
int currentterm = miniterm;
miniterm = 0;
networkinfo();
miniterm = currentterm;
if (currentmode == DIGIMODE)
trxmode = DIGIMODE;
break;
}
case 36: /* CLOFF */
{
cluster = FREQWINDOW;
break;
}
case 37: /* RECONNECT */
{
if ((nopacket == 0) && (packetinterface > 0)) {
cleanup_telnet();
init_packet();
packet();
}
break;
}
case 38: /* EXIT=QUIT */
{
writeparas();
cleanup_telnet();
endwin();
puts("\n\nThanks for using TLF.. 73\n");
exit(0);
break;
}
case 39: /* CQDELAY */
{
mvprintw(12, 29, "CQD: pgup/dwn", cqdelay);
refreshp();
x = 1;
while (x) {
x = onechar();
switch (x) {
case 156:{
if (cqdelay <= 60) {
cqdelay++;
attron(COLOR_PAIR(C_HEADER) | A_STANDOUT);
mvprintw(0, 19, " ");
mvprintw(0, 19, "%i", cqdelay);
break;
}
}
case 157:{
if (cqdelay >= 1) {
cqdelay--;
attron(COLOR_PAIR(C_HEADER) | A_STANDOUT);
mvprintw(0, 19, " ");
mvprintw(0, 19, "%i", cqdelay);
break;
}
default:
x = 0;
}
}
}
if (use_rxvt == 0)
attron(COLOR_PAIR(NORMCOLOR) | A_BOLD);
else
attron(COLOR_PAIR(NORMCOLOR));
mvprintw(12, 29 + strlen(hiscall), "");
break;
}
case 40: /* ADIF */
{
write_adif();
break;
}
case 41: /* SYNC */
{
if (strlen(synclogfile) > 0)
synclog(synclogfile);
scroll_log();
/** \todo register return value */
readcalls();
clear_display();
break;
}
case 42: /* RESCORE */
{
/** \todo register return value */
readcalls();
clear_display();
break;
}
case 43: /* SCVOLUME - set soundcard volume */
{
volumebuffer = atoi(sc_volume);
mvprintw(12, 29, "Vol: pgup/dwn");
refreshp();
usleep(500000);
mvprintw(12, 29, "Vol: ");
mvprintw(12, 29, "Vol: %d", volumebuffer);
x = 1;
while (x) {
x = onechar();
switch (x) {
case 156:{
if (volumebuffer < 95)
volumebuffer += 5;
break;
}
case 157:{
if (volumebuffer >= 5)
volumebuffer -= 5;
break;
}
default:
x = 0;
}
attron(COLOR_PAIR(COLOR_GREEN) | A_STANDOUT);
mvprintw(12, 34, " ");
mvprintw(12, 34, "%d", volumebuffer);
if (volumebuffer >= 0 && volumebuffer <= 99)
sprintf(sc_volume, "%d", volumebuffer);
netkeyer(K_STVOLUME, sc_volume);
}
clear_display();
break;
}
case 44: /* SCAN */
{
int currentterm = miniterm;
miniterm = 0;
testaudio();
clear_display();
miniterm = currentterm;
break;
}
case 45: /* DEBUG */
{
debug_tty();
clear_display();
break;
}
case 46: /* MINITERM ON/OFF */
{
if (miniterm == 1)
miniterm = 0;
else
miniterm = 1;
break;
}
case 47: /* RTTY Initialize mode (MFJ1278B controller) */
{
strcpy(buffer, "MODE VB");
buffer[7] = '\015';
buffer[8] = 'K';
buffer[9] = '\015';
buffer[10] = '\0';
sendbuf();
trxmode = DIGIMODE;
break;
}
case 48: /* SOUND */
{
clear_display();
record();
clear_display();
break;
}
case 50: /* CHARS */
{
mvprintw(13, 29, "Autosend: (0 (off), 2..5 chars) ?");
refreshp();
x = 1;
/* wait for correct input or ESC */
while ((x != 0) && ((x < 2) || (x > 5)) ) {
x = onechar();
if (x == 27)
break;
x = x - '0';
}
/* remember new setting */
if (x != 27)
cwstart = x;
if (cwstart)
mvprintw(13,29, "Autosend now: %1d ",
cwstart);
else
mvprintw(13,29, "Autosend now: OFF ");
refreshp();
break;
}
default:
{
nopar = 1;
}
}
if (nopar != 1) {
mvprintw(12, 29, "OK ! ");
writeparas();
} else {
if ((nopacket ==0) && (packetinterface > 0))
packet();
}
refreshp();
if (use_rxvt == 0)
attron(COLOR_PAIR(NORMCOLOR) | A_BOLD);
else
attron(COLOR_PAIR(NORMCOLOR));
mvprintw(12, 29, " ");
mvprintw(12, 29, "");
refreshp();
hiscall[0] = '\0';
return (0);
}
void Foam::kineticTheoryModel::solve(const volTensorField& gradUat)
{
if (!kineticTheory_)
{
return;
}
const scalar sqrtPi = sqrt(constant::mathematical::pi);
surfaceScalarField phi(1.5*rhoa_*phia_*fvc::interpolate(alpha_));
volTensorField dU(gradUat.T()); //fvc::grad(Ua_);
volSymmTensorField D(symm(dU));
// NB, drag = K*alpha*beta,
// (the alpha and beta has been extracted from the drag function for
// numerical reasons)
volScalarField Ur(mag(Ua_ - Ub_));
volScalarField betaPrim(alpha_*(1.0 - alpha_)*draga_.K(Ur));
// Calculating the radial distribution function (solid volume fraction is
// limited close to the packing limit, but this needs improvements)
// The solution is higly unstable close to the packing limit.
gs0_ = radialModel_->g0
(
min(max(alpha_, scalar(1e-6)), alphaMax_ - 0.01),
alphaMax_
);
// particle pressure - coefficient in front of Theta (Eq. 3.22, p. 45)
volScalarField PsCoeff
(
granularPressureModel_->granularPressureCoeff
(
alpha_,
gs0_,
rhoa_,
e_
)
);
// 'thermal' conductivity (Table 3.3, p. 49)
kappa_ = conductivityModel_->kappa(alpha_, Theta_, gs0_, rhoa_, da_, e_);
// particle viscosity (Table 3.2, p.47)
mua_ = viscosityModel_->mua(alpha_, Theta_, gs0_, rhoa_, da_, e_);
dimensionedScalar Tsmall
(
"small",
dimensionSet(0 , 2 ,-2 ,0 , 0, 0, 0),
1.0e-6
);
dimensionedScalar TsmallSqrt = sqrt(Tsmall);
volScalarField ThetaSqrt(sqrt(Theta_));
// dissipation (Eq. 3.24, p.50)
volScalarField gammaCoeff
(
12.0*(1.0 - sqr(e_))*sqr(alpha_)*rhoa_*gs0_*(1.0/da_)*ThetaSqrt/sqrtPi
);
// Eq. 3.25, p. 50 Js = J1 - J2
volScalarField J1(3.0*betaPrim);
volScalarField J2
(
0.25*sqr(betaPrim)*da_*sqr(Ur)
/(max(alpha_, scalar(1e-6))*rhoa_*sqrtPi*(ThetaSqrt + TsmallSqrt))
);
// bulk viscosity p. 45 (Lun et al. 1984).
lambda_ = (4.0/3.0)*sqr(alpha_)*rhoa_*da_*gs0_*(1.0+e_)*ThetaSqrt/sqrtPi;
// stress tensor, Definitions, Table 3.1, p. 43
volSymmTensorField tau(2.0*mua_*D + (lambda_ - (2.0/3.0)*mua_)*tr(D)*I);
if (!equilibrium_)
{
// construct the granular temperature equation (Eq. 3.20, p. 44)
// NB. note that there are two typos in Eq. 3.20
// no grad infront of Ps
// wrong sign infront of laplacian
fvScalarMatrix ThetaEqn
(
fvm::ddt(1.5*alpha_*rhoa_, Theta_)
+ fvm::div(phi, Theta_, "div(phi,Theta)")
==
fvm::SuSp(-((PsCoeff*I) && dU), Theta_)
+ (tau && dU)
+ fvm::laplacian(kappa_, Theta_, "laplacian(kappa,Theta)")
+ fvm::Sp(-gammaCoeff, Theta_)
+ fvm::Sp(-J1, Theta_)
+ fvm::Sp(J2/(Theta_ + Tsmall), Theta_)
);
ThetaEqn.relax();
ThetaEqn.solve();
}
else
{
// equilibrium => dissipation == production
// Eq. 4.14, p.82
volScalarField K1(2.0*(1.0 + e_)*rhoa_*gs0_);
volScalarField K3
(
0.5*da_*rhoa_*
(
(sqrtPi/(3.0*(3.0-e_)))
*(1.0 + 0.4*(1.0 + e_)*(3.0*e_ - 1.0)*alpha_*gs0_)
+1.6*alpha_*gs0_*(1.0 + e_)/sqrtPi
)
);
volScalarField K2
(
4.0*da_*rhoa_*(1.0 + e_)*alpha_*gs0_/(3.0*sqrtPi) - 2.0*K3/3.0
);
volScalarField K4(12.0*(1.0 - sqr(e_))*rhoa_*gs0_/(da_*sqrtPi));
volScalarField trD(tr(D));
volScalarField tr2D(sqr(trD));
volScalarField trD2(tr(D & D));
volScalarField t1(K1*alpha_ + rhoa_);
volScalarField l1(-t1*trD);
volScalarField l2(sqr(t1)*tr2D);
volScalarField l3
(
4.0
*K4
*max(alpha_, scalar(1e-6))
*(2.0*K3*trD2 + K2*tr2D)
);
Theta_ = sqr((l1 + sqrt(l2 + l3))/(2.0*(alpha_ + 1.0e-4)*K4));
}
Theta_.max(1.0e-15);
Theta_.min(1.0e+3);
volScalarField pf
(
frictionalStressModel_->frictionalPressure
(
alpha_,
alphaMinFriction_,
alphaMax_,
Fr_,
eta_,
p_
)
);
PsCoeff += pf/(Theta_+Tsmall);
PsCoeff.min(1.0e+10);
PsCoeff.max(-1.0e+10);
// update particle pressure
pa_ = PsCoeff*Theta_;
// frictional shear stress, Eq. 3.30, p. 52
volScalarField muf
(
frictionalStressModel_->muf
(
alpha_,
alphaMax_,
pf,
D,
phi_
)
);
// add frictional stress
mua_ += muf;
mua_.min(1.0e+2);
mua_.max(0.0);
Info<< "kinTheory: max(Theta) = " << max(Theta_).value() << endl;
volScalarField ktn(mua_/rhoa_);
Info<< "kinTheory: min(nua) = " << min(ktn).value()
<< ", max(nua) = " << max(ktn).value() << endl;
Info<< "kinTheory: min(pa) = " << min(pa_).value()
<< ", max(pa) = " << max(pa_).value() << endl;
}
/**
* Generates post of @p state, by constructing posts of lesser level and
* doing the union of these states with projection over the prefix
*
* @param aut: base automaton
* @param state: initial state we are generating post for
* @param level: level of inception
* @param prefix: list of variables for projection
* @return MTBDD representing the post of the state @p state
*/
MacroTransMTBDD GetMTBDDForPost(Automaton & aut, TStateSet* state, unsigned level, PrefixListType & prefix) {
// Convert MTBDD from VATA to MacroStateRepresentation
if (level == 0) {
// Is Leaf State set
LeafStateSet* lState = reinterpret_cast<LeafStateSet*>(state);
StateType stateValue = lState->state;
TransMTBDD *stateTransition = getMTBDDForStateTuple(aut, Automaton::StateTuple({stateValue}));
int projecting = getProjectionVariable(level, prefix);
StateDeterminizatorFunctor sdf;
if (projecting > 0) {
AdditionApplyFunctor adder;
TransMTBDD projected = stateTransition->Project(
[stateTransition, projecting](size_t var) {return var < projecting;}, adder);
return sdf(projected);
} else {
// Convert to TStateSet representation
return sdf(*stateTransition);
}
} else {
MacroStateSet* mState = reinterpret_cast<MacroStateSet*>(state);
// Look into cache
#if (USE_BDDCACHE == true)
if(BDDCache.inCache(mState, level)) {
return BDDCache.lookUp(mState, level);
}
#endif
StateSetList states = mState->getMacroStates();
// get post for all states under lower level
TStateSet* front;
MacroStateDeterminizatorFunctor msdf;
#if USE_PRUNED_UNION_FUNCTOR
MacroPrunedUnionFunctor muf(level-1);
#else
MacroUnionFunctor muf;
#endif
MacroTransMTBDD detResultMtbdd(new MacroStateSet());
// get first and determinize it
//const MacroTransMTBDD & frontPost = GetMTBDDForPost(aut, front, level-1, prefix);
int projecting = getProjectionVariable(level-1, prefix);
/*MacroTransMTBDD detResultMtbdd = (level == 1) ? frontPost : (msdf(frontPost)).Project(
[&frontPost, projecting](size_t var) {return var < projecting;}, muf);*/
// do the union of posts represented as mtbdd
while(!states.empty()) {
front = states.back();
states.pop_back();
if(front->isEmpty()) {
continue;
}
const MacroTransMTBDD & nextPost = GetMTBDDForPost(aut, front, level-1, prefix);
detResultMtbdd = muf(detResultMtbdd, (level == 1) ? nextPost : (msdf(nextPost)).Project(
[&nextPost, projecting](size_t var) {return var < projecting;}, muf));
}
// cache the results
#if (USE_BDDCACHE == true)
BDDCache.storeIn(mState, detResultMtbdd, level);
#endif
// do projection and return;
return detResultMtbdd;
}
}
