Mercurial > repos > public > sbplib

annotate vandermonde.m @ 755:14f0058356f2 feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with feature/grids
author Martin Almquist <malmquist@stanford.edu>
date Fri, 15 Jun 2018 18:10:26 -0700
parents 184833fe4c0e
children 0090a86d8b72
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
660
184833fe4c0e Add monomial, vandermonde and stencil cretation tools
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1 % Create vandermonde matrix for points x and polynomials of order p
184833fe4c0e Add monomial, vandermonde and stencil cretation tools
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
2 % x and p are vectors
184833fe4c0e Add monomial, vandermonde and stencil cretation tools
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
3 % v is a length(x) by length(p) matrix
184833fe4c0e Add monomial, vandermonde and stencil cretation tools
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4 function V = vandermonde(x, p)
184833fe4c0e Add monomial, vandermonde and stencil cretation tools
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5 V = sym(zeros(length(x), length(p))); % Is there a way to make this work for both double and sym
184833fe4c0e Add monomial, vandermonde and stencil cretation tools
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6
184833fe4c0e Add monomial, vandermonde and stencil cretation tools
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7 for i = 1:length(p)
184833fe4c0e Add monomial, vandermonde and stencil cretation tools
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
8 V(:, i) = mononomial(x,p(i));
184833fe4c0e Add monomial, vandermonde and stencil cretation tools
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
9 end
184833fe4c0e Add monomial, vandermonde and stencil cretation tools
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10 end