Rev | Line | |
---|
[37] | 1 | function monom = genmonoms(newton_m,x); |
---|
| 2 | %GENMONOMS Internal function used in SOS programs |
---|
| 3 | |
---|
| 4 | % Author Johan Löfberg |
---|
| 5 | % $Id: genmonoms.m,v 1.1 2006/03/30 13:56:54 joloef Exp $ |
---|
| 6 | |
---|
| 7 | |
---|
| 8 | precalc = []; |
---|
| 9 | monom = []; |
---|
| 10 | for i = 1:size(newton_m,2) |
---|
| 11 | powers = unique(newton_m(:,i)); |
---|
| 12 | for j = 1:length(powers) |
---|
| 13 | precalc{1+powers(j),i} = x(i)^powers(j); |
---|
| 14 | end |
---|
| 15 | end |
---|
| 16 | |
---|
| 17 | monom = []; |
---|
| 18 | for i = 1:size(newton_m,1) |
---|
| 19 | temp = 1; |
---|
| 20 | for j = 1:size(newton_m,2) |
---|
| 21 | temp = temp*precalc{newton_m(i,j)+1,j};%x(j)^newton_m(i,j); |
---|
| 22 | end |
---|
| 23 | monom = [monom;temp]; |
---|
| 24 | end |
---|
Note: See
TracBrowser
for help on using the repository browser.