ReadMeshFromGambit1D
所属分类:matlab编程
开发工具:matlab
文件大小:1KB
下载次数:8
上传日期:2012-07-02 22:29:44
上 传 者:
lasido
说明: 从网格划分软件Gambit读取1维网格的Matlab程序,来自《节点间断迦辽金方法》一书
(Matlab program read from the meshing software Gambit one-dimensional grid, from the node discontinuous Galerkin method, " a book)
function [Nv, VX, K, EToV] = MeshReader1DGambit(FileName)
% function [Nv, VX, K, EToV] = MeshReader1DGambit(FileName)
% Purpose : Read in basic grid information to build grid
%
% NOTE : gambit(Fluent, Inc) *.neu format is assumed
% Open file and read intro
Fid = fopen(FileName, 'rt');
for i=1:6; line = fgetl(Fid); end
% Find number of nodes and number of elements
dims = fscanf(Fid, '%d');
Nv = dims(1); K = dims(2);
for i=1:2; line = fgetl(Fid); end
% read node coordinates
VX = (1:Nv);
for i = 1:Nv
line = fgetl(Fid);
tmpx = sscanf(line, '%lf');
VX(i) = tmpx(2);
end
for i=1:2; line = fgetl(Fid); end
% read element to node connectivity
EToV = zeros(K, 2);
for k = 1:K
line = fgetl(Fid);
tmpcon = sscanf(line, '%lf');
EToV(k,1) = tmpcon(4); EToV(k,2) = tmpcon(5);
end
% Close file
st = fclose(Fid);
return
近期下载者:
相关文件:
收藏者: