int main ( int argc, char **argv ) {
std::shared_ptr<Vipper::IRouter> nextRouter;
std::shared_ptr<Vipper::IPresenter> presenter;
std::shared_ptr<Vipper::IRouter> router;
auto renderer = std::make_shared<BlockyFalls::CSDLRenderer>();
auto titleScreenRouter = std::make_shared<BlockyFalls::CTitleScreenRouter>(renderer);
titleScreenRouter->initWithDefaults();
nextRouter = titleScreenRouter;
presenter = titleScreenRouter->getPresenter();
std::vector<std::shared_ptr<Vipper::IRouter>> routers;
do {
if ( nextRouter != nullptr ) {
if ( router != nullptr ) {
router->onRelinquishFocus();
if ( !router->isFinished() ) {
routers.push_back( router );
}
}
router = nextRouter;
router->onFocus();
nextRouter = nullptr;
}
presenter = router->getPresenter();
presenter->present();
renderer->update();
renderer->render();
SDL_Delay(33);
nextRouter = router->route();
if ( router->isFinished() && nextRouter == nullptr ) {
nextRouter = routers.back();
routers.pop_back();
}
} while ( nextRouter != router );
renderer->shutdown();
return 0;
}
/**
* Save an MD workspace to a vts/vtu file.
* @param workspace: the workspace which is to be saved.
* @param filename: the name of the file to which the workspace is to be saved.
* @param normalization: the visual normalization option
* @param recursionDepth: the recursion depth for MDEvent Workspaces determines
* @param compressorType: the compression type used by VTK
* from which level data should be displayed
*/
void SaveMDWorkspaceToVTKImpl::saveMDWorkspace(
const Mantid::API::IMDWorkspace_sptr &workspace,
const std::string &filename, VisualNormalization normalization,
int recursionDepth, const std::string &compressorType) {
auto isHistoWorkspace =
boost::dynamic_pointer_cast<Mantid::API::IMDHistoWorkspace>(workspace) !=
nullptr;
auto fullFilename = getFullFilename(filename, isHistoWorkspace);
const vtkXMLWriter::CompressorType compressor = [&compressorType] {
if (compressorType == "NONE") {
return vtkXMLWriter::NONE;
} else if (compressorType == "ZLIB") {
return vtkXMLWriter::ZLIB;
} else {
// This should never happen.
Mantid::Kernel::Logger g_log("SaveMDWorkspaceToVTK");
g_log.warning("Incorrect CompressorType: " + compressorType +
". Using CompressorType=NONE.");
return vtkXMLWriter::NONE;
}
}();
// Define a time slice.
auto time = selectTimeSliceValue(*workspace);
// Get presenter and data set factory set up
auto factoryChain =
getDataSetFactoryChain(isHistoWorkspace, normalization, time);
auto presenter = getPresenter(isHistoWorkspace, workspace, recursionDepth);
// Create the vtk data
NullProgressAction nullProgressA;
NullProgressAction nullProgressB;
auto dataSet =
presenter->execute(factoryChain.get(), nullProgressA, nullProgressB);
// Do an orthogonal correction
dataSet = getDataSetWithOrthogonalCorrection(dataSet, presenter.get(),
workspace, isHistoWorkspace);
// ParaView 5.1 checks the range of the entire signal array, including blank
// cells.
if (isHistoWorkspace) {
auto structuredGrid = vtkStructuredGrid::SafeDownCast(dataSet);
vtkIdType imageSize = structuredGrid->GetNumberOfCells();
vtkNew<vtkFloatArray> signal;
signal->SetNumberOfComponents(1);
signal->SetNumberOfTuples(imageSize);
auto oldSignal = structuredGrid->GetCellData()->GetScalars();
for (vtkIdType index = 0; index < imageSize; ++index) {
if (structuredGrid->IsCellVisible(index)) {
signal->SetComponent(index, 0, oldSignal->GetTuple1(index));
} else {
signal->SetComponent(index, 0, std::numeric_limits<float>::quiet_NaN());
}
}
structuredGrid->GetCellData()->SetScalars(signal.GetPointer());
}
// Write the data to the file
vtkSmartPointer<vtkXMLWriter> writer = getXMLWriter(isHistoWorkspace);
auto writeSuccessFlag =
writeDataSetToVTKFile(writer, dataSet, fullFilename, compressor);
if (!writeSuccessFlag) {
throw std::runtime_error("SaveMDWorkspaceToVTK: VTK could not write "
"your data set to a file.");
}
}
