Mercurial > repos > public > sbplib

annotate vandermonde.m @ 768:dec0447cbf2c feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Make assertSize more general, allow picking dimensions
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 26 Jun 2018 14:07:35 +0200
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