main.cpp

#include <iostream>
#include <sycl/sycl.hpp>

#include "ConfigReader.h"
#include "InterpolationFiller.h"
#include "MobileMappingReader.h"
#include "PngWriter.h"
#include "RadialFiller.h"
#include "Rasterizer.h"
#include "Sorter.h"
#include "SYCLState.h"

int main(const int argc, const char *argv[]) {
	if (argc < 2) {
		std::cerr << "Provide a path to a config file!\n";
		return 1;
	}

	const auto config = EFTDEM::ConfigReader::readConfig(argv[1]);

	int debug;
	if (config.numDebugLines != 0) debug = std::max(config.numDebugLines, EFTDEM_DEBUG);
	else debug = EFTDEM_DEBUG * 25;

	EFTDEM::PointCloud pointCloud{config.width, config.height};

	EFTDEM::MobileMappingReader::readPointsFromFile(pointCloud, config.inputPath, debug);

	pointCloud.heights = std::vector<float>(pointCloud.width * pointCloud.height, 0);

	const auto startTime = std::chrono::high_resolution_clock::now();

	// Trigger write back to host on buffer destruction
	{
		sycl::queue queue;
		if (config.useGPU) queue = sycl::queue{sycl::gpu_selector_v, sycl::property::queue::in_order{}};
		else queue = sycl::queue{sycl::cpu_selector_v, sycl::property::queue::in_order{}};

		EFTDEM::SYCLState syclState{
			queue,
			std::nullopt,
			std::nullopt,
			std::nullopt
		};
		std::cout << "Selected " << syclState.queue.get_device().get_info<sycl::info::device::name>() << " as the SYCL device.\n";

		EFTDEM::Sorter::sortPointCloud(pointCloud, syclState, debug);

		EFTDEM::Rasterizer::rasterizePointCloud(pointCloud, syclState, debug);

		EFTDEM::RadialFiller::fill(pointCloud, syclState, config.numClosingIterations, debug);
		//EFTDEM::InterpolationFiller::fill(pointCloud, syclState, config.numClosingIterations, debug);
	}

	const auto endTime = std::chrono::high_resolution_clock::now();

	std::cout << "Elapsed time for SYCL operations: " << std::chrono::duration_cast<std::chrono::milliseconds>(endTime - startTime).count() << " ms\n";

	EFTDEM::PngWriter::exportPointCloud(pointCloud, config.outputPath);

	return 0;
}