Mercurial > repos > public > sbplib

annotate +noname/calculateErrors.m @ 657:b59345f905f0 feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Make noname.calculateErrors parallel
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 24 Oct 2017 16:32:00 +0200
parents 818d52d4928f
children 1201eb16557e
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
616
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1 % [discr, trueSolution] = schemeFactory(m)
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
2 % where trueSolution should be a timeSnapshot of the true solution a time T
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
3 % T is the end time
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4 % m are grid size parameters.
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5 % N are number of timesteps to use for each gird size
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6 % timeOpt are options for the timeStepper
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7 function e = calculateErrors(schemeFactory, T, m, N, errorFun, timeOpt)
657
b59345f905f0 Make noname.calculateErrors parallel
Jonatan Werpers <jonatan@werpers.com>
parents: 616
diff changeset
8 %TODO: Ability to choose paralell or not
616
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
9 assertType(schemeFactory, 'function_handle');
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10 assertNumberOfArguments(schemeFactory, 1);
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
11 assertScalar(T);
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
12 assert(length(m) == length(N), 'Vectors m and N must have the same length');
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
13 assertType(errorFun, 'function_handle');
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
14 assertNumberOfArguments(errorFun, 2);
657
b59345f905f0 Make noname.calculateErrors parallel
Jonatan Werpers <jonatan@werpers.com>
parents: 616
diff changeset
15 default_arg('timeOpt', struct());
b59345f905f0 Make noname.calculateErrors parallel
Jonatan Werpers <jonatan@werpers.com>
parents: 616
diff changeset
16
616
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
17
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
18 e = [];
657
b59345f905f0 Make noname.calculateErrors parallel
Jonatan Werpers <jonatan@werpers.com>
parents: 616
diff changeset
19 parfor i = 1:length(m)
616
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
20 done = timeTask('m = %3d ', m(i));
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
21
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
22 [discr, trueSolution] = schemeFactory(m(i));
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
23
657
b59345f905f0 Make noname.calculateErrors parallel
Jonatan Werpers <jonatan@werpers.com>
parents: 616
diff changeset
24 timeOptTemp = timeOpt;
b59345f905f0 Make noname.calculateErrors parallel
Jonatan Werpers <jonatan@werpers.com>
parents: 616
diff changeset
25 timeOptTemp.k = T/N(i);
b59345f905f0 Make noname.calculateErrors parallel
Jonatan Werpers <jonatan@werpers.com>
parents: 616
diff changeset
26 ts = discr.getTimestepper(timeOptTemp);
616
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
27 ts.stepTo(N(i), true);
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
28 approxSolution = discr.getTimeSnapshot(ts);
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
29
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
30 e(i) = errorFun(trueSolution, approxSolution);
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
31
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
32 fprintf('e = %.4e', e(i))
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
33 done()
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
34 end
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
35 fprintf('\n')
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
36 end
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
37
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
38
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
39 %% Example error function
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
40 % u_true = grid.evalOn(dr.grid, @(x,y)trueSolution(T,x,y));
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
41 % err = u_true-u_false;
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
42 % e(i) = norm(err)/norm(u_true);
818d52d4928f Add helper function for convergence runs on discretizations
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
43 % % e(i) = sqrt(err'*d.H*d.J*err/(u_true'*d.H*d.J*u_true));