Mercurial > repos > public > sbplib

annotate +noname/calcSol.m @ 970:23d9ca6755be feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator.
author Martin Almquist <malmquist@stanford.edu>
date Tue, 25 Dec 2018 07:23:38 +0100
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