home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 2.ddi / MUTOOLS2.DI$ / INITT.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  1.5 KB  |  73 lines
  1. function [tl,tm,tr,cl,cm,cr,nop] = initt(mat,blk,ts)
  2.  [nr,nc] = size(mat);
  3.  inc = max([abs(ts) 0.3]);
  4.  nop = 0;
  5.  if ts < 0
  6.    tl = 2*ts;
  7.    tm = ts;
  8.    tr = 0;
  9.  elseif ts > 0
  10.    tl = 0;
  11.    tm = ts;
  12.    tr = 2*ts;
  13.  else
  14.    tl = -0.001;
  15.    tm = 0;
  16.    tr = 0.001;
  17.    inc = 0.3;
  18.  end
  19.  mt = mat;
  20.  ml = mat;
  21.  mm = mat;
  22.  mr = mat;
  23.  ur = mat(1:blk(1,2),blk(1,1)+1:nc);
  24.  ll = mat(blk(1,2)+1:nr,1:blk(1,1));
  25.  ml(1:blk(1,2),blk(1,1)+1:nc) = exp(tl)*ur;
  26.  ml(blk(1,2)+1:nr,1:blk(1,1)) = exp(-tl)*ll;
  27.  mr(1:blk(1,2),blk(1,1)+1:nc) = exp(tr)*ur;
  28.  mr(blk(1,2)+1:nr,1:blk(1,1)) = exp(-tr)*ll;
  29.  mm(1:blk(1,2),blk(1,1)+1:nc) = exp(tm)*ur;
  30.  mm(blk(1,2)+1:nr,1:blk(1,1)) = exp(-tm)*ll;
  31.  cl = norm(ml);
  32.  cm = norm(mm);
  33.  cr = norm(mr);
  34.  done = 0;
  35.  while done == 0
  36.    if cl >= cm
  37.      if cm <= cr
  38.        done = 1;
  39.      else
  40.        tl = tm;
  41.        cl = cm;
  42.        tm = tr;
  43.        cm = cr;
  44.        bestt = tr;
  45.        tr = tr+inc;
  46.        mr(1:blk(1,2),blk(1,1)+1:nc) = exp(tr)*ur;
  47.        mr(blk(1,2)+1:nr,1:blk(1,1)) = exp(-tr)*ll;
  48.        cr = norm(mr);
  49.      end
  50.    else
  51.      tr = tm;
  52.      cr = cm;
  53.      tm = tl;
  54.      cm = cl;
  55.      bestt = tl;
  56.      tl = tl-inc;
  57.      ml(1:blk(1,2),blk(1,1)+1:nc) = exp(tl)*ur;
  58.      ml(blk(1,2)+1:nr,1:blk(1,1)) = exp(-tl)*ll;
  59.      cl = norm(ml);
  60.    end
  61.    if abs(bestt) > 8
  62.      done = 1;
  63.      tl = bestt;
  64.      tm = bestt;
  65.      tr = bestt;
  66.      cl = cm;
  67.      cr = cm;
  68.      nop = -1;
  69.    end
  70.  end
  71. %
  72. % Copyright MUSYN INC 1991,  All Rights Reserved
  73.