home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 5.ddi / SPLINES.DI$ / PPMAK.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  1.8 KB  |  50 lines
  1. function pp=ppmak(breaks,coefs,d)
  2. % PPMAK    Make a pp.
  3. %
  4. %     pp = ppmak(breaks,coefs[,d])
  5. %
  6. %  puts together a pp function from the breaks and coefficients input or
  7. %  requested. 
  8. % 
  9. %    If  d  is not explicitly given, program expects  coefs(d,(k,l)) ,
  10. %  and infers  l := length(breaks)-1  and then  k := (# cols(coefs))/l .
  11. % 
  12. %    If  d  is explicitly given, program expects  coefs((d,l),k) , and
  13. %  infers  k := # cols(coefs), and  l := (# rows(coefs))/d .
  14.  
  15. % C. de Boor / latest change: June 1, 1989
  16. % Copyright (c) 1990-92 by Carl de Boor and The MathWorks, Inc.
  17.  
  18.  
  19. if (nargin==0);
  20.    breaks=input('Give the (l+1)-vector of breaks  >');  
  21.    coefs=input('Give the (d by (k*l)) matrix of local pol. coefficients  >');
  22. end
  23. if (nargin<3),
  24.    [d,kl]=size(coefs);
  25.    if (d==0),
  26.       error('The coefficient sequence is empty!')
  27.    end
  28.    l=length(breaks)-1;k=fix(kl/l);
  29.    if (k<=0)|(k*l~=kl);
  30.       fprintf('The given number %.0f of polynomial pieces is incompatible',l)
  31.       fprintf(' with the total number %.0f of coefficients supplied!\n',kl)
  32.       error('')
  33.    elseif (~isempty(find(diff(breaks)<0))),
  34.       error('The breakpoint sequence should be nondecreasing!')
  35.    else
  36.       % the pp-format expects coefs in array  (d*l) by k, while the standard
  37.       % input supplies them in an array d by (k*l) . This requires the
  38.       % following shuffling, from  D+d(-1+K + k(-1+L))=D-d +(K-k)d + dkL
  39.       % to  D+d(-1+L + l(-1+K)=D-d +(L-l)d + dlK .
  40.       c=coefs(:);([1-k:0]'*ones(1,l)+k*ones(k,1)*[1:l])';
  41.       coefs=[1-d:0]'*ones(1,kl)+d*ones(d,1)*(ans(:)');
  42.       coefs(:)=c(coefs);
  43.    end
  44. else
  45.    [dl,k]=size(coefs); l = dl/d;
  46.    if (l+1~=length(breaks)|dl~=l*d), 
  47.       error('input arguments have incompatible sizes.'),return,end
  48. end
  49. pp=[10 d l breaks(:)' k coefs(:)'];
  50.