Function main

postprocessing/validation/XDMFCubeCheck/src/main.cpp:57–211

}

int main(int argc, char* argv[])
{
	utils::Args args;
	args.addOption("diff", 'd', "Maximum allowed difference (Default 0.001)", utils::Args::Required, false);
	args.addOption("verbose", 'V', "Enable verbose output", utils::Args::No, false);
	args.addAdditionalOption("output.h5", "The SeisSol wave field output file");

	switch (args.parse(argc, argv)) {
	case utils::Args::Error:
		return 1;
	case utils::Args::Help:
		return 127;
	default:
		break;
	}

	bool verbose = args.getArgument("verbose", false);

	hid_t hdfFile = H5Fopen(args.getAdditionalArgument<const char*>("output.h5"),
			H5F_ACC_RDONLY, H5P_DEFAULT);

	// Read vertices
	hid_t vertData = H5Dopen(hdfFile, "/geometry", H5P_DEFAULT);
	hid_t space = H5Dget_space(vertData);
	hsize_t dims[2]; // Assuming to dimensions
	H5Sget_simple_extent_dims(space, dims, 0L);

	double* vertices = new double[dims[0]*dims[1]];
	H5Dread(vertData, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, vertices);

	H5Sclose(space);
	H5Dclose(vertData);

	// Read cells
	hid_t cellData = H5Dopen(hdfFile, "/connect", H5P_DEFAULT);
	space = H5Dget_space(cellData);
	H5Sget_simple_extent_dims(space, dims, 0L);

	unsigned long* cells = new unsigned long[dims[0]*dims[1]];
	H5Dread(cellData, H5T_NATIVE_ULONG, H5S_ALL, H5S_ALL, H5P_DEFAULT, cells);

	hsize_t numCells = dims[0];

	H5Sclose(space);
	H5Dclose(cellData);

	// Read the 9 variables
	double* var[9];
	const char* varNames[9] = {"/sigma_xx", "/sigma_yy", "/sigma_zz", "/sigma_xy", "/sigma_yz", "/sigma_xz", "/u", "/v", "/w"};

	hid_t memSpace = H5Screate_simple(1, &numCells, &numCells);
	H5Sselect_all(memSpace);

	for (unsigned int i = 0; i < 9; i++) {
		hid_t varData = H5Dopen(hdfFile, varNames[i], H5P_DEFAULT);
		space = H5Dget_space(varData);
		H5Sget_simple_extent_dims(space, dims, 0L);