convergence_analysis.m File Reference

Functions
	errors (i, j)
	fprintf ('h=%.6f, final L2_rho=%.3e\n', h, data(end, 7))

Variables
	folders = {'tri', 'quad', 'hexa'}
	Nfolders = 3
	Nfiles = 4
	errors = zeros(Nfolders,Nfiles)
for	i
end end	filename = 'quad/timehistory_GuermondQuartapelle_k0_cart10x10_dt=0.000500.txt'
Extract mesh size	text = fileread(filename)
	h = str2double(regexp(text, 'MeshSize:\s*([0-9.e+-]+)', 'tokens', 'once'))
Load data matrix	data = readmatrix(filename, 'CommentStyle', '#')

Function Documentation

errors()

errors	(	i	,
		j
	)

fprintf()

fprintf	(	'	h = %.6f,
		final	L2_rho = %.3e\n',
		h	,
		data(end, 7)
	)

Variable Documentation

data

Load data matrix data = readmatrix(filename, 'CommentStyle', '#')

errors

errors = zeros(Nfolders,Nfiles)

filename

end end filename = 'quad/timehistory_GuermondQuartapelle_k0_cart10x10_dt=0.000500.txt'

folders

folders = {'tri', 'quad', 'hexa'}

h

h = str2double(regexp(text, 'MeshSize:\s*([0-9.e+-]+)', 'tokens', 'once'))

i

for i

Initial value:

=1:Nfolders
    for j=1:Nfiles
        data = extract_data_from_file(file_name)

Nfiles

Nfiles = 4

Nfolders

Nfolders = 3

text

Extract mesh size text = fileread(filename)