Mercurial > repos > public > sbplib
annotate vandermonde.m @ 721:c29eca357028 feature/grids
Change name of grid parameter in multiblock.DiffOp to avoid conflict with package
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Mon, 19 Mar 2018 17:11:15 +0100 |
parents | 184833fe4c0e |
children | 0090a86d8b72 |
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 |