简介这份资源是面向CFD研究者与工程师的CGNS MATLAB接口库用于在MATLAB环境中直接读写、分析CGNS格式的流体仿真数据解决跨软件数据交换与后处理效率低的问题。压缩包共514个文件约1.33MB以476个m脚本为主体辅以13个c与12个h源码、mexw64/mexmaci/mexglx等多平台MEX二进制文件及lib、inc、readme等覆盖从底层C库封装到MATLAB调用层的完整链路。已有777人学习下载。借助它用户可读取网格、流场与边界条件将MATLAB计算结果写回CGNS格式并利用内置函数完成积分、微分方程求解与统计分析再通过MATLAB绘图快速生成流线、速度矢量与压力分布图配合脚本还能批量处理多区域、多时间序列文件实现自动化后处理与算法迭代是扩展MATLAB在流体模拟领域应用范围的实用工具。1. CGNS 在 MATLAB 里到底能干什么别被“版本”两个字带偏很多人搜“cgns matlab版本”脑子里想的是找一个安装包双击装完就能在 MATLAB 里读写 CGNS 文件。实际干过 CFD 后处理的人都知道CGNS 本身是一套跨平台的数据标准不是某个软件的版本号。MATLAB 官方并没有内置 CGNS 读写函数所谓“MATLAB 版本”通常指两条路一是用 MEX 把 C 语言的 CGNS 库包一层给 MATLAB 调二是用纯 MATLAB 脚本解析 CGNS 的 ADF/HDF5 底层结构。前者性能好、功能全后者部署简单但只适合小规模网格。如果你手头有几十万节点的结构网格还要带边界条件和流场解那基本只能走 MEX 路线。这篇笔记就按这个现实来拆先讲清 CGNS 在 MATLAB 里的数据映射逻辑再给一套能跑通的最小读写流程最后把参数设置和踩坑点摊开说。2. 用 MEX 把 CGNS 库接进 MATLAB编译链路与最小读写闭环2.1 为什么绕不开 MEX纯 MATLAB 解析的边界在哪CGNS 文件内部是一棵节点树每个节点有类型、维度、数据数组和子节点。纯 MATLAB 可以用h5read读 HDF5 底座但 CGNS 的 ADF 层做了大量元数据封装比如Zone_t下面挂GridCoordinates_t再下面才是DataArray_t。你如果只用h5read拿到的是扁平数据集得自己重建父子关系。我试过用纯脚本读一个 12 万节点的四面体网格光解析节点树就花了 40 多秒内存峰值到 1.8 GB。MEX 方案直接调cg_open、cg_zone_read、cg_coord_read同样的网格 0.3 秒读完内存只多了一个 C 结构体的开销。所以选型逻辑很简单网格节点数超过 5 万或者需要写边界条件、流场解就上 MEX只是偶尔看一眼小网格的坐标纯脚本凑合能用。2.2 编译环境怎么配头文件、库路径和 mex 命令假设你已经从 CGNS 官方源码编好了动态库Windows 下得到cgns.dll和cgns.libLinux 下是libcgns.so。MATLAB 里用mex编译时关键是把头文件目录和库目录喂进去。下面是一个 Linux 下的编译命令示例Windows 把-I和-L换成对应路径即可。# 在 MATLAB 命令行里执行不是系统终端 mex -v ... -I/usr/local/include ... -L/usr/local/lib ... -lcgns ... cgns_mex.c-v是打开详细输出第一次编译一定要加能看到链接器到底找了哪些库。-I后面跟cgns.h所在目录-L跟libcgns.so所在目录-lcgns告诉链接器库名。编译成功后当前目录会出现cgns_mex.mexa64。如果报undefined reference to cg_open九成是-lcgns位置不对把它放到源文件后面。Windows 下还要注意位数匹配MATLAB 是 64 位CGNS 库也必须是 64 位混用 32 位会直接崩。2.3 读一个 Zone 的坐标和流场MEX 函数参数怎么传MEX 函数的入口是mexFunction里面用mxGetPr拿 MATLAB 数组指针调 CGNS C 接口再把结果塞回plhs。下面这段代码实现一个最小功能打开 CGNS 文件读第一个 Base 下第一个 Zone 的顶点坐标和名为Pressure的流场变量。#include mex.h #include cgnslib.h void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { char *filename mxArrayToString(prhs[0]); // 输入文件名 int fn, base, zone, coord; cgsize_t size[3]; double *x, *y, *z, *p; cg_open(filename, CG_MODE_READ, fn); // 打开文件 cg_base_read(fn, 1, NULL, NULL, NULL); // 读第一个 Base cg_zone_read(fn, 1, 1, NULL, size); // 读第一个 Zone 的尺寸 int nnode size[0]; // 顶点数 plhs[0] mxCreateDoubleMatrix(nnode, 1, mxREAL); plhs[1] mxCreateDoubleMatrix(nnode, 1, mxREAL); plhs[2] mxCreateDoubleMatrix(nnode, 1, mxREAL); x mxGetPr(plhs[0]); y mxGetPr(plhs[1]); z mxGetPr(plhs[2]); cg_coord_read(fn, 1, 1, CoordinateX, RealDouble, 1, nnode, x); cg_coord_read(fn, 1, 1, CoordinateY, RealDouble, 1, nnode, y); cg_coord_read(fn, 1, 1, CoordinateZ, RealDouble, 1, nnode, z); cg_field_read(fn, 1, 1, 1, Pressure, RealDouble, 1, nnode, p); cg_close(fn); }逻辑说明cg_open的第三个参数是模式读用CG_MODE_READ。cg_zone_read的size数组前三个元素是顶点数后面是单元数。cg_coord_read的第四个参数是坐标名CGNS 标准里就是CoordinateX/Y/Z大小写敏感。cg_field_read的第四个参数是解索引通常从 1 开始。参数坑点cgsize_t在 64 位编译下是long long如果你用int接会截断大网格的节点数。另外mxArrayToString拿到的字符串要手动mxFree否则反复调用会漏内存。2.4 写回 CGNS边界条件怎么挂到 Zone 上写文件比读多一步先建 Base、建 Zone、写坐标、写解最后挂边界条件。下面这段把 MATLAB 里的坐标和压力写成一个新 CGNS 文件并给第一个面加一个BCWall边界。cg_open(out.cgns, CG_MODE_WRITE, fn); cg_base_write(fn, Base, 3, 3, base); // 3 维3 维 cg_zone_write(fn, base, Zone1, size, Structured, zone); cg_coord_write(fn, base, zone, RealDouble, CoordinateX, x, coord); cg_coord_write(fn, base, zone, RealDouble, CoordinateY, y, coord); cg_coord_write(fn, base, zone, RealDouble, CoordinateZ, z, coord); cg_sol_write(fn, base, zone, FlowSolution, Vertex, sol); cg_field_write(fn, base, zone, sol, RealDouble, Pressure, p, field); cg_bc_write(fn, base, zone, 1, Wall, BCWall, bc); // 第一个面 cg_close(fn);cg_base_write的第二个参数是 Base 名第三、四个参数是物理维度和拓扑维度三维网格都填 3。cg_zone_write的size数组要按[nnode, ncell, 0]填。cg_sol_write的Vertex表示解存在顶点上如果存在单元中心就换CellCenter。cg_bc_write的第四个参数是面索引结构化网格从 1 开始。注意边界条件必须在解写完之后挂顺序反了 CGNS 库会报CG_ERROR。3. 参数设置与数据映射把 MATLAB 数组对齐到 CGNS 节点树3.1 维度顺序陷阱MATLAB 列优先 vs CGNS 行优先MATLAB 数组是列优先CGNS 底层按行优先存。一个nx * ny * nz的结构化网格在 MATLAB 里习惯写成X(i,j,k)但 CGNS 的cg_coord_write期望的是一维数组顺序是k变化最快。如果你直接把X(:)传进去读出来会发现坐标在 Z 方向错位。正确做法是用permute调整维度顺序再展平。% 假设 X 是 nx-by-ny-by-nz X_cgns permute(X, [3, 2, 1]); % 变成 nz-by-ny-by-nx X_flat X_cgns(:); % 展平后 k 变化最快参数说明permute的第二个参数是维度顺序[3,2,1]把原来的第三维换到第一维。展平后传给 MEX 函数。这个坑我踩过两次第一次读出来的网格像被拧过的麻花排查了一下午才发现是维度顺序反了。3.2 边界条件类型映射BCWall、BCFarfield 怎么选CGNS 标准定义了几十种边界条件类型常用的就几个。下面这张表是 MATLAB 侧字符串到 CGNS 枚举的映射写 MEX 时直接查表。MATLAB 侧字符串CGNS 枚举适用场景wallBCWall无滑移固壁farfieldBCFarfield远场边界inletBCInflow速度入口outletBCOutflow压力出口symmetryBCSymmetryPlane对称面选型逻辑如果做外流远场用BCFarfield物面用BCWall如果做内流入口用BCInflow并指定速度方向出口用BCOutflow。注意BCInflow需要额外写InflowBC_t结构体指定速度分量比BCWall多一层嵌套MEX 里要调cg_bc_dataset_write。3.3 流场变量命名Pressure、VelocityX 这些名字不能乱写CGNS 对解变量的名字没有强制标准但主流求解器都遵循一套约定。你写Pressure读的时候也得用Pressure写成pressure或P就找不到。下面是我常用的变量名对照。物理量推荐 CGNS 名数据类型压力PressureRealDoubleX 速度VelocityXRealDoubleY 速度VelocityYRealDoubleZ 速度VelocityZRealDouble温度TemperatureRealDouble密度DensityRealDouble如果读一个别人给的 CGNS 文件不确定变量名先用cg_nsols和cg_sol_info把解的名字列出来再决定读哪个。别硬猜猜错了返回全零你还以为数据没写进去。4. 避坑与排查CGNS 在 MATLAB 里翻车的五个典型场景4.1 现象cg_open返回错误码 0但文件明明存在原因路径里有中文或空格CGNS C 库对非 ASCII 路径支持不好。解决把文件挪到纯英文无空格目录或者用getShortPathName转短路径。我一般直接在tempdir下操作省事。4.2 现象读出来的坐标全是零但cg_coord_read没报错原因坐标名大小写写错了比如写了Coordinatex。CGNS 是大小写敏感的CoordinateX和Coordinatex是两个不同节点。解决用cg_coord_info把 Zone 里所有坐标名打印出来照着抄。4.3 现象MEX 编译通过一运行 MATLAB 就崩原因MEX 文件里用了mxGetPr拿复数数组的实部或者cgsize_t和int混用导致栈溢出。解决编译时加-largeArrayDims所有 CGNS 尺寸变量统一用cgsize_tMATLAB 数组索引用mwIndex。4.4 现象写边界条件时cg_bc_write返回CG_ERROR原因面索引超出 Zone 的面数或者边界条件类型和网格拓扑不匹配。比如二维网格写了三维的BCWall。解决先调cg_nbocos确认面数再检查cg_zone_type返回的是Structured还是Unstructured。4.5 现象读大网格时内存暴涨MATLAB 被系统杀掉原因MEX 里一次性mxCreateDoubleMatrix分配了所有节点数组加上 CGNS 库内部的缓冲区内存翻倍。解决分块读每次读 10 万节点用cg_coord_read的range_min和range_max参数控制范围MATLAB 侧用matfile追加写盘。5. 进阶技巧用 MATLAB 做 CGNS 网格质量检查与批量转换5.1 网格质量指标怎么在 MATLAB 里算读完坐标后用delaunay或delaunayn重建四面体连接关系然后算每个单元的雅可比行列式。下面这段算最小雅可比和长宽比。% 假设 X, Y, Z 是顶点坐标T 是四面体索引矩阵 dt delaunayTriangulation(X, Y, Z); T dt.ConnectivityList; J zeros(size(T,1), 1); for i 1:size(T,1) v1 [X(T(i,1)), Y(T(i,1)), Z(T(i,1))]; v2 [X(T(i,2)), Y(T(i,2)), Z(T(i,2))]; v3 [X(T(i,3)), Y(T(i,3)), Z(T(i,3))]; v4 [X(T(i,4)), Y(T(i,4)), Z(T(i,4))]; J(i) abs(det([v2-v1; v3-v1; v4-v1])); end minJ min(J);delaunayTriangulation对大规模点集比较慢超过 50 万点建议用 MEX 调 CGNS 自带的cg_conn_read直接拿连接关系省去重建。det算的是平行六面体体积除以 6 才是四面体体积但比较相对质量时不用除。5.2 批量转换把多个 CGNS 文件转成 MAT 格式如果手头有几十个 CGNS 文件要统一处理写个循环调 MEX 函数每个文件存一个.mat。注意每次循环后clear掉大变量否则内存会累积。files dir(*.cgns); for k 1:length(files) [x, y, z, p] cgns_mex(files(k).name); save([files(k).name, .mat], x, y, z, p, -v7.3); clear x y z p; end-v7.3是必须的普通-v7格式单个变量不能超过 2 GB。clear放在save之后确保写盘完成再释放。5.3 一个我常用的验证习惯每次写完 CGNS 文件不要直接信。用 MEX 读回来和原始 MATLAB 数组做isequal或norm差。我一般会算坐标的max(abs(x_read - x_orig))超过 1e-12 就说明写入或读取环节有精度损失。这个习惯帮我逮到过好几次维度顺序错误和数据类型不匹配。希望帮到你。本文还有配套的精品资源点击获取
简介这份资源是面向CFD研究者与工程师的CGNS MATLAB接口库用于在MATLAB环境中直接读写、分析CGNS格式的流体仿真数据解决跨软件数据交换与后处理效率低的问题。压缩包共514个文件约1.33MB以476个m脚本为主体辅以13个c与12个h源码、mexw64/mexmaci/mexglx等多平台MEX二进制文件及lib、inc、readme等覆盖从底层C库封装到MATLAB调用层的完整链路。已有777人学习下载。借助它用户可读取网格、流场与边界条件将MATLAB计算结果写回CGNS格式并利用内置函数完成积分、微分方程求解与统计分析再通过MATLAB绘图快速生成流线、速度矢量与压力分布图配合脚本还能批量处理多区域、多时间序列文件实现自动化后处理与算法迭代是扩展MATLAB在流体模拟领域应用范围的实用工具。1. CGNS 在 MATLAB 里到底能干什么别被“版本”两个字带偏很多人搜“cgns matlab版本”脑子里想的是找一个安装包双击装完就能在 MATLAB 里读写 CGNS 文件。实际干过 CFD 后处理的人都知道CGNS 本身是一套跨平台的数据标准不是某个软件的版本号。MATLAB 官方并没有内置 CGNS 读写函数所谓“MATLAB 版本”通常指两条路一是用 MEX 把 C 语言的 CGNS 库包一层给 MATLAB 调二是用纯 MATLAB 脚本解析 CGNS 的 ADF/HDF5 底层结构。前者性能好、功能全后者部署简单但只适合小规模网格。如果你手头有几十万节点的结构网格还要带边界条件和流场解那基本只能走 MEX 路线。这篇笔记就按这个现实来拆先讲清 CGNS 在 MATLAB 里的数据映射逻辑再给一套能跑通的最小读写流程最后把参数设置和踩坑点摊开说。2. 用 MEX 把 CGNS 库接进 MATLAB编译链路与最小读写闭环2.1 为什么绕不开 MEX纯 MATLAB 解析的边界在哪CGNS 文件内部是一棵节点树每个节点有类型、维度、数据数组和子节点。纯 MATLAB 可以用h5read读 HDF5 底座但 CGNS 的 ADF 层做了大量元数据封装比如Zone_t下面挂GridCoordinates_t再下面才是DataArray_t。你如果只用h5read拿到的是扁平数据集得自己重建父子关系。我试过用纯脚本读一个 12 万节点的四面体网格光解析节点树就花了 40 多秒内存峰值到 1.8 GB。MEX 方案直接调cg_open、cg_zone_read、cg_coord_read同样的网格 0.3 秒读完内存只多了一个 C 结构体的开销。所以选型逻辑很简单网格节点数超过 5 万或者需要写边界条件、流场解就上 MEX只是偶尔看一眼小网格的坐标纯脚本凑合能用。2.2 编译环境怎么配头文件、库路径和 mex 命令假设你已经从 CGNS 官方源码编好了动态库Windows 下得到cgns.dll和cgns.libLinux 下是libcgns.so。MATLAB 里用mex编译时关键是把头文件目录和库目录喂进去。下面是一个 Linux 下的编译命令示例Windows 把-I和-L换成对应路径即可。# 在 MATLAB 命令行里执行不是系统终端 mex -v ... -I/usr/local/include ... -L/usr/local/lib ... -lcgns ... cgns_mex.c-v是打开详细输出第一次编译一定要加能看到链接器到底找了哪些库。-I后面跟cgns.h所在目录-L跟libcgns.so所在目录-lcgns告诉链接器库名。编译成功后当前目录会出现cgns_mex.mexa64。如果报undefined reference to cg_open九成是-lcgns位置不对把它放到源文件后面。Windows 下还要注意位数匹配MATLAB 是 64 位CGNS 库也必须是 64 位混用 32 位会直接崩。2.3 读一个 Zone 的坐标和流场MEX 函数参数怎么传MEX 函数的入口是mexFunction里面用mxGetPr拿 MATLAB 数组指针调 CGNS C 接口再把结果塞回plhs。下面这段代码实现一个最小功能打开 CGNS 文件读第一个 Base 下第一个 Zone 的顶点坐标和名为Pressure的流场变量。#include mex.h #include cgnslib.h void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { char *filename mxArrayToString(prhs[0]); // 输入文件名 int fn, base, zone, coord; cgsize_t size[3]; double *x, *y, *z, *p; cg_open(filename, CG_MODE_READ, fn); // 打开文件 cg_base_read(fn, 1, NULL, NULL, NULL); // 读第一个 Base cg_zone_read(fn, 1, 1, NULL, size); // 读第一个 Zone 的尺寸 int nnode size[0]; // 顶点数 plhs[0] mxCreateDoubleMatrix(nnode, 1, mxREAL); plhs[1] mxCreateDoubleMatrix(nnode, 1, mxREAL); plhs[2] mxCreateDoubleMatrix(nnode, 1, mxREAL); x mxGetPr(plhs[0]); y mxGetPr(plhs[1]); z mxGetPr(plhs[2]); cg_coord_read(fn, 1, 1, CoordinateX, RealDouble, 1, nnode, x); cg_coord_read(fn, 1, 1, CoordinateY, RealDouble, 1, nnode, y); cg_coord_read(fn, 1, 1, CoordinateZ, RealDouble, 1, nnode, z); cg_field_read(fn, 1, 1, 1, Pressure, RealDouble, 1, nnode, p); cg_close(fn); }逻辑说明cg_open的第三个参数是模式读用CG_MODE_READ。cg_zone_read的size数组前三个元素是顶点数后面是单元数。cg_coord_read的第四个参数是坐标名CGNS 标准里就是CoordinateX/Y/Z大小写敏感。cg_field_read的第四个参数是解索引通常从 1 开始。参数坑点cgsize_t在 64 位编译下是long long如果你用int接会截断大网格的节点数。另外mxArrayToString拿到的字符串要手动mxFree否则反复调用会漏内存。2.4 写回 CGNS边界条件怎么挂到 Zone 上写文件比读多一步先建 Base、建 Zone、写坐标、写解最后挂边界条件。下面这段把 MATLAB 里的坐标和压力写成一个新 CGNS 文件并给第一个面加一个BCWall边界。cg_open(out.cgns, CG_MODE_WRITE, fn); cg_base_write(fn, Base, 3, 3, base); // 3 维3 维 cg_zone_write(fn, base, Zone1, size, Structured, zone); cg_coord_write(fn, base, zone, RealDouble, CoordinateX, x, coord); cg_coord_write(fn, base, zone, RealDouble, CoordinateY, y, coord); cg_coord_write(fn, base, zone, RealDouble, CoordinateZ, z, coord); cg_sol_write(fn, base, zone, FlowSolution, Vertex, sol); cg_field_write(fn, base, zone, sol, RealDouble, Pressure, p, field); cg_bc_write(fn, base, zone, 1, Wall, BCWall, bc); // 第一个面 cg_close(fn);cg_base_write的第二个参数是 Base 名第三、四个参数是物理维度和拓扑维度三维网格都填 3。cg_zone_write的size数组要按[nnode, ncell, 0]填。cg_sol_write的Vertex表示解存在顶点上如果存在单元中心就换CellCenter。cg_bc_write的第四个参数是面索引结构化网格从 1 开始。注意边界条件必须在解写完之后挂顺序反了 CGNS 库会报CG_ERROR。3. 参数设置与数据映射把 MATLAB 数组对齐到 CGNS 节点树3.1 维度顺序陷阱MATLAB 列优先 vs CGNS 行优先MATLAB 数组是列优先CGNS 底层按行优先存。一个nx * ny * nz的结构化网格在 MATLAB 里习惯写成X(i,j,k)但 CGNS 的cg_coord_write期望的是一维数组顺序是k变化最快。如果你直接把X(:)传进去读出来会发现坐标在 Z 方向错位。正确做法是用permute调整维度顺序再展平。% 假设 X 是 nx-by-ny-by-nz X_cgns permute(X, [3, 2, 1]); % 变成 nz-by-ny-by-nx X_flat X_cgns(:); % 展平后 k 变化最快参数说明permute的第二个参数是维度顺序[3,2,1]把原来的第三维换到第一维。展平后传给 MEX 函数。这个坑我踩过两次第一次读出来的网格像被拧过的麻花排查了一下午才发现是维度顺序反了。3.2 边界条件类型映射BCWall、BCFarfield 怎么选CGNS 标准定义了几十种边界条件类型常用的就几个。下面这张表是 MATLAB 侧字符串到 CGNS 枚举的映射写 MEX 时直接查表。MATLAB 侧字符串CGNS 枚举适用场景wallBCWall无滑移固壁farfieldBCFarfield远场边界inletBCInflow速度入口outletBCOutflow压力出口symmetryBCSymmetryPlane对称面选型逻辑如果做外流远场用BCFarfield物面用BCWall如果做内流入口用BCInflow并指定速度方向出口用BCOutflow。注意BCInflow需要额外写InflowBC_t结构体指定速度分量比BCWall多一层嵌套MEX 里要调cg_bc_dataset_write。3.3 流场变量命名Pressure、VelocityX 这些名字不能乱写CGNS 对解变量的名字没有强制标准但主流求解器都遵循一套约定。你写Pressure读的时候也得用Pressure写成pressure或P就找不到。下面是我常用的变量名对照。物理量推荐 CGNS 名数据类型压力PressureRealDoubleX 速度VelocityXRealDoubleY 速度VelocityYRealDoubleZ 速度VelocityZRealDouble温度TemperatureRealDouble密度DensityRealDouble如果读一个别人给的 CGNS 文件不确定变量名先用cg_nsols和cg_sol_info把解的名字列出来再决定读哪个。别硬猜猜错了返回全零你还以为数据没写进去。4. 避坑与排查CGNS 在 MATLAB 里翻车的五个典型场景4.1 现象cg_open返回错误码 0但文件明明存在原因路径里有中文或空格CGNS C 库对非 ASCII 路径支持不好。解决把文件挪到纯英文无空格目录或者用getShortPathName转短路径。我一般直接在tempdir下操作省事。4.2 现象读出来的坐标全是零但cg_coord_read没报错原因坐标名大小写写错了比如写了Coordinatex。CGNS 是大小写敏感的CoordinateX和Coordinatex是两个不同节点。解决用cg_coord_info把 Zone 里所有坐标名打印出来照着抄。4.3 现象MEX 编译通过一运行 MATLAB 就崩原因MEX 文件里用了mxGetPr拿复数数组的实部或者cgsize_t和int混用导致栈溢出。解决编译时加-largeArrayDims所有 CGNS 尺寸变量统一用cgsize_tMATLAB 数组索引用mwIndex。4.4 现象写边界条件时cg_bc_write返回CG_ERROR原因面索引超出 Zone 的面数或者边界条件类型和网格拓扑不匹配。比如二维网格写了三维的BCWall。解决先调cg_nbocos确认面数再检查cg_zone_type返回的是Structured还是Unstructured。4.5 现象读大网格时内存暴涨MATLAB 被系统杀掉原因MEX 里一次性mxCreateDoubleMatrix分配了所有节点数组加上 CGNS 库内部的缓冲区内存翻倍。解决分块读每次读 10 万节点用cg_coord_read的range_min和range_max参数控制范围MATLAB 侧用matfile追加写盘。5. 进阶技巧用 MATLAB 做 CGNS 网格质量检查与批量转换5.1 网格质量指标怎么在 MATLAB 里算读完坐标后用delaunay或delaunayn重建四面体连接关系然后算每个单元的雅可比行列式。下面这段算最小雅可比和长宽比。% 假设 X, Y, Z 是顶点坐标T 是四面体索引矩阵 dt delaunayTriangulation(X, Y, Z); T dt.ConnectivityList; J zeros(size(T,1), 1); for i 1:size(T,1) v1 [X(T(i,1)), Y(T(i,1)), Z(T(i,1))]; v2 [X(T(i,2)), Y(T(i,2)), Z(T(i,2))]; v3 [X(T(i,3)), Y(T(i,3)), Z(T(i,3))]; v4 [X(T(i,4)), Y(T(i,4)), Z(T(i,4))]; J(i) abs(det([v2-v1; v3-v1; v4-v1])); end minJ min(J);delaunayTriangulation对大规模点集比较慢超过 50 万点建议用 MEX 调 CGNS 自带的cg_conn_read直接拿连接关系省去重建。det算的是平行六面体体积除以 6 才是四面体体积但比较相对质量时不用除。5.2 批量转换把多个 CGNS 文件转成 MAT 格式如果手头有几十个 CGNS 文件要统一处理写个循环调 MEX 函数每个文件存一个.mat。注意每次循环后clear掉大变量否则内存会累积。files dir(*.cgns); for k 1:length(files) [x, y, z, p] cgns_mex(files(k).name); save([files(k).name, .mat], x, y, z, p, -v7.3); clear x y z p; end-v7.3是必须的普通-v7格式单个变量不能超过 2 GB。clear放在save之后确保写盘完成再释放。5.3 一个我常用的验证习惯每次写完 CGNS 文件不要直接信。用 MEX 读回来和原始 MATLAB 数组做isequal或norm差。我一般会算坐标的max(abs(x_read - x_orig))超过 1e-12 就说明写入或读取环节有精度损失。这个习惯帮我逮到过好几次维度顺序错误和数据类型不匹配。希望帮到你。本文还有配套的精品资源点击获取
