Mercurial > repos > public > sbplib

annotate +noname/calcSol.m @ 1299:73e52c74baac feature/boundary_optimized_grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Closed branch feature/boundary_optimized_grids
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 08 Jul 2020 18:20:57 +0200
parents 127e562bd6d3
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
243
127e562bd6d3 Added function to calculate a solution at a given time.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1 % Calculates the solution of a discr at a given time using aligned timesteps.
127e562bd6d3 Added function to calculate a solution at a given time.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
2 % Returns the solution as a grid function as defined in +grid
127e562bd6d3 Added function to calculate a solution at a given time.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
3 function gf = calculateSolution(discr, T, tsOpt)
127e562bd6d3 Added function to calculate a solution at a given time.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4 k_max = discr.getTimestep(tsOpt);
127e562bd6d3 Added function to calculate a solution at a given time.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5 [k,N] = alignedTimestep(k_max,T);
127e562bd6d3 Added function to calculate a solution at a given time.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6 tsOpt.k = k;
127e562bd6d3 Added function to calculate a solution at a given time.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7 ts = discr.getTimestepper(tsOpt);
127e562bd6d3 Added function to calculate a solution at a given time.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
8
127e562bd6d3 Added function to calculate a solution at a given time.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
9 gf = ts.stepN(N-ts.n);
127e562bd6d3 Added function to calculate a solution at a given time.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10 end