读取ugrid格式文件

#include <iostream>
#include <vector>
#include <string>
#include <fstream>

using namespace std;

int main()
{
	int node_num, tri_bd_num, quad_bd_num, tet_num, pyramid_num, prism_num, hex_num;
	node_num = tri_bd_num = quad_bd_num = tet_num = pyramid_num = prism_num = hex_num = 0;
	ofstream fout("test.ugrid");
	ifstream fin("ugrid_bin.ugrid", ios::binary | ios::in);

	if (!fin.is_open())
	{
		cout << "文件未正常读入" << endl;
		return 0;
	}
	if (!fout.is_open())
	{
		cout << "文件输出有误" << endl;
	}

	fin.read((char*)&node_num, sizeof(int));
	fin.read((char*)&tri_bd_num, sizeof(int));
	fin.read((char*)&quad_bd_num, sizeof(int));
	fin.read((char*)&tet_num, sizeof(int));
	fin.read((char*)&pyramid_num, sizeof(int));
	fin.read((char*)&prism_num, sizeof(int));
	fin.read((char*)&hex_num, sizeof(int));

	fout << node_num << " " << tri_bd_num << " " << quad_bd_num << " " << tet_num << " " << pyramid_num << " " << prism_num << " " << hex_num << endl;

	// 边界面数目和实单元数目
	int bd_fac_num = tri_bd_num + quad_bd_num;
	int ele_num = tet_num + pyramid_num + prism_num + hex_num;
	
	// 点坐标
	double* Node = new double[node_num * 3];

	// 三角形边界面的点索引
	int *tri_bd_fac = new int[tri_bd_num*3];
	// 四边形边界面的点索引
	int* quad_bd_fac = new int[quad_bd_num*4];
	// 边界面元对应的surface ID
	int* fac_type = new int[bd_fac_num];
	// 四面体单元的点索引
	int* tet = new int[tet_num * 4];
	// 金字塔单元的点索引
	int* pyramid = new int[pyramid_num * 5];
	// 三棱柱单元的点索引
	int* prism = new int[prism_num * 6];
	// 六面体单元的点索引
	int* hex = new int[hex_num * 8];

	// 读取点坐标
	fin.read((char*)Node, sizeof(double) * 3 * node_num);
	for (auto i = 0; i < node_num; i++)
	{
		fout << Node[i * 3 + 0] << " " << Node[i * 3 + 1] << " " << Node[i * 3 + 2] << endl;
	}

	// 读取三角形边界面的点索引
	fin.read((char*)tri_bd_fac, sizeof(int) * 3 * tri_bd_num);
	for (auto i = 0; i < tri_bd_num; i++)
	{
		fout << tri_bd_fac[i * 3 + 0] << " " << tri_bd_fac[i * 3 + 1] << " " << tri_bd_fac[i * 3 + 2] << endl;
	}
	
	// 读取四角形边界面的点索引
	fin.read((char*)quad_bd_fac, sizeof(int) * 4 * quad_bd_num);
	for (auto i = 0; i < quad_bd_num; i++)
	{
		fout << quad_bd_fac[i * 4 + 0] << " " << quad_bd_fac[i * 4 + 1] << " " 
			<< quad_bd_fac[i * 4 + 2]<<" " << quad_bd_fac[i * 4 + 3] << endl;
	}
	
	// 读取边界面元对应的 surface ID
	fin.read((char*)fac_type, sizeof(int)*bd_fac_num);
	for (auto i = 0; i < bd_fac_num; i++)
	{
		fout << fac_type[i] << endl;
	}
	
	// 读取四面体单元的点索引
	fin.read((char*)tet, sizeof(int) * 4 * tet_num);
	for (auto i = 0; i < tet_num; i++)
	{
		fout << tet[4 * i + 0] << " " << tet[4 * i + 1] << " "
			<< tet[4 * i + 2] << " " << tet[4 * i + 3] << endl;
	}

	// 读取金字塔单元的点索引
	fin.read((char*)pyramid, sizeof(int) * 5 * pyramid_num);
	for (auto i = 0; i < pyramid_num; i++)
	{
		fout << pyramid[4 * i + 0] << " " << pyramid[4 * i + 1] << " "
			<< pyramid[4 * i + 2] << " " << pyramid[4 * i + 3] << " "
			<< pyramid[4 * i + 4] << endl;
	}

	// 读取三棱柱单元的点索引
	fin.read((char*)prism, sizeof(int) * 6 * prism_num);
	for (auto i = 0; i < prism_num; i++)
	{
		fout << prism[4 * i + 0] << " " << prism[4 * i + 1] << " "
			<< prism[4 * i + 2] << " " << prism[4 * i + 3] << " "
			<< prism[4 * i + 4] << " " << prism[4 * i + 5] << endl;
	}

	// 读取四面体单元的点索引
	fin.read((char*)hex, sizeof(int) * 8 * hex_num);
	for (auto i = 0; i < hex_num; i++)
	{
		fout << hex[4 * i + 0] << " " << hex[4 * i + 1] << " "
			<< hex[4 * i + 2] << " " << hex[4 * i + 3] << " "
			<< hex[4 * i + 4] << " " << hex[4 * i + 5] << " "
			<< hex[4 * i + 6] << " " << hex[4 * i + 7] << endl;
	}

	fin.close();
	fout.close();

	// 点坐标
	delete[] Node;
	delete[] tri_bd_fac;
	delete[] quad_bd_fac;
	delete[] fac_type;
	delete[] tet;
	delete[] pyramid;
	delete[] prism;
	delete[] hex;

	return 0;
}