void RCSSeparateResponse::set()
{
assert(m_request.getOCRequest());
auto resObj = m_request.getResourceObject().lock();
if (!resObj)
{
throw RCSBadRequestException{ "ResourceObject is unspecified(or destroyed)!" };
}
if (m_done) throw RCSBadRequestException{ "The response is already set!" };
auto ocRequest = m_request.getOCRequest();
auto response = std::make_shared< OC::OCResourceResponse >();
response->setRequestHandle(ocRequest->getRequestHandle());
response->setResourceHandle(ocRequest->getResourceHandle());
response->setResponseResult(OC_EH_OK);
// TODO the response should be different by the request interface.
response->setResourceRepresentation(
RCSRepresentation::toOCRepresentation(resObj->toRepresentation()));
invokeOCFunc(OC::OCPlatform::sendResponse, response);
m_done = true;
}
SEXP stochasticEnumeration(SEXP nVertices_sexp, SEXP budget_sexp, SEXP seed_sexp, SEXP options_sexp)
{
BEGIN_RCPP
int seed = Rcpp::as<int>(seed_sexp);
int nVertices = Rcpp::as<int>(nVertices_sexp);
int budget = Rcpp::as<int>(budget_sexp);
Rcpp::List options = options_sexp;
bool outputSamples = false;
if(!options.containsElementNamed("reduceChains")) throw std::runtime_error("Unable to find option named reduceChains");
if(!options.containsElementNamed("graphRepresentation")) throw std::runtime_error("Unable to find option named graphRepresentation");
bool reduceChains = Rcpp::as<bool>(options("reduceChains"));
std::string graphRepresentationString = Rcpp::as<std::string>(options("graphRepresentation"));
graphRepresentation representation = toRepresentation(graphRepresentationString);
if (options.containsElementNamed("outputSamples"))
{
outputSamples = Rcpp::as<bool>(options("outputSamples"));
}
boost::mt19937 randomSource;
randomSource.seed(seed);
chordalGraph::stochasticEnumerationArgs args(randomSource);
args.nVertices = nVertices;
args.budget = budget;
args.outputSamples = outputSamples;
int maxEdges = ((nVertices - 1)*nVertices / 2);
//Set up text progress bar
Rcpp::Function txtProgressBar("txtProgressBar");
Rcpp::Function close("close");
Rcpp::Function setTxtProgressBar("setTxtProgressBar");
Rcpp::RObject barHandle = txtProgressBar(Rcpp::Named("style") = 3, Rcpp::Named("min") = 0, Rcpp::Named("max") = maxEdges, Rcpp::Named("initial") = 0);
//List for samples
Rcpp::List samples(maxEdges+1);
std::vector<bool> exactVector;
std::vector<int> minimumSizeForExactVector;
std::vector<std::string> estimatesAsStrings;
for (int nEdges = 0; nEdges < maxEdges+1; nEdges++)
{
setTxtProgressBar(barHandle, nEdges);
args.nEdges = nEdges;
args.samples.clear();
if(reduceChains)
{
if(representation == listRepresentation)
{
chordalGraph::stochasticEnumerationReduceChains<chordalGraph::cliqueTree>(args);
}
else
{
chordalGraph::stochasticEnumerationReduceChains<chordalGraph::cliqueTreeAdjacencyMatrix>(args);
}
}
else
{
if(representation == listRepresentation)
{
chordalGraph::stochasticEnumeration<chordalGraph::cliqueTree>(args);
}
else
{
chordalGraph::stochasticEnumeration<chordalGraph::cliqueTreeAdjacencyMatrix>(args);
}
}
estimatesAsStrings.push_back(args.estimate.str());
exactVector.push_back(args.exact);
if (args.outputSamples)
{
Rcpp::List currentEdgesSamples;
convert(currentEdgesSamples, args.samples);
samples[nEdges] = currentEdgesSamples;
}
if(args.minimumSizeForExact == -1)
{
minimumSizeForExactVector.push_back(NA_INTEGER);
}
else
{
minimumSizeForExactVector.push_back(args.minimumSizeForExact);
}
}
close(barHandle);
SEXP estimatesAsStrings_sexp = Rcpp::wrap(estimatesAsStrings);
SEXP minimumSizeForExactVector_sexp = Rcpp::wrap(minimumSizeForExactVector);
Rcpp::Function mpfrFunction("mpfr");
Rcpp::List listResult = Rcpp::List::create(Rcpp::Named("data") = mpfrFunction(estimatesAsStrings_sexp, Rcpp::Named("prec", 50)), Rcpp::Named("exact") = Rcpp::wrap(exactVector), Rcpp::Named("minimumSizeForExact") = minimumSizeForExactVector_sexp);
if (args.outputSamples)
{
listResult("samples") = samples;
}
return listResult;
END_RCPP
}
SEXP stochasticEnumerationSpecificEdges(SEXP nVertices_sexp, SEXP nEdges_sexp, SEXP budget_sexp, SEXP seed_sexp, SEXP options_sexp)
{
BEGIN_RCPP
int seed = Rcpp::as<int>(seed_sexp);
int nVertices = Rcpp::as<int>(nVertices_sexp);
int nEdges = Rcpp::as<int>(nEdges_sexp);
int budget = Rcpp::as<int>(budget_sexp);
Rcpp::List options = options_sexp;
bool outputSamples = false;
if(!options.containsElementNamed("reduceChains")) throw std::runtime_error("Unable to find option named reduceChains");
if(!options.containsElementNamed("graphRepresentation")) throw std::runtime_error("Unable to find option named graphRepresentation");
bool reduceChains = Rcpp::as<bool>(options("reduceChains"));
std::string graphRepresentationString = Rcpp::as<std::string>(options("graphRepresentation"));
graphRepresentation representation = toRepresentation(graphRepresentationString);
if (options.containsElementNamed("outputSamples"))
{
outputSamples = Rcpp::as<bool>(options("outputSamples"));
}
boost::mt19937 randomSource;
randomSource.seed(seed);
chordalGraph::stochasticEnumerationArgs args(randomSource);
args.nEdges = nEdges;
args.nVertices = nVertices;
args.budget = budget;
args.outputSamples = outputSamples;
if(reduceChains)
{
if(representation == listRepresentation)
{
chordalGraph::stochasticEnumerationReduceChains<chordalGraph::cliqueTree>(args);
}
else
{
chordalGraph::stochasticEnumerationReduceChains<chordalGraph::cliqueTreeAdjacencyMatrix>(args);
}
}
else
{
if(representation == listRepresentation)
{
chordalGraph::stochasticEnumeration<chordalGraph::cliqueTree>(args);
}
else
{
chordalGraph::stochasticEnumeration<chordalGraph::cliqueTreeAdjacencyMatrix>(args);
}
}
std::string estimateAsString = args.estimate.str();
SEXP estimateAsString_sexp = Rcpp::wrap(estimateAsString);
SEXP minimumSizeForExact_sexp;
if(args.minimumSizeForExact != -1) minimumSizeForExact_sexp = Rcpp::wrap(args.minimumSizeForExact);
else minimumSizeForExact_sexp = Rcpp::wrap(NA_INTEGER);
Rcpp::Function mpfrFunction("mpfr");
Rcpp::List listResult = Rcpp::List::create(Rcpp::Named("data") = mpfrFunction(estimateAsString_sexp, Rcpp::Named("prec", 50)), Rcpp::Named("exact") = Rcpp::wrap(args.exact), Rcpp::Named("minimumSizeForExact") = minimumSizeForExact_sexp);
if (args.outputSamples)
{
Rcpp::List samples;
convert(samples, args.samples);
listResult("samples") = samples;
}
return listResult;
END_RCPP
}
