%%%%%%%%%%%%%%%%%%%%% Word driven tesselation  %%%%%%%%%%%%%%%%%%%
%
% TTESS
% Usage: TTESS(p,q,r,word,ti,drawIT,T0)
%
% Inputs:  p,q,r  - real numbers with 1/p +1/q +1/r < 1 
%                   see SCHWARTZ for the interpretation of p,q,r
%          word   - the sequence of reflections to drive the tesellation 
%          ti     - title of tesselation
%          drawIT - draw initial triangle
%          T0     - optional starting triangle  
%
%          the sequence of reflections given by word are used to tile the disc
%          with a p-q-r schwartz triangle

function T = ttess(p,q,r,word,ti,drawIT, T0)

% colour plotting options
vc1 = 'r';
vc2 = 'y';
vc3 = 'g';
cc  = 'y';
icc = 'k';
invisible = get(gcf,'Color');
verttxt =  [inttostr(p)  inttostr(q)  inttostr(r)]; 

% initializations
I = sqrt(-1);  

if (nargin >= 6)
     DIT = drawIT;
    else 
     DIT = 1;
   end; 
   
if (nargin >= 7)
     T = T0;
    else 
     T =schwartz(p,q,r);
   end; 


boundary(ti,'','w'); 
if DIT plottri(T,'w'); end;   
%plotvert(T,vc1,vc2,vc3);    % remove comment character to plot vertices in colour 
%plotvrtt(T,verttxt);        % remove comment character to plot vertex orders 

word = word(:);
n = length(word);
for i = 1:n
    side = word(i);
    T2 = trirefl(T, side);
    if (side >= 'a')&(side <= 'c')
        plottri(T2,'w');
        %plotvert(T2,vc1,vc2,vc3);  % remove comment character to plot vertices in colour 
        %plotvrtt(T2,verttxt);      % remove comment character to plot vertex orders  
        end;
        T = T2;
    end;