Mercurial > repos > public > sbplib

annotate +parametrization/equal_step_size.m @ 724:a9e8c9d71307 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Modify penalty for Neumann in Head2d so that data is for normal derivative and not u_x or u_y.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 14 Mar 2018 13:04:25 -0700
parents 3a3cf386bb7e
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1 % Calculates M so that the stepsize m/M is as close to n/M as possible
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
2 function M = equal_step_size(n,N,m)
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
3 M = round((m*(N-1)+n)/n);
48b6fb693025 Initial commit.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4 end