GetMesh.m

Go to the documentation of this file.

function Mesh=GetMesh(cFileName)
% function Mesh=GetMesh(cFileName)
% Initialization of the Mesh structure from a FreeFEM++ mesh file
%
% Description plus longue...
%
% Parameters:
%  cFileName: file name (string)
%
% Return values:
%  Mesh: mesh structure
%
% Generated fields of Mesh:
%  q: @code $\int_a^b f(x) dx$ @endcode
%  me:
%
% usage:  Th=GetMesh('carre.msh')
% This program is open source.  For license terms, see the COPYING file.
%
% --------------------------------------------------------------------
% ATTRIBUTION NOTICE:
% This product includes software developed for the RBmatlab project at
% (C) Universities of Stuttgart and Münster, Germany.
%
% RBmatlab is a MATLAB software package for model reduction with an
% emphasis on Reduced Basis Methods. The project is maintained by
% M. Dihlmann, M. Drohmann, B. Haasdonk, M. Ohlberger and M. Schaefer.
% For Online Documentation and Download we refer to www.morepas.org.
% --------------------------------------------------------------------
[fid,message]=fopen(cFileName,'r');
if ( fid == -1 )
  error([message,' : ',cFileName]);
end
n=fscanf(fid,'%d',3); ; % n(1) -> number of vertices    
                      % n(2) -> number of triangles     
                      % n(3) -> number of boundary edges
me=zeros(3,n(2));
ql=zeros(1,n(1));
q=zeros(2,n(1));
mel=zeros(1,n(2));
for i=1:n(1)
  d_tmp=fscanf(fid,'%g',2);
  q(1,i)=d_tmp(1);
  q(2,i)=d_tmp(2);
  ql(i)=fscanf(fid,'%d',1);
end
for i=1:n(2)
  i_tmp=fscanf(fid,'%d',4);
  me(1,i)=i_tmp(1);
  me(2,i)=i_tmp(2);
  me(3,i)=i_tmp(3);
  mel(i)=i_tmp(4);
end

for i=1:n(3)
  i_tmp=fscanf(fid,'%d',3);
  be(1,i)=i_tmp(1);
  be(2,i)=i_tmp(2);
  bel(i) =i_tmp(3);
end

fclose(fid);

Mesh=struct('q',q,'me',me,'ql',ql,'mel',mel,'be',be,'bel',bel, ...
            'nq',size(q,2), ...
            'nme',size(me,2), ...
            'nbe',size(be,2), ...
            'areas',ComputeArea(q,me),...
            'lbe',EdgeLength(be,q));